Sequence motifs and free energies of selected natural and non-natural nucleosome positioning DNA sequences.

Our laboratories recently completed SELEX experiments to isolate DNA sequences that most-strongly favor or disfavor nucleosome formation and positioning, from the entire mouse genome or from even more diverse pools of chemically synthetic random sequence DNA. Here we directly compare these selected natural and non-natural sequences. We find that the strongest natural positioning sequences have affinities for histone binding and nucleosome formation that are sixfold or more lower than those possessed by many of the selected non-natural sequences. We conclude that even the highest-affinity sequence regions of eukaryotic genomes are not evolved for the highest affinity or nucleosome positioning power. Fourier transform calculations on the selected natural sequences reveal a special significance for nucleosome positioning of a motif consisting of approximately 10 bp periodic placement of TA dinucleotide steps. Contributions to histone binding and nucleosome formation from periodic TA steps are more significant than those from other periodic steps such as AA (=TT), CC (=GG) and more important than those from the other YR steps (CA (=TG) and CG), which are reported to have greater conformational flexibility in protein-DNA complexes even than TA. We report the development of improved procedures for measuring the free energies of even stronger positioning sequences that may be isolated in the future, and show that when the favorable free energy of histone-DNA interactions becomes sufficiently large, measurements based on the widely used exchange method become unreliable.     

Principles of sequence-dependent flexure of DNA

The curvature of a bent rod may be defined in several different, but equivalent ways. The best way of describing the curvature of double-helical DNA is by an angle of turning per base step. Curvature comes mainly from the angle of roll between successive base-pairs, and this is defined as positive when the angle opens up on the minor groove side of the bases. DNA forms a plane curve if the roll angle values along the molecule alternate periodically between positive and negative, with a complete period equal to the helical repeat. It is known from studies of crystallized oligomers that the roll angles for particular dinucleotide steps have preferred values, or lie in preferred ranges of values. Therefore the formation of a plane curve will be easier with some base sequences of DNA than with others. We set up a computer algorithm for determining the ease with which DNA of given sequence will adopt a curved form. The algorithm has two different sets of constants: in model 1 the base step parameters come from an inspection of crystallized oligomers, and in model 2 data from a statistical survey of the incidence of dinucleotide steps in a large number of samples of chicken erythrocyte core DNA is incorporated. Both forms of the algorithm successfully locate the dyad of the nucleosome sequence (modulo 10) in a frog gene, and suggest strongly that sequence-dependent flexural properties of DNA play a part in the recognition of binding sites by nucleosome cores.     

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.     

Correlation between the flexibility and periodic dinucleotide patterns in yeast nucleosomal DNA sequences

Nucleosome formation and positioning, which play important roles in a number of biological processes, are thought to be related to the distinctive periodic dinucleotide patterns observed in the DNA sequence wrapped around the protein octamer. Previous research shows that flexibility is a key structural property of a nucleosomal DNA sequence. However, the relationship between the flexibility and the periodic dinucleotide patterns has received little attention in research in the past. In this study, we propose the use of three different models to measure the flexibility of yeast DNA sequences. Although the three models involve different parameters, they deliver consistent results showing that yeast nucleosomal DNA sequences are more flexible than non-nucleosomal ones. In contrast to random flexibility values along non-nucleosomal DNA sequences, the flexibility of nucleosomal DNA sequences shows a clear periodicity of 10.14 base pairs, which is consistent with the periodicity of dinucleotide distributions. We also demonstrate that there is a strong relationship between the peak positions of the flexibility and the dinucleotide frequencies. Correlation between the flexibility and the dinucleotide patterns of CA/TG, CG, GC, GG/CC, AG/CT, AC/GT and GA/TC are positive with an average value of 0.5946. The highest correlation is shown by CA/TG with a value of 0.7438 and the lowest correlation is shown by AA/TT with a value of −0.7424. The source codes and data sets are available for downloading on http://www.hy8.com/bioinformatics.htm.     

Cheminformatics methods for novel nanopore analysis of HIV DNA termini

BACKGROUND: Channel current feature extraction methods, using Hidden Markov Models (HMMs) have been designed for tracking individual-molecule conformational changes. This information is derived from observation of changes in ionic channel current blockade "signal" upon that molecule's interaction with (and occlusion of) a single nanometer-scale channel in a "nanopore detector". In effect, a nanopore detector transduces single molecule events into channel current blockades. HMM analysis tools described are used to help systematically explore DNA dinucleotide flexibility, with particular focus on HIV's highly conserved (and highly flexible/reactive) viral DNA termini. One of the most critical stages in HIV's attack is the binding between viral DNA and the retroviral integrase, which is influenced by the dynamic-coupling induced high flexibility of a CA/TG dinucleotide positioned precisely two base-pairs from the blunt terminus of the duplex viral DNA. This suggests the study of a family of such CA/TG dinucleotide molecules via nanopore measurement and cheminformatics analysis. RESULTS: HMMs are used for level identification on the current blockades, HMM/EM with boosted variance emissions are used for level projection pre-processing, and time-domain FSAs are used to parse the level-projected waveform for kinetic information. The observed state kinetics of the DNA hairpins containing the CA/TG dinucleotide provides clear evidence for HIV's selection of a peculiarly flexible/interactive DNA terminus.     

Structural properties of replication origins in yeast DNA sequences

Sequence-dependent DNA flexibility is an important structural property originating from the DNA 3D structure. In this paper, we investigate the DNA flexibility of the budding yeast (S. Cerevisiae) replication origins on a genome-wide scale using flexibility parameters from two different models, the trinucleotide and the tetranucleotide models. Based on analyzing average flexibility profiles of 270 replication origins, we find that yeast replication origins are significantly rigid compared with their surrounding genomic regions. To further understand the highly distinctive property of replication origins, we compare the flexibility patterns between yeast replication origins and promoters, and find that they both contain significantly rigid DNAs. Our results suggest that DNA flexibility is an important factor that helps proteins recognize and bind the target sites in order to initiate DNA replication. Inspired by the role of the rigid region in promoters, we speculate that the rigid replication origins may facilitate binding of proteins, including the origin recognition complex (ORC), Cdc6, Cdt1 and the MCM2-7 complex.     

Sequence-dependent DNA structure: dinucleotide conformational maps

We have used a computational model to calculate the potential energy surface for dinucleotide steps in double helical DNA as a function of the two principal degrees of freedom, slide and shift. By using a virtual bond to model the constraints imposed by the sugar-phosphate backbone, twist, roll, tilt and rise can be simultaneously optimised for any given values of slide and shift. Thus we have been able to construct complete conformational maps for all step types. For some steps, the maps agree well with experimental data from X-ray crystal structures, but other steps appear to be strongly perturbed by the effects of context (conformational coupling with the neighbouring steps). The optimised values of twist and roll show sequence-dependent variations consistent with the crystal structure data. The conformational maps allow us to construct adiabatic paths, and hence calculate the flexibility of each step with respect to slide and shift. Again the results agree well with the available experimental assignments of flexibility: YR steps, CA/TG and CG, are the most flexible and RR steps, such as AA, the least flexible.     

Sequence-dependent DNA structure: tetranucleotide conformational maps

A database of X-ray crystal structures of double helical DNA oligomers has been used to analyse the role of the sugar-phosphate backbone in coupling the conformational properties of neighbouring dinucleotide steps. The base step parameters which are most strongly coupled to the backbone degrees of freedom are slide and shift, and these are the two dinucleotide step parameters which show strong correlations along a sequence: the value of slide follows the values in the neighbouring steps, whereas shift tends to alternate. This conformational coupling is mediated by the shared furanose rings at the step junctions: a change in the value of slide causes a change in the mean value of the same strand 3' and 5'-chi torsion angle, and a change in the mean value of the 3' and 5' sugar pseudo-rotation phase angle, P; a change in the value of shift causes a difference between the same strand 3' and 5'-chi in A-DNA and a difference between the 3' and 5'-P in B-DNA. We have used a database of tetranucleotide X-ray crystal structures to parameterise a simple model for the coupling of slide and shift. Using this junction model together with our dinucleotide step potential energy maps described previously, we can in principle calculate the structure of any DNA oligomer. The parameterisation indicates that the rotational step parameters are accurate to within 5 degrees, and the translational step parameters are accurate to within 0.5 A. The model has been used to study the potential energy surfaces of all possible tetranucleotide sequences, and the calculations agree well with the experimental data from X-ray crystal structures. Some dinucleotide steps are context independent (AA/TT, AT and TA), because the conformational properties of all possible neighbouring steps are compatible. When the conformational properties of the neighbours are not compatible, the behaviour of a step cannot be understood at the dinucleotide level. Thus the conformations of CG, GC and GG/CC are all strongly context dependent. The remaining mixed sequence steps show weakly context-dependent behaviour. The approach allows the calculation of the relative stability and flexibility of tetranucleotide sequences, and the results indicate why TATA is used as an origin of replication. Clear predictions are made about sequences which have not yet been characterised crystallographically. In particular, poly(CCA).poly(TGG) is predicted to have an unusual structure which lies between the C and D-DNA polymorphs.     

Sequences that adopt non-B-DNA conformation in form V DNA as probed by enzymic methylation

pBR322 form V DNA is a highly torsionally strained molecule with a linking number of zero. We have used sequence-specific DNA methylases as probes for B-DNA in this molecule, exploiting the inability of methylases to methylate single-stranded DNA and Z-DNA, both of which are known to occur in form V DNA. Some sequences in form V DNA were shown to be totally in the B-form, others were totally in an altered, unmethylatable conformation, while still other sites appeared to exist partly in altered and partly in normal B-conformation. Some potential Z-forming sequences (alternating pyrimidine/purine) of less than seven base-pairs were not in the Z conformation in form V DNA, whereas others did adopt an altered structure, indicating a modulating influence of flanking sequences. Furthermore, regions of imperfect alternating pyrimidine/purine structure were sometimes capable of adopting an altered structure. In addition, some regions of altered structure had no apparent Z-forming sequences, nor were they in polypurine stretches, which have also been proposed to form left-handed DNA. These non-B-DNA conformations may represent novel left-handed helical structures or sequences that become single stranded under torsional strain. Long regions of either altered (unmethylatable) DNA or B-DNA were not always observed. In fact, one region showed three transitions between B-like DNA and altered structure within 26 base-pairs.     

Sequence-dependent dynamics of duplex DNA: the applicability of a dinucleotide model

The short-time (submicrosecond) bending dynamics of duplex DNA were measured to determine the effect of sequence on dynamics. All measurements were obtained from a single site on duplex DNA, using a single, site-specific modified base containing a rigidly tethered, electron paramagnetic resonance active spin probe. The observed dynamics are interpreted in terms of single-step sequence-dependent bending force constants, determined from the mean squared amplitude of bending relative to the end-to-end vector using the modified weakly bending rod model. The bending dynamics at a single site are a function of the sequence of the nucleotides constituting the duplex DNA. We developed and examined several dinucleotide-based models for flexibility. The models indicate that the dominant feature of the dynamics is best explained in terms of purine- and pyrimidine-type steps, although distinction is made among all 10 unique steps: It was found that purine-purine steps (which are the same as pyrimidine-pyrimidine steps) were near average in flexibility, but the pyrimidine-purine steps (5' to 3') were nearly twice as flexible, whereas purine-pyrimidine steps were more than half as flexible as average DNA. Therefore, the range of stepwise flexibility is approximately fourfold and is characterized by both the type of base pair step (pyrimidine/purine combination) and the identity of the bases within the pair (G, A, T, or C). All of the four models considered here underscore the complexity of the dependence of dynamics on DNA sequence with certain sequences not satisfactorily explainable in terms of any dinucleotide model. These findings provide a quantitative basis for interpreting the dynamics and kinetics of DNA-sequence-dependent biological processes, including protein recognition and chromatin packaging.     

Sequence dependence of DNA conformational flexibility

By using conformational free energy calculations, we have studied the sequence dependence of flexibility and its anisotropy along various conformational variables of DNA base pairs. The results show the AT base step to be very flexible along the twist coordinate. On the other hand, homonucleotide steps, GG(CC) and AA(TT), are among the most rigid sequences. For the roll motion that would correspond to a bend, the TA step is most flexible, while the GG(CC) step is least flexible. The flexibility of roll is quite anisotropic; the ratio of fluctuations toward the major and minor grooves is the largest for the GC step and the smallest for the AA(TT) and CG steps. Propeller twisting of base pairs is quite flexible, especially of A.T base pairs; propeller twist can reach 19 degrees by thermal fluctuation. We discuss the effect of electrostatic parameters, comparison with available experimental results, and biological relevance of these results.     

Unzipping mechanism of the double-stranded DNA unwinding by a hexameric helicase: quantitative analysis of the rate of the dsDNA unwinding, processivity and kinetic step-size of the Escherichia coli DnaB helicase using rapid quench-flow method

Kinetics of the double-stranded (ds) DNA unwinding by the Escherichia coli replicative helicase DnaB protein has been examined under single-turnover conditions using the chemical quench-flow technique. The unwinding reaction proceeds through an initial conformational transition followed by the unwinding catalytic steps and the release of the single-stranded (ss) DNA. Analyses of the reaction as a function of the number of base-pairs in the dsDNA reveal that the number of catalytic steps is not strictly proportional to the length of the dsDNA. As the helicase approaches the end of the substrate, the remaining approximately 11 bp of the DNA melts without catalytic participation of the enzyme. The kinetic step-size of the DnaB helicase, i.e. the number of the base-pairs unwound in a single catalytic step is only 1.4(+/- 0.2). The low value of the step-size indicates that the helicase unwinds a single base-pair in a single catalytic step. Thus, the DnaB helicase unzips the dsDNA in a reverse process to the zipping mechanism of the non-enzymatic double helix formation. The protein is a fast helicase that at 25 degrees C unwinds approximately 291 bp/s, much faster than previously thought, and the unwinding rate can be much higher at higher temperatures. However, the ATP-state of the enzyme has an increased dissociation rate, resulting in only a moderate unwinding processivity, P = 0.89(+/- 0.03), little dependent on the temperature. The conformational transition of the DnaB helicase-DNA complex, preceding the unwinding, is an intrinsic transition of the enzyme from the stationary conformation to the ATP-state of the helicase.     

Intrinsic bending and deformability at the T-A step of CCTTTAAAGG: a comparative analysis of T-A and A-T steps within A-tracts

Introduction of a T-A or pyrimidine-purine step into a straight and rigid A-tract can cause a positive roll deformation that kinks the DNA helix at that step. In CCTTTAAAGG, the central T-A step has an 8.6 degrees bend toward the major groove. We report the structural analysis of CCTTTAAAGG and a comparison with 25 other representative crystal structures from the NDB containing at least four consecutive A or T bases. On average, more local bending occurs at the disruptive T-A step (8.21 degrees ) than at an A-T step (5.71 degrees ). In addition, A-tracts containing an A-T step are more bent than are pure A-tracts, and hence A-A and A-T steps are not equivalent. All T-A steps examined exhibit positive roll, bending towards the major groove, while A-T steps display negative roll and bend slightly towards the minor groove. This illustrates how inherent negative and positive roll are, respectively, at A-T and T-A steps within A-tracts. T-A steps are more deformable, showing larger and more variable deformations of minor groove width, rise, cup, twist, and buckle. Standard deviations of twist, rise, and cup for T-A steps are 6.66 degrees, 0.55 A, and 15.90 degrees, versus 2.28 degrees, 0.21 A, and 2.99 degrees for A-T steps. Packing constraints determine which local values of these helical parameters an individual T-A step will adopt. For instance, with CCTTTAAAGG and three isomorphous structures, CGATTAATCG, CGATATATCG, and CGATCGATCG, crystal packing forces lead to a series of correlated changes: widened minor groove, large slide, low twist, and large rise. The difference in helical parameters between A-T steps lying within A-tracts, versus A-T steps within alternating AT sequences, demonstrates the importance of neighboring steps on the conformation of a given dinucleotide step.     

Sequence-dependent flexibility in promoter sequences

The non-neighbor interactions between base-pairs were taken into account to calculate the angular parameters (Omega, rho and tau) describing the orientation of successive base-pair planes and the translation parameters (D(y)) along the long axis of base-pair steps for 36 independent tetramers. A statistical mechanical model was proposed to predict the DNA flexibility that is mainly related to the thermal fluctuations at individual base-pair steps. The DNA flexibility can be described by the root-mean-square deviation of the end-to-end distance of DNA helical structure. The present model was then used to investigate the extreme flexible pattern in prokaryotic and eukaryotic promoter sequences. The results demonstrated several extreme flexible regions related to functionally important elements exist both in prokaryotic promoters and in eukaryotic promoters, DNA flexibility and AT content are highly correlated. The probabilities finding flexibility pattern in promoter sequences were also estimated statistically. The biological implications were discussed briefly.     

A-form conformational motifs in ligand-bound DNA structures

Recognition and biochemical processing of DNA requires that proteins and other ligands are able to distinguish their DNA binding sites from other parts of the molecule. In addition to the direct recognition elements embedded in the linear sequence of bases (i.e. hydrogen bonding sites), these molecular agents seemingly sense and/or induce an "indirect" conformational response in the DNA base-pairs that facilitates close intermolecular fitting. As part of an effort to decipher this sequence-dependent structural code, we have analyzed the extent of B-->A conformational conversion at individual base-pair steps in protein and drug-bound DNA crystal complexes. We take advantage of a novel structural parameter, the position of the phosphorus atom in the dimer reference frame, as well as other documented measures of local helical structure, e.g. torsion angles, base-pair step parameters. Our analysis pinpoints ligand-induced conformational changes that are difficult to detect from the global perspective used in other studies of DNA structure. The collective data provide new structural details on the conformational pathway connecting A and B-form DNA and illustrate how both proteins and drugs take advantage of the intrinsic conformational mechanics of the double helix. Significantly, the base-pair steps which exhibit pure A-DNA conformations in the crystal complexes follow the scale of A-forming tendencies exhibited by synthetic oligonucleotides in solution and the known polymorphism of synthetic DNA fibers. Moreover, most crystallographic examples of complete B-to-A deformations occur in complexes of DNA with enzymes that perform cutting or sealing operations at the (O3'-P) phosphodiester linkage. The B-->A transformation selectively exposes sugar-phosphate atoms, such as the 3'-oxygen atom, ordinarily buried within the chain backbone for enzymatic attack. The forced remodeling of DNA to the A-form also provides a mechanism for smoothly bending the double helix, for controlling the widths of the major and minor grooves, and for accessing the minor groove edges of individual base-pairs.