Local DNA stretching mimics the distortion caused by the TATA box-binding protein

The origin of different bending stiffness between double-stranded RNA and DNA revealed by magnetic tweezers and simulations

The subtle differences in the chemical structures of double-stranded (ds) RNA and DNA lead to significant variations in their biological roles and medical implications, largely due to their distinct biophysical properties, such as bending stiffness. Although it is well known that A-form dsRNA is stiffer than B-form dsDNA under physiological salt conditions, the underlying cause of this difference remains unclear. In this study, we employ high-precision magnetic-tweezer experiments along with molecular dynamics simulations and reveal that the relative bending stiffness between dsRNA and dsDNA is primarily determined by the structure- and salt-concentration-dependent ion distribution around their helical structures. At near-physiological salt conditions, dsRNA shows a sparser ion distribution surrounding its phosphate groups compared to dsDNA, causing its greater stiffness. However, at very high monovalent salt concentrations, phosphate groups in both dsRNA and dsDNA become fully neutralized by excess ions, resulting in a similar intrinsic bending persistence length of approximately 39 nm. This similarity in intrinsic bending stiffness of dsRNA and dsDNA is coupled to the analogous fluctuations in their total groove widths and further coupled to the similar fluctuation of base-pair inclination, despite their distinct A-form and B-form helical structures.

Modeling the Homologous Recombination Process: Methods, Successes and Challenges

Homologous recombination (HR) is a fundamental process common to all species. HR aims to faithfully repair DNA double strand breaks. HR involves the formation of nucleoprotein filaments on DNA single strands (ssDNA) resected from the break. The nucleoprotein filaments search for homologous regions in the genome and promote strand exchange with the ssDNA homologous region in an unbroken copy of the genome. HR has been the object of intensive studies for decades. Because multi-scale dynamics is a fundamental aspect of this process, studying HR is highly challenging, both experimentally and using computational approaches. Nevertheless, knowledge has built up over the years and has recently progressed at an accelerated pace, borne by increasingly focused investigations using new techniques such as single molecule approaches. Linking this knowledge to the atomic structure of the nucleoprotein filament systems and the succession of unstable, transient intermediate steps that takes place during the HR process remains a challenge; modeling retains a very strong role in bridging the gap between structures that are stable enough to be observed and in exploring transition paths between these structures. However, working on ever-changing long filament systems submitted to kinetic processes is full of pitfalls. This review presents the modeling tools that are used in such studies, their possibilities and limitations, and reviews the advances in the knowledge of the HR process that have been obtained through modeling. Notably, we will emphasize how cooperative behavior in the HR nucleoprotein filament enables modeling to produce reliable information.

Atomistic Molecular Dynamics Simulations of DNA in Complex 3D Arrangements for Comparison with Lower Resolution Structural Experiments

Atomic-level computer simulations are a very useful tool for describing the structure and dynamics of complex biomolecules such as DNA and for providing detail at a resolution where experimental techniques cannot arrive. Molecular dynamics (MD) simulations of mechanically distorted DNA caused by agents like supercoiling and protein binding are computationally challenging due to the large size of the associated systems and timescales. However, nowadays they are achievable thanks to the efficient usage of GPU and to the improvements of continuum solvation models. This together with the concurrent improvements in the resolution of single-molecule experiments, such as atomic force microscopy (AFM), makes possible the convergence between the two. Here we present detailed protocols for doing so: for performing molecular dynamics (MD) simulations of DNA adopting complex three-dimensional arrangements and for comparing the outcome of the calculations with single-molecule experimental data with a lower resolution than atomic.

The emergence of sequence-dependent structural motifs in stretched, torsionally constrained DNA

The double-helical structure of DNA results from canonical base pairing and stacking interactions. However, variations from steady-state conformations resulting from mechanical perturbations in cells have physiological relevance but their dependence on sequence remains unclear. Here, we use molecular dynamics simulations showing sequence differences result in markedly different structural motifs upon physiological twisting and stretching. We simulate overextension on different sequences of DNA ((AA)12, (AT)12, (CC)12 and (CG)12) with supercoiling densities at 200 and 50 mM salt concentrations. We find that DNA denatures in the majority of stretching simulations, surprisingly including those with over-twisted DNA. GC-rich sequences are observed to be more stable than AT-rich ones, with the specific response dependent on the base pair order. Furthermore, we find that (AT)12 forms stable periodic structures with non-canonical hydrogen bonds in some regions and non-canonical stacking in others, whereas (CG)12 forms a stacking motif of four base pairs independent of supercoiling density. Our results demonstrate that 20-30% DNA extension is sufficient for breaking B-DNA around and significantly above cellular supercoiling, and that the DNA sequence is crucial for understanding structural changes under mechanical stress. Our findings have important implications for the activities of protein machinery interacting with DNA in all cells.

DNA-protein interaction: identification, prediction and data analysis

Life in living organisms is dependent on specific and purposeful interaction between other molecules. Such purposeful interactions make the various processes inside the cells and the bodies of living organisms possible. DNA-protein interactions, among all the types of interactions between different molecules, are of considerable importance. Currently, with the development of numerous experimental techniques, diverse methods are convenient for recognition and investigating such interactions. While the traditional experimental techniques to identify DNA-protein complexes are time-consuming and are unsuitable for genome-scale studies, the current high throughput approaches are more efficient in determining such interaction at a large-scale, but they are clearly too costly to be practice for daily applications. Hence, according to the availability of much information related to different biological sequences and clearing different dimensions of conditions in which such interactions are formed, with the developments related to the computer, mathematics, and statistics motivate scientists to develop bioinformatics tools for prediction the interaction site(s). Until now, there has been much progress in this field. In this review, the factors and conditions governing the interaction and the laboratory techniques for examining such interactions are addressed. In addition, developed bioinformatics tools are introduced and compared for this reason and, in the end, several suggestions are offered for the promotion of such tools in prediction with much more precision.

Integrating multi-scale data on homologous recombination into a new recognition mechanism based on simulations of the RecA-ssDNA/dsDNA structure

RecA protein is the prototypical recombinase. Members of the recombinase family can accurately repair double strand breaks in DNA. They also provide crucial links between pairs of sister chromatids in eukaryotic meiosis. A very broad outline of how these proteins align homologous sequences and promote DNA strand exchange has long been known, as are the crystal structures of the RecA-DNA pre- and postsynaptic complexes; however, little is known about the homology searching conformations and the details of how DNA in bacterial genomes is rapidly searched until homologous alignment is achieved. By integrating a physical model of recognition to new modeling work based on docking exploration and molecular dynamics simulation, we present a detailed structure/function model of homology recognition that reconciles extremely quick searching with the efficient and stringent formation of stable strand exchange products and which is consistent with a vast body of previously unexplained experimental results.

Molecular modeling of nucleic Acid structure: electrostatics and solvation

This unit presents an overview of computer simulation techniques as applied to nucleic acid systems, ranging from simple in vacuo molecular modeling techniques to more complete all-atom molecular dynamics treatments that include an explicit representation of the environment. The third in a series of four units, this unit focuses on critical issues in solvation and the treatment of electrostatics. UNITS 7.5 & 7.8 introduced the modeling of nucleic acid structure at the molecular level. This included a discussion of how to generate an initial model, how to evaluate the utility or reliability of a given model, and ultimately how to manipulate this model to better understand its structure, dynamics, and interactions. Subject to an appropriate representation of the energy, such as a specifically parameterized empirical force field, the techniques of minimization and Monte Carlo simulation, as well as molecular dynamics (MD) methods, were introduced as a way of sampling conformational space for a better understanding of the relevance of a given model. This discussion highlighted the major limitations with modeling in general. When sampling conformational space effectively, difficult issues are encountered, such as multiple minima or conformational sampling problems, and accurately representing the underlying energy of interaction. In order to provide a realistic model of the underlying energetics for nucleic acids in their native environments, it is crucial to include some representation of solvation (by water) and also to properly treat the electrostatic interactions. These subjects are discussed in detail in this unit.

Automatic workflow for the classification of local DNA conformations

Background

A growing number of crystal and NMR structures reveals a considerable structural polymorphism of DNA architecture going well beyond the usual image of a double helical molecule. DNA is highly variable with dinucleotide steps exhibiting a substantial flexibility in a sequence-dependent manner. An analysis of the conformational space of the DNA backbone and the enhancement of our understanding of the conformational dependencies in DNA are therefore important for full comprehension of DNA structural polymorphism.

Results

A detailed classification of local DNA conformations based on the technique of Fourier averaging was published in our previous work. However, this procedure requires a considerable amount of manual work. To overcome this limitation we developed an automatic classification method consisting of the combination of supervised and unsupervised approaches. A proposed workflow is composed of k-NN method followed by a non-hierarchical single-pass clustering algorithm. We applied this workflow to analyze 816 X-ray and 664 NMR DNA structures released till February 2013. We identified and annotated six new conformers, and we assigned four of these conformers to two structurally important DNA families: guanine quadruplexes and Holliday (four-way) junctions. We also compared populations of the assigned conformers in the dataset of X-ray and NMR structures.

Conclusions

In the present work we developed a machine learning workflow for the automatic classification of dinucleotide conformations. Dinucleotides with unassigned conformations can be either classified into one of already known 24 classes or they can be flagged as unclassifiable. The proposed machine learning workflow permits identification of new classes among so far unclassifiable data, and we identified and annotated six new conformations in the X-ray structures released since our previous analysis. The results illustrate the utility of machine learning approaches in the classification of local DNA conformations.

Tight associations between transcription promoter type and epigenetic variation in histone positioning and modification

Background

Transcription promoters are fundamental genomic cis-elements controlling gene expression. They can be classified into two types by the degree of imprecision of their transcription start sites: peak promoters, which initiate transcription from a narrow genomic region; and broad promoters, which initiate transcription from a wide-ranging region. Eukaryotic transcription initiation is suggested to be associated with the genomic positions and modifications of nucleosomes. For instance, it has been recently shown that histone with H3K9 acetylation (H3K9ac) is more likely to be distributed around broad promoters rather than peak promoters; it can thus be inferred that there is an association between histone H3K9 and promoter architecture.

Results

Here, we performed a systematic analysis of transcription promoters and gene expression, as well as of epigenetic histone behaviors, including genomic position, stability within the chromatin, and several modifications. We found that, in humans, broad promoters, but not peak promoters, generally had significant associations with nucleosome positioning and modification. Specifically, around broad promoters histones were highly distributed and aligned in an orderly fashion. This feature was more evident with histones that were methylated or acetylated; moreover, the nucleosome positions around the broad promoters were more stable than those around the peak ones. More strikingly, the overall expression levels of genes associated with broad promoters (but not peak promoters) with modified histones were significantly higher than the levels of genes associated with broad promoters with unmodified histones.

Conclusion

These results shed light on how epigenetic regulatory networks of histone modifications are associated with promoter architecture.

Modeling the early stage of DNA sequence recognition within RecA nucleoprotein filaments

Homologous recombination is a fundamental process enabling the repair of double-strand breaks with a high degree of fidelity. In prokaryotes, it is carried out by RecA nucleofilaments formed on single-stranded DNA (ssDNA). These filaments incorporate genomic sequences that are homologous to the ssDNA and exchange the homologous strands. Due to the highly dynamic character of this process and its rapid propagation along the filament, the sequence recognition and strand exchange mechanism remains unknown at the structural level. The recently published structure of the RecA/DNA filament active for recombination (Chen et al., Mechanism of homologous recombination from the RecA-ssDNA/dsDNA structure, Nature 2008, 453, 489) provides a starting point for new exploration of the system. Here, we investigate the possible geometries of association of the early encounter complex between RecA/ssDNA filament and double-stranded DNA (dsDNA). Due to the huge size of the system and its dense packing, we use a reduced representation for protein and DNA together with state-of-the-art molecular modeling methods, including systematic docking and virtual reality simulations. The results indicate that it is possible for the double-stranded DNA to access the RecA-bound ssDNA while initially retaining its Watson–Crick pairing. They emphasize the importance of RecA L2 loop mobility for both recognition and strand exchange.

Overstretching of a 30 bp DNA duplex studied with steered molecular dynamics simulation: effects of structural defects on structure and force-extension relation

Single-molecule experiments on polymeric DNA show that the molecule can be overstretched at nearly constant force by about 70% beyond its relaxed contour length. In this publication we use steered molecular dynamics (MD) simulation to study the effect of structural defects on force-extension curves and structures at high elongation in a 30 base pair duplex pulled by its torsionally unconstrained 5' -5' ends. The defect-free duplex shows a plateau in the force-extension curve at 120 pN in which large segments with inclined and paired bases ("S-DNA") near both ends of the duplex coexist with a central B-type segment separated from the former by small denaturation bubbles. In the presence of a base mismatch or a nick, force-extension curves are very similar to the ones of the defect-free duplex. For the duplex with a base mismatch, S-type segments with highly inclined base pairs are not observed; rather, the overstretched duplex consists of B-type segments separated by denaturation bubbles. The nicked duplex evolves, via a two-step transition, into a two-domain structure characterized by a large S-type segment coexisting with several short S-type segments which are separated by short denaturation bubbles. Our results suggest that in the presence of nicks the force-extension curve of highly elongated duplex DNA might reflect locally highly inhomogeneous stretching.

Deforming DNA: from physics to biology

The DNA double helix has become a modern icon which symbolizes our understanding of the molecular basis of life. It is less widely recognized that the double helix proposed by Watson and Crick more than half a century ago is a remarkably adaptable molecule that can undergo major conformational rearrangements without being irreversibly damaged. Indeed, DNA deformation is an intrinsic feature of many of the biological processes in which it is involved. Over the last two decades, single-molecule experiments coupled with molecular modeling have transformed our understanding of DNA flexibility, while the accumulation of high-resolution structures of DNA-protein complexes have demonstrated how organisms can exploit this property as a useful feature for preserving, reading, replicating, and packaging the genetic message. In this Minireview we summarize the information now available on the extreme--and the less extreme--deformations of the double helix.

Ets-1 p51 and p42 isoforms differentially modulate Stromelysin-1 promoter according to induced DNA bend orientation

The Stromelysin-1 gene promoter contains a palindrome of two Ets-binding sites (EBS) that bind the p51 and p42 isoforms of the human Ets-1-transcription factor. A previous study established that full gene transactivation is associated with a ternary complex consisting of two p51 bound to the two EBS on the promoter. p42, only able to bind one of the two EBS, induces only very weak activity. Here, we investigate the mechanism by which the Stromelysin-1 promoter discriminates between p51 and p42. The differential stoichiometry of the two Ets-1 isoforms arises from the Stromelysin-1 EBS palindrome. The ternary complex requires the presence of two inhibitory domains flanking the DNA-binding domain and the ability to form an intramolecular autoinhibition module. Most importantly, the p51-ternary and the p42-binary complexes induce DNA curvatures with opposite orientations. These results establish that differential DNA bending, via p51 and p42 differential binding, is correlated with the Stromelysin-1 promoter activation process.

Protein mechanics: a route from structure to function

In order to better understand the mechanical properties of proteins, we have developed simulation tools which enable these properties to be analysed on a residue-by-residue basis. Although these calculations are relatively expensive with all-atom protein models, good results can be obtained much faster using coarse-grained approaches. The results show that proteins are surprisingly heterogeneous from a mechanical point of view and that functionally important residues often exhibit unusual mechanical behaviour. This finding offers a novel means for detecting functional sites and also potentially provides a route for understanding the links between structure and function in more general terms.

NMR analysis of duplex d(CGCGATCGCG)2 modified by Λ- and Δ-[Ru(bpy)2(m-GHK)]Cl2 and DNA photocleavage study

The interaction of the diastereomeric complexes Lambda-[Ru(bpy)2(m-GHK)]Cl2 and Delta-[Ru(bpy)2(m-GHK)]Cl2 (bpy is 2,2'-bipyridine, GHK is glycine-L-histidine-L-lysine) with the deoxynucleotide duplex d(5'-CGCGATCGCG)2 was studied by means of 1H NMR spectroscopy. At a Delta-isomer to DNA ratio of 1:1, significant shifts for the metal complex are observed, whereas there is negligible effect on the oligonucleotide protons and only one intermolecular nuclear Overhauser effect (NOE) is present at the 2D nuclear Overhauser enhancement spectroscopy spectrum. The 1Eta NMR spectrum at ratio 2:1 is characterized by a slight shift for the Delta-isomer's bpy aromatic protons as well as significant shifts for the decanucleotide G4 H1' and Eta2'', A5 H2, G10 H1', T6 NH and G2 NH protons. Furthermore, at ratio 2:1, 11 intermolecular NOEs are observed. The majority of the NOEs involve the sugar Eta2' and Eta2'' protons sited in the major groove of the decanucleotide. Increasing the Delta-isomer to d(CGCGATCGCG)2 ratio to 5:1 results in noteworthy spectral changes. The Delta-isomer's proton shifts are reduced, whereas significant shifts are observed for the decanucleotide protons, especially the sugar protons, as well as for the exchangeable protons. Interaction is characterized by the presence of only one intermolecular NOE. Furthermore, there is significant broadening of the imino proton signals as the ratio of the Delta-isomer to DNuAlpha increases, which is attributed to the opening of the two strands of the duplex. The Lambda-isomer, on the other hand, approaches the minor groove of the oligonucleotide and interacts only weakly, possibly by electrostatic interactions. Photocleavage studies were also conducted with the plasmid pUC19 and a 158-bp restriction fragment, showing that both diastereomers cleave DNA with similar efficiency, attacking mainly the guanines of the sequence probably by generating active oxygen species.

Structure, dynamics, and elasticity of free 16s rRNA helix 44 studied by molecular dynamics simulations

Molecular dynamics (MD) simulations were employed to investigate the structure, dynamics, and local base-pair step deformability of the free 16S ribosomal helix 44 from Thermus thermophilus and of a canonical A-RNA double helix. While helix 44 is bent in the crystal structure of the small ribosomal subunit, the simulated helix 44 is intrinsically straight. It shows, however, substantial instantaneous bends that are isotropic. The spontaneous motions seen in simulations achieve large degrees of bending seen in the X-ray structure and would be entirely sufficient to allow the dynamics of the upper part of helix 44 evidenced by cryo-electron microscopic studies. Analysis of local base-pair step deformability reveals a patch of flexible steps in the upper part of helix 44 and in the area proximal to the bulge bases, suggesting that the upper part of helix 44 has enhanced flexibility. The simulations identify two conformational substates of the second bulge area (bottom part of the helix) with distinct base pairing. In agreement with nuclear magnetic resonance (NMR) and X-ray studies, a flipped out conformational substate of conserved 1492A is seen in the first bulge area. Molecular dynamics (MD) simulations reveal a number of reversible alpha-gamma backbone flips that correspond to transitions between two known A-RNA backbone families. The flipped substates do not cumulate along the trajectory and lead to a modest transient reduction of helical twist with no significant influence on the overall geometry of the duplexes. Despite their considerable flexibility, the simulated structures are very stable with no indication of substantial force field inaccuracies.

Minor groove deformability of DNA: a molecular dynamics free energy simulation study

The conformational deformability of nucleic acids can influence their function and recognition by proteins. A class of DNA binding proteins including the TATA box binding protein binds to the DNA minor groove, resulting in an opening of the minor groove and DNA bending toward the major groove. Explicit solvent molecular dynamics simulations in combination with the umbrella sampling approach have been performed to investigate the molecular mechanism of DNA minor groove deformations and the indirect energetic contribution to protein binding. As a reaction coordinate, the distance between backbone segments on opposite strands was used. The resulting deformed structures showed close agreement with experimental DNA structures in complex with minor groove-binding proteins. The calculated free energy of minor groove deformation was approximately 4-6 kcal mol(-1) in the case of a central TATATA sequence. A smaller equilibrium minor groove width and more restricted minor groove mobility was found for the central AAATTT and also a significantly ( approximately 2 times) larger free energy change for opening the minor groove. The helical parameter analysis of trajectories indicates that an easier partial unstacking of a central TA versus AT basepair step is a likely reason for the larger groove flexibility of the central TATATA case.

Probing protein mechanics: residue-level properties and their use in defining domains

It is becoming clear that, in addition to structural properties, the mechanical properties of proteins can play an important role in their biological activity. It nevertheless remains difficult to probe these properties experimentally. Whereas single-molecule experiments give access to overall mechanical behavior, notably the impact of end-to-end stretching, it is currently impossible to directly obtain data on more local properties. We propose a theoretical method for probing the mechanical properties of protein structures at the single-amino acid level. This approach can be applied to both all-atom and simplified protein representations. The probing leads to force constants for local deformations and to deformation vectors indicating the paths of least mechanical resistance. It also reveals the mechanical coupling that exists between residues. Results obtained for a variety of proteins show that the calculated force constants vary over a wide range. An analysis of the induced deformations provides information that is distinct from that obtained with measures of atomic fluctuations and is more easily linked to residue-level properties than normal mode analyses or dynamic trajectories. It is also shown that the mechanical information obtained by residue-level probing opens a new route for defining so-called dynamical domains within protein structures.

The relative flexibility of B-DNA and A-RNA duplexes: database analysis

An extensive analysis of structural databases is carried out to investigate the relative flexibility of B-DNA and A-RNA duplexes in crystal form. Our results show that the general anisotropic concept of flexibility is not very useful to compare the deformability of B-DNA and A-RNA duplexes, since the flexibility patterns of B-DNA and A-RNA are quite different. In other words, 'flexibility' is a dangerous word for describing macromolecules, unless it is clearly defined. A few soft essential movements explain most of the natural flexibility of A-RNA, whereas many are necessary for B-DNA. Essential movements occurring in naked B-DNAs are identical to those necessary to deform DNA in DNA-protein complexes, which suggest that evolution has designed DNA-protein complexes so that B-DNA is deformed according to its natural tendency. DNA is generally more flexible, but for some distortions A-RNA is easier to deform. Local stiffness constants obtained for naked B-DNAs and DNA complexes are very close, demonstrating that global distortions in DNA necessary for binding to proteins are the result of the addition of small concerted deformations at the base-pair level. Finally, it is worth noting that in general the picture of the relative deformability of A-RNA and DNA derived from database analysis agrees very well with that derived from state of the art molecular dynamics (MD) simulations.

Elucidating a key intermediate in homologous DNA strand exchange: structural characterization of the RecA-triple-stranded DNA complex using fluorescence resonance energy transfer

The RecA protein of Escherichia coli plays essential roles in homologous recombination and restarting stalled DNA replication forks. In vitro, the protein mediates DNA strand exchange between single-stranded (ssDNA) and homologous double-stranded DNA (dsDNA) molecules that serves as a model system for the in vivo processes. To date, no high-resolution structure of the key intermediate, comprised of three DNA strands simultaneously bound to a RecA filament (RecA-tsDNA complex), has been reported. We present a systematic characterization of the helical geometries of the three DNA strands of the RecA-tsDNA complex using fluorescence resonance energy transfer (FRET) under physiologically relevant solution conditions. FRET donor and acceptor dyes were used to label different DNA strands, and the interfluorophore distances were inferred from energy transfer efficiencies measured as a function of the base-pair separation between the two dyes. The energy transfer efficiencies were first measured on a control RecA-dsDNA complex, and the calculated helical parameters (h approximately 5 A, Omega(h) approximately 20 degrees ) were consistent with structural conclusions derived from electron microscopy (EM) and other classic biochemical methods. Measurements of the helical parameters for the RecA-tsDNA complex revealed that all three DNA strands adopt extended and unwound conformations similar to those of RecA-bound dsDNA. The structural data are consistent with the hypothesis that this complex is a late, post-strand-exchange intermediate with the outgoing strand shifted by about three base-pairs with respect to its registry with the incoming and complementary strands. Furthermore, the bases of the incoming and complementary strands are displaced away from the helix axis toward the minor groove of the heteroduplex, and the bases of the outgoing strand lie in the major groove of the heteroduplex. We present a model for the strand exchange intermediate in which homologous contacts preceding strand exchange arise in the minor groove of the substrate dsDNA.

The stretched DNA geometry of recombination and repair nucleoprotein filaments

The RecA protein of Escherichia coli plays essential roles in homologous recombination and restarting stalled DNA replication forks. In vitro, the protein mediates DNA strand exchange between single-stranded (ssDNA) and homologous double-stranded DNA (dsDNA) molecules that serves as a model system for the in vivo processes. To date, no high-resolution structure of the key intermediate, comprised of three DNA strands simultaneously bound to a RecA filament (RecA x tsDNA complex), has been elucidated by classical methods. Here we review the systematic characterization of the helical geometries of the three DNA strands of the RecA x tsDNA complex using fluorescence resonance energy transfer (FRET) under physiologically relevant solution conditions. Measurements of the helical parameters for the RecA x tsDNA complex are consistent with the hypothesis that this complex is a late, poststrand-exchange intermediate with the outgoing strand shifted by about three base pairs with respect to its registry with the incoming and complementary strands. All three strands in the RecA x tsDNA complex adopt extended and unwound conformations similar to those of RecA-bound ssDNA and dsDNA.

Molecular modeling of nucleic acid structure: electrostatics and solvation

This unit presents an overview of computer simulation techniques as applied to nucleic acid systems, ranging from simple in vacuo molecular modeling techniques to more complete all-atom molecular dynamics treatments that include an explicit representation of the environment. The third in a series of four units, this unit focuses on critical issues in solvation and the treatment of electrostatics.

Dynamic simulations of 13 TATA variants refine kinetic hypotheses of sequence/activity relationships 1 1Edited by B. Honig

The fundamental relationship between DNA sequence/deformability and biological function has attracted numerous experimental and theoretical studies. A classic prototype system used for such studies in eukaryotes is the complex between the TATA element transcriptional regulator and the TATA-box binding protein (TBP). The recent crystallographic study by Burley and co-workers demonstrated the remarkable structural similarity contrasted to different transcriptional activity of 11 TBP/DNA complexes in which the DNAs differed by single base-pairs. By simulating these TATA variants and two other single base-pair variants that were not crystallizable, we uncover sequence-dependent structural, energetic, and flexibility properties that tailor TATA elements to TBP interactions, complementing many previous studies by refining kinetic hypotheses on sequence/activity correlations. The factors that combine to produce favorable elements for TBP activity include overall flexibility; minor groove widening, as well as roll, rise, and shift increases at the ends of the TATA element; untwisting within the TATA element accompanied by large roll at the TATA element ends; and relatively low maximal water densities around the DNA. These features accompany the severe deformation induced by the minor-groove binding protein, which kinks the TATA element at the ends and displaces local water molecules to form stabilizing hydrophobic contacts. Interestingly, the preferred bending direction itself is not a significant predictor of activity disposition, although certain variants (such as wild-type AdMLP, 5'-TATA4G-3', and inactive A29, 5'-TA6G-3') exhibit large preferred bends in directions consistent with their activity or inactivity (major groove and minor groove bends, respectively). These structural, flexibility, and hydration preferences, identified here and connected to a new crystallographic study of a larger group of DNA variants than reported to date, highlight the profound influence of single base-pair DNA variations on DNA motion. Our refined kinetic hypothesis suggests the functional implications of these motions in a kinetic model of TATA/TBP recognition, inviting further theoretical and experimental research.

Modeling multi-component protein-DNA complexes: the role of bending and dimerization in the complex of p53 dimers with DNA

We used molecular modeling to study the optimal conformation of the complex between two p53 DNA-binding domain monomers and a 12 base-pair target DNA sequence. The complex was constructed using experimental data on the monomer binding conformation and a new approach to deform the target DNA sequence. Combined with an internal/helicoidal coordinate model of DNA, this approach enables us to bend the target sequence in a controlled way while respecting the contacts formed with each p53 monomer. The results show that the dimeric complex favors DNA bending towards the major groove at the dimer junction by a value close to experimental findings. In contrast to inferences from earlier models, the calculation of key contributions to the free energy of the complexes indicates a determinant role for DNA in the formation of the complex with the dimer of the p53 DNA-binding domains.

Twisting and stretching single DNA molecules

The elastic properties of DNA are essential for its biological function. They control its bending and twisting as well as the induction of structural modifications in the molecule. These can affect its interaction with the cell machinery. The response of a single DNA molecule to a mechanical stress can be precisely determined in single-molecule experiments which give access to an accurate measurement of the elastic parameters of DNA.

Optimization of nucleic acid sequences

Base sequence influences the structure, mechanics, dynamics, and interactions of nucleic acids. However, studying all possible sequences for a given fragment leads to a number of base combinations that increases exponentially with length. We present here a novel methodology based on a multi-copy approach enabling us to determine which base sequence favors a given structural change or interaction via a single energy minimization. This methodology, termed ADAPT, has been implemented starting from the JUMNA molecular mechanics program by adding special nucleotides, "lexides," containing all four bases, whose contribution to the energy of the system is weighted by continuously variable coefficients. We illustrate the application of this approach in the case of double-stranded DNA by determining the optimal sequences satisfying structural (B-Z transition), mechanical (intrinsic curvature), and interaction (ligand-binding) properties.

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.

DNA structure control by polycationic species: polyamine, cobalt ammines, and di-metallo transition metal chelates

Many polycationic species bind to DNA and induce structural changes. The work reported here is the first phase of a program whose long-term aim is to create a class of simple and inexpensive sequence-selective compounds that will enable enhanced DNA structure control for a wide range of applications. Three classes of molecule have been included in this work: the polyamine spermine (charge: 4(+)) and spermidine (charge: 3(+)) (which are known to induce a wide range of DNA conformational changes but whose binding modes are still not well understood); cobalt (III) amine transition metal complexes as potential polyamine mimics and [Fe(H(2)O)(6)](3+); and the first member of a new class of di-metallo tris-chelated cylinders of helical structure (charge 4(+)). Temperature-dependent absorption, circular dichroism, linear dichroism, gel electrophoresis, and molecular modeling data are presented. The cobalt amines prove to be effective polyamine mimics, although their binding appears to be restricted to backbone and major groove. All the ligands stabilize the DNA, but the 4(+) di-iron tris-chelate does so comparatively weakly and seems to have a preference for single-stranded DNA. All the molecules studied bend the DNA, with the di-iron tris-chelate having a particularly dramatic effect even at very low drug load.

Binding mechanisms of TATA box-binding proteins: DNA kinking is stabilized by specific hydrogen bonds

One of the common mechanisms of DNA bending by minor groove-binding proteins is the insertion of protein side chains between basepair steps, exemplified in TBP (TATA box-binding protein)/DNA complexes. At the central basepair step of the TATA box TBP produces a noticeable decrease in twist and an increase in roll, while engaging in hydrogen bonds with the bases and sugars. This suggests a mechanism for the stabilization of DNA kinks that was explored here with ab initio quantum mechanical calculations and molecular dynamics/potential of mean force calculations. The hydrogen bonds are found to contribute the energy necessary to drive the conformational transition at the central basepair step. The Asn, Thr, and Gly residues involved in hydrogen bonding to the DNA bases and sugar oxygens form a relatively rigid motif in TBP. The interaction of this motif with DNA is found to be responsible for inducing the untwisting and rolling of the central basepair step. Notably, direct readout is shown not to be capable of discriminating between AA and AT steps, as the strength of the hydrogen bonds between TBP and the DNA are the same for both sequences. Rather, the calculated free energy cost for an equivalent conformational transition is found to be sequence-dependent, and is calculated to be higher for AA steps than for AT steps.

TATA box DNA deformation with and without the TATA box-binding protein

DNA ring closure methods have been applied to TATA box DNA and its complex with the TATA box-binding protein (TBP). The J factors for cyclization (effective concentrations of one DNA end about the other) have been measured using cyclization kinetics, with and without bound TBP, for 18 DNA constructs containing the adenovirus major late promoter TATA box (TATAAAAG) separated by a variable helical phasing adapter from sequence-induced A-tract DNA bends. Six phasing lengths were used at three overall DNA lengths each. Cyclization kinetics were also measured in the absence of protein for the same set of molecules bearing a mutant TATA box (TACAAAAG). The results suggest that the TATA box DNA itself is strongly bent and anisotropically flexible, in a direction opposite to the bend induced by TBP, and that the mutant TACA box is much less bent/flexible. The bending and flexibility of the free DNA may govern the energetics of recognition of different DNA sequences by TBP, and the intrinsic bend may act to repress transcription complex assembly in the absence of TBP. The cyclization kinetics of TBP-DNA complexes in solution predict a geometry generally consistent with crystal structures, which show dramatic bending and unwinding. The novel observation of TBP-induced topoisomers suggests that this minicircle approach is able to distinguish TBP-induced unwinding from writhe (these cancel out in larger DNA), and this in turn suggests that changes in supercoiling in small topological domains can control TBP binding.

DNA stretching and compression: large-scale simulations of double helical structures

Computer-simulated elongation and compression of A - and B -DNA structures beyond the range of thermal fluctuations provide new insights into high energy "activated" forms of DNA implicated in biochemical processes, such as recombination and transcription. All-atom potential energy studies of regular poly(dG).poly(dC) and poly(dA).poly(dT) double helices, stretched from compressed states of 2.0 A per base-pair step to highly extended forms of 7.0 A per residue, uncover four different hyperfamilies of right-handed structures that differ in mutual base-pair orientation and sugar-phosphate backbone conformation. The optimized structures embrace all currently known right-handed forms of double-helical DNA identified in single crystals as well as non-canonical forms, such as the original "Watson-Crick" duplex with trans conformations about the P-O5' and C5'-C4' backbone bonds. The lowest energy minima correspond to canonical A and B -form duplexes. The calculations further reveal a number of unusual helical conformations that are energetically disfavored under equilibrium conditions but become favored when DNA is highly stretched or compressed. The variation of potential energy versus stretching provides a detailed picture of dramatic conformational changes that accompany the transitions between various families of double-helical forms. In particular, the interchanges between extended canonical and non-canonical states are reminiscent of the cooperative transitions identified by direct stretching experiments. The large-scale, concerted changes in base-pair inclination, brought about by changes in backbone and glycosyl torsion angles, could easily give rise to the observed sharp increase in force required to stretch single DNA molecules more than 1.6-1.65 times their canonical extension. Our extended duplexes also help to tie together a number of previously known structural features of the RecA-DNA complex and offer a self-consistent stereochemical model for the single-stranded/duplex DNA recognition brought in register by recombination proteins. The compression of model duplexes, by contrast, yields non-canonical structures resembling the deformed steps in crystal complexes of DNA with the TATA-box binding protein (TBP). The crystalline TBP-bound DNA steps follow the calculated compression-elongation pattern of an unusual "vertical" duplex with base planes highly inclined with respect to the helical axis, exposed into the minor groove, and accordingly accessible for recognition.Significantly, the double helix can be stretched by a factor of two and compressed roughly in half before its computed internal energy rises sharply. The energy profiles show that DNA extension-compression is related not only to the variation of base-pair Rise but also to concerted changes of Twist, Roll, and Slide. We suggest that the high energy "activated" forms calculated here are critical for DNA processing, e.g. nucleo-protein recognition, DNA/RNA synthesis, and strand exchange.

Unraveling proteins: a molecular mechanics study

An internal coordinate molecular mechanics study of unfolding peptide chains by external stretching has been carried out to predict the type of force spectra that may be expected from single-molecule manipulation experiments currently being prepared. Rather than modeling the stretching of a given protein, we have looked at the behavior of simple secondary structure elements (alpha-helix, beta-ribbon, and interacting alpha-helices) to estimate the magnitude of the forces involved in their unfolding or separation and the dependence of these forces on the way pulling is carried out as well as on the length of the structural elements. The results point to a hierarchy of forces covering a surprisingly large range and to important orientational effects in the response to external stress.

Modeling DNA deformations induced by minor groove binding proteins

Molecular modeling is used to demonstrate that the major structural deformations of DNA caused by four different minor groove binding proteins, TBP, SRY, LEF-1, and PurR, can all be mimicked by stretching the double helix between two 3'-phosphate groups flanking the binding region. This deformation reproduces the widening of the minor groove and the overall bending and unwinding of DNA caused by protein binding. It also reproduces the principal kinks associated with partially intercalated amino acid side chains, observed with such interactions. In addition, when protein binding involves a local transition to an A-like conformation, phosphate neutralization, via the formation of protein-DNA salt bridges, appears to favor the resulting deformation.

Modeling the mechanics of a DNA oligomer

DNA stretching and strand separation have been studied by molecular mechanics using an oligomer which has been the subject of nanomanipulation experiments (Noy et al., Chem. Biol. 4, 519, 1997). Adiabatic mapping of conformational energy carried out as a function of stretching leads to force/extension curves in good correlation with the experimental results. Other types of deformation are also modeled and compared with the experimental results obtained on polymeric DNA. The results highlight overall similarities, but point to thermodynamic differences and also to local base sequence effects which can be expected to play an important role at the level of biologically induced structural deformations.

A model for parallel triple helix formation by RecA: single-single association with a homologous duplex via the minor groove

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.

Selective binding of the TATA box-binding protein to the TATA box-containing promoter: analysis of structural and energetic factors

We report the results of an energy-based exploration of the components of selective recognition of the TATA box-binding protein (TBP) to a TATA box sequence that includes 1) the interaction between the hydrophobic Leu, Pro, and Phe residues of TBP with the TA, AT, AA, TT, and CG steps, by ab initio quantum mechanical calculations; and 2) the free energy penalty, calculated from molecular dynamics/potential of mean force simulations, for the conformational transition from A-DNA and B-DNA into the TA-DNA form of DNA observed in a complex with TBP. The GTAT, GATT, GAAT, and GTTT tetramers were explored. The results show that 1) the discrimination of TA, AT, AA, TT, or CG steps by TBP cannot rest on their interaction with the inserting Phe side chains; 2) the steric clash between the bulky and hydrophobic Pro and Leu residues and the protruding -NH2 group of guanine is responsible for the observed selectivity against any Gua-containing basepair; 3) the Pro and Leu residues cannot selectively discriminate among TA, AT, AA, or TT steps; and 4) the calculated energy required to achieve the TA-DNA conformation of DNA that is observed in the complex with TBP appears to be a key determinant for the observed selectivity against the AT, AA, and TT steps. The simulations also indicate that only the TA step can form a very efficient interbase hydrogen bond network in the TA-DNA conformation. Such an energetically stabilizing network is not achievable in the AA and TT steps. While it is viable in the AT step, structural constraints render the hydrogen bonding network energetically ineffective there.

Inherent DNA curvature and flexibility correlate with TATA box functionality

Four 1.5 ns molecular dynamics (MD) simulations were performed on the d(GCTATAAAAGGG).d(CCCTTTTATAGC) double helix dodecamer bearing the Adenovirus major late promoter TATA element and three iso-composition mutants for which physical and biochemical data are available from the same laboratory. Three of these DNA sequences experimentally induce tight binding with the TATA box binding protein (TBP) and induce high transcription rates; the other DNA sequence induces much lower TBP binding and transcription. The x-ray crystal structures have previously shown that the duplex DNA in DNA-TBP complexes are highly bent. We performed and analyzed MD simulations for these four DNAs, whose experimental structures are not available, in order to address the issue of whether inherent DNA structure and flexibility play a role in establishing these observed preferences. A comparison of the experimental and simulated results demonstrated that DNA duplex sequence-dependent curvature and flexibility play a significant role in TBP recognition, binding, and transcriptional activation.

RecA binding to a single double-stranded DNA molecule: a possible role of DNA conformational fluctuations

Most genetic regulatory mechanisms involve protein-DNA interactions. In these processes, the classical Watson-Crick DNA structure sometimes is distorted severely, which in turn enables the precise recognition of the specific sites by the protein. Despite its key importance, very little is known about such deformation processes. To address this general question, we have studied a model system, namely, RecA binding to double-stranded DNA. Results from micromanipulation experiments indicate that RecA binds strongly to stretched DNA; based on this observation, we propose that spontaneous thermal stretching fluctuations may play a role in the binding of RecA to DNA. This has fundamental implications for the protein-DNA binding mechanism, which must therefore rely in part on a combination of flexibility and thermal fluctuations of the DNA structure. We also show that this mechanism is sequence sensitive. Theoretical simulations support this interpretation of our experimental results, and it is argued that this is of broad relevance to DNA-protein interactions.

DNA sequence-dependent deformability deduced from protein-DNA crystal complexes

The deformability of double helical DNA is critical for its packaging in the cell, recognition by other molecules, and transient opening during biochemically important processes. Here, a complete set of sequence-dependent empirical energy functions suitable for describing such behavior is extracted from the fluctuations and correlations of structural parameters in DNA-protein crystal complexes. These elastic functions provide useful stereochemical measures of the local base step movements operative in sequence-specific recognition and protein-induced deformations. In particular, the pyrimidine-purine dimers stand out as the most variable steps in the DNA-protein complexes, apparently acting as flexible "hinges" fitting the duplex to the protein surface. In addition to the angular parameters widely used to describe DNA deformations (i.e., the bend and twist angles), the translational parameters describing the displacements of base pairs along and across the helical axis are analyzed. The observed correlations of base pair bending and shearing motions are important for nonplanar folding of DNA in nucleosomes and other nucleoprotein complexes. The knowledge-based energies also offer realistic three-dimensional models for the study of long DNA polymers at the global level, incorporating structural features beyond the scope of conventional elastic rod treatments and adding a new dimension to literal analyses of genomic sequences.

Excess counterion binding and ionic stability of kinked and branched DNA

We compute the excess number of counterions associated with kinked and branched DNA, and the ionic stabilities of these structures as a function of chain length and both sodium and magnesium salt concentration, using numerical counterion condensation theory. The DNA structures are modeled as two or more finite lines of phosphate charges radiating from the kink or junction center. The number of excess counterions around the (40-90 degrees) kinked duplex is very small (at most four). The geometries of large three- and four-way DNA junctions (with > 50 base pairs per arm) in solutions containing low to moderate NaCl concentrations, by contrast, accumulate a substantial number of excess sodium ions (> 20) but no more than 15 magnesium counterions. The excess number of counterions surrounding the kinked linear chain and the branched DNA structures either remains invariant or increases with chain length, tending to reach a plateau value. Open configurations, such as the planar Y-shaped three-way junction (with three 120 degrees inter-arm angles) and the 90 degrees cross-shaped four-way junction, are ionically more stable than compact geometries, such as pyramidal three-way junctions and X-shaped four-way junctions, over the entire range of salt concentration considered (10(-5)-10(-1) M NaCl or MgCl2). The ionic stabilities of the compact forms increase with increasing salt concentration and become comparable to those of the extended geometries at high salt (especially when magnesium is the supporting salt).

Sequence-dependent dynamics of TATA-Box binding sites

We have carried out two nanosecond-length molecular dynamics simulations on a DNA oligomer, d(GCGTAAAAAAAACGC)2, which contains a weak binding site for the TATA-box binding protein. An analysis of the resulting trajectories shows that this oligomer behaves differently from a related oligomer [d(GCGTATATAAAACGC)2] studied earlier using the same protocol (Flatters, D., M. Young, D. L. Beveridge, and R. Lavery. 1997. Conformational properties of the TATA-box binding sequence of DNA. J. Biomol. Struct. & Dyn. 14:757-765), and which contains a strong binding site for the same protein. The two basepair mutations that relate these oligomers lead to significant changes in time-averaged structure and in dynamic behavior, which extend over entire length of the oligomer and appear to be compatible with the experimentally observed decrease of binding and functional activity. These results suggest that molecular dynamics simulations, taking into account explicit solvent and counterions, and avoiding the truncation of electrostatic interactions, are a powerful tool for investigating the indirect aspects of protein-nucleic acid recognition.

Interaction of human SRY protein with DNA: a molecular dynamics study

Molecular dynamics simulations have been conducted to study the interaction of human sex-determining region Y (hSRY) protein with DNA. For this purpose, simulations of the hSRY high mobility group (HMG) domain (hSRY-HMG) with and without its DNA target site, a DNA octamer, and the DNA octamer alone have been carried out, employing the NMR solution structure of hSRY-HMG-DNA complex as a starting model. Analyses of the simulation results demonstrated that the interaction between hSRY and DNA was hydrophobic, just a few hydrogen bonds and only one water molecule as hydrogen-bonding bridge were observed at the protein-DNA interface. These two hydrophobic cores in the hSRY-HMG domain were the physical basis of hSRY-HMG-DNA specific interaction. They not only maintained the stability of the complex, but also primarily caused the DNA deformation. The salt bridges formed between the positive-charged residues of hSRY and phosphate groups of DNA made the phosphate electroneutral, which was advantageous for the deformation of DNA and the formation of a stable complex. We predicted the structure of hSRY-HMG domain in the free state and found that both hSRY and DNA changed their conformations to achieve greater complementarity of geometries and properties during the binding process; that is, the protein increased the angle between its long and short arms to accommodate the DNA, and the DNA became bent severely to adapt to the protein, although the conformational change of DNA was more severe than that of the hSRY-HMG domain. The sequence specificity and the role of residue Met9 are also discussed.

DNA bending by small, mobile multivalent cations

We propose a purely electrostatic mechanism by which small, mobile, multivalent cations can induce DNA bending. A multivalent cation binds at the entrance to the B-DNA major groove, between the two phosphate strands, electrostatically repelling sodium counterions from the neighboring phosphates. The unscreened phosphates on both strands are strongly attracted to the groove-bound cation. This leads to groove closure, accompanied by DNA bending toward the cationic ligand. We explicitly treat the dynamic character of the cation-DNA interaction using an adiabatic approximation, noting that DNA bending is much slower than the diffusion of nonspecifically bound, mobile cations. We make semiquantitative estimates of the free energy components of bending-electrostatic (with a sigmoidal distance-dependent dielectric function), elastic, and entropic cation localization-and find that the equilibrium state is bent B-DNA stabilized with a self-localized cation. This is a bending polaron, formation of which should be critically dependent on the strength of electrostatic interaction and the concentration of highly mobile cations available for self-localization. We predict that the resultant bend will be large (approximately 20-40 degrees), smooth (because it is spread over 6 bp), and infrequent. The stability of such a bend can be variable, from transient to highly stable (static) bending, observable with standard curvature-measuring techniques. We further predict that this bending mechanism will have an unusual sequence dependence: sequences with less binding specificity will be more bent, unless the specific binding site is in the major groove.

Progressive DNA bending is made possible by gradual changes in the torsion angle of the glycosyl bond

Structural comparisons have led to the suggestion that the conformational rearrangement that would be required to change A-DNA into the TA-DNA form of DNA observed in the complex with the TATA box binding protein (TBP) could be completed by modifying only the value of the glycosyl bond chi by approximately 45 degrees. The lack of a high number of crystal structures of this type makes it difficult to conclude whether a smooth transition from A-DNA to TA-DNA can occur without disrupting at any point either the Watson-Crick base pairing or the A-DNA conformation of the backbone. To explore the possibility of such a smooth transition, constrained molecular dynamics simulations were carried out for the double-stranded dodecamer d(GGTATATAAAAC), in which a transition from A-DNA to TA-DNA was induced by modifying only the chi angle values. The results demonstrate the feasibility of a continuous path in the A-DNA to TA-DNA transition. Varying extents of DNA curvature are also attainable, by maintaining the A-DNA backbone structure and Watson-Crick hydrogen bonding while changing the chi angle value smoothly from that in A-DNA to one corresponding to B-DNA.

Use of a 3D structure data base for understanding sequence-dependent conformational aspects of DNA

The roll-twist-slide correlation in the DNA crystal structures that are collected in the Nucleic Acid Data Base is analyzed in order to obtain a general understanding of the effects of the nucleotide sequence on the 3D structure of a dinucleotide step. It is concluded that the differences between the pyrimidine bases and the purine bases in terms of their physical shapes are the major factors that determine the stereochemical characteristics of the steps through base to backbone and base to base interactions. The characteristics are further modulated by the differences between the A:T and G:C base-pairs, which can be explained by enhancement of the purine-pyrimidine asymmetry in the A:T base-pair.