A Refined Prediction Method for Gel Retardation of DNA Oligonucleotides from Dinucleotide Step Parameters: Reconciliation of DNA Bending Models with Crystal Structure Data

Abstract

The development and assessment of a prediction method for gel retardation and sequence dependent curvature of DNA based on dinucleotide step parameters are described. The method is formulated using the Babcock-Olson equations for base pair step geometry (1) and employs Monte Carlo simulated annealing for parameter optimization against experimental data. The refined base pair step parameters define a structural construct which, when the width of observed parameter distributions is taken into account, is consistent with the results of DNA oligonucleotide crystal structures. The predictive power of the method is demonstrated and tested via comparisons with DNA bending data on sets of sequences not included in the training set, including A-tracts with and without periodic helix phasing, phased A₄T₄ and T₄A₄ motifs, a sequence with a phased GGGCCC motif, some “unconventional” helix phasing sequences, and three short fragments of kinetoplast DNA from Crithidia fasiculata that exhibit significantly different behavior on non-denaturing polyacrylamide gels. The nature of the structural construct produced by the methodology is discussed with respect to static and dynamic models of structure and representations of bending and bendability. An independent theoretical account of sequence dependent chemical footprinting results is provided. Detailed analysis of sequences with A-tract induced axis bending forms the basis for a critical discussion of the applicability of wedge models, junction models and non A-tract, general sequence models for understanding the origin of DNA curvature at the molecular level.     

Molecular Dynamics Study of Intrinsic Stability in Six RNA Terminal Loop Motifs

Abstract

Single stranded RNA molecules can assume a wide range of tertiary structures beyond the canonical A-form double helix. Certain sequences, termed motifs, are more common than a random distribution would suggest. The existence of such motifs can be rationalized in structural terms. In this study, we have investigated the intrinsic structural stability of RNA terminal loop motifs using multiple MD simulations in explicit water. Representative loops were chosen from the major tetraloop motifs, including also the U-turn motif. Not all loops retain their folded starting structure, but lowering the temperature to 277 K, or adding adjacent base pairs from the stem to which the motif is attached, helps stabilizing the folded loop structure.     

On the mechanism of interaction between histone II B1 and DNA and histone II B2 and DNA

A mechanism is suggested at the molecular level whereby histone IIB₂ can act as a cross-link between two (or possibly three) adjacent and parallel strands of DNA double helix some 40 Å apart. Application of Prothero's rule and the Lewis probability functions indicate the probable locations of three a-helices and a number of β-turns. This, coupled with the requirement that the tertiary conformation of the histone be complementary to the DNA molecules and for as many basic groups as possible to bind to phosphate oxygens, allows us to suggest, on the basis of model building using accurate space-filling (CPK) models, a complex conformation that achieves this.A similar process applied to histone IIB₁, whose complete amino acid sequence is also known, shows the location of five probable a-helices, a number of β-turns, and a segment of β-pleated sheet. The basic amino acids are gathered in four groupings. Model building experiments suggest that histone IIB₁ forms a complex strut joining four parallel strands of DNA double helix that form a diamond with diameters 100 and 40 Å. In both these models the purpose and function of a fair proportion of the individual amino acids can be specified.This paper is the third and last of a series in this Journal in which models are presented for the tertiary conformation and function of all five histones of known (in whole or in part) amino acid sequence. This suggests that all five are concerned in packing the long DNA double helix, which may be in a “square helix” form, into the confined space of the chromosome. The hypotheses may be tested by a direct investigation of nucleoprotein in situ to see if these 40, 70, and 100 Å interhelical distances can be detected by biophysical methods.     

The flexibility of DNA double crossover molecules

Double crossover molecules are DNA structures containing two Holliday junctions connected by two double helical arms. There are several types of double crossover molecules, differentiated by the relative orientations of their helix axes, parallel or antiparallel, and by the number of double helical half-turns (even or odd) between the two crossovers. They are found as intermediates in meiosis and they have been used extensively in structural DNA nanotechnology for the construction of one-dimensional and two-dimensional arrays and in a DNA nanomechanical device. Whereas the parallel double helical molecules are usually not well behaved, we have focused on the antiparallel molecules; antiparallel molecules with an even number of half-turns between crossovers (termed DAE molecules) produce a reporter strand when ligated, facilitating their characterization in a ligation cyclization assay. Hence, we have estimated the flexibility of antiparallel DNA double crossover molecules by means of ligation-closure experiments. We are able to show that these molecules are approximately twice as rigid as linear duplex DNA.     

An Overview of Structural DNA Nanotechnology

Structural DNA Nanotechnology uses unusual DNA motifs to build target shapes and arrangements. These unusual motifs are generated by reciprocal exchange of DNA backbones, leading to branched systems with many strands and multiple helical domains. The motifs may be combined by sticky ended cohesion, involving hydrogen bonding or covalent interactions. Other forms of cohesion involve edge-sharing or paranemic interactions of double helices. A large number of individual species have been developed by this approach, including polyhedral catenanes, a variety of single-stranded knots, and Borromean rings. In addition to these static species, DNA-based nanomechanical devices have been produced that are ultimately targeted to lead to nanorobotics. Many of the key goals of structural DNA nanotechnology entail the use of periodic arrays. A variety of 2D DNA arrays have been produced with tunable features, such as patterns and cavities. DNA molecules have be used successfully in DNA-based computation as molecular representations of Wang tiles, whose self-assembly can be programmed to perform a calculation. About 4 years ago, on the fiftieth anniversary of the double helix, the area appeared to be at the cusp of a truly exciting explosion of applications; this was a correct assessment, and much progress has been made in the intervening period.     

Common fold in helix-hairpin-helix proteins

Helix–hairpin–helix (HhH) is a widespread motif involved in non-sequence-specific DNA binding. The majority of HhH motifs function as DNA-binding modules, however, some of them are used to mediate protein–protein interactions or have acquired enzymatic activity by incorporating catalytic residues (DNA glycosylases). From sequence and structural analysis of HhH-containing proteins we conclude that most HhH motifs are integrated as a part of a five-helical domain, termed (HhH)₂ domain here. It typically consists of two consecutive HhH motifs that are linked by a connector helix and displays pseudo-2-fold symmetry. (HhH)₂ domains show clear structural integrity and a conserved hydrophobic core composed of seven residues, one residue from each α-helix and each hairpin, and deserves recognition as a distinct protein fold. In addition to known HhH in the structures of RuvA, RadA, MutY and DNA-polymerases, we have detected new HhH motifs in sterile alpha motif and barrier-to-autointegration factor domains, the α-subunit of Escherichia coli RNA-polymerase, DNA-helicase PcrA and DNA glyco­s­y­lases. Statistically significant sequence similarity of HhH motifs and pronounced structural conservation argue for homology between (HhH)₂ domains in different protein families. Our analysis helps to clarify how non-symmetric protein motifs bind to the double helix of DNA through the formation of a pseudo-2-fold symmetric (HhH)₂ functional unit.     

Interaction of Two Peptide-Acridine Conjugates Containing the SPKK Peptide Motif with DNA and Chromatin

Abstract

The interaction between DNA and two peptide-acridine conjugates containing one (1) or two (2) moieties of the Ser-Pro-Lys-Lys (SPKK) minor groove-binding peptide motif has been studied by a combination of hydrodynamic, biochemical and spectroscopic methods including diffusion-enhanced luminescence energy transfer (DELET) measurements with a Tb(III) lanthanide chelate as donor. Viscometric titrations do not reveal any significant difference between the two hybrid molecules which both unwind (by about 15°) and extend the DNA similarly. DELET measurements show that the acridinyl chromophore of compounds 1 and 2 is much more accessible than that of a simple monointercalating drug such as acridine orange or ethidium. The accessibility factor increases proportionally with the peptide length, reflecting the extent of perturbation imposed upon the intercalating chromophore by the binding to DNA of the peptide moiety of the hybrids. Experiments with the osmium tetroxide- bispyridine reagent indicate that the two hybrid compounds both affect the local conformation of DNA rendering certain thymine residues conspicuously accessible to the probe. The drug-induced sites of hyperreactivity towards OsO₄ in DNA are very similar with the exception of a short run of three T residues which is attacked more strongly in the presence of tetrapeptide-acridine conjugate 1 than with the octapeptide-acridine conjugate 2. These results are fully in agreement with previous footprinting studies and support the view that a minimum of two SPKK motifs is required to mimic the AT-specific minor groove binding antibiotic netropsin. On the basis of the DNA-binding properties of these two peptide-acridine hybrids, we present DNA-binding models in which the acridinyl moiety of compound 1 protrudes slightly outside the double helix but remains more or less parallel to the plane of the base-pairs. In contrast, with compound 2, where the octapeptide SPKKSPKK is bound to the minor groove, we postulate that the chromophore lies only partially overlapped with the base pairs in the intercalation site and, in addition, the heterocyclic chromophore is significantly tilted with respect to the double helix axis. Electric linear dichroism and DELET measurements with chromatin reveal that the presence of histone proteins affects the intercalative binding of compound 2 while it has practically no effect on the binding of compound 1.     

Flexibility Of DNA In 2:1 Drug-DNA Complexes - Simultaneous Binding Of Two DAPI Molecules To DNA

Abstract

Simultaneous binding of two DAPI molecules in the minor groove of (dA)₁₅.(dT)₁₅ B-DNA helix has been simulated by molecular mechanics calculations. The energy minimised structure shows some novel features in relation to binding of DAPI molecules as well as the flexibility of the grooves of DNA helices. The minor groove of the helix expands locally considerably (to 15 Å) to accommodate the two DAPI molecules and is achieved by positive propeller twisting of base pairs at the binding site concomitant with small variations in the local nucleotide stereochemistry. The expansion also brings forth simultaneously a contraction in the width of the major groove spread over to a few phosphates. These findings demonstrate another facet of the flexible stereochemistry of DNA helices in which the local features are significantly altered without being propagated beyond a few base pairs, and with the rest of the regions retaining the normal structure. Both the DAPI molecules are engaged in specific hydrogen bonds with the bases and non specific interactions with phosphates. Stacking interactions of DAPI molecules between themselves as well as with sugar-phosphate backbone contribute to the stability of the complex. The studies provide a stereochemical support to the experimental findings that under high drug-DNA ratio DAPI could bind in the 2:1 ratio.     

A Self-Consistent Formulation for Analysis and Generation of Non-Uniform DNA Structures

Abstract

A fully self-consistent formulation is described here for the analysis and generation of base-pairs in non-uniform DNA structures, in terms of various local parameters. It is shown that the internal “wedge parameters” are mathematically related to the parameters describing the base-pair orientation with respect to an external helix axis. Hence any one set of three translation and three rotation parameters are necessary and sufficient to completely describe the relative orientation of the base-pairs comprising a step (or doublet). A general procedure is outlined for obtaining an average or global helix axis from the local helix axes for each step. A graphical representation of the local helix axes in the form of a polar plot is also shown and its application for estimating the curvature of oligonucleotide structures is illustrated, with examples of both A and B type structures.     

Left-Handed Intercalated DNA Double Helix: Rendezvous of Ethidium and Actinomycin D in the Z-Helical Conformation Space

Abstract

It is now very well recognized that the DNA double helix is conformationally pluralistic and that this flexibility is derived from internal motions due to backbone torsions. But what is less apparent is that such internal motions can occur in a correlated fashion and express themselves in a wide variety of structural motifs and phenomena. For example, flexibility inherent in the DNA molecule can lead to a family of Z-DNA, LZ1 and LZ2 being the two extremes and correlated internal motion can cause LZ1?LZ2 transition. More interestingly, such motions manifest themselves as breathing modes on the DNA lattice resulting in the sequence specific intercalation sites. Following a detailed stereochemical analyses we observed that the intercalation site for ethidium is located at the dCpdG sequence of the intercalated LZ1 helix (LZ1*) while that for actinomycin D is located at the dGpdC sequence of the intercalated LZ2 helix (LZ2*). From the stereochemistry of the drug binding we make experimentally testable predictions which are in fact supported by a few recent experimental studies. These studies also show that a left-handed intercalated B-DNA model is a viable intermediate in the Z to B transition which can hold the drug with binding energy comparable to that of the intercalated right-handed B-DNA.     

Melting of Cross-Linked DNA IV. Methods for Computer Modeling of Total Influence on DNA Melting of Monofunctional Adducts,Intrastrand and Interstrand Cross-Links Formed by Molecules of an Antitumor Drug

Abstract

A theoretical method is developed for calculation of melting curves of covalent complexes of DNA with antitumor drugs. The method takes into account all the types of chemical modifications of the double helix caused by platinum compounds and DNA alkylating agents: 1) monofunctional adducts bound to one nucleotide; 2) intrastrand cross-links which appear due to bidentate binding of a drug molecule to two nucleotides that are included into the same DNA strand; 3) interstrand cross-links caused by bidentate binding of a molecule to two nucleotides of different strands. The developed calculation method takes into account the following double helix alterations at sites of chemical modifications: 1) a change in stability of chemically modified base pairs and neighboring ones, that is caused by all the types of chemical modifications; 2) a change in the energy of boundaries between helical and melted regions at sites of chemical modification (local alteration of the factor of cooperativity of DNA melting), that is caused by all the types of chemical modifications, too; 3) a change in the loop entropy factor of melted regions that include interstrand cross-links; 4) the prohibition of divergence of DNA strands in completely melted DNA molecules, which is caused by interstrand cross-links only. General equations are derived, and three calculation methods are proposed to calculate DNA melting curves and the parameters that characterize the helix-coil transition.     

5-Deazaflavins: New Very Efficient DNA Photosensitisers,Synthesis of Oligonucleotide Conjugates

Abstract

In order to study electron transfer processes through the DNA double helix, we have synthesised a series of 5-deazaflavin derivatives 1–4 and demonstrated their ability to induce very efficiently 2′-deoxyguanosine and DNA oxidations by electron transfer from guanine. 5-Deazaflavin-oligonucleotide adducts were also synthesised for the study of electron transfers through double or triple helix formed with their complementary sequence.     

DNA Models for A,B, C and D Conformations Related to Fiber X-ray,Infrared and NMR Measurements

Abstract

A conformational analysis of the A, B, C and D DNA forms was made in order to establish molecular models presenting a good agreement with experimental data obtained from fiber X-ray, infrared linear dichroism and ³¹P NMR. The proposed models have been refined and do present good stereochemistry and optimized H-bond distances between bases associated with the Watson-Crick pairing. The DNA conformations proposed are a left handed double helix for the C form and right handed helices for A, B and D. Relations to conformational transitions between these forms are discussed.     

Recovery of the poisoned topoisomerase II for DNA religation: coordinated motion of the cleavage core revealed with the microsecond atomistic simulation

Type II topoisomerases resolve topological problems of DNA double helices by passing one duplex through the reversible double-stranded break they generated on another duplex. Despite the wealth of information in the cleaving operation, molecular understanding of the enzymatic DNA ligation remains elusive. Topoisomerase poisons are widely used in anti-cancer and anti-bacterial therapy and have been employed to entrap the intermediates of topoisomerase IIβ with religatable DNA substrate. We removed drug molecules from the structure and conducted molecular dynamics simulations to investigate the enzyme-mediated DNA religation. The drug-unbound intermediate displayed transitions toward the resealing-compliant configuration: closing distance between the cleaved DNA termini, B-to-A transformation of the double helix, and restoration of the metal-binding motif. By mapping the contact configurations and the correlated motions between enzyme and DNA, we identified the indispensable role of the linker preceding winged helix domain (WHD) in coordinating the movements of TOPRIM, the nucleotide-binding motifs, and the bound DNA substrate during gate closure. We observed a nearly vectorial transition in the recovery of the enzyme and identified the previously uncharacterized roles of Asn508 and Arg677 in DNA rejoining. Our findings delineate the dynamic mechanism of the DNA religation conducted by type II topoisomerases.     

A Model for Parallel Triple Helix Formation by Ree A: Single-Strand Association with a Homologous Duple via the Minor Groove

Abstract

The nucleoproteic filaments of RecA polymerized on single stranded DNA are able to integrate double stranded DNA in a coaxial arrangement (with DNA stretched by a factor 1.5), to recognize homologous sequences in the duplex and to perform strand exchange between the single stranded and double stranded molecules. While experimental results favor the hypothesis of an invasion of the minor groove of the duplex by the single strand, parallel minor groove triple helices have never been isolated or even modeled, the minor groove offering little space for a third strand to interact. Based on an internal coordinate modeling study, we show here that such a structure is perfectly conceivable when the two interacting oligomers are stretched by a factor 1.5, in order to open the minor groove of the duplex. The model helix presents characteristics that coincide with known experimental data on unwinding, base pair inclination and inter-proton distances. Moreover, we show that extension and unwinding stabilize the triple helix. New patterns of triplet interaction via the minor groove are presented.     

MPD and DNA Bending in Crystals and in Solution

Bending of 15 to 24° is observed within crystal structures ofB-DNA duplexes, is strongly sequence-dependent, and exhibits no correlation with the concentration of MPD (2-methyl-2,4-pentanediol) in the crystallizing solution. Two types of bends are observed: facultative bends or flexible hinges at junctions between regions of G·C and A·T base-pairs, and a persistent and almost obligatory bend at the center of the sequence R-G-C-Y. Only A-tracts are characteristically straight and unbent in every crystal structure examined to date. A detailed examination of normal vector plots for individual strands of a double helix provides an explanation, in terms of the stacking properties of guanine and adenine bases. The effect of high MPD concentrations, in both solution and crystal, is to decrease local bending somewhat without removing it altogether. MPD gel retardation experiments provide no basis for choosing among the three models that seek to explain macroscopic curvature of DNA by means of microscopic bending: junction bending, bent A-tracts, or bent general- sequence DNA. Crystallographic data on the straightness of A-tracts, the bendability of non-A sequences, and the identity of inclination angles in A-tract and non-A-tractB-DNA support only the general-sequence bending model. The pre-melting transition observed in A-tract DNA probably represents a relaxation of stiff adenine stacks to a flexible conformation more typical of general-sequence DNA.     

HELANAL: A Program to Characterize Helix Geometry in Proteins

Abstract

A detailed analysis of structural and position dependent characteristic features of helices will give a better understanding of the secondary structure formation in globular proteins. Here we describe an algorithm that quantifies the geometry of helices in proteins on the basis of their C^α atoms alone. The Fortran program HELANAL can extract the helices from the PDB files and then characterises the overall geometry of each helix as being linear, curved or kinked, in terms of its local structural features, viz. local helical twist and rise, virtual torsion angle, local helix origins and bending angles between successive local helix axes. Even helices with large radius of curvature are unambiguously identified as being linear or curved. The program can also be used to differentiate a kinked helix and other motifs, such as helix-loop-helix or a helix-turn-helix (with a single residue linker) with the help of local bending angles. In addition to these, the program can also be used to characterise the helix start and end as well as other types of secondary structures.     

FOUR-STRANDED DNA AS DETERMINED BY ELECTRON MICROSCOPY

Pneumococcus DNA, of weight-average molecular weight 1.6 million by light scattering, had a weight-average length of 4300 A by electron microscopy. Thus, the average mass per unit length was 370 molecular-weight units per A, or approximately two times that expected (208) for a Watson-Crick double helix. This corresponds to an average of 3.6 strands per molecule, which is close to that obtained by other methods. Morphologically, all the particles in the micrographs were relatively stiff, and had a cross-sectional height of 20 to 30 A. Some divided into two stiff branches of the same height, apparently double helical. Where the branches combined into one (minimally four-stranded) structure they apparently lay side by side in close association.     

Is DNA really a double helix?

Do the two chains of the DNA molecule coil round one another plectonemically ? If so, what is the approximate value of Lk (the linking number) for any closed, circular DNA molecule? Experiments using gel electrophoresis have shown that supercoiled DNA molecules usually migrate in a series of discrete bands. The only tenable explanation for this quantized behavior is that the molecules in one band all have the same value of Lk and that this value differs by unity from that of the adjacent bands. Various experiments in which circular DNA is unwound by known amounts show that (given this assumption) Lk for relaxed DNA is very roughly equal to $\text{N}\text{10}$ (where N is the number of base-pairs), as expected from the classical double helix.The original model for the double helix was right-handed. The experimental evidence for this feature is suggestive but not yet completely compelling.     

Index to Authors,Volume 3

Abstract

Results of calculations using various empirical potentials suggest that base pair buckling, which commonly occurs in DNA crystal structures, is sufficient to eliminate the steric clash at CpG steps in B-DNA, originating from the base pair propeller twisting. The buckling is formed by an inclination of cytosines while deviations of guanines from a plane perpendicular to the double helix axis are unfavorable. The buckling is accompanied by an increased vertical separation of the base pair centers but the buckled arrangement of base pairs is at least as stable as when the vertical separation is normal and buckle zero. In addition, room is created by the increased vertical separation for the bases to propeller twist as is observed in DNA crystal structures. Further stabilization of base stacking is introduced into the buckled base pair arrangement by roll opening the base pairs into the double helix minor groove. The roll may lead to the double helix bending and liberation of guanines from the strictly perpendicular orientation to the double helix axis. The liberated guanines further contribute to the base pair buckling and stacking improvement. This work also suggests a characteristic very stable DNA structure promoted by nucleotide sequences in which runs of purines follow runs of pyrimidine bases.