Dissecting protein-induced DNA looping dynamics in real time

Many proteins that interact with DNA perform or enhance their specific functions by binding simultaneously to multiple target sites, thereby inducing a loop in the DNA. The dynamics and energies involved in this loop formation influence the reaction mechanism. Tethered particle motion has proven a powerful technique to study in real time protein-induced DNA looping dynamics while minimally perturbing the DNA–protein interactions. In addition, it permits many single-molecule experiments to be performed in parallel. Using as a model system the tetrameric Type II restriction enzyme SfiI, that binds two copies of its recognition site, we show here that we can determine the DNA–protein association and dissociation steps as well as the actual process of protein-induced loop capture and release on a single DNA molecule. The result of these experiments is a quantitative reaction scheme for DNA looping by SfiI that is rigorously compared to detailed biochemical studies of SfiI looping dynamics. We also present novel methods for data analysis and compare and discuss these with existing methods. The general applicability of the introduced techniques will further enhance tethered particle motion as a tool to follow DNA–protein dynamics in real time.     

Solution parameters modulating DNA binding specificity of the restriction endonuclease EcoRV

The DNA binding stringency of restriction endonucleases is crucial for their proper function. The X-ray structures of the specific and non-cognate complexes of the restriction nuclease EcoRV are considerably different suggesting significant differences in the hydration and binding free energies. Nonetheless, the majority of studies performed at pH 7.5, optimal for enzymatic activity, have found a < 10-fold difference between EcoRV binding constants to the specific and nonspecific sequences in the absence of divalent ions. We used a recently developed self-cleavage assay to measure EcoRV-DNA competitive binding and to evaluate the influence of water activity, pH and salt concentration on the binding stringency of the enzyme in the absence of divalent ions. We find the enzyme can readily distinguish specific and nonspecific sequences. The relative specific-nonspecific binding constant increases strongly with increasing neutral solute concentration and with decreasing pH. The difference in number of associated waters between specific and nonspecific DNA-EcoRV complexes is consistent with the differences in the crystal structures. Despite the large pH dependence of the sequence specificity, the osmotic pressure dependence indicates little change in structure with pH. The large osmotic pressure dependence means that measurement of protein-DNA specificity in dilute solution cannot be directly applied to binding in the crowded environment of the cell. In addition to divalent ions, water activity and pH are key parameters that strongly modulate binding specificity of EcoRV.     

Participation of water in Hin recombinase--DNA recognition.

The participation of water molecules in the interaction between the Hin recombinase and its operator DNA has been detected by analysis of the dissociation constant in the presence of varying concentrations of neutral solutes and cosolvents. The dissociation constant as measured by gel mobility shift assays increased as the concentration of dimethyl sulfoxide, glycerol, sucrose, or polyethylene glycol was increased. Osmotic pressure is the only property that correlates with the change in the dissociation constant for all compounds. This data indicates that binding of a small population of water molecules accompanies formation of the Hin-DNA complex, and points to a novel role for solvent molecules in assisting site specific interaction between DNA-binding proteins and their cognate DNA sequence.     

Inferring primase-DNA specific recognition using a data driven approach

DNA–protein interactions play essential roles in all living cells. Understanding of how features embedded in the DNA sequence affect specific interactions with proteins is both challenging and important, since it may contribute to finding the means to regulate metabolic pathways involving DNA–protein interactions. Using a massive experimental benchmark dataset of binding scores for DNA sequences and a machine learning workflow, we describe the binding to DNA of T7 primase, as a model system for specific DNA–protein interactions. Effective binding of T7 primase to its specific DNA recognition sequences triggers the formation of RNA primers that serve as Okazaki fragment start sites during DNA replication.     

From Nonspecific DNA–Protein Encounter Complexes to the Prediction of DNA–Protein Interactions

DNA–protein interactions are involved in many essential biological activities. Because there is no simple mapping code between DNA base pairs and protein amino acids, the prediction of DNA–protein interactions is a challenging problem. Here, we present a novel computational approach for predicting DNA-binding protein residues and DNA–protein interaction modes without knowing its specific DNA target sequence. Given the structure of a DNA-binding protein, the method first generates an ensemble of complex structures obtained by rigid-body docking with a nonspecific canonical B-DNA. Representative models are subsequently selected through clustering and ranking by their DNA–protein interfacial energy. Analysis of these encounter complex models suggests that the recognition sites for specific DNA binding are usually favorable interaction sites for the nonspecific DNA probe and that nonspecific DNA–protein interaction modes exhibit some similarity to specific DNA–protein binding modes. Although the method requires as input the knowledge that the protein binds DNA, in benchmark tests, it achieves better performance in identifying DNA-binding sites than three previously established methods, which are based on sophisticated machine-learning techniques. We further apply our method to protein structures predicted through modeling and demonstrate that our method performs satisfactorily on protein models whose root-mean-square Cα deviation from native is up to 5 Å from their native structures. This study provides valuable structural insights into how a specific DNA-binding protein interacts with a nonspecific DNA sequence. The similarity between the specific DNA–protein interaction mode and nonspecific interaction modes may reflect an important sampling step in search of its specific DNA targets by a DNA-binding protein.     

Ultrasensitive isolation,identification and quantification of DNA–protein adducts by ELISA-based RADAR assay

Enzymes that form transient DNA–protein covalent complexes are targets for several potent classes of drugs used to treat infectious disease and cancer, making it important to establish robust and rapid procedures for analysis of these complexes. We report a method for isolation of DNA–protein adducts and their identification and quantification, using techniques compatible with high-throughput screening. This method is based on the RADAR assay for DNA adducts that we previously developed (Kiianitsa and Maizels (2013) A rapid and sensitive assay for DNA–protein covalent complexes in living cells. Nucleic Acids Res., 41:e104), but incorporates three key new steps of broad applicability. (i) Silica-assisted ethanol/isopropanol precipitation ensures reproducible and efficient recovery of DNA and DNA–protein adducts at low centrifugal forces, enabling cell culture and DNA precipitation to be carried out in a single microtiter plate. (ii) Rigorous purification of DNA–protein adducts by a procedure that eliminates free proteins and free nucleic acids, generating samples suitable for detection of novel protein adducts (e.g. by mass spectroscopy). (iii) Identification and quantification of DNA–protein adducts by direct ELISA assay. The ELISA-based RADAR assay can detect Top1–DNA and Top2a–DNA adducts in human cells, and gyrase–DNA adducts in Escherichia coli. This approach will be useful for discovery and characterization of new drugs to treat infectious disease and cancer, and for development of companion diagnostics assays for individualized medicine.     

A capture approach for supercoiled plasmid DNA using a triplex-forming oligonucleotide

Proteins that recognize and bind specific sites in DNA are essential for regulation of numerous biological functions. Such proteins often require a negative supercoiled DNA topology to function correctly. In current research, short linear DNA is often used to study DNA–protein interactions. Although linear DNA can easily be modified, for capture on a surface, its relaxed topology does not accurately resemble the natural situation in which DNA is generally negatively supercoiled. Moreover, specific binding sequences are flanked by large stretches of non-target sequence in vivo. Here, we present a straightforward method for capturing negatively supercoiled plasmid DNA on a streptavidin surface. It relies on the formation of a temporary parallel triplex, using a triple helix forming oligonucleotide containing locked nucleic acid nucleotides. All materials required for this method are commercially available. Lac repressor binding to its operator was used as model system. Although the dissociation constants for both the linear and plasmid-based operator are in the range of 4 nM, the association and dissociation rates of Lac repressor binding to the plasmid-based operator are ∼18 times slower than on a linear fragment. This difference underscores the importance of using a physiologically relevant DNA topology for studying DNA–protein interactions.     

Selective Glutaraldehyde-Mediated Coupling of Proteins to the 3′-Adenine Terminus of Polymerase Chain Reaction Products

Attachment of proteins to the 3′ end of DNA increases stability of the DNA in serum and retards clearance of DNA by major organs, thereby enhancing in vivo half-life and therapeutic potential of DNA. Unfortunately, the length of DNA molecules that can be produced with 3 ′ modifications by solid-phase synthesis for protein attachment is limited to 45–60 nucleotides due to uncertainties about sequence fidelity for longer oligonucleotides. Here we describe selective covalent coupling of proteins or other molecules to the 3′-adenine overhang of unlabeled and fluorophore-labeled double-stranded polymerase chain reaction products putatively at the N⁶ position of adenine using 2.5% glutaraldehyde at pH 6.0 and 4°C for at least 16 h. Gel mobility shift analyses and fluorescence analyses of the shifted bands supported conjugate formation between double-stranded polymerase chain reaction products and β2-microglobulin. In addition, blunt-ended DNA ladder fragments treated with glutaraldehyde at 4°C showed no evidence of DNA–DNA or DNA–protein conjugate formation. With the present cold glutaraldehyde technique, longer DNA–3′-protein conjugates might be easily mass-produced. The protein portion of a DNA–3′-protein conjugate could possess functionality as well, such as receptor binding for cell entry, cytotoxicity, or opsonization.     

Bacteriophage T4 endonuclease II,a promiscuous GIY-YIG nuclease,binds as a tetramer to two DNA substrates

The oligomerization state and mode of binding to DNA of the GIY-YIG endonuclease II (EndoII) from bacteriophage T4 was studied using gel filtration and electrophoretic mobility shift assays with a set of mutants previously found to have altered enzyme activity. At low enzyme/DNA ratios all mutants except one bound to DNA only as tetramers to two DNA substrates. The putatively catalytic E118 residue actually interfered with DNA binding (possibly due to steric hindrance or repulsion between the glutamate side chain and DNA), as shown by the ability of E118A to bind stably also as monomer or dimer to a single substrate. The tetrameric structure of EndoII in the DNA–protein complex is surprising considering the asymmetry of the recognized sequence and the predominantly single-stranded nicking. Combining the results obtained here with those from our previous in vivo studies and the recently obtained crystal structure of EndoII E118A, we suggest a model where EndoII translocates DNA between two adjacent binding sites and either nicks one strand of one or both substrates bound by the tetramer, or nicks both strands of one substrate. Thus, only one or two of the four active sites in the tetramer is catalytically active at any time.     

Mechanism of the formation of DNA-protein cross-links by antitumor cisplatin

DNA–protein cross-links are formed by various DNA-damaging agents including antitumor platinum drugs. The natures of these ternary DNA–Pt–protein complexes (DPCLs) can be inferred, yet much remains to be learned about their structures and mechanisms of formation. We investigated the origin of these DPCLs and their cellular processing on molecular level using gel electrophoresis shift assay. We show that in cell-free media cisplatin [cis-diamminedichloridoplatinum(II)] forms DPCLs more effectively than ineffective transplatin [trans-diamminedichloridoplatinum(II)]. Mechanisms of transformation of individual types of plain DNA adducts of the platinum complexes into the DPCLs in the presence of several DNA-binding proteins have been also investigated. The DPCLs are formed by the transformation of DNA monofunctional and intrastrand cross-links of cisplatin. In contrast, interstrand cross-links of cisplatin and monofunctional adducts of transplatin are stable in presence of the proteins. The DPCLs formed by cisplatin inhibit DNA polymerization or removal of these ternary lesions from DNA by nucleotide excision repair system more effectively than plain DNA intrastrand or monofunctional adducts. Thus, the bulky DNA–protein cross-links formed by cisplatin represent a more distinct and persisting structural motif recognized by the components of downstream cellular systems processing DNA damage considerably differently than the plain DNA adducts of this metallodrug.     

DNA-tension dependence of restriction enzyme activity reveals mechanochemical properties of the reaction pathway

Type II restriction endonucleases protect bacteria against phage infections by cleaving recognition sites on foreign double-stranded DNA (dsDNA) with extraordinary specificity. This capability arises primarily from large conformational changes in enzyme and/or DNA upon target sequence recognition. In order to elucidate the connection between the mechanics and the chemistry of DNA recognition and cleavage, we used a single-molecule approach to measure rate changes in the reaction pathway of EcoRV and BamHI as a function of DNA tension. We show that the induced-fit rate of EcoRV is strongly reduced by such tension. In contrast, BamHI is found to be insensitive, providing evidence that both substrate binding and hydrolysis are not influenced by this force. Based on these results, we propose a mechanochemical model of induced-fit reactions on DNA, allowing determination of induced-fit rates and DNA bend angles. Finally, for both enzymes a strongly decreased association rate is obtained on stretched DNA, presumably due to the absence of intradomain dissociation/re-association between non-specific sites (jumping). The obtained results should apply to many other DNA-associated proteins.     

Massively parallel characterization of restriction endonucleases

Restriction endonucleases are highly specific in recognizing the particular DNA sequence they act on. However, their activity is affected by sequence context, enzyme concentration and buffer composition. Changes in these factors may lead to either ineffective cleavage at the cognate restriction site or relaxed specificity allowing cleavage of degenerate ‘star’ sites. Additionally, uncharacterized restriction endonucleases and engineered variants present novel activities. Traditionally, restriction endonuclease activity is assayed on simple substrates such as plasmids and synthesized oligonucleotides. We present and use high-throughput Illumina sequencing-based strategies to assay the sequence specificity and flanking sequence preference of restriction endonucleases. The techniques use fragmented DNA from sequenced genomes to quantify restriction endonuclease cleavage on a complex genomic DNA substrate in a single reaction. By mapping millions of restriction site–flanking reads back to the Escherichia coli and Drosophila melanogaster genomes we were able to quantitatively characterize the cognate and star site activity of EcoRI and MfeI and demonstrate genome-wide decreases in star activity with engineered high-fidelity variants EcoRI-HF and MfeI-HF, as well as quantify the influence on MfeI cleavage conferred by flanking nucleotides. The methods presented are readily applicable to all type II restriction endonucleases that cleave both strands of double-stranded DNA.     

Analytical biochemistry of DNA–protein assemblies from crude cell extracts

Purification of specific DNA–protein complexes is a challenging task, as the involved interactions can be both electrostatic/H-bond and hydrophobic. The chromatographic stringency needed to obtain reasonable purifications uses salts and detergents. However, these components elicit the removal of proteins unspecifically bound to the chromatographic support itself, thus contaminating the purification products. In this work, a photocleavable linker connected the target oligonucleotidic sequence to the chromatographic beads so as to allow the irradiation-based release of the purified DNA–protein complexes off the beads. Our bioanalytical conditions were validated by purifying the tetracycline repressor protein onto a specific oligonucleotide. The purification factor was unprecedented, with a single contaminant. The robustness of our method was challenged by applying it to the purification of multiprotein assemblies forming onto DNA damage-mimicking oligonucleotides. The purified components were identified as well-known DNA repair proteins, and were shown to retain their enzymatic activities, as seen by monitoring DNA ligation products. Remarkably, kinase activities, also monitored, were found to be distinct on the beads and on the purified DNA–protein complexes, showing the benefits to uncouple the DNA–protein assemblies from the beads for a proper understanding of biochemical regulatory mechanisms involved in the DNA–protein assemblies.     

Analyzing the forces binding a restriction endonuclease to DNA using a synthetic nanopore

Restriction endonucleases are used prevalently in recombinant DNA technology because they bind so stably to a specific target sequence and, in the presence of cofactors, cleave double-helical DNA specifically at a target sequence at a high rate. Using synthetic nanopores along with molecular dynamics (MD), we have analyzed with atomic resolution how a prototypical restriction endonuclease, EcoRI, binds to the DNA target sequence—GAATTC—in the absence of a Mg²⁺ ion cofactor. We have previously shown that there is a voltage threshold for permeation of DNA bound to restriction enzymes through a nanopore that is associated with a nanonewton force required to rupture the complex. By introducing mutations in the DNA, we now show that this threshold depends on the recognition sequence and scales linearly with the dissociation energy, independent of the pore geometry. To predict the effect of mutation in a base pair on the free energy of dissociation, MD is used to qualitatively rank the stability of bonds in the EcoRI–DNA complex. We find that the second base in the target sequence exhibits the strongest binding to the protein, followed by the third and first bases, with even the flanking sequence affecting the binding, corroborating our experiments.     

Mutational and in vitro protein-binding studies on centromere DNA from Saccharomyces cerevisiae.

Centromeres on chromosomes in the yeast Saccharomyces cerevisiae contain approximately 140 base pairs (bp) of DNA. The functional centromere (CEN) region contains three important sequence elements (I, PuTCACPuTG; II, 78 to 86 bp of high-AT DNA; and III, a conserved 25-bp sequence with internal bilateral symmetry). Various point mutations or deletions in the element III region have a profound effect on CEN function in vivo, indicating that this DNA region is a key protein-binding site. This has been confirmed by the use of two in vitro assays to detect binding of yeast proteins to DNA fragments containing wild-type or mutationally altered CEN3 sequences. An exonuclease III protection assay was used to demonstrate specific binding of proteins to the element III region of CEN3. In addition, a gel DNA fragment mobility shift assay was used to characterize the binding reaction parameters. Sequence element III mutations that inactivate CEN function in vivo also prevent binding of proteins in the in vitro assays. The mobility shift assay indicates that double-stranded DNAs containing sequence element III efficiently bind proteins in the absence of sequence elements I and II, although the latter sequences are essential for optimal CEN function in vivo.     

Differences in hydration coupled to specific and nonspecific competitive binding and to specific DNA Binding of the restriction endonuclease BamHI

Using the osmotic stress technique together with a self-cleavage assay we measure directly differences in sequestered water between specific and nonspecific DNA-BamHI complexes as well as the numbers of water molecules released coupled to specific complex formation. The difference between specific and nonspecific binding free energy of the BamHI scales linearly with solute osmolal concentration for seven neutral solutes used to set water activity. The observed osmotic dependence indicates that the nonspecific DNA-BamHI complex sequesters some 120-150 more water molecules than the specific complex. The weak sensitivity of the difference in number of waters to the solute identity suggests that these waters are sterically inaccessible to solutes. This result is in close agreement with differences in the structures determined by x-ray crystallography. We demonstrate additionally that when the same solutes that were used in competition experiments are used to probe changes accompanying the binding of free BamHI to its specific DNA sequence, the measured number of water molecules released in the binding process is strikingly solute-dependent (with up to 10-fold difference between solutes). This result is expected for reactions resulting in a large change in a surface exposed area.