Changes in solvation during DNA binding and cleavage are critical to altered specificity of the EcoRI endonuclease

The biophysics of water in cell biology: perspectives on a keystone for both marine sciences and cancer research

The biophysics of water, has been debated over more than a century. Although its importance is still underestimated, significant breakthroughs occurred in recent years. The influence of protein condensation on water availability control was documented, new findings on water-transport proteins emerged, and the way water molecules rearrange to minimize free energy at interfaces was deciphered, influencing membrane thermodynamics. The state of knowledge continued to progress in the field of deep-sea marine biology, highlighting unknown effects of high hydrostatic pressure and/or temperature on interactions between proteins and ligands in extreme environments, and membrane structure adaptations. The role of osmolytes in protein stability control under stress is also discussed here in relation to fish egg hydration/buoyancy. The complexity of water movements within the cell is updated, all these findings leading to a better view of their impact on many cellular processes. The way water flow and osmotic gradients generated by ion transport work together to produce the driving force behind cell migration is also relevant to both marine biology and cancer research. Additional common points concern water dynamic changes during the neoplastic transformation of cells and tissues, or embryo development. This could improve imaging techniques, early cancer diagnosis, and understanding of the molecular and physiological basis of buoyancy for many marine species.

In-silicon studies on hydration in EcoRI-cognate DNA complex

Restriction endonucleases (REs) cleave DNA at specific site in presence of Mg2+ ion. Experiments further emphasize the role of hydration in metal ion specificity and sequence specificity of DNA cleavage. However, the relation between hydration and specificity has not been understood till date. This leads us to study via all-atom molecular dynamics (MD) simulations how the hydration around the scissile phosphate group changes in presence of Mg2+ and Ca2+ and depend on the DNA sequence. We observe the least number of hydrogen bonds around the scissile phosphate group in presence of Mg2+ ion. We further find that the hydrogen bonds decrease at the scissile phosphate on mutating one base pair in the cleavage region of the DNA in Mg2+ loaded EcoRI-DNA complex. We also perform steered MD simulations and observe that the rate of decrease of fraction of hydrogen bonds is slower in the mutated complex than the unmutated complex.

Effect of osmolytes on the EcoRI endonuclease: Insights into hydration and protein dynamics from molecular dynamics simulations

Osmolytes play an important role in cellular physiology by modulating the properties of proteins, including their molecular specificity. EcoRI is a model restriction enzyme whose specificity to DNA is altered in the presence of osmolytes. Here, we investigate the effect of two different osmolytes, glycerol and DMSO, on the dynamics and hydration of the EcoRI enzyme using molecular dynamics simulations. Our results show that the osmolytes, alter the essential dynamics of EcoRI. Particularly, we observe that the dynamics of the arm region of EcoRI which is involved in DNA binding is significantly altered. In addition, conformational free energy analyses reveals that the osmolytes bring about a change in the landscape similar to that of EcoRI bound to cognate DNA. We further observe that the hydration of the enzyme for each of the osmolyte is different, indicating that the mechanism of action of each of these osmolytes could be different. Further analyses of interfacial water dynamics using rotational autocorrelation function reveals that while the protein surface contributes to a slower tumbling motion of water, osmolytes, additionally contribute to the slowing of the angular motion of the water molecules. Entropy analysis also corroborates with this finding. We also find that the slowed rotational motion of interfacial waters in the presence of osmolytes contributes to a slowed relaxation of the hydrogen bonds between the interfacial waters and the functionally important residues in the protein. Taken together, our results show that osmolytes alter the dynamics of the protein by altering the dynamics of water. This altered dynamics, mediated by the changes in the water dynamics and hydrogen bonds with functionally important residues, may contribute to the altered specificity of EcoRI in the presence of osmolytes.

Relaxed Cleavage Specificity of Hyperactive Variants of Escherichia coli RNase E on RNA I

RNase E is an essential enzyme in Escherichia coli. The cleavage site of this single-stranded specific endoribonuclease is well-characterized in many RNA substrates. Here, we report that the upregulation of RNase E cleavage activity by a mutation that affects either RNA binding (Q36R) or enzyme multimerization (E429G) was accompanied by relaxed cleavage specificity. Both mutations led to enhanced RNase E cleavage in RNA I, an antisense RNA of ColE1-type plasmid replication, at a major site and other cryptic sites. Expression of a truncated RNA I with a major RNase E cleavage site deletion at the 5'-end (RNA I_-5) resulted in an approximately twofold increase in the steady-state levels of RNA I_-5 and the copy number of ColE1-type plasmid in E. coli cells expressing wild-type or variant RNase E compared to those expressing RNA I. These results indicate that RNA I_-5 does not efficiently function as an antisense RNA despite having a triphosphate group at the 5'-end, which protects the RNA from ribonuclease attack. Our study suggests that increased cleavage rates of RNase E lead to relaxed cleavage specificity on RNA I and the inability of the cleavage product of RNA I as an antisense regulator in vivo does not stem from its instability by having 5'-monophosphorylated end.

Effects of Hydrostatic Pressure on the Thermodynamics of CspB-Bs Interactions with the ssDNA Template

Understanding the thermodynamic mechanisms of adaptation of biomacromolecules to high hydrostatic pressure can help shed light on how piezophilic organisms can survive at pressures reaching over 1000 atmospheres. Interaction of proteins with nucleic acids is one of the central processes that allow information flow encoded in the sequence of DNA. Here, we report the results of a study on the interaction of cold shock protein B from Bacillus subtilis (CspB-Bs) with heptadeoxythymine template (pDT7) as a function of temperature and hydrostatic pressure. Experimental data collected at different CspB-Bs:pDT7 ratios were analyzed using a thermodynamic linkage model that accounts for both protein unfolding and CspB-Bs:pDT7 binding. The global fit to the model provided estimates of the stability of CspB-Bs, ΔG_Prot^o, the volume change upon CspB-Bs unfolding, ΔV_Prot, the association constant for CspB-Bs:pDT7 complex, K_a^o, and the volume changes upon pDT7 single-stranded DNA (ssDNA) template binding, ΔV_Bind. The protein, CspB-Bs, unfolds with an increase in hydrostatic pressure (ΔV_Prot < 0). Surprisingly, our study showed that ΔV_Bind < 0, which means that the binding of CspB-Bs to ssDNA is stabilized by an increase in hydrostatic pressure. Thus, CspB-Bs binding to pDT7 represents a case of linked equilibrium in which folding and binding react differently upon an increase in hydrostatic pressure: protein folding/unfolding equilibrium favors the unfolded state, while protein-ligand binding equilibrium favors the bound state. These opposing effects set a "maximum attainable" pressure tolerance to the protein-ssDNA complex under given conditions.

Probing DNA-protein interactions using single-molecule diffusivity contrast

Single-molecule fluorescence investigations of protein-nucleic acid interactions require robust means to identify the binding state of individual substrate molecules in real time. Here, we show that diffusivity contrast, widely used in fluorescence correlation spectroscopy at the ensemble level and in single-particle tracking on individual (but slowly diffusing) species, can be used as a general readout to determine the binding state of single DNA molecules with unlabeled proteins in solution. We first describe the technical basis of drift-free single-molecule diffusivity measurements in an anti-Brownian electrokinetic trap. We then cross-validate our method with protein-induced fluorescence enhancement, a popular technique to detect protein binding on nucleic acid substrates with single-molecule sensitivity. We extend an existing hydrodynamic modeling framework to link measured diffusivity to particular DNA-protein structures and obtain good agreement between the measured and predicted diffusivity values. Finally, we show that combining diffusivity contrast with protein-induced fluorescence enhancement allows simultaneous mapping of binding stoichiometry and location on individual DNA-protein complexes, potentially enhancing single-molecule views of relevant biophysical processes.

Structural and dynamical heterogeneity of water trapped inside Na+-pumping KR2 rhodopsin in the dark state

Photoisomerization in the retinal leads to a channel opening in rhodopsins that triggers translocation or pumping of ions/protons. Crystal structures of rhodopsins contain several structurally conserved water molecules. It has been suggested that water plays an active role in facilitating the ion pumping/translocation process by acting as a lubricant in these systems. In this paper, we systematically investigate the localization, structure, dynamics, and energetics of the water molecules along the channel for the resting/dark state of KR2 rhodopsin. By employing several microseconds long atomistic molecular dynamics simulation of this trans-membrane protein system, we demonstrate the presence of five distinct water containing pockets/cavities separated by gateways controlled by protein side-chains. There exists a strong hydrogen bonded network involving these buried water molecules and functionally important key residues. We present evidence of significant structural and dynamical heterogeneity in the water molecules present in these cavities, with very rare exchange between them. The exchange time scale of such buried water with the bulk has an extremely wide range, from tens of nanoseconds to >1.5 µs. The translational and rotational dynamics of buried water are found to be strongly dependent on the protein cavity size and local interactions with a classic signature of trapped diffusion and rotational anisotropy.

Mapping interfacial hydration in ETS-family transcription factor complexes with DNA: a chimeric approach

Hydration of interfaces is a major determinant of target specificity in protein/DNA interactions. Interfacial hydration is a highly variable feature in DNA recognition by ETS transcription factors and functionally relates to cellular responses to osmotic stress. To understand how hydration is mediated in the conserved ETS/DNA binding interface, secondary structures comprising the DNA contact surface of the strongly hydrated ETS member PU.1 were substituted, one at a time, with corresponding elements from its sparsely hydrated relative Ets-1. The resultant PU.1/Ets-1 chimeras exhibited variably reduced sensitivity to osmotic pressure, indicative of a distributed pattern of interfacial hydration in wildt-ype PU.1. With the exception of the recognition helix H3, the chimeras retained substantially high affinities. Ets-1 residues could therefore offset the loss of favorable hydration contributions in PU.1 via low-water interactions, but at the cost of decreased selectivity at base positions flanking the 5'-GGA-3' core consensus. Substitutions within H3 alone, which contacts the core consensus, impaired binding affinity and PU.1 transactivation in accordance with the evolutionary separation of the chimeric residues involved. The combined biophysical, bioinformatics and functional data therefore supports hydration as an evolved specificity determinant that endows PU.1 with more stringent sequence selection over its ancestral relative Ets-1.

Temperature and osmotic stress dependence of the thermodynamics for binding linker histone H10, Its carboxyl domain (H10-C) or globular domain (H10-G) to B-DNA

Linker histones (H1) are the basic proteins in higher eukaryotes that are responsible for the final condensation of chromatin. In contrast to the nucleosome core histone proteins, the role of H1 in compacting DNA is not clearly understood. In this study ITC was used to measure the binding constant, enthalpy change, and binding site size for the interactions of H10, or its C-terminal (H10-C) and globular (H10-G) domains to highly polymerized calf-thymus DNA at temperatures from 288 K to 308 K. Heat capacity changes, ΔCp, for these same H10 binding interactions were estimated from the temperature dependence of the enthalpy changes. The enthalpy changes for binding H10, H10-C, or H10-G to CT-DNA are all endothermic at 298 K, becoming more exothermic as the temperature is increased. The ΔH for binding H10-G to CT-DNA is exothermic at temperatures above approximately 300 K. Osmotic stress experiments indicate that the binding of H10 is accompanied by the release of approximately 35 water molecules. We estimate from our naked DNA titration results that the binding of the H10 to the nucleosome places the H10 protein in close contact with approximately 41 DNA bp. The breakdown is that the H10 carboxyl terminus interacts with 28 bp of linker DNA on one side of the nucleosome, the H10 globular domain binds directly to 7 bp of core DNA, and shields another 6 linker DNA bases, 3 bp on either side of the nucleosome where the linker DNA exits the nucleosome core.

Mechanisms of Protein Translocation on DNA Are Differentially Responsive to Water Activity

Water plays important but poorly understood roles in the functions of most biomolecules. We are interested in understanding how proteins use diverse search mechanisms to locate specific sites on DNA; here we present a study of the role of closely associated waters in diverse translocation mechanisms. The bacterial DNA adenine methyltransferase, Dam, moves across large segments of DNA using an intersegmental hopping mechanism, relying in part on movement through bulk water. In contrast, other proteins, such as the bacterial restriction endonuclease EcoRI, rely on a sliding mechanism, requiring the protein to stay closely associated with DNA. Here we probed how these two mechanistically distinct proteins respond to well-characterized osmolytes, dimethyl sulfoxide (DMSO), and glycerol. The ability of Dam to move over large segments of DNA is not impacted by either osmolyte, consistent with its minimal reliance on a sliding mechanism. In contrast, EcoRI endonuclease translocation is significantly enhanced by DMSO and inhibited by glycerol, providing further corroboration that these proteins rely on distinct translocation mechanisms. The well-established similar effects of these osmolytes on bulk water, and their differential effects on macromolecule-associated waters, support our results and provide further evidence of the importance of water in interactions between macromolecules and their ligands.

Relevance of the Drag Force during Controlled Translocation of a DNA-Protein Complex through a Glass Nanocapillary

Combination of glass nanocapillaries with optical tweezers allowed us to detect DNA-protein complexes in physiological conditions. In this system, a protein bound to DNA is characterized by a simultaneous change of the force and ionic current signals from the level observed for the bare DNA. Controlled displacement of the protein away from the nanocapillary opening revealed decay in the values of the force and ionic current. Negatively charged proteins EcoRI, RecA, and RNA polymerase formed complexes with DNA that experienced electrophoretic force lower than the bare DNA inside nanocapillaries. Force profiles obtained for DNA-RecA in our system were different than those in the system with nanopores in membranes and optical tweezers. We suggest that such behavior is due to the dominant impact of the drag force comparing to the electrostatic force acting on a DNA-protein complex inside nanocapillaries. We explained our results using a stochastic model taking into account the conical shape of glass nanocapillaries.

Relaxed specificity of prokaryotic DNA methyltransferases results in DNA site-specific modification of RNA/DNA heteroduplexes

RNA/DNA hybrid duplexes regularly occur in nature, for example in transcriptional R loops. Their susceptibility to modification by DNA-specific or RNA-specific enzymes is, thus, a biologically relevant question, which, in addition, has possible biotechnological implications. In this study, we investigated the activity of four isospecific DNA methyltransferases (M.EcoVIII, M.LlaCI, M.HindIII, M.BstZ1II) toward an RNA/DNA duplex carrying one 5'-AAGCUU-3'/3'-TTCGAA-5' target sequence. The analyzed enzymes belong to the β-group of adenine N6-methyltransferases and recognize the palindromic DNA sequence 5'-AAGCTT-3'/3'-TTCGAA-5'. Under standard conditions, none of these isospecific enzymes could detectibly methylate the RNA/DNA duplex. However, the addition of agents that generally relax specificity, such as dimethyl sulfoxide (DMSO) and glycerol, resulted in substantial methylation of the RNA/DNA duplex by M.EcoVIII and M.LlaCI. Only the DNA strand of the RNA/DNA duplex was methylated. The same was not observed for M.HindIII or M.BstZ1II. This is, to our knowledge, the first report that demonstrates such activity by prokaryotic DNA methyltransferases. Possible applications of these findings in a laboratory practice are also discussed.

Enzyme-free and label-free signal amplification for monitoring endonuclease activity and inhibition via hybridization chain reaction

A facile and sensitive method for studying EcoRI activity and inhibition had been developed based on EcoRI-triggered hybridization chain reaction.

Using single-turnover kinetics with osmotic stress to characterize the EcoRV cleavage reaction

Type II restriction endonucleases require metal ions to specifically cleave DNA at canonical sites. Despite the wealth of structural and biochemical information, the number of Mg(2+) ions used for cleavage by EcoRV, in particular, at physiological divalent ion concentrations has not been established. In this work, we employ a single-turnover technique that uses osmotic stress to probe reaction kinetics between an initial specific EcoRV-DNA complex formed in the absence of Mg(2+) and the final cleavage step. With osmotic stress, complex dissociation before cleavage is minimized and the reaction rates are slowed to a convenient time scale of minutes to hours. We find that cleavage occurs by a two-step mechanism that can be characterized by two rate constants. The dependence of these rate constants on Mg(2+) concentration and osmotic pressure gives the number of Mg(2+) ions and water molecules coupled to each kinetic step of the EcoRV cleavage reaction. Each kinetic step is coupled to the binding 1.5-2.5 Mg(2+) ions, the uptake of ∼30 water molecules, and the cleavage of a DNA single strand. We suggest that each kinetic step reflects an independent, rate-limiting conformational change of each monomer of the dimeric enzyme that allows Mg(2+) ion binding. This modified single-turnover protocol has general applicability for metalloenzymes.

Dynamics of water molecules and sodium ions in solid hydrates of nucleotides

The dynamics of the co-ordinating water and metal cations in solid hydrates of nucleotide salts is explored with solid-state NMR spectroscopy and DFT calculations.

Effects of 2'-O-methyl nucleotide substitution on EcoRI endonuclease cleavage activities

To investigate the effect of sugar pucker conformation on DNA-protein interactions, we used 2′-O-methyl nucleotide (2′-OMeN) to modify the EcoRI recognition sequence -TGAATTCT-, and monitored the enzymatic cleavage process using FRET method. The 2′-O-methyl nucleotide has a C3′-endo sugar pucker conformation different from the C2′-endo sugar pucker conformation of native DNA nucleotides. The initial reaction velocities were measured and the kinetic parameters, K_m and V_max were derived using Michaelis-Menten equation. Experimental results showed that 2′-OMeN substitutions for the EcoRI recognition sequence decreased the cleavage efficiency for A2, A3 and T4 substitutions significantly, and 2′-OMeN substitution for T5 residue inhibited the enzymatic activity completely. In contrast, substitutions for G1 and C6 could maintain the original activity. 2′-fluoro nucleic acid (2′-FNA) and locked nucleic acid (LNA) having similar C3′-endo sugar pucker conformation also demonstrated similar enzymatic results. This position-dependent enzymatic cleavage property might be attributed to the phosphate backbone distortion caused by the switch from C2′-endo to C3′-endo sugar pucker conformation, and was interpreted on the basis of the DNA-EcoRI structure. These 2′-modified nucleotides could behave as a regulatory element to modulate the enzymatic activity in vitro, and this property will have potential applications in genetic engineering and biomedicine.

Quantifying the molecular origins of opposite solvent effects on protein-protein interactions

Although the nature of solvent-protein interactions is generally weak and non-specific, addition of cosolvents such as denaturants and osmolytes strengthens protein-protein interactions for some proteins, whereas it weakens protein-protein interactions for others. This is exemplified by the puzzling observation that addition of glycerol oppositely affects the association constants of two antibodies, D1.3 and D44.1, with lysozyme. To resolve this conundrum, we develop a methodology based on the thermodynamic principles of preferential interaction theory and the quantitative characterization of local protein solvation from molecular dynamics simulations. We find that changes of preferential solvent interactions at the protein-protein interface quantitatively account for the opposite effects of glycerol on the antibody-antigen association constants. Detailed characterization of local protein solvation in the free and associated protein states reveals how opposite solvent effects on protein-protein interactions depend on the extent of dewetting of the protein-protein contact region and on structural changes that alter cooperative solvent-protein interactions at the periphery of the protein-protein interface. These results demonstrate the direct relationship between macroscopic solvent effects on protein-protein interactions and atom-scale solvent-protein interactions, and establish a general methodology for predicting and understanding solvent effects on protein-protein interactions in diverse biological environments.

Luminescent and colorimetric strategies for the label-free DNA-based detection of enzyme activity

Enzymes are critically involved in maintaining normal cellular physiology through the catalysis of highly specific and tightly regulated chemical reactions. The inhibition or undesired activation of particular enzymatic functions has been associated with the pathogenesis of a number of diseases. Consequently, the aberrant activity of certain enzymes can be regarded as biomarkers for the diagnosis or monitoring of particular diseases. With rapid technological advances in the field of DNA nanotechnology, oligonucleotides have recently emerged as attractive recognition units for monitoring the activity of enzymes compared with organic small molecules or protein antibodies. In this review article, we present an overview of advantages and versatility of the 'label-free' approach for the fabrication of DNA-based sensing platforms using colorimetric and luminescent molecules as signal transducing units and highlight recent examples of label-free strategies that have been employed for monitoring enzyme activity.

Label-free fluorescent assay for real-time monitoring site-specific DNA cleavage by EcoRI endonuclease

DNA cleavage reaction catalyzed by nucleases is essential in many important biological processes and medicinal chemistry. Therefore, it is important to develop reliable and facile methods to assay nuclease activity. With this goal in mind, we report a fluorescent assay for label-free, facile, and real-time monitoring of DNA cleavage by EcoRI endonuclease using SYBR Green I (SGI) as a signal probe. The fluorescence of SGI dramatically increased when the free SGI was mixed with double-stranded DNA (dsDNA) substrate. Upon interacting with EcoRI, which cleaves the dsDNA into small fragments, the weakened interaction between SGI and the shortened DNA fragments caused a decrease in fluorescence of SGI. EcoRI-DNA interaction was real-time studied by monitoring fluorescence change with the prolonging of interaction time. The important kinetic parameters, including Michaelis-Menten constant (K(M)) and maximum initial velocity (V(max)), were accurately calculated, which is consistent with previously reported studies. Site-specific DNA cleavage by EcoRI endonuclease has also been verified by gel electrophoresis analysis, which indicated that this method is a simple and effective approach to assay DNA cleavage reaction. Specificity investigation demonstrated that EcoRI-DNA interactions can be studied with high selectivity. Compared with previously reported methods, this approach is selective, simple, convenient and cost-efficient without any labeling of the probe or of the target.

Correlated structural kinetics and retarded solvent dynamics at the metalloprotease active site

Solvent dynamics can play a major role in enzyme activity, but obtaining an accurate, quantitative picture of solvent activity during catalysis is quite challenging. Here, we combine terahertz spectroscopy and X-ray absorption analyses to measure changes in the coupled water-protein motions during peptide hydrolysis by a zinc-dependent human metalloprotease. These changes were tightly correlated with rearrangements at the active site during the formation of productive enzyme-substrate intermediates and were different from those in an enzyme-inhibitor complex. Molecular dynamics simulations showed a steep gradient of fast-to-slow coupled protein-water motions around the protein, active site and substrate. Our results show that water retardation occurs before formation of the functional Michaelis complex. We propose that the observed gradient of coupled protein-water motions may assist enzyme-substrate interactions through water-polarizing mechanisms that are remotely mediated by the catalytic metal ion and the enzyme active site.

Modeling of the water network at protein-RNA interfaces

Water plays an important role in the mediation of biomolecular interactions. Thus, accurate prediction and evaluation of water-mediated interactions is an important element in the computational design of interfaces involving proteins, RNA, and DNA. Here, we use an algorithm (WATGEN) to predict the locations of interfacial water molecules for a data set of 224 protein-RNA interfaces. The accuracy of the prediction is validated against water molecules present in the X-ray structures of 105 of these complexes. The complexity of the water networks is deconvoluted through definition of the characteristics of each water molecule based on its bridging properties between the protein and RNA and on its depth in the interface with respect to the bulk solvent. This approach has the potential for scoring the water network for incorporation into the computational design of protein-RNA complexes.

Stabilizing labile DNA-protein complexes in polyacrylamide gels

The electrophoretic mobility-shift assay (EMSA) is one of the most popular tools in molecular biology for measuring DNA-protein interactions. EMSA, as standardly practiced today, works well for complexes with association binding constants K(a)>10(9) M(-1) under normal conditions of salt and pH. Many DNA-protein complexes are not stable enough so that they dissociate while moving through the gel matrix giving smeared bands that are difficult to quantitate reliably. In this work we demonstrate that the addition of the osmolyte triethylene glycol to polyacrylamide gels dramatically stabilizes labile restriction endonuclease EcoRI complexes with nonspecific DNA sequences enabling quantitation of binding using EMSA. The significant improvement of the technique resulting from the addition of osmolytes to the gel matrix greatly extends the range of binding constants of protein-DNA complexes that can be investigated using this widely used assay. Extension of this approach to other techniques used for separating bound and free components such as gel chromatography and CE is straightforward.

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(2+) 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.

Hydration of Watson-Crick base pairs and dehydration of Hoogsteen base pairs inducing structural polymorphism under molecular crowding conditions

It has been revealed recently that molecular crowding, which is one of the largest differences between in vivo and in vitro conditions, is a critical factor determining the structure, stability, and function of nucleic acids. However, the effects of molecular crowding on Watson-Crick and Hoogsteen base pairs remain unclear. In order to investigate directly and quantitatively the molecular crowding effects on base pair types in nucleic acids, we designed intramolecular parallel- and antiparallel-stranded DNA duplexes consisting of Hoogsteen and Watson-Crick base pairs, respectively, as well as an intramolecular parallel-stranded triplex containing both types of base pairs. Thermodynamic analyses demonstrated that the values of free energy change at 25 degrees C for Hoogsteen base-pair formations decreased from +1.45 +/- 0.15 to +1.09 +/- 0.13 kcal mol(-1), and from -1.89 +/- 0.13 to -2.71 +/- 0.11 kcal mol(-1) in the intramolecular duplex and triplex, respectively, when the concentration of PEG 200 (polyethylene glycol with average molecular weight 200) increased from 0 to 20 wt %. However, corresponding values for Watson-Crick formation in the duplex and triplex increased from -10.2 +/- 0.2 to -8.7 +/- 0.1 kcal mol(-1), and from -10.8 +/- 0.2 to -9.2 +/- 0.2 kcal mol(-1), respectively. Furthermore, it was revealed that the opposing effects of molecular crowding on the Hoogsteen and Watson-Crick base pairs were due to different behaviors of water molecules binding to the DNA strands.

Sequence-dependence and prediction of nucleotide solvent accessibility in double stranded DNA

Solvent accessibility of amino acid residues in proteins has been widely studied and many methods for its prediction from sequence and evolutionary information are available. Some of the advantages of studying amino acid solvent accessibility also apply to DNA. However, currently there are no methods to estimate the solvent accessibility of nucleotides, as most works on DNA structures have focused on elastic deformations and other structural attributes. In this work, an attempt has been made to analyze the distribution of different nucleotides in various accessibility ranges. Effect of neighboring nucleotides on the predictability of exposure has been evaluated by developing a linear perceptron model that takes sequence information as the input. Five different types of solvent accessibility (overall nucleotide, side chain, main chain, polar and non-polar) have been predicted. From the analysis, it is observed that Thymine stands out in terms of its higher exposed surface area, particularly its side chain and non-polar atoms. It is also concluded that the solvent accessibility of a nucleotide strongly depends on its sequence neighbors and can be predicted with fair success using this information.

Specific versus nonspecific isothermal DNA amplification through thermophilic polymerase and nicking enzyme activities

Rapid isothermal nucleic acid amplification technologies can enable diagnosis of human pathogens and genetic variations in a simple, inexpensive, user-friendly format. The isothermal exponential amplification reaction (EXPAR) efficiently amplifies short oligonucleotides called triggers in less than 10 min by means of thermostable polymerase and nicking endonuclease activities. We recently demonstrated that this reaction can be coupled with upstream generation of trigger oligonucleotides from a genomic target sequence, and with downstream visual detection using DNA-functionalized gold nanospheres. The utility of EXPAR in clinical diagnostics is, however, limited by a nonspecific background amplification phenomenon, which is further investigated in this report. We found that nonspecific background amplification includes an early phase and a late phase. Observations related to late phase background amplification are in general agreement with literature reports of ab initio DNA synthesis. Early phase background amplification, which limits the sensitivity of EXPAR, differs however from previous reports of nonspecific DNA synthesis. It is observable in the presence of single-stranded oligonucleotides following the EXPAR template design rules and generates the trigger sequence expected for the EXPAR template present in the reaction. It appears to require interaction between the DNA polymerase and the single-stranded EXPAR template. Early phase background amplification can be suppressed or eliminated by physically separating the template and polymerase until the final reaction temperature has been reached, thereby enabling detection of attomolar starting trigger concentrations.

WATGEN: an algorithm for modeling water networks at protein-protein interfaces

Water molecules at protein-protein interfaces contribute to the close packing of atoms and ensure complementarity between the protein surfaces, as well as mediating polar interactions. Therefore, modeling of interface water is of importance in understanding the structural basis of biomolecular association. We present an algorithm, WATGEN, which predicts locations for water molecules at a protein-protein or protein-peptide interface, given the atomic coordinates of the protein and peptide. A key element of the WATGEN algorithm is the prediction of water sites that can form multiple hydrogen bonds that bridge the binding interface. Trial calculations were performed on water networks predicted by WATGEN at 126 protein-peptide interfaces (X-ray resolutions <or= 2.0 A), using different criteria for water placement. The energies of the predicted water networks were evaluated in AMBER8 and used in the choice of parameters for WATGEN. The 126 interfaces include 1264 experimentally determined bridging water sites, and the WATGEN algorithm predicts 72 and 88% of these sites within 1.5 and 2.0 A, respectively. The predicted number of water molecules at each interface was much higher than the number of water molecules identified experimentally. Therefore, random placement of the same number of water molecules as that predicted at each interface was performed as a control, and resulted in only 22 and 40% of water sites placed within 1.5 and 2.0 A of experimental sites, respectively. Based on these data, we conclude that WATGEN can accurately predict the location of water molecules at a protein-peptide interface, and this may be of value for understanding the energetics and specificity of biomolecular association.

Electrophoretic mobility shift assay (EMSA) for detecting protein-nucleic acid interactions

The gel electrophoresis mobility shift assay (EMSA) is used to detect protein complexes with nucleic acids. It is the core technology underlying a wide range of qualitative and quantitative analyses for the characterization of interacting systems. In the classical assay, solutions of protein and nucleic acid are combined and the resulting mixtures are subjected to electrophoresis under native conditions through polyacrylamide or agarose gel. After electrophoresis, the distribution of species containing nucleic acid is determined, usually by autoradiography of 32P-labeled nucleic acid. In general, protein-nucleic acid complexes migrate more slowly than the corresponding free nucleic acid. In this protocol, we identify the most important factors that determine the stabilities and electrophoretic mobilities of complexes under assay conditions. A representative protocol is provided and commonly used variants are discussed. Expected outcomes are briefly described. References to extensions of the method and a troubleshooting guide are provided.

Effects of ethylene glycol on the torsion elastic constant and hydrodynamic radius of p30δ DNA

Upon increasing the concentration of ethylene glycol (EG) at 37 degrees C, the twist energy parameter, E(T), which governs the supercoiling free energy, was recently found to undergo a decreasing (or reverse) sigmoidal transition with a midpoint near 20 w/v % EG. In this study, the effects of adding 20 w/v % EG on the torsion elastic constant (alpha) of linear p30delta DNA and on the hydrodynamic radius (R(H)) of a synthetic 24 bp duplex DNA were examined at both 40 and 20 degrees C. The time-resolved fluorescence intensity and fluorescence polarization anisotropy (FPA) of intercalated ethidium were measured in order to assess the effects of 20 w/v % EG on: (1) alpha; (2) R(H); (3) the lifetimes of intercalated and non-intercalated dye; (4) the amplitude of dye wobble in its binding site; and (5) the binding constant for intercalation. The effects of 20 w/v % EG on the circular dichroism (CD) spectrum of the DNA and on the emission spectrum of the free dye were also measured. At 40 degrees C, addition of 20 w/v % EG caused a substantial (1.27- to 1.35-fold) increase in alpha, a significant change in the CD spectrum, and a very small, marginally significant increase in R(H), but little or no change in the amplitude of dye wobble in its binding site or the lifetime of intercalated dye. Together with previously reported measurements of E(T), these results imply that the bending elastic constant of DNA is significantly decreased by 20 w/v % EG at 40 degrees C. At 20 degrees C, addition of 20 w/v % EG caused a marginally significant decrease in alpha and very little change in any other measured properties. Also at 20 degrees C, addition of 30 w/v % betaine caused a marginally significant increase in alpha and significant but modest change in the CD spectrum, but very little change in any other properties.

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.

Sequence-dependent enhancement of hydrolytic deamination of cytosines in DNA by the restriction enzyme PspGI

Hydrolytic deamination of cytosines in DNA creates uracil and, if unrepaired, these lesions result in C to T mutations. We have suggested previously that a possible way in which cells may prevent or reduce this chemical reaction is through the binding of proteins to DNA. We use a genetic reversion assay to show that a restriction enzyme, PspGI, protects cytosines within its cognate site, 5'-CCWGG (W is A or T), against deamination under conditions where no DNA cleavage can occur. It decreases the rate of cytosine deamination to uracil by 7-fold. However, the same protein dramatically increases the rate of deaminations within the site 5'-CCSGG (S is C or G) by approximately 15-fold. Furthermore, a similar increase in cytosine deaminations is also seen with a catalytically inactive mutant of the enzyme showing that endonucleolytic ability of the protein is dispensable for its mutagenic action. The sequences of the mutants generated in the presence of PspGI show that only one of the cytosines in CCSGG is predominantly converted to thymine. Our results are consistent with PspGI 'sensitizing' the cytosine in the central base pair in CCSGG for deamination. Remarkably, PspGI sensitizes this base for damage despite its inability to form stable complexes at CCSGG sites. These results can be explained if the enzyme has a transient interaction with this sequence during which it flips the central cytosine out of the helix. This prediction was validated by modeling the structure of PspGI-DNA complex based on the structure of the related enzyme Ecl18kI which is known to cause base-flipping.

Sequestered water and binding energy are coupled in complexes of lambda Cro repressor with non-consensus binding sequences

We use the osmotic pressure dependence of dissociation rates and relative binding constants to infer differences in sequestered water among complexes of lambda Cro repressor with varied DNA recognition sequences. For over a 1000-fold change in association constant, the number of water molecules sequestered by non-cognate complexes varies linearly with binding free energy. One extra bound water molecule is coupled with the loss of approximately 150 cal/mol complex in binding free energy. Equivalently, every tenfold decrease in binding constant at constant salt and temperature is associated with eight to nine additional water molecules sequestered in the non-cognate complex. The relative insensitivity of the difference in water molecules to the nature of the osmolyte used to probe the reaction suggests that the water is sterically sequestered. If the previously measured changes in heat capacity for lambda Cro binding to different non-cognate sequences are attributed solely to this change in water, then the heat capacity change per incorporated water is almost the same as the difference between ice and water. The associated changes in enthalpies and entropies, however, indicate that the change in complex structure involves more than a simple incorporation of fixed water molecules that act as adaptors between non-complementary surfaces.

Solvent-induced changes in photochemical activity and conformation of photosystem I particles by glycerol

It has been shown that a large number of water molecules coordinate with the pigments and subunits of photosystem I (PSI); however, the function of these water molecules remains to be clarified. In this study, the photosynthetic properties of PSI from spinach were investigated using different spectroscopic and activity measurements under conditions of decreasing water content caused by increasing concentrations of glycerol. The results show that glycerol addition caused pronounced changes in the photochemical activity of PSI particles. At low concentrations (<60%, v/v), glycerol stimulated the rate of oxygen uptake in PSI particles, while higher concentrations of glycerol cause inhibition of PSI activity. The capacity of P700 photooxidation also increased with glycerol concentrations lower than 60%. In contrast, this capacity decreased at higher glycerol concentrations. On the other hand, glycerol addition considerably affected the distribution of the bulk and red antenna chlorophyll (Chl) forms or states, with the population of red-shifted Chl forms augmented with increasing glycerol. In addition, glycerol-treated PSI particles showed a blue shift of the tryptophan fluorescence emission maximum and an increase in their capacity to bind the hydrophobic probe 1-anilino-8-naphthalene sulfonate, indicating a more non-polar environment for tryptophan residues and increased exposure of hydrophobic surfaces.

Trapping DNA-protein binding reactions with neutral osmolytes for the analysis by gel mobility shift and self-cleavage assays

We take advantage of our previous observation that neutral osmolytes can strongly slow down the rate of DNA–protein complex dissociation to develop a method that uses osmotic stress to ‘freeze’ mixtures of DNA–protein complexes and prevent further reaction enabling analysis of the products. We apply this approach to the gel mobility shift assay and use it to modify a self-cleavage assay that uses the nuclease activity of the restriction endonucleases to measure sensitively their specific binding to DNA. At sufficiently high concentrations of neutral osmolytes the cleavage reaction can be triggered at only those DNA fragments with initially bound enzyme. The self-cleavage assay allows measurement of binding equilibrium and kinetics directly in solution avoiding the intrinsic problems of gel mobility shift and filter binding assays while providing the same sensitivity level. Here we compare the self-cleavage and gel mobility shift assays applied to the DNA binding of EcoRI and BamHI restriction endonucleases. Initial results indicate that BamHI dissociation from its specific DNA sequence is strongly linked to water activity with the half-life time of the specific complex increasing ∼20-fold from 0 to 1 osmolal betaine.

Interfacial water as a "hydration fingerprint" in the noncognate complex of BamHI

The molecular code of specific DNA recognition by proteins as a paradigm in molecular biology remains an unsolved puzzle primarily because of the subtle interplay between direct protein-DNA interaction and the indirect contribution from water and ions. Transformation of the nonspecific, low affinity complex to a specific, high affinity complex is accompanied by the release of interfacial water molecules. To provide insight into the conversion from the loose to the tight form, we characterized the structure and energetics of water at the protein-DNA interface of the BamHI complex with a noncognate sequence and in the specific complex. The fully hydrated models were produced with Grand Canonical Monte Carlo simulations. Proximity analysis shows that water distributions exhibit sequence dependent variations in both complexes and, in particular, in the noncognate complex they discriminate between the correct and the star site. Variations in water distributions control the number of water molecules released from a given sequence upon transformation from the loose to the tight complex as well as the local entropy contribution to the binding free energy. We propose that interfacial waters can serve as a "hydration fingerprint" of a given DNA sequence.

Factors influencing resistance of UV-irradiated DNA to the restriction endonuclease cleavage

DNA molecules of pUC19, pBR322 and PhiX174 were irradiated by various doses of UV light and the irradiated molecules were cleaved by about two dozen type II restrictases. The irradiation generally blocked the cleavage in a dose-dependent way. In accordance with previous studies, the (A + T)-richness and the (PyPy) dimer content of the restriction site belongs among the factors that on average, cause an increase in the resistance of UV damaged DNA to the restrictase cleavage. However, we observed strong effects of UV irradiation even with (G + C)-rich and (PyPy)-poor sites. In addition, sequences flanking the restriction site influenced the protection in some cases (e.g. HindIII), but not in others (e.g. SalI), whereas neoschizomer couples SmaI and AvaI, or SacI and Ecl136II, cleaved the UV-irradiated DNA similarly. Hence the intrastrand thymine dimers located in the recognition site are not the only photoproduct blocking the restrictases. UV irradiation of the A-form generally made the irradiated DNA less resistant to restrictase cleavage than irradiation in the B-form and in some cases, the A-form completely protected the UV-irradiated DNA against the damage recognized by the restrictases. The present results also demonstrate that the UV irradiation approach used to generate partial digests in genomic DNA studies, can be extended to the (G + C)-rich and (PyPy)-poor restriction sites. The present extensive and quantitative data can be used in genomic applications of UV damage probing by restrictases.