Sequence-structure relationships in DNA oligomers: a computational approach

Vibrational dynamics of DNA. I. Vibrational basis modes and couplings

Carrying out density functional theory calculations of four DNA bases, base derivatives, Watson-Crick (WC) base pairs, and multiple-layer base pair stacks, we studied vibrational dynamics of delocalized modes with frequency ranging from 1400 to 1800 cm(-1). These modes have been found to be highly sensitive to structure fluctuation and base pair conformation of DNA. By identifying eight fundamental basis modes, it is shown that the normal modes of base pairs and multilayer base pair stacks can be described by linear combinations of these vibrational basis modes. By using the Hessian matrix reconstruction method, vibrational coupling constants between the basis modes are determined for WC base pairs and multilayer systems and are found to be most strongly affected by the hydrogen bonding interaction between bases. It is also found that the propeller twist and buckle motions do not strongly affect vibrational couplings and basis mode frequencies. Numerically simulated IR spectra of guanine-cytosine and adenine-thymine bases pairs as well as of multilayer base pair stacks are presented and described in terms of coupled basis modes. It turns out that, due to the small interlayer base-base vibrational interactions, the IR absorption spectrum of multilayer base pair system does not strongly depend on the number of base pairs.

Structure of the alpha-homo-DNA:RNA duplex and the function of twist and slide to catalogue nucleic acid duplexes

High-resolution NMR studies of an alpha-homo-DNA:RNA duplex reveal the formation of a right-handed parallel-oriented helix. It differs significantly from a standard A- or B-type helix by a small twist value (26.2 degrees ), which leads to a helical pitch of 13.7 base pairs per helical turn, a negative inclination (-1.78 Angstrom) and a large x displacement (5.90 Angstrom). The rise (3.4 Angstrom) is similar to that found in B-DNA. The solution of this new helix structure has stimulated us to develop a mathematical and geometrical model based on slide and twist parameters to describe nucleic acid duplexes. All existing duplexes can be positioned within this landscape, which can be used to understand the helicalization process.

Genomic data analysis using DNA structure: an analysis of conserved nongenic sequences and ultraconserved elements

Recent comparative studies of the human and mouse genomes have revealed sets of conserved nongenic sequences (CNGs) and sets of ultraconserved elements (UCEs). Both sets of sequences, which exhibit extremely high levels of conservation, extend over hundreds of bases and have no known function. Since there is no detectable sequence homology between paralogous CNGs or UCEs in either of the species, an alignment-free technique is needed for their analysis. We have previously compiled a database of the structural properties of all 32,896 unique DNA octamers, including information on stability, the minimum energy conformation, and flexibility. We have used Fourier techniques to analyze the UCEs and CNGs in terms of their octamer structural properties, to reveal structural correlations which may indicate possible functions for some of these sequences.

Model identification for DNA sequence-structure relationships

We investigate the use of algebraic state-space models for the sequence dependent properties of DNA. By considering the DNA sequence as an input signal, rather than using an all atom physical model, computational efficiency is achieved. A challenge in deriving this type of model is obtaining its structure and estimating its parameters. Here we present two candidate model structures for the sequence dependent structural property Slide and a method of encoding the models so that a recursive least squares algorithm can be applied for parameter estimation. These models are based on the assumption that the value of Slide at a base-step is determined by the surrounding tetranucleotide sequence. The first model takes the four bases individually as inputs and has a median root mean square deviation of 0.90 A. The second model takes the four bases pairwise and has a median root mean square deviation of 0.88 A. These values indicate that the accuracy of these models is within the useful range for structure prediction. Performance is comparable to published predictions of a more physically derived model, at significantly less computational cost.

Structural DNA Profiles: Single Sequence Queries

Structural DNA profiles use the structural properties of the constituent octamers either to observe any characteristics of a single sequence that are unusual (a single sequence query) or to visualize a pattern common to a set of sequences (a multiple sequence query). They are an aid in understanding structural reasons for functional DNA activity. Profiles that answer single sequence queries are introduced and Profile Manager (a software application developed to automate profile generation) is presented. Two sequences that are similar by their nucleotide composition but are known to be very different by structure are analyzed, resulting in useful illustrations that agree with the experimental nuclear magnetic resonance structures.

Prediction of atomic structure from sequence for double helical DNA oligomers

DNA can adopt different conformations depending on the base sequence, solvent, electrolyte composition and concentration, pH, temperature, and interaction with proteins. Here we present a model for calculating the three-dimensional atomic structure of double-stranded DNA oligomers. A theoretical energy function is used for calculating the interactions within the base steps and an empirical backbone function is used to restrict the conformational space accessible to the bases and to account for the conformational coupling of neighboring steps in a sequence. Conformational searching on large structures or a large number of structures is possible, because each base step can be described by just two primary degrees of freedom (slide and shift). A genetic algorithm is used to search for low-energy structures in slide-shift space, and this allows very rapid optimization of DNA oligomers. The other base step parameters have been previously optimized for all possible slide-shift sequence combinations, and a heuristic algorithm is used to add the atomic details of the backbone conformation in the final step of the calculation. The structures obtained by this method are very similar to the corresponding X-ray crystal structures observed experimentally. The average RMSD is 2.24 Angstroms for a set of 20 oligomer structures. For 15 of these sequences, the X-ray crystal structure is the global energy minimum. The other 5 are bistable sequences that have B-form global energy minima but crystallize as A-DNA.

A structural similarity analysis of double-helical DNA

A database of the structural properties of all 32,896 unique DNA octamer sequences has been calculated, including information on stability, the minimum energy conformation and flexibility. The contents of the database have been analysed using a variety of Euclidean distance similarity measures. A global comparison of sequence similarity with structural similarity shows that the structural properties of DNA are much less diverse than the sequences, and that DNA sequence space is larger and more diverse than DNA structure space. Thus, there are many very different sequences that have very similar structural properties, and this may be useful for identifying DNA motifs that have similar functional properties that are not apparent from the sequences. On the other hand, there are also small numbers of almost identical sequences that have very different structural properties, and these could give rise to false-positives in methods used to identify function based on sequence alignment. A simple validation test demonstrates that structural similarity can differentiate between promoter and non-promoter DNA. Combining structural and sequence similarity improves promoter recall beyond that possible using either similarity measure alone, demonstrating that there is indeed information available in the structure of double-helical DNA that is not readily apparent from the sequence.

Understanding nucleic acids using synthetic chemistry

This Account describes work done in these laboratories that has used synthetic, physical organic, and biological chemistry to understand the roles played by the nucleobases, sugars, and phosphates of DNA in the molecular recognition processes central to genetics. The number of nucleobases has been increased from 4 to 12, generating an artificially expanded genetic information system. This system is used today in the clinic to monitor the levels of HIV and hepatitis C viruses in patients, helping to manage patient care. Work with uncharged phosphate replacements suggests that a repeating charge is a universal feature of genetic molecules operating in water and will be found in extraterrestrial life (if it is ever encountered). The use of ribose may reflect prebiotic processes in the presence of borate-containing minerals, which stabilize ribose formed from simple organic precursors. A new field, synthetic biology, is emerging on the basis of these experiments, where chemistry mimics biological processes as complicated as Darwinian evolution.

DNA basepair step deformability inferred from molecular dynamics simulations

The sequence-dependent DNA deformability at the basepair step level was investigated using large-scale atomic resolution molecular dynamics simulation of two 18-bp DNA oligomers: d(GCCTATAAACGCCTATAA) and d(CTAGGTGGATGACTCATT). From an analysis of the structural fluctuations, the harmonic potential energy functions for all 10 unique steps with respect to the six step parameters have been evaluated. In the case of roll, three distinct groups of steps have been identified: the flexible pyrimidine-purine (YR) steps, intermediate purine-purine (RR), and stiff purine-pyrimidine (RY). The YR steps appear to be the most flexible in tilt and partially in twist. Increasing stiffness from YR through RR to RY was observed for rise, whereas shift and slide lack simple trends. A proposed measure of the relative importance of couplings identifies the slide-rise, twist-roll, and twist-slide couplings to play a major role. The force constants obtained are of similar magnitudes to those based on a crystallographic ensemble. However, the current data have a less complicated and less pronounced sequence dependence. A correlation analysis reveals concerted motions of neighboring steps and thus exposes limitations in the dinucleotide model. The comparison of DNA deformability from this and other studies with recent quantum-chemical stacking energy calculations suggests poor correlation between the stacking and flexibility.

Natural selection and the emergence of a mutation phenotype: an update of the evolutionary synthesis considering mechanisms that affect genome variation

Most descriptions of evolution assume that all mutations are completely random with respect to their potential effects on survival. However, much like other phenotypic variations that affect the survival of the descendants, intrinsic variations in the probability, type, and location of genetic change can feel the pressure of natural selection. From site-specific recombination to changes in polymerase fidelity and repair of DNA damage, an organism's gene products affect what genetic changes occur in its genome. Through the action of natural selection on these gene products, potentially favorable mutations can become more probable than random. With examples from variation in bacterial surface proteins to the vertebrate immune response, it is clear that a great deal of genetic change is better than "random" with respect to its potential effect on survival. Indeed, some potentially useful mutations are so probable that they can be viewed as being encoded implicitly in the genome. An updated evolutionary theory includes emergence, under selective pressure, of genomic information that affects the probability of different classes of mutation, with consequences for genome survival.

DNA sequence-dependent deformability?insights from computer simulations

The article reviews some recent developments in studying DNA sequence-dependent deformability, with emphasis on computer modeling. After a brief outline of available experimental techniques, we proceed to computational methods and focus on atomic-resolution molecular dynamics (MD) simulations. A sequence-dependent local (base-pair step) force field inferred from MD is compared with force fields obtained by other techniques. Various methods for establishing global (flexible-rod) DNA elastic constants are reviewed, including an approach based on atomic resolution MD. The problem of defining the global deformation variables, as well as the question of anisotropy and nonlocal effects, are discussed. As an example, both local and global deformability calculations from atomic-resolution MD of EcoRI dodecamer are presented.

Sequence-dependent DNA structure: a database of octamer structural parameters

We have constructed the potential energy surfaces for all unique tetramers, hexamers and octamers in double helical DNA, as a function of the two principal degrees of freedom, slide and shift at the central step. From these potential energy maps, we have calculated a database of structural and flexibility properties for each of these sequences. These properties include: the values of each of the six step parameters (twist roll, tilt, rise, slide and shift), for each step of the sequence; flexibility measures for both decrease and increase in each property value from the minimum energy conformation for the central step; and the deviation from the path of a hypothetical straight octamer. In an analysis of structural change as a function of sequence length, we observe that almost all DNA tends to B-DNA and becomes less flexible. A more detailed analysis of octamer properties has allowed us to determine the structural preferences of particular sequence elements. GGC and GCC sequences tend to confer bistability, low stability and a predisposition to A-form DNA, whereas AA steps strongly prefer B-DNA and inhibit A-structures. There is no correlation between flexibility and intrinsic curvature, but bent DNA is less stable than straight. The most difficult deformation is undertwisting. The TA step stands out as the most flexible sequence element with respect to decreasing twist and increasing roll. However, as with the structural properties, this behavior is highly context-dependent and some TA steps are very straight.

Cross-modulation of the pKa of nucleobases in a single-stranded hexameric-RNA due to tandem electrostatic nearest-neighbor interactions

The pH titration studies (pH 6.7-12.1) in a series of dimeric, trimeric, tetrameric, pentameric, and hexameric oligo-RNA molecules [GpA (2a), GpC (3a), GpApC (5), GpA(1)pA(2)pC (6), GpA(1)pA(2)pA(3)pC (7), GpA(1)pA(2)pA(3)pA(4)pC (8)] have shown that the pK(a) of N(1)-H of 9-guaninyl could be measured not only from its own deltaH8G, but also from the aromatic marker protons of other constituent nucleobases. The relative chemical shift differences [Deltadelta((N)(-)(D))] between the protons in various nucleotide residues in the oligo-RNAs at the neutral (N) and deprotonated (D) states of the guanine moiety show that the generation of the 5'-(9-guanylate ion) in oligo-RNAs 2-8 reduces the stability of the stacked helical RNA conformation owing to the destabilizing anion(G(-))-pi/dipole(Im(delta)(-)) interaction. This destabilizing effect in the deprotonated RNA is, however, opposed by the electrostatically attractive atom-pisigma (major) as well as the anion(G(-))-pi/dipole(Py(delta)(+)) (minor) interactions. Our studies have demonstrated that the electrostatically repulsive anion(G(-))-pi/dipole(Im(delta)(-)) interaction propagates from the first to the third nucleobase quite strongly in the oligo-RNAs 6-8, causing destacking of the helix, and then its effect is gradually reduced, although it is clearly NMR detectable along the RNA chain. Thus, such specific generation of a charge at a single nucleobase moiety allows us to explore the relative strength of stacking within a single-stranded helix. The pK(a) of 5'-Gp residue from its own deltaH8G in the hexameric RNA 8 is found to be 9.76 +/- 0.01; it, however, varies from 9.65 +/- 0.01 to 10.5 +/- 0.07 along the RNA chain as measured from the other marker protons (H2, H8, H5, and H6) of 9-adeninyl and 1-cytosinyl residues. This nucleobase-dependent modulation of pK(a)s (DeltapK(a) +/- 0.9) of 9-guaninyl obtained from other nucleobases in the hexameric RNA 8 represents a difference of ca. 5.1 kJ mol(-)(1), which has been attributed to the variable strength of electrostatic interactions between the electron densities of the involved atoms in the offset stacked nucleobases as well as with that of the phosphates. The chemical implication of this variable pK(a) for guanin-9-yl deprotonation as obtained from all other marker protons of each nucleotide residue within a ssRNA molecule is that it enables us to experimentally understand the variation of the electronic microenvironment around each constituent nucleobase along the RNA chain in a stepwise manner with very high accuracy without having to make any assumption. This means that the pseudoaromaticity of neighboring 9-adeninyl and next-neighbor nucleobases within a polyanionic sugar-phosphate backbone of a ssRNA can vary from one case to another due to cross-modulation of an electronically coupled pi system by a neighboring nucleobase. This modulation may depend on the sequence context, spatial proximity of the negatively charged phosphates, as well as whether the offset stacking is ON or OFF. The net outcome of this electrostatic interaction between the neighbors is creation of new sequence-dependent hybrid nucleobases in an oligo- or polynucleotide whose properties are unlike the monomeric counterpart, which may have considerable biological implications.

The structure of DNA in the nucleosome core

The 1.9-A-resolution crystal structure of the nucleosome core particle containing 147 DNA base pairs reveals the conformation of nucleosomal DNA with unprecedented accuracy. The DNA structure is remarkably different from that in oligonucleotides and non-histone protein-DNA complexes. The DNA base-pair-step geometry has, overall, twice the curvature necessary to accommodate the DNA superhelical path in the nucleosome. DNA segments bent into the minor groove are either kinked or alternately shifted. The unusual DNA conformational parameters induced by the binding of histone protein have implications for sequence-dependent protein recognition and nucleosome positioning and mobility. Comparison of the 147-base-pair structure with two 146-base-pair structures reveals alterations in DNA twist that are evidently common in bulk chromatin, and which are of probable importance for chromatin fibre formation and chromatin remodelling.