Dynamic DNA contacts observed in the NMR structure of winged helix protein-DNA complex.

Genesis is an HNF-3/fkh homologous protein. By using multi-dimensional NMR techniques, we have obtained the solution structure and backbone dynamics of Genesis complexed with a 17 base-pair DNA. Our results indicate that both the local folding and dynamic properties of Genesis are perturbed when it binds to the DNA site. Our data show that a conserved flexible amino acid sequence (wing 1) makes dynamic contacts to DNA in the complex and a short helix is induced by Genesis-DNA interactions. Our data indicate that, unlike the HNF-3gamma/DNA complex, a magnesium ion is not required in forming the stable Genesis-DNA complex.     

Amino acid substitutions in a long flexible sequence influence thermodynamics and internal dynamic properties of winged helix protein genesis and its DNA complex

The hepatocyte nuclear factor (HNF)-3 homologous DNA binding domain is a highly conserved motif that contains a well-folded helix-turn-helix motif and two highly dynamic wings. Although the function and the properties of this motif have been intensively studied, the role of the internal wing (wing 1) is not well understood. In this study, amino acid substitutions were introduced into wing 1 of a conserved HNF-3 homologous protein, Genesis, and heteronuclear NMR, circular dichroism, DNA gel-shift assay, and fluorescent methods were employed to study and compare the properties of both wild-type and variant Genesis proteins. The data indicate that even though the substitutions are located on a dynamic wing outside the hydrophobic core sequences, they still globally influence biophysical properties of DNA-free Genesis and its DNA complex.     

Local dynamics of DNA probed with optical absorption spectroscopy of bound ethidium bromide.

We have studied the local dynamics of calf thymus double-helical DNA by means of an "optical labeling" technique. The study has been performed by measuring the visible absorption band of the cationic dye ethidium bromide, both free in solution and bound to DNA, in the temperature interval 360-30 K and in two different solvent conditions. The temperature dependence of the absorption line shape has been analyzed within the framework of the vibronic coupling theory, to extract information on the dynamic properties of the system; comparison of the thermal behavior of the absorption band of free and DNA-bound ethidium bromide gave information on the local dynamics of the double helix in the proximity of the chromophore. For the dye free in solution, large spectral heterogeneity and coupling to a "bath" of low-frequency (soft) modes is observed; moreover, anharmonic motions become evident at suitably high temperatures. The average frequency of the soft modes and the amplitude of anharmonic motions depend upon solvent composition. For the DNA-bound dye, at low temperatures, heterogeneity is decreased, the average frequency of the soft modes is increased, and anharmonic motions are hindered. However, a new dynamic regime characterized by a large increase in anharmonic motions is observed at temperatures higher than approximately 280 K. The DNA double helix therefore appears to provide, at low temperatures, a rather rigid environment for the bound chromophore, in which conformational heterogeneity is reduced and low-frequency motions (both harmonic vibrations and anharmonic contributions) are hindered. The system becomes anharmonic at approximately 180 K; however, above approximately 280 K, anharmonicity starts to increase much more rapidly than for the dye free in solution; this can be attributed to the onset of wobbling of the dye in its intercalation site, which is likely connected with the onset of (functionally relevant) DNA motions, involving local opening/unwinding of the double helix. As shown by parallel measurements of the melting curves, these motions precede the melting of the double helix and depend upon solvent composition much more than does the melting itself.     

Nontarget DNA binding shapes the dynamic landscape for enzymatic recognition of DNA damage

The DNA repair enzyme human uracil DNA glycosylase (UNG) scans short stretches of genomic DNA and captures rare uracil bases as they transiently emerge from the DNA duplex via spontaneous base pair breathing motions. The process of DNA scanning requires that the enzyme transiently loosen its grip on DNA to allow stochastic movement along the DNA contour, while engaging extrahelical bases requires motions on a more rapid timescale. Here, we use NMR dynamic measurements to show that free UNG has no intrinsic dynamic properties in the millisecond to microsecond and subnanosecond time regimes, and that the act of binding to nontarget DNA reshapes the dynamic landscape to allow productive millisecond motions for scanning and damage recognition. These results suggest that DNA structure and the spontaneous dynamics of base pairs may drive the evolution of a protein sequence that is tuned to respond to this dynamic regime.     

Structure and dynamics of the complex of single stranded DNA binding protein of Escherichia coli with circular single stranded DNA of filamentous phages

We have analyzed the equilibrium and nonequilibrium properties of the complex of the single stranded DNA binding protein of Escherichia coli (EcoSSB) and circular single stranded DNA of filamentous phages M13mp8 and F1 using static and dynamic light scattering, analytical ultracentrifugation and electron microscopy. Upon binding to the single stranded DNA the EcoSSB tetramer replaces an equivalent volume of water trapped within the coiled single stranded DNA and hinders the folding of the single stranded DNA into secondary structures at all salt concentrations. The salt dependent compaction of the stoichiometric complex can be described assuming a flexible polyelectrolyte chain. The solution structure of the macromolecular complex is a random coil and in the electron microscope a beaded flexible structure of the complex with a bead diameter of 6 nm appears at all salt concentrations used. The internal motions of the stoichiometric complex can be described by the Rouse-Zimm model of polymer dynamics. The segmental mobility of the complex can be correlated with changes in the binding site size of the EcoSSB tetramer; it indicates the presence of interactions between EcoSSB tetramers bound to single stranded DNA.     

Strand opening-deficient Escherichia coli RNA polymerase facilitates investigation of closed complexes with promoter DNA: effects of DNA sequence and temperature

Formation of the strand-separated, open complex between RNA polymerase and a promoter involves several intermediates, the first being the closed complex in which the DNA is fully base-paired. This normally short lived complex has been difficult to study. We have used a mutant Escherichia coli RNA polymerase, deficient in promoter DNA melting, and variants of the P(R) promoter of bacteriophage lambda to model the closed complex intermediate at physiologically relevant temperatures. Our results indicate that in the closed complex, RNA polymerase recognizes base pairs as double-stranded DNA even in the region that becomes single-stranded in the open complex. Additionally, a particular base pair in the -35 region engages in an important interaction with the RNA polymerase, and a DNase I-hypersensitive site, pronounced in the promoter DNA of the open complex, was not present. The effect of temperature on closed complex formation was found to be small over the temperature range from 15 to 37 degrees C. This suggests that low temperature complexes of wild type RNA polymerase and promoter DNA may adequately model the closed complex.     

Comparison of backbone dynamics of the type III antifreeze protein and antifreeze-like domain of human sialic acid synthase

Antifreeze proteins (AFPs) are found in a variety of cold-adapted (psychrophilic) organisms to promote survival at subzero temperatures by binding to ice crystals and decreasing the freezing temperature of body fluids. The type III AFPs are small globular proteins that consist of one α-helix, three 3₁₀-helices, and two β-strands. Sialic acids play important roles in a variety of biological functions, such as development, recognition, and cell adhesion and are synthesized by conserved enzymatic pathways that include sialic acid synthase (SAS). SAS consists of an N-terminal catalytic domain and a C-terminal antifreeze-like (AFL) domain, which is similar to the type III AFPs. Despite having very similar structures, AFL and the type III AFPs exhibit very different temperature-dependent stability and activity. In this study, we have performed backbone dynamics analyses of a type III AFP (HPLC12 isoform) and the AFL domain of human SAS (hAFL) at various temperatures. We also characterized the structural/dynamic properties of the ice-binding surfaces by analyzing the temperature gradient of the amide proton chemical shift and its correlation with chemical shift deviation from random coil. The dynamic properties of the two proteins were very different from each other. While HPLC12 was mostly rigid with a few residues exhibiting slow motions, hAFL showed fast internal motions at low temperature. Our results provide insight into the molecular basis of thermostability and structural flexibility in homologous psychrophilic HPLC12 and mesophilic hAFL proteins.     

The 'glass transition' in protein dynamics: what it is, why it occurs, and how to exploit it

All proteins undergo a dramatic change in their dynamical properties at approximately 200 K. Above this temperature, their dynamic behavior is dominated by large-scale collective motions of bonded and nonbonded groups of atoms. At lower temperatures, simple harmonic vibrations predominate. The transition has been described as a 'glass transition' to emphasize certain similarities between the change in dynamic behavior of individual protein molecules and the changes in viscosity and other properties of liquids when they form a glass. The glass transition may reflect the intrinsic temperature dependence of the motions of atoms in the protein itself, in the bound solvent on the surface of the protein, or it may reflect contributions from both. Protein function is significantly altered below this transition temperature; a fact that can be exploited to trap normally unstable intermediates in enzyme-catalyzed reactions and stabilize them for periods long enough to permit their characterization by high-resolution protein crystallography.     

Molecular dynamics studies on free and bound targets of the bovine papillomavirus type I e2 protein: the protein binding effect on DNA and the recognition mechanism

Molecular dynamics simulations of a total duration of 30 ns in explicit solvent were carried out on the BPV-1-E2 protein complexed to a high-affinity DNA target containing the two hydrogen-bonded ACCG.CGGT half-sites separated by the noncontacted ACGT sequence. The analysis of the trajectories focuses on the DNA structure and on the dynamics. The data are compared to those issued from recent simulations made on three free targets that recognize E2 with different affinities. E2 does not drastically perturb the mechanic properties of the free DNA: the structural relationships between the BI/BII backbone substates and some helical parameters are preserved in the complex despite a severe slowing down of the phosphate group motions. The structures of both free and bound half-sites are very close to each other although the conformational space explored by these regions is narrowed when they are contacted by the protein. The enhanced plasticity found in the best free target spacers, mainly manifested through the backbone motions, allows a clear overlap between several free and bound global DNA features such as the base displacement. Furthermore, this flexibility is preserved in the complex. Our results support the hypothesis that E2 takes advantage of free predistorted structures that may minimize the DNA deformation cost. In addition, we observe that E2 is far from totally stiffening the DNA, suggesting that the entropic penalty inherent in the complex formation could be limited.     

Solving traveling salesman problems with DNA molecules encoding numerical values

We introduce a DNA encoding method to represent numerical values and a biased molecular algorithm based on the thermodynamic properties of DNA. DNA strands are designed to encode real values by variation of their melting temperatures. The thermodynamic properties of DNA are used for effective local search of optimal solutions using biochemical techniques, such as denaturation temperature gradient polymerase chain reaction and temperature gradient gel electrophoresis. The proposed method was successfully applied to the traveling salesman problem, an instance of optimization problems on weighted graphs. This work extends the capability of DNA computing to solving numerical optimization problems, which is contrasted with other DNA computing methods focusing on logical problem solving.     

Slow motions in chicken villin headpiece subdomain probed by cross-correlated NMR relaxation of amide NH bonds in successive residues

The villin headpiece subdomain (HP36) is a widely used system for protein-folding studies. Nuclear magnetic resonance cross-correlated relaxation rates arising from correlated fluctuations of two N-H^N dipole-dipole interactions involving successive residues were measured at two temperatures at which HP36 is at least 99% folded. The experiment revealed the presence of motions slower than overall tumbling of the molecule. Based on the theoretical analysis of the spectral densities we show that the structural and dynamic contributions to the experimental cross-correlated relaxation rate can be separated under certain conditions. As a result, dynamic cross-correlated order parameters describing slow microsecond-to-millisecond motions of N-H bonds in neighboring residues can be introduced for any extent of correlations in the fluctuations of the two bond vectors. These dynamic cross-correlated order parameters have been extracted for HP36. The comparison of their values at two different temperatures indicates that when the temperature is raised, slow motions increase in amplitude. The increased amplitude of these fluctuations may reflect the presence of processes directly preceding the unfolding of the protein.     

Dynamics of the Mrf-2 DNA-binding domain free and in complex with DNA

Mrf-2 is a member of a new class of DNA-binding proteins known as the AT-rich interaction domain family or ARID. Chemical shift indices and characteristic NOE values indicate that the three-dimensional structure of the Mrf-2 ARID in complex with DNA is nearly identical to that of the free protein. The backbone dynamics of the Mrf-2 domain free and in complex with DNA have been characterized by (15)N NMR relaxation measurements and model-free analysis. Chemical shift perturbations and dynamic studies suggest that two flexible interhelical loops, the flexible C-terminal tail, and one alpha-helix are involved in DNA recognition, indicating the importance of protein dynamics in DNA binding. Some well-structured regions, in particular the putative DNA-contacting helix, in Mrf-2 show a decrease in the order parameters (S(2)) upon complex formation. The less well-structured loops and the unstructured C-terminus show reduced flexibility upon DNA binding. In addition, the model-free analysis indicates motions on the picosecond to nanosecond and micro- to millisecond time scales at the DNA-binding surface of the bound Mrf-2 ARID, suggesting a model where interactions between the protein and DNA are highly dynamic.     

DNA binding of a non-sequence-specific HMG-D protein is entropy driven with a substantial non-electrostatic contribution

The thermal properties of two forms of the Drosophila melanogaster HMG-D protein, with and without its highly basic 26 residue C-terminal tail (D100 and D74) and the thermodynamics of their non-sequence-specific interaction with linear DNA duplexes were studied using scanning and titration microcalorimetry, spectropolarimetry, fluorescence anisotropy and FRET techniques at different temperatures and salt concentrations. It was shown that the C-terminal tail of D100 is unfolded at all temperatures, whilst the state of the globular part depends on temperature in a rather complex way, being completely folded only at temperatures close to 0 degrees C and unfolding with significant heat absorption at temperatures below those of the gross denaturational changes. The association constant and thus Gibbs energy of binding for D100 is much greater than for D74 but the enthalpies of their association are similar and are large and positive, i.e. DNA binding is a completely entropy-driven process. The positive entropy of association is due to release of counterions and dehydration upon forming the protein/DNA complex. Ionic strength variation showed that electrostatic interactions play an important but not exclusive role in the DNA binding of the globular part of this non-sequence-specific protein, whilst binding of the positively charged C-terminal tail of D100 is almost completely electrostatic in origin. This interaction with the negative charges of the DNA phosphate groups significantly enhances the DNA bending. An important feature of the non-sequence-specific association of these HMG boxes with DNA is that the binding enthalpy is significantly more positive than for the sequence-specific association of the HMG box from Sox-5, despite the fact that these proteins bend the DNA duplex to a similar extent. This difference shows that the enthalpy of dehydration of apolar groups at the HMG-D/DNA interface is not fully compensated by the energy of van der Waals interactions between these groups, i.e. the packing density at the interface must be lower than for the sequence-specific Sox-5 HMG box.     

Backbone dynamics of receptor binding and antigenic regions of a Pseudomonas aeruginosa pilin monomer

Pilin is the major structural protein that forms type IV pili of various pathogenic bacteria, including Pseudomonas aeruginosa. Pilin is involved in attachment of the bacterium to host cells during infection, in the initiation of immune response, and serves as a receptor for a variety of bacteriophage. We have used (15)N nuclear magnetic resonance relaxation measurements to probe the backbone dynamics of an N-terminally truncated monomeric pilin from P. aeruginosa strain K122-4. (15)N-T(1), -T(2), and [(1)H]-(15)N nuclear Overhauser enhancement measurements were carried out at three magnetic field strengths. The measurements were interpreted using the Lipari-Szabo model-free analysis, which reveals the amplitude of spatial restriction for backbone N-NH bond vectors with respect to nano- to picosecond time-scale motions. Regions of well-defined secondary structure exhibited consistently low-amplitude spatial fluctuations, while the terminal and loop regions showed larger amplitude motions in the subnano- to picosecond time-scale. Interestingly, the C-terminal disulfide loop region that contains the receptor binding domain was found to be relatively rigid on the pico- to nanosecond time-scale but exhibited motion in the micro- to millisecond time-scale. It is notable that this disulfide loop displays a conserved antigenic epitope and mediates binding to the asialo-GM(1) cell surface receptor. The present study suggests that a rigid backbone scaffold mediates attachment to the host cell receptor, and also maintains the conformation of the conserved antigenic epitope for antibody recognition. In addition, slower millisecond time-scale motions are likely to be crucial for conferring a range of specificity for these interactions. Characterization of pilin dynamics will aid in developing a detailed understanding of infection, and will facilitate the design of more efficient anti-adhesin synthetic vaccines and therapeutics against pathogenic bacteria containing type IV pili.     

Sequence-dependent motions of DNA: a normal mode analysis at the base-pair level

Computer simulation of the dynamic structure of DNA can be carried out at various levels of resolution. Detailed high resolution information about the motions of DNA is typically collected for the atoms in a few turns of double helix. At low resolution, by contrast, the sequence-dependence features of DNA are usually neglected and molecules with thousands of base pairs are treated as ideal elastic rods. The present normal mode analysis of DNA in terms of six base-pair "step" parameters per chain residue addresses the dynamic structure of the double helix at intermediate resolution, i.e., the mesoscopic level of a few hundred base pairs. Sequence-dependent effects are incorporated into the calculations by taking advantage of "knowledge-based" harmonic energy functions deduced from the mean values and dispersion of the base-pair "step" parameters in high-resolution DNA crystal structures. Spatial arrangements sampled along the dominant low frequency modes have end-to-end distances comparable to those of exact polymer models which incorporate all possible chain configurations. The normal mode analysis accounts for the overall bending, i.e., persistence length, of the double helix and shows how known discrepancies in the measured twisting constants of long DNA molecules could originate in the deformability of neighboring base-pair steps. The calculations also reveal how the natural coupling of local conformational variables affects the global motions of DNA. Successful correspondence of the computed stretching modulus with experimental data requires that the DNA base pairs be inclined with respect to the direction of stretching, with chain extension effected by low energy transverse motions that preserve the strong van der Waals' attractions of neighboring base-pair planes. The calculations further show how one can "engineer" the macroscopic properties of DNA in terms of dimer deformability so that polymers which are intrinsically straight in the equilibrium state exhibit the mesoscopic bending anisotropy essential to DNA curvature and loop formation.     

Local deformations revealed by dynamics simulations of DNA polymerase Beta with DNA mismatches at the primer terminus

Nanosecond dynamics simulations for DNA polymerase beta (pol beta)/DNA complexes with three mismatched base-pairs, namely GG, CA, or CC (primer/template) at the DNA polymerase active site, are performed to investigate the mechanism of polymerase opening and how the mispairs may affect the DNA extension step; these trajectories are compared to the behavior of a pol beta/DNA complex with the correct GC base-pair, and assessed with the aid of targeted molecular dynamics (TMD) simulations of all systems from the closed to the open enzyme state. DNA polymerase conformational changes (subdomain closing and opening) have been suggested to play a critical role in DNA synthesis fidelity, since these changes are associated with the formation of the substrate-binding pocket for the nascent base-pair. Here we observe different large C-terminal subdomain (thumb) opening motions in the simulations of pol beta with GC versus GG base-pairs. Whereas the conformation of pol beta in the former approaches the observed open state in the crystal structures, the enzyme in the latter does not. Analyses of the motions of active-site protein/DNA residues help explain these differences. Interestingly, rotation of Arg258 toward Asp192, which coordinates both active-site metal ions in the closed "active" complex, occurs rapidly in the GG simulation. We have previously suggested that this rotation is a key slow step in the closed to open transition. TMD simulations also point to a unique pathway for Arg258 rotation in the GG mispair complex. Simulations of the mismatched systems also reveal distorted geometries in the active site of all the mispair complexes examined. The hierarchy of the distortions (GG>CC>CA) parallels the experimentally deduced inability of pol beta to extend these mispairs. Such local distortions would be expected to cause inefficient DNA extension and polymerase dissociation and thereby might lead to proofreading by an extrinsic exonuclease. Thus, our studies on the dynamics of pol beta opening in mismatch systems provide structural and dynamic insights to explain experimental results regarding inefficient DNA extension following misincorporation; these details shed light on how proofreading may be invoked by the abnormal active-site geometry.     

Diffusion-controlled protein–DNA association: Influence of segemental diffusion of the DNA

The intrachain reaction theory of Wilemki and Fixman [(1974) J. Chem. Phys. 60 , 866–877] is used to assess the influence of internal DNA motions on various protein–DNA association schemes. It is found that, for large proteins, the diffusional association rate can be totally dominated by these motions rather than by the free-trans-lational diffusion rates. Also, the time required for the diffusion together of two DNA segments is estimated. This estimate can be used to provide an upper limit for the rate of intrachain cyclization, and also for the effective intrachain transfer rate of a protein bound to a DNA chain.