Effects of DNA sequence and histone-histone interactions on nucleosome placement

Using competitive reconstitution, we have refined the parameters for the binding of histone octamers to artificial nucleosome-positioning sequences of the form: (A/T3nn(G/C)3nn. We find that the optimal period between flexible segments is approximately 10.1 base-pairs, supporting the view that the DNA on the nucleosome surface is overwound. The strongest requirement for flexible DNA is near the protein dyad. However, we see no indication of changes in DNA helical repeat in this region. Using a series of repetitive sequences, we confirm that neither all A/T-rich nor all G/C-rich regions are identical in promoting nucleosome formation. Surprisingly, A/T-rich segments containing the TpA step, subject to purine-purine clash in the minor groove, favor nucleosome formation over sequences lacking this step. Short tracts of adenine residues are found to position on the histone surface like other A/T-rich regions, in the manner predicted by the direction of their sequence-directed bends as determined by electrophoretic methods. Tracts containing five adenine residues are extremely aniostropic in their flexibility and are strongly detrimental to nucleosome formation when positioned for major groove compression. Longer adenine tracts are found to position near the ends of the nucleosomal DNA. However, other positions may be occupied by an A12 tract, with only a minor penalty in the free energy of nucleosome formation. Overall, reconstituted nucleosome positions are translationally degenerate, suggesting a weak dependence on DNA flexibility for nucleosome positioning. Dinucleosomal reconstitutions on tandem dimers of the 5 S RNA gene of Lytechinus variegatus demonstrate a weak phasing dependence for the interaction between nucleosomes. This interaction is maximal for the 202 base-pair repeat and suggests a co-operative mechanism for the formation of ordered nucleosomal arrays based on a combination of DNA flexibility and nucleosome-nucleosome interactions.     

Recognition and processing of the origin of transfer DNA by conjugative relaxase TrwC

Relaxases are DNA strand transferases that catalyze the initial and final stages of DNA processing during conjugative cell-to-cell DNA transfer. Upon binding to the origin of transfer (oriT) DNA, relaxase TrwC melts the double helix. The three-dimensional structure of the relaxase domain of TrwC in complex with its cognate DNA at oriT shows a fold built on a two-layer alpha/beta sandwich, with a deep narrow cleft that houses the active site. The DNA includes one arm of an extruded cruciform, an essential feature for specific recognition. This arm is firmly embraced by the protein through a beta-ribbon positioned in the DNA major groove and a loop occupying the minor groove. It is followed by a single-stranded DNA segment that enters the active site, after a sharp U-turn forming a hydrophobic cage that traps the N-terminal methionine. Structural analysis combined with site-directed mutagenesis defines the architecture of the active site.     

Interaction of nucleic acids with electrically charged surfaces. VI. A comparative study on the electrochemical behaviour of native and denatured DNAs at graphite electrodes.

Adsorption and electrochemical oxidation of deoxyribonucleic acid (DNA) at a pyrolytic graphite electrode (PGE) and a paraffin wax-impregnated spectroscopic graphite electrode (WISGE) were studied using differential pulse voltammetry. DNA is adsorbed at the surface of the graphite electrodes in a broad range of potentials including the potentials of electrochemical oxidation of DNA. Both native and denatured DNAs yield two single, well-defined and separated peaks, G and A, on the differential pulse voltammograms at the PGE and WISGE. The more negative peak, G, corresponds to electrochemical oxidation of adenine residues. Peaks G and A of native DNA occur at the same potentials as peaks G and A of denatured DNA. However, electrochemical oxidation of adenine and guanine residues at graphite electrodes is markedly suppressed in native DNA. The heights of the peaks G and A represent a sensitive indicator of the helix-coil transition of DNA. An analysis of the product of interaction of a sample of native DNA with a large pyrolytic graphite electrode in the presence of formaldehyde at approximately neutral pH did not prove changes in the secondary structure of native DNA due to its interaction with the graphite electrode. It is suggested that the decreased differential pulse-voltammetric activity of native DNA is connected with its decreased flexibility.     

Interaction of berenil with the tyrT DNA sequence studied by footprinting and molecular modelling. Implications for the design of sequence-specific DNA recognition agents.

We have developed a technique of partially-restrained molecular mechanics enthalpy minimisation which enables the sequence-dependence of the DNA binding of a non-intercalating ligand to be studied for arbitrary sequences of considerable length (greater than = 60 base-pairs). The technique has been applied to analyse the binding of berenil to the minor groove of a 60 base-pair sequence derived from the tyrT promoter; the results are compared with those obtained by DNAse I and hydroxyl radical footprinting on the same sequence. The calculated and experimentally observed patterns of binding are in good agreement. Analysis of the modelling data highlights the importance of DNA flexibility in ligand binding. Further, the electrostatic component of the interaction tends to favour binding to AT-rich regions, whilst the van der Waals interaction energy term favours GC-rich ones. The results also suggest that an important contribution to the observed preference for binding in AT-rich regions arises from lower DNA perturbation energies and is not accompanied by reduced DNA structural perturbations in such sequences. It is therefore concluded that those modes of DNA distortion favourable to binding are probably more flexible in AT-rich regions. The structure of the modelled DNA sequence has also been analysed in terms of helical parameters. For the DNA energy-minimised in the absence of berenil, certain helical parameters show marked sequence-dependence. For example, purine-pyrimidine (R-Y) base pairs show a consistent positive buckle whereas this feature is consistently negative for Y-R pairs. Further, CG steps show lower than average values of slide while GC steps show lower than average values of rise. Similar analysis of the modelling data from the calculations including berenil highlights the importance of DNA flexibility in ligand binding. We observe that the binding of berenil induces characteristic responses in different helical parameters for the base-pairs around the binding site. For example, buckle and tilt tend to become more negative to the 5'-side of the binding site and more positive to the 3'-side, while the base steps at either side of the centre of the site show increased twist and decreased roll.     

Structural and mutagenic analysis of foot-and-mouth disease virus 3C protease reveals the role of the beta-ribbon in proteolysis

The 3C protease (3C(pro)) from foot-and-mouth disease virus (FMDV), the causative agent of a widespread and economically devastating disease of domestic livestock, is a potential target for antiviral drug design. We have determined the structure of a new crystal form of FMDV 3C(pro), a chymotrypsin-like cysteine protease, which reveals features that are important for catalytic activity. In particular, we show that a surface loop which was disordered in previous structures adopts a beta-ribbon structure that is conformationally similar to equivalent regions on other picornaviral 3C proteases and some serine proteases. This beta-ribbon folds over the peptide binding cleft and clearly contributes to substrate recognition. Replacement of Cys142 at the tip of the beta-ribbon with different amino acids has a significant impact on enzyme activity and shows that higher activity is obtained with more hydrophobic side chains. Comparison of the structure of FMDV 3C(pro) with homologous enzyme-peptide complexes suggests that this correlation arises because the side chain of Cys142 contacts the hydrophobic portions of the P2 and P4 residues in the peptide substrate. Collectively, these findings provide compelling evidence for the role of the beta-ribbon in catalytic activity and provide valuable insights for the design of FMDV 3C(pro) inhibitors.     

Interaction of nucleic acids with electrically charged surfaces. VI. A comparative study on the electrochemical behaviour of native and denatured DNAs at graphite electrodes

Adsorption and electrochemical oxidation of deoxyribonucleic acid (DNA) at a pyrolytic graphite electrode (PGE) and a paraffin wax-impregnated spectroscopic graphite electrode (WISGE) were studied using differential pulse voltammetry. DNA is adsorbed at the surface of the graphite electrodes in a broad range of potentials including the potentials of electrochemical oxidation of DNA. Both native and denatured DNAs yield two single, well-defined and separated peaks, G and A, on the differential pulse voltammograms at the PGE and WISGE. The more negative peak, G, corresponds to electrochemical oxidation of guanine residues, whereas the more positive peak, A, corresponds to electrochemical oxidation of adenine residues. Peaks G and A of native DNA occur at the same potentials as peaks G and A of denatured DNA. However, electrochemical oxidation of adenine and guanine residues at graphite electrodes is markedly suppressed in native DNA. The heights of the peaks G and A represent a sensitive indicator of the helix-coil transition of DNA. An analysis of the product of interaction of a sample of native DNA with a large pyrolytic graphite electrode in the presence of formaldehyde at approximately neutral pH did not prove changes in the secondary structure of native DNA due to its interaction with the graphite electrode. It is suggested that the decreased differential pulse-voltammetric activity of native DNA is connected with its decreased flexibility.     

Counterion dependent variation of DNA secondary structure in (A . T) clusters: evidence by use of netropsin as a structural probe.

The interaction of the oligopeptide antibiotic netropsin (Nt) with (A . T) regions of DNA is characterized by a spectrum of discrete modes. This has been revealed by viscometric analysis, at 20 degrees C and 0.2 M "counterions", for NaDNA in a preceding and for NH4DNA in this paper. The increase of DNA contour length as induced by one Nt molecule was found to depend on the special mode only, while the respective stiffening is generally higher for NH4DNA. The latter property is interpreted in terms of an enhanced flexibility, relative to that of NaDNA, of the (A . T) cluster segments before complex formation. For some of the interaction modes of the DNA-Nt systems a difference in the number of corresponding binding sites has been observed. This phenomenon is understood by assuming an influence of the counterion species upon existing equilibria between different forms of the (A . T) cluster secondary structure. Not less than 5 to 10% of the total DNA are effected in this manner. Upper limits for the local differences in the axial rise per base pair are 0.04 nm and 0.02 nm.     

Relative flexibility of DNA and RNA: a molecular dynamics study

State of the art molecular dynamics simulations are used to study the structure, dynamics, molecular interaction properties and flexibility of DNA and RNA duplexes in aqueous solution. Special attention is paid to the deformability of both types of structures, revisiting concepts on the relative flexibility of DNA and RNA duplexes. Our simulations strongly suggest that the concepts of flexibility, rigidity and deformability are much more complex than usually believed, and that it is not always true that DNA is more flexible than RNA.     

Amino acid substitution in the C-terminal arm domain of HU-2 results in an enhanced affinity for DNA.

Three mutants of the Escherichia coli hupA gene, encoding the HU-2 protein, were constructed by synthetic oligodeoxyribonucleotide-directed, site-specific mutagenesis on M13mp18 vectors. The resulting HupAN10, HupAN11 and HupAN12 proteins contained Thr59-->Lys, Gln64-->Lys and Asn53-->Arg substitutions, respectively. These amino acid (aa) changes increased the positive charge of the N-terminal half of the two-strand, antiparallel beta-ribbon of the arm structure, which is believed to be a domain for DNA binding. The three mutant proteins bound to DNA more tightly than wild-type HU-2, and their affinities for DNA increased in the order of HupAN10, HupAN11, HupAN12. The mutant proteins showed a slightly increased HU activity for supporting Mu phage development. A mutant HU-2 protein with increased basicity, but with an altered aa sequence in the arm region due to a frameshift mutation, was also constructed. This mutant protein showed a reduced affinity to DNA and was unable to support Mu growth, suggesting that a unique aa sequence of the arm domain, rather than mere basicity of this domain, is required for efficient binding to DNA.     

The influence of terminal 3'', 5'' phosphates on conformations of dApdA.

Addition of 3' and 5' terminal phosphates to dApdA causes a decrease in conformational flexibility. pdApdAp has much fewer conformers with energies below 2.5 kcal./mole than dApdA. THE A, B and Watson-Crick (34) helices are the most preferred forms. Other important conformations are in the trans domain of psi. Thus, flexibility in psi as well as in omega and omega, and in the sugar pucker is indicated. The transformation from the B helix to the Watson-Crick helix follows a low energy path. This is significant since Watson-Crick conformations may be important for intercalation into nucleic acid polymers (40-42) above the dimer level. The B helix is preferred over the A form in these large DNA subunits.     

The bacterial histone-like protein HU specifically recognizes similar structures in all nucleic acids. DNA,RNA, and their hybrids

HU, a major component of the bacterial nucleoid, shares properties with histones, high mobility group proteins (HMGs), and other eukaryotic proteins. HU, which participates in many major pathways of the bacterial cell, binds without sequence specificity to duplex DNA but recognizes with high affinity DNA repair intermediates. Here we demonstrate that HU binds to double-stranded DNA, double-stranded RNA, and linear DNA-RNA duplexes with a similar low affinity. In contrast to this nonspecific binding to total cellular RNA and to supercoiled DNA, HU specifically recognizes defined structures common to both DNA and RNA. In particular HU binds specifically to nicked or gapped DNA-RNA hybrids and to composite RNA molecules such as DsrA, a small non-coding RNA. HU, which modulates DNA architecture, may play additional key functions in the bacterial machinery via its RNA binding capacity. The simple, straightforward structure of its binding domain with two highly flexible beta-ribbon arms and an alpha-helical platform is an alternative model for the elaborate binding domains of the eukaryotic proteins that display dual DNA- and RNA-specific binding capacities.     

Structure of a complex of catabolite gene activator protein and cyclic AMP refined at 2.5 A resolution

The structure of a dimer of the Escherichia coli catabolite gene activator protein has been refined at 2.5 A resolution to a crystallographic R-factor of 20.7% starting with coordinates fitted to the map at 2.9 A resolution. The two subunits are in different conformations and each contains one bound molecule of the allosteric activator, cyclic AMP. The amino-terminal domain is linked to the smaller carboxy-terminal domain by a nine-residue hinge region that exists in different conformations in the two subunits, giving rise to approximately a 30 degree rotation between the positions of the small domains relative to the larger domains. The amino-terminal domain contains an antiparallel beta-roll structure in which the interstrand hydrogen bonding is well-determined. The beta-roll can be described as a long antiparallel beta-ribbon that folds into a right-handed supercoil and forms part of the cyclic AMP binding site. Each cyclic AMP molecule is in an anti conformation and has ionic and hydrogen bond interactions with both subunits.     

Ultrafast interfacial solvation dynamics in specific protein DNA recognition

An overwhelming number of structural and functional studies on specific protein–DNA complexes reveal the existence of water molecules at the interaction interface. What role does the interfacial water molecules play in determining the specificity of association is thus a critical question. Herein, we have explored the dynamical role of minor groove water molecules and DNA side chain flexibility in lambda repressor–operator DNA interaction using well-characterized DNA minor groove binder dye, Hoechst 33258. The most striking finding of our studies reveals that the solvation time scale corresponding to the minor groove water molecules (∼50 ps) and DNA side chain flexibility (∼10 ns) remain unaltered even in protein–DNA complex in comparison to unbound operator DNA. The temperature dependent study further reveals the slower exchange of minor grove water molecules with bulk water in DNA–protein complex in comparison to the unbound DNA. Detailed structural studies including circular dichroism (CD) and Förster resonance energy transfer (FRET) have also been performed to elucidate the interaction between protein and DNA.     

A theoretical model for the prediction of sequence-dependent nucleosome thermodynamic stability

A theoretical model for predicting nucleosome thermodynamic stability in terms of DNA sequence is advanced. The model is based on a statistical mechanical approach, which allows the calculation of the canonical ensemble free energy involved in the competitive nucleosome reconstitution. It is based on the hypothesis that nucleosome stability mainly depends on the bending and twisting elastic energy to transform the DNA intrinsic superstructure into the nucleosomal structure. The ensemble average free energy is calculated starting from the intrinsic curvature, obtained by integrating the dinucleotide step deviations from the canonical B-DNA and expressed in terms of a Fourier series, in the framework of first-order elasticity. The sequence-dependent DNA flexibility is evaluated from the differential double helix thermodynamic stability. A large number of free-energy experimental data, obtained in different laboratories by competitive nucleosome reconstitution assays, are successfully compared to the theoretical results. They support the hypothesis that the stacking energies are the major factor in DNA rigidity and could be a measure of DNA stiffness. A dual role of DNA intrinsic curvature and flexibility emerges in the determination of nucleosome stability. The difference between the experimental and theoretical (elastic) nucleosome-reconstitution free energy for the whole pool of investigated DNAs suggests a significant role for the curvature-dependent DNA hydration and counterion interactions, which appear to destabilize nucleosomes in highly curved DNAs. This model represents an attempt to clarify the main features of the nucleosome thermodynamic stability in terms of physical-chemical parameters and suggests that in molecular systems with a large degree of complexity, the average molecular properties dominate over the local features, as in a statistical ensemble.     

Refined structure of the gene 5 DNA binding protein from bacteriophage fd

The three-dimensional structure of the gene 5 DNA binding protein (G5BP) from bacteriophage fd has been determined from a combination of multiple isomorphous replacement techniques, partial refinements and deleted fragment difference Fourier syntheses. The structure was refined using restrained parameter least-squares and difference Fourier methods to a final residual of R = 0.217 for the 3528 statistically significant reflections present to 2.3 A resolution. In addition to the 682 atoms of the protein, 12 solvent molecules were included. We describe here the dispositions and orientations of the amino acid side-chains and their interactions as visualized in the G5BP structure. The G5BP monomer of 87 peptide units is almost entirely in the beta-conformation, organized as a three-stranded sheet, a two-stranded beta-ribbon and a broad connecting loop. There is no alpha-helix present in the molecule. Two G5BP monomers are tightly interlocked about an intermolecular dyad axis to form a compact dimer unit of about 55 A X 45 A X 36 A. The dimer is characterized by two symmetry-related antiparallel clefts that traverse the monomer surfaces essentially perpendicular to the dyad axis. From the three-stranded antiparallel beta-sheet, formed from the first two-thirds of the sequence, extend three tyrosine residues (26, 34, 41), a lysine (46) and two arginine residues (16, 21) that, as indicated by other physical and chemical experiments, are directly involved in DNA binding. Other residues likely to share binding responsibility are arginine 80 extending from the beta-ribbon and phenylalanine 73 from the tip of this loop, but as provided, however, by the opposite monomer within each G5BP dimer pair. Thus, both symmetry-related DNA binding sites have a composite nature and include contributions from both elements of the dimer. The gene 5 dimer is clearly the active binding species, and the two monomers within the dyad-related pair are so structurally contiguous that one cannot be certain whether the isolated monomer would maintain its observed crystal structure. This linkage is manifested primarily as a skeletal core of hydrophobic residues that extends from the center of each monomer continuously through an intermolecular beta-barrel that joins the pair. Protruding from the major area of density of each monomer is an elongated wing of tenuous structure comprising residues 15 through 32, which is, we believe, intimately involved in DNA binding. This wing appears to be dynamic and mobile, even in the crystal and, therefore, is likely to undergo conformational change in the presence of the ligand.     

Infrared linear dichroism studies of DNA-drug complexes: quantitative determination of the drug-induced restriction of the B-A transition.

The B-A transition of films or fibers of NaDNA occurs at a relative humidity of 75-85%. The fraction of DNA that changed the conformation from B to A form can be determined quantitatively by infrared linear dichroism. DNA-binding drugs can 'freeze' a fraction of DNA in the B form. This fraction of DNA is in the B form and cannot be converted to A-DNA even at a reduced relative humidity of 54%. The 'freezing' potentiality of various drugs can be described by the 'freezing' index, FI, expressed in base pairs per added drug. Drugs with a high value of FI (more than eight base pairs per drug) were observed among both intercalating and groove-binding drugs. High values of FI imply restriction of the conformational flexibility of DNA significantly going beyond the binding site of the drug. This long-range effect of drugs on the conformational flexibility of DNA may be connected with the molecular mechanism of drug action. The freezing index FI is a new quantitative parameter of drug-DNA interaction that should be considered as a valuable tool for drug design.