DNA-bending on ligand binding; change of DNA persistence length.

This paper simulates the helix-characteristic changes of apparent DNA persistence length caused by randomly distributed helix bends as induced, e.g., by DNA-bound ligand molecules. The parameters varied are the constant angle gamma of helix bending and the size alpha of the DNA drug binding site, but also the degree of DNA-ligand binding cooperativity and the helix-unwinding angle. If the size of the binding site is comparable with the helix pitch, the influence of phasing between helix bends and helix screw upon the apparent persistence length is obvious. In the accompanying paper experimental data are analyzed in terms of this theoretical background.     

Human flap endonuclease-1: conformational change upon binding to the flap DNA substrate and location of the Mg2+ binding site

Human flap endonuclease-1 (FEN-1) is a member of the structure-specific endonuclease family and is a key enzyme in DNA replication and repair. FEN-1 recognizes the 5'-flap DNA structure and cleaves it, a specialized endonuclease function essential for the processing of Okazaki fragments during DNA replication and for the repair of 5'-end single-stranded tails from nicked double-stranded DNA substrates. Magnesium is a cofactor required for nuclease activity. We have used Fourier transform infrared (FTIR) spectroscopy to better understand how Mg2+ and flap DNA interact with human FEN-1. FTIR spectroscopy provides three fundamentally new insights into the structural changes induced by the interaction of FEN-1 with substrate DNA and Mg2+. First, FTIR difference spectra in the amide I vibrational band (1600-1700 cm(-1)) reveal a change in the secondary structure of FEN-1 induced by substrate DNA binding. Quantitative analysis of the FTIR spectra indicates a 4% increase in helicity upon DNA binding or about 14 residues converted from disordered to helical conformations. The observation that the residues are disordered without DNA strongly implicates the flexible loop region. The conversion to helix also suggests a mechanism for locking the flexible loop region around the bound DNA. This is the first direct experimental evidence for a binding mechanism that involves a secondary structural change of the protein. Second, in contrast with DNA binding, no change is observed in the secondary structure of FEN-1 upon Mg2+ binding to the wild type or to the noncleaving D181A mutant. Third, the FTIR results provide direct evidence (via the carboxylate ligand band at 1535 cm(-1)) that not only is D181 a ligand to Mg2+ in the human enzyme but Mg2+ binding does not occur in the D181A mutant which lacks this ligand.     

Novel allosteric conformation of human HB revealed by the hydration and anion effects on O(2) binding.

The effect of anions on the stability of different functional conformations of Hb is examined through the determination of the dependence of O(2) affinity on water activity (a(w)). The control of a(w) is effected by varying the sucrose osmolal concentration in the bathing solution according to the "osmotic stress" method. Thus, the hydration change following Hb oxygenation is determined as a function of Cl(-) and of DPG concentration. We find that only approximately 25 additional water molecules bind to human Hb during the deoxy-to-oxy conformation transition in the absence of anions, in contrast with approximately 72 that bind in the presence of more than 50 mM Cl(-) or more than 15 microM DPG. We demonstrate that the increase in the hydration change linked with oxygenation is coupled with anion binding to the deoxy-Hb. Hence, we propose that the deoxy-Hb coexists in two allosteric conformations which depend on whether anion is bound or not: the tense T-state, with low oxygen affinity and anion bound, or a new allosteric P-state, with intermediate oxygen affinity and free of bound anions. The intrinsic oxygen affinity of this unforeseen P-state and the differential binding of Cl(-), DPG, and H(2)O between states P and T and P and R are characteristics which are consistent with those expected for a putative intermediate allosteric state of Hb. These findings represent a new opportunity to explore the structure-function relationships of hemoglobin regulation.     

The role of hydration on the mechanism of allosteric regulation: in situ measurements of the oxygen-linked kinetics of water binding to hemoglobin

We report here the first direct measurements of changes in protein hydration triggered by a functional binding. This task is achieved by weighing hemoglobin (Hb) and myoglobin films exposed to an atmosphere of 98% relative humidity during oxygenation. The binding of the first oxygen molecules to Hb tetramer triggers a change in protein conformation, which increases binding affinity to the remaining empty sites giving rise to the appearance of cooperative phenomena. Although crystallographic data have evidenced that this structural change increases the protein water-accessible surface area, isobaric osmotic stress experiments in aqueous cosolutions have shown that water binding is linked to Hb oxygenation. Now we show that the differential hydration between fully oxygenated and fully deoxygenated states of these proteins, determined by weighing protein films with a quartz crystal microbalance, agree with the ones determined by osmotic stress in aqueous cosolutions, from the linkage between protein oxygen affinity and water activity. The agreements prove that the changes in water activity brought about by adding osmolytes to the buffer solution shift biochemical equilibrium in proportion to the number of water molecules associated with the reaction. The concomitant kinetics of oxygen and of water binding to Hb have been also determined. The data show that the binding of water molecules to the extra protein surface exposed on the transition from the low-affinity T to the high-affinity R conformations of hemoglobin is the rate-limiting step of Hb cooperative reaction. This evidences that water binding is a crucial step on the allosteric mechanism regulating cooperative interactions, and suggests the possibility that environmental water activity might be engaged in the kinetic control of some important reactions in vivo.     

Subdividing repressor function: DNA binding affinity, selectivity, and allostery can be altered by amino acid substitution of nonconserved residues in a LacI/GalR homologue

Many mutations that impact protein function occur at residues that do not directly contact ligand. To understand the functional contributions from the sequence that links the DNA-binding and regulatory domains of the LacI/GalR homologues, we have created a chimeric protein (LLhP), which comprises the LacI DNA-binding domain, the LacI linker, and the PurR regulatory domain. Although DNA binding site residues are identical in LLhP and LacI, thermodynamic measurements of DNA binding affinity show that LLhP does not discriminate between alternative DNA ligands as well as LacI. In addition, small-angle scattering experiments show that LLhP is more compact than LacI. When DNA is released, LacI shows a 20 A increase in length that was previously attributed to unfolding of the linker. This change is not seen in apo-LLhP, even though the linker sequences of the two proteins are identical. Together, results indicate that long-range functional and structural changes are propagated across the interface that forms between the linker and regulatory domain. These changes could be mediated via the side chains of several linker residues that contact the regulatory domains of the naturally occurring proteins, LacI and PurR. Substitution of these residues in LLhP leads to a range of functional effects. Four variants exhibit altered affinity for DNA, with no changes in selectivity or allosteric response. Another two result in proteins that bind operator DNA with very low affinity and no allosteric response, similar to LacI binding nonspecific DNA sequences. Two more substitutions simultaneously diminish affinity, enhance allostery, and profoundly alter DNA ligand selectivity. Thus, positions within the linker can be varied to modulate different aspects of repressor function.     

Sequence Dependencies of DNA Deformability and Hydration in the Minor Groove

DNA deformability and hydration are both sequence-dependent and are essential in specific DNA sequence recognition by proteins. However, the relationship between the two is not well understood. Here, systematic molecular dynamics simulations of 136 DNA sequences that differ from each other in their central tetramer revealed that sequence dependence of hydration is clearly correlated with that of deformability. We show that this correlation can be illustrated by four typical cases. Most rigid basepair steps are highly likely to form an ordered hydration pattern composed of one water molecule forming a bridge between the bases of distinct strands, but a few exceptions favor another ordered hydration composed of two water molecules forming such a bridge. Steps with medium deformability can display both of these hydration patterns with frequent transition. Highly flexible steps do not have any stable hydration pattern. A detailed picture of this correlation demonstrates that motions of hydration water molecules and DNA bases are tightly coupled with each other at the atomic level. These results contribute to our understanding of the entropic contribution from water molecules in protein or drug binding and could be applied for the purpose of predicting binding sites.     

Hydration of counterions interacting with DNA double helix: a molecular dynamics study

In the present work, molecular dynamics simulations have been carried out to study the dependence of counterion distribution around the DNA double helix on the character of ion hydration. The simulated systems consisted of DNA fragment d(CGCGAATTCGCG) in water solution with the counterions Na⁺, K⁺, Cs⁺ or Mg²⁺. The characteristic binding sites of the counterions with DNA and the changes in their hydration shell have been determined. The results show that due to the interaction with DNA at least two hydration shells of the counterions undergo changes. The first hydration shell of Na⁺, K⁺, Cs⁺, and Mg²⁺ counterions in the bulk consists of six, seven, ten, and six water molecules, respectively, while the second one has several times higher values. The Mg²⁺ and Na⁺ counterions, constraining water molecules of the first hydration shell, mostly form with DNA water-mediated contacts. In this case the coordination numbers of the first hydration shell do not change, while the coordination numbers of the second one decrease about twofold. The Cs⁺ and K⁺ counterions that do not constrain surrounding water molecules may be easily dehydrated, and when interacting with DNA their first hydration shell may be decreased by three and five water molecules, respectively. Due to the dehydration effect, these counterions can squeeze through the hydration shell of DNA to the bottom of the double helix grooves. The character of ion hydration establishes the correlation between the coordination numbers of the first and the second hydration shells.

Open image in new window

Graphical Abstract Hydration of counterions interacting with DNA double helix

     

Double threading through DNA: NMR structural study of a bis-naphthalene macrocycle bound to a thymine-thymine mismatch

The macrocyclic bis-naphthalene macrocycle (2,7-BisNP), belonging to the cyclobisintercalator family of DNA ligands, recognizes T-T mismatch sites in duplex DNA with high affinity and selectivity, as evidenced by thermal denaturation experiments and NMR titrations. The binding of this macrocycle to an 11-mer DNA oligonucleotide containing a T-T mismatch was studied using NMR spectroscopy and NMR-restrained molecular modeling. The ligand forms a single type of complex with the DNA, in which one of the naphthalene rings of the ligand occupies the place of one of the mismatched thymines, which is flipped out of the duplex. The second naphthalene unit of the ligand intercalates at the A-T base pair flanking the mismatch site, leading to encapsulation of its thymine residue via double stacking. The polyammonium linking chains of the macrocycle are located in the minor and the major grooves of the oligonucleotide and participate in the stabilization of the complex by formation of hydrogen bonds with the encapsulated thymine base and the mismatched thymine remaining inside the helix. The study highlights the uniqueness of this cyclobisintercalation binding mode and its importance for recognition of DNA lesion sites by small molecules.     

The osmotic sensitivity of netropsin analogue binding to DNA

The binding of a netropsin analogue to random sequence DNA, monitored by CD, is seen dependent on the concentration of neutral solutes. The binding free energy decreases linearly with solute osmolal concentration and the magnitude of the effect is insensitive to the chemical identity of the solute fur betaine, sorbitol, and triethylene glycol. These solutes appear to modulate binding through their effect on water activity and changes in the hydration of the drug and DNA in the complex reaction, not through a direct interaction with the reactants or the product. The dependence of binding constant on solute concentration can be interpreted as an additional binding of some 50–60 extra solute excluding water molecules by the complex. A water sensitivity of drug binding is further seen from the dependence of binding constants on the type of anion in solution. Anions in the Hofmeister series strongly affect bulk water free energies and entropies. The differences in netropsin analogue binding to DNA with Cl⁻, F⁻, and CIO are consistent with the effect observed with neutral solutes. The ability to measure changes in water binding associated with a specific DNA interaction is a first step toward correlating changes in hydration with the strength and specificity of binding. © 1995 John Wiley & Sons, Inc.     

Conformational flexibility of B-DNA at 0.74 A resolution: d(CCAGTACTGG)(2)

The affinity and specificity of a ligand for its DNA site is a function of the conformational changes between the isolated and complexed states. Although the structures of a hydroxypyrrole-imidazole-pyrrole polyamide dimer with 5'-CCAGTACTGG-3' and the trp repressor recognizing the sequence 5'-GTACT-3' are known, the baseline conformation of the DNA site would contribute to our understanding of DNA recognition by these ligands. The 0.74 A resolution structure of a B-DNA double helix, 5'-CCAGTACTGG-3', has been determined by X-ray crystallography. Six of the nine phosphates, two of four bound calcium ions and networks of water molecules hydrating the oligonucleotide have alternate conformations. By contrast, nine of the ten bases have a single, unique conformation with hydrogen atoms visible in most cases. The polyamide molecules alter the geometry of the phosphodiester backbone, and the water molecules mediating contacts in the trp repressor/operator complex are conserved in the unliganded DNA. Furthermore, the multiple conformational states, ions and hydration revealed by this ultrahigh resolution structure of a B-form oligonucleotide are potentially general considerations for understanding DNA-binding affinity and specificity by ligands.     

Quantitative measurement of water diffusion lifetimes at a protein/DNA interface by NMR

Hydration site lifetimes of slowly diffusing water molecules at the protein/DNA interface of the vnd/NK-2 homeodomain DNA complex were determined using novel three-dimensional NMR techniques. The lifetimes were calculated using the ratios of ROE and NOE cross-relaxation rates between the water and the protein backbone and side chain amides. This calculation of the lifetimes is based on a model of the spectral density function of the water-protein interaction consisting of three timescales of motion: fast vibrational/rotational motion, diffusion into/out of the hydration site, and overall macromolecular tumbling. The lifetimes measured ranged from approximately 400 ps to more than 5 ns, and nearly all the slowly diffusing water molecules detected lie at the protein/DNA interface. A quantitative analysis of relayed water cross-relaxation indicated that even at very short mixing times, 5 ms for ROESY and 12 ms for NOESY, relay of magnetization can make a small but detectable contribution to the measured rates. The temperature dependences of the NOE rates were measured to help discriminate direct dipolar cross-relaxation from chemical exchange. Comparison with several X-ray structures of homeodomain/DNA complexes reveals a strong correspondence between water molecules in conserved locations and the slowly diffusing water molecules detected by NMR. A homology model based on the X-ray structures was created to visualize the conserved water molecules detected at the vnd/NK-2 homeodomain DNA interface. Two chains of water molecules are seen at the right and left sides of the major groove, adjacent to the third helix of the homeodomain. Two water-mediated hydrogen bond bridges spanning the protein/DNA interface are present in the model, one between the backbone of Phe8 and a DNA phosphate, and one between the side chain of Asn51 and a DNA phosphate. The hydrogen bond bridge between Asn51 and the DNA might be especially important since the DNA contact made by the invariant Asn51 residue, seen in all known homeodomain/DNA structures, is critical for binding affinity and specificity.     

The effect of structural water molecules on the normal mode spectrum of dTn . dAn x dTn DNA

In this paper we present a theoretical treatment of triplex B type DNA hydration using normal mode calculation techniques. Discrete solvent is added as spines of hydration in the Watson-Crick and Crick-Hoogsteen grooves as well as water bridges between the Phosphate groups. The effect of binding the discrete structural waters on the normal mode of vibration of the system was studied by introducing a parameter, Xw, that is proportional to the degree of water binding and inversely proportional to the relative humidity (RH) of the system. We examined the variation of the dipole moments of characteristic modes with Xw. The results show that there is a direct relationship between the degree of binding of the water molecules to the atoms in the triple helix, the relative humidity of the system and the conformation and stability of the triple helix. At high RH and Xw = 0:0 the triple helix has mostly B type conformation characteristics, with C'2 -endo sugars. The emergence of normal modes of vibration characteristic to the A type conformation (C'3 - endo sugars) at Xw = 0:4 and 60% RH indicates a conformational shift towards A-type for some of the sugars between Xw = 0.2 (80% RH) and Xw = 0.4 (60% RH). These results are in agreement with the "economy of hydration hypothesis" of Saenger (Saenger et al., 1986) which maintains that the main difference in the hydration of A- and B- forms of DNA is the presence of water bridges between adjacent Phosphate groups in the low-hydration A-form but not in the B- form. Free energy calculations for the triplex DNA with structural waters show that there is a minimum of the free energy at Xw = 0.2 and the free energy increases with Xw and becomes larger than the free energy of the B conformation without structural waters for Xw equal to and larger than 0.4. This result indicates that the B conformation is more stable with bound structural water molecules (for degrees of water binding that are not over 20% higher than the degree of binding between bulk water molecules). The structural water molecules are bound much tighter in the A conformation than in the B conformation. The model predicts that the B to A transition occurs at higher relative humidities in D2O than in H2O. Part of these results (Dadarlat, 1997) have been subsequently confirmed by the experimental work and MD simulations of Ouali (Ouali et al., 1997). The experimental results showed that the N-type sugars corresponding to the A conformation are clearly detected below 75% RH.     

Binding of histone H1 to DNA is described by an allosteric model

Equilibrium binding data were analyzed to characterize the interaction of the linker histone H1 degrees with unmodified T4 phage DNA. Data were cast into the Scatchard-type plot described by McGhee and von Hippel and fit to their eponymous model for nonspecific binding of ligand to DNA. The data were not fit by the simple McGhee-von Hippel model, nor fit satisfactorily by the inclusion of a cooperativity parameter. Instead, the interaction appeared to be well described by Crothers' allosteric model, in which the higher affinity of the protein for one conformational form of the DNA drives an allosteric transition of the DNA to the conformational form with higher affinity (form 2). At 214 mM Na(+), the observed affinity K for an isolated site on unmodified T4 bacteriophage DNA in the form 2 conformation is 4.5 x 10(7) M(-1). The binding constant for an isolated site on DNA in the conformation with lower affinity, form 1, appears to be about 10-fold lower. Binding affinity is dependent on ion concentration: the magnitude of K is about 10-fold higher at 14 mM (5.9 x 10(8) M(-1) for form 2 DNA) than at 214 mM Na(+) concentration.     

DNA recognition of a 24-mer peptide derived from RecA protein

A novel 24-residue peptide (L2-G), Ile-Arg-Met-Lys-Ile-Gly-Val-Met-Phe-Gly-Asn-Pro-Glu-Thr-Thr-Thr-Gly-Gly-Asn-Ala-Leu-Lys-Phe-Tyr, derived from RecA can discriminate a single-stranded DNA (ssDNA) from a double-stranded DNA (dsDNA) and a new developed support with this peptide recognizes not dsDNA but ssDNA. The 24-mer peptide with L2 and helix G amino acids of Escherichia coli RecA protein showed the ssDNA binding property with more than 1000 times affinity difference for the dsDNA. However, truncated 15-mer peptide showed no ssDNA binding activity. In the ssDNA binding, L2-G changed its conformation with the perturbation of an alpha-helix structure. The ssDNA binding and the DNA discrimination property of this peptide were due to almost all L2 and helix G amino acids, respectively. This result is useful to design synthetic peptides as functional materials for DNA recognition.     

Use of the parmbsc0 force field and trajectory analysis to study the binding of netropsin to the DNA fragment (5'CCAATTGG)(2) in the presence of excess NaCl salt in aqueous solution

The parmbsc0 force field was applied to study in detail the binding of netropsin, at a salt concentration of 0.28M Na(+), to the minor groove of an 8-mer (5'CCAATTGG)(2) DNA duplex forming a netropsin·DNA complex which previously has been characterized by X-ray crystallography, albeit with the use of closely related DNA duplexes. The X-ray structure revealed that the terminal guanidinium and amidinium groups of netropsin interact with the extreme ends of the palindromic AATT sequence of the receptor DNA. The parmbsc0 parameters of B-DNA and AMBER v9 parameters of netropsin generated a stable 6ns molecular dynamics (MD) trajectory for a 1:1 class I binding motif of this complex. Trajectory analysis for the salt and hydration effects on the binding of netropsin to the 8-mer DNA duplex revealed that 18 water molecules and 2 Na(+) are displaced from the DNA upon netropsin binding. A hydration density map of the complex parallels the X-ray data showing that two structured water molecules are localized near the netropsin guanidinium and amidinium groups forming H-bond bridges between the receptor and the ligand.