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.     

Spatial translational motions of base pairs in DNA molecules: Application of the extended matrix generator method

We have used the elementary generator matrices outlined in the preceding paper to examine the conformational plasticity of the nucleic acid double helix. Here we investigate kinked DNA structures made up of alternating B- and A-type helices and intrinsically curved duplexes perturbed by the intercalation of ligands. We model the B-to-A transition by the lateral translation of adjacent base pairs, and the intercalation of ligands by the vertical displacement of neighboring residues. We report a complete set of average configuration-dependent parameters, ranging from scalars (i.e., persistence lengths) to first- and second-order tensor parameters (i.e., average second moments of inertia), as well as approximations of the associated spatial distributions of the DNA and their angular correlations. The average structures of short chains (of lengths less than 100 base pairs) with local kinks or intrinsically curved sequences are essentially rigid rods. At the smallest chain lengths (10 base pairs), the kinked and curved chains exhibit similar average properties, although they are structurally perturbed compared to the standard B-DNA duplex. In contrast, at lengths of 200 base pairs, the curved and kinked chains are more compact on average and are located in a different space from the standard B- or A-DNA helix. While A-DNA is shorter and thicker than B-DNA in x-ray models, the long flexible A-DNA helix is thinner and more extended on average than its B-DNA counterpart because of more limited fluctuations in local structure. Curved polymers of 50 base pairs or longer also show significantly greater asymmetry than other DNAs (in terms of the distribution of base pairs with respect to the center of gravity of the chain). The intercalation of drugs in the curved DNA straightens and extends the smoothly deformed template. The dimensions of the average ellipsoidal boundaries defining the configurations of the intercalated polymers are roughly double those of the intrinsically curved chain. The altered proportions and orientations of these density functions reflect the changing shape and flexibility of the double helix. The calculations shed new light on the possible structural role of short A-DNA fragments in long B-type duplexes and also offer a model for understanding how GC-specific intercalative ligands can straighten naturally curved DNA. The mechanism is not immediately obvious from current models of DNA curvature, which attribute the bending of the chain to a perturbed structure in repeating tracts of A · T base pairs. © 1994 John Wiley & Sons, Inc.     

G . T base-pairs in a DNA helix: the crystal structure of d(G-G-G-G-T-C-C-C)

The synthetic deoxyoctanucleotide d(G-G-G-G-T-C-C-C) crystallizes as an A-type DNA double helix containing two adjacent G . T base-pair mismatches. The structure has been refined to an R-factor of 14% at 2.1 A resolution with 104 solvent molecules located. The two G . T mismatches adopt the "wobble" form of base-pairing. The mismatched bases are linked by a network of water molecules interacting with the exposed functional groups in both the major and minor grooves. The presence of two mispaired bases in the octamer has surprisingly little effect on the global structure of the helix or the backbone and glycosidic torsional angles. Base stacking around the mismatch is perturbed, but the central G-T step shows particularly good base overlap, which may contribute to the relatively high stability of this oligomer.     

The crystal structure of d(CCCCGGGG): a new A-form variant with an extended backbone conformation

The crystal structure of d(CCCCGGGG) has been determined at a resolution of 2.25 A. The oligomers crystallize as A-DNA duplexes occupying crystallographic two-fold axes. The backbone conformation is, in general, similar to that observed in previously reported crystal structures of A-DNA fragments, except for the central linkage, where it adopts an extended structure resulting from all trans conformation at the P-O5'-C5'-C4' bonds. This type of conformation facilitates interstrand stacking between the guanines at the C-G site. The local helix twist at this step is very small (25 degrees) compared to an overall average of 33.5 degrees. The unique structure of the C-G base-pair step, namely the extended backbone and the distinct stacking geometry, may be an important feature in the recognition mechanism between double-stranded DNA molecules and restriction endonucleases such as Msp I, which cuts the sequence CCGG very specifically with a rate unaffected by neighboring base pairs.     

Electric dichroism and bending amplitudes of DNA fragments according to a simple orientation function for weakly bent rods

The linear dichroism is calculated for DNA fragments in their thermal bending equilibrium. These calculations are given for relatively short fragments, where bent molecules can be described by an arc model. Using the measured value of 350 A for the persistence length, the limit dichroism (corresponding to complete alignment) decreases due to thermal bending, e.g., for a fragment with 100 base pairs to 80% of the value expected for straight molecules. Thermal bending should lead to a strong continuous decrease of the dichroism with increasing chain length, which is not observed, however, in electric dichroism experiments due to electric stretching. The influence of the electric field on the bending equilibrium is described by a contribution to the bending energy, which is calculated from the movement of charge equivalents against the potential gradient upon bending. The charge equivalents, which are assigned to the helix ends, are derived from the dipole moments causing the stationary degree of orientation. By this procedure the energy term inducing DNA stretching is given for induced, permanent, and saturating induced dipole models without introduction of any additional parameter. The stationary dichroism at a given electric field strength is then calculated according to an arc model by integration over all angles of orientation of helix axes or chords with respect to the field vector, and at each of these angles the contribution to the dichroism is calculated by integration over all helices with different degrees of bending. Orientation functions obtained by this procedure are fitted to dichroism data measured for various restriction fragments. Optimal fits are found for an induced dipole model with saturation of the polarizability. The difference between orientation functions with and without electric stretching is used to evaluate dichroism bending amplitudes. Both chain length and field strength dependence of bending amplitudes are consistent with experimental amplitudes derived from the dichroism decay in low salt buffers containing multivalent ions like Mg2+, spermine, or [CoNH3)6]3+. Bending amplitudes can be used to evaluate the persistence length from electrooptical data obtained for a single DNA restriction fragment. Bending and stretching effects are considerable already at relatively low chain length, and thus should not be neglected in any quantitative evaluation of experimental data.     

Conformational effects of DNA stretching

DNA is an extensible molecule, and an extended conformation of DNA is involved in some biological processes. We have examined the effect of elongation stress on the conformational properties of DNA base pairs by conformational analysis. The calculations show that stretching does significantly affect the conformational properties and flexibilities of base pairs. In particular, we have found that the propeller twist in base pairs reverses its sign upon stretching. The energy profile analysis indicates that electrostatic interactions make a major contribution to the stabilization of the positive-propeller-twist configuration in stretched DNA. This stretching also results in a monotonic decrease in the helical twist angle, tending to unwind the double helix. Fluctuations in most variables initially increase upon stretching, because of unstacking of base pairs, but then the fluctuations decrease as DNA is stretched further, owing to the formation of specific interactions between base pairs induced by the positive propeller twist. Thus, the stretching of DNA has particularly significant effects upon DNA flexibility. These changes in both the conformation and flexibility of base pairs probably have a role in functional interactions with proteins.     

Conformational Effects of DNA Stretching

DNA is an extensible molecule, and an extended conformation of DNA is involved in some biological processes. We have examined the effect of elongation stress on the conformational properties of DNA base pairs by conformational analysis. The calculations show that stretching does significantly affect the conformational properties and flexibilities of base pairs. In particular, we have found that the propeller twist in base pairs reverses its sign upon stretching. The energy profile analysis indicates that electrostatic interactions make a major contribution to the stabilization of the positive-propeller-twist configuration in stretched DNA. This stretching also results in a monotonic decrease in the helical twist angle, tending to unwind the double helix. Fluctuations in most variables initially increase upon stretching, because of unstacking of base pairs, but then the fluctuations decrease as DNA is stretched further, owing to the formation of specific interactions between base pairs induced by the positive propeller twist. Thus, the stretching of DNA has particularly significant effects upon DNA flexibility. These changes in both the conformation and flexibility of base pairs probably have a role in functional interactions with proteins.     

Molecular structure of an A-DNA decamer d(ACCGGCCGGT)

The molecular structure of the DNA decamer d(ACCGGCCGGT) has been solved and refined by single-crystal X-ray-diffraction analysis at 0.20 nm to a final R-factor of 18.0%. The decamer crystallizes as an A-DNA double helical fragment with unit-cell dimensions of a = b = 3.923 nm and c = 7.80 nm in the space group P6(1)22. The overall conformation of this A-DNA decamer is very similar to that of the fiber model for A-DNA which has a large average base-pair tilt and hence a wide and shallow minor groove. This structure is in contrast to that of several A-DNA octamers in which the molecules all have low base-pair-tilt angles (8-12 degrees) resulting in an appearance intermediate between B-DNA and A-DNA. The average helical parameters of this decamer are typical of A-DNA with 10.9 base pairs/turn of helix, an average helical twist angle of 33.1 degrees, and a base-pair-tilt angle of 18.2 degrees. However, the CpG step in this molecule has a low local-twist angle of 24.5 degrees, similar to that seen in other A-DNA oligomers, and therefore appears to be an intrinsic stacking pattern for this step. The molecules pack in the crystal using a recurring binding motif, namely, the terminal base pair of one helix abuts the surface of the shallow minor groove of another helix. In addition, the GC base pairs have large propeller-twist angles, unlike those found most other A-DNA structures.     

The Crystal Structure of d(CCCCGGGG): A New A-Form Variant with an Extended Backbone Conformation

Abstract

The crystal structure of d(CCCCGGGG) has been determined at a resolution of 2.25Å. The oligomers crystallize as A-DNA duplexes occupying crystallographic two-fold axes. The backbone conformation is, in general, similar to that observed in previously reported crystal structures of A-DNA fragments, except for the central linkage, where it adopts an extended structure resulting from all trans conformation at the P-05′-C5′-C4′ bonds. This type of conformation facilitates interstrand stacking between the guanines at the C-G site. The local helix twist at this step is very small (25°) compared to an overall average of 33.5°. The unique structure of the C-G base-pair step, namely the extended backbone and the distinct stacking geometry, may be an important feature in the recognition mechanism between double- stranded DNA molecules and restriction endonucleases such as Msp I, which cuts the sequence CCGG very specifically with a rate unaffected by neighboring base pairs.     

The reaction of aromatic peptides with double helical DNA. Quantitative characterisation of a two step reaction scheme

The binding of LysTrpLys and LysTyrLys to calf thymus DNA has been investigated by the field jump method using fluorescence detection. Two separate relaxation processes, clearly distinguished on the time scale and by opposite ampli- tudes, are observed for the binding of LysTrpLys to DNA with ~ 30000 base pairs. The concentration dependence of the relaxation time constants demonstrates a mechanism with a bimolecular step followed by a slow intramolecular transition with a forward rate of 6.4 X 103 s^?1 and an equilibrium constant of 11. Measurements at various degrees of peptide binding demonstrate that the binding mechanism associated with low binding rates is restricted to a rather low number of binding sites (roughly one site in 15 base pairs). The binding of LysTyrLys to the same DNA is not associated with relaxation pro- cesses of opposite amplitudes; nevertheless two processes could be identified and assigned to a two step mechanism corre- sponding to that observed in the case of LysTrpLys. In the presence of sonicated DNA both peptides show a single relaxa- tion process with characteristics similar to those observed for the slow process in the binding to high molecular DNA. The data indicate that the intramolecular step is faster for low than for high molecular DNA. These results suggest an assignment of the intramolecular step to an insertion of the aromatic residues into the DNA associated with bending of the helix. The increase in the rate of the intramolecular step with decreasing chain length of the DNA may then be explained by a higher flexibility of the double helix at lower chain lengths.     

Crystallographic study of one turn of G/C-rich B-DNA

The DNA decamer d(CCAGGCCTGG) has been studied by X-ray crystallography. At a nominal resolution of 1.6 A, the structure was refined to R = 16.9% using stereochemical restraints. The oligodeoxyribonucleotide forms a straight B-DNA double helix with crystallographic dyad symmetry and ten base-pairs per turn. In the crystal lattice, DNA fragments stack end-to-end along the c-axis to form continuous double helices. The overall helical structure and, notably, the groove dimensions of the decamer are more similar to standard, fiber diffraction-determined B-DNA than A-tract DNA. A unique stacking geometry is observed at the CA/TG base-pair step, where an increased rotation about the helix axis and a sliding motion of the base-pairs along their long axes leads to a superposition of the base rings with neighboring carbonyl and amino functions. Three-center (bifurcated) hydrogen bonds are possible at the CC/GG base-pair steps of the decamer. In their common sequence elements, d(CCAGGCCTGG) and the related G.A mismatch decamer d(CCAAGATTGG) show very similar three-dimensional structures, except that d(CCAGGCCTGG) appears to have a less regularly hydrated minor groove. The paucity of minor groove hydration in the center of the decamer may be a general feature of G/C-rich DNA and explain its relative instability in the B-form of DNA.     

Modelling extreme stretching of DNA.

Molecular modelling with Jumna is used to study extreme stretching of the DNA double helix. The results, which correlate well with recent nanomanipulation experiments, show how the double helix can be extended to twice its normal length before its base pairs break. Depending on the way the duplex is stretched two types of conformation can occur, either an unwound flat ribbon or a narrow fibre with negatively inclined base pairs. The energetics of both types of deformation are similar and existing structures show that at least the flat ribbon form can exist locally under biological conditions.     

Base pair buckling can eliminate the interstrand purine clash at the CpG steps in B-DNA caused by the base pair propeller twisting

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

Stretching single-stranded DNA: interplay of electrostatic, base-pairing, and base-pair stacking interactions.

Recent single-macromolecule observations revealed that the force/extension characteristics of single-stranded DNA (ssDNA) are closely related to solution ionic concentration and DNA sequence composition. To understand this, we studied the elastic property of ssDNA through the Monte Carlo implementation of a modified freely jointed chain (FJC), with electrostatic, base-pairing, and base-pair stacking interactions all incorporated. The simulated force-extension profiles for both random and designed sequences have attained quantitative agreements with the experimental data. In low-salt solution, electrostatic interaction dominates, and at low forces, the molecule can be more easily aligned than an unmodified FJC. In high-salt solution, secondary hairpin structure appears in ssDNA by the formation of base pairs between complementary bases, and external stretching causes a hairpin-coil structural transition, which is continuous for ssDNA made of random sequences. In designed sequences such as poly(dA-dT) and poly(dG-dC), the stacking potential between base pairs encourages the aggregation of base pairs into bulk hairpins and makes the hairpin-coil transition a discontinuous (first-order) process. The sensitivity of elongation to the base-pairing rule is also investigated. The comparison of modeling calculations and the experimental data suggests that the base pairing of single-stranded polynucleotide molecules tends to form a nested and independent planar hairpin structure rather than a random intersecting pattern.     

An effective mesoscopic model of double-stranded DNA

Watson and Crick’s epochal presentation of the double helix structure in 1953 has paved the way to intense exploration of DNA’s vital functions in cells. Also, recent advances of single molecule techniques have made it possible to probe structures and mechanics of constrained DNA at length scales ranging from nanometers to microns. There have been a number of atomistic scale quantum chemical calculations or molecular level simulations, but they are too computationally demanding or analytically unfeasible to describe the DNA conformation and mechanics at mesoscopic levels. At micron scales, on the other hand, the wormlike chain model has been very instrumental in describing analytically the DNA mechanics but lacks certain molecular details that are essential in describing the hybridization, nano-scale confinement, and local denaturation. To fill this fundamental gap, we present a workable and predictive mesoscopic model of double-stranded DNA where the nucleotides beads constitute the basic degrees of freedom. With the inter-strand stacking given by an interaction between diagonally opposed monomers, the model explains with analytical simplicity the helix formation and produces a generalized wormlike chain model with the concomitant large bending modulus given in terms of the helical structure and stiffness. It also explains how the helical conformation undergoes overstretch transition to the ladder-like conformation at a force plateau, in agreement with the experiment.