The stretched intermediate model of B-Z DNA transition

There have been numerous attempts to describe the mechanism of B-Z transition. Our simulations based on the stochastic difference equation with length algorithm show that a short DNA oligomer will tend to unwind and overstretch during the transition. The overstretching of DNA is then understood from the Zhou, Zhang, and Ou-Yang model. Unlike the Harvey model, the stretched intermediate model does not pose any steric dilemma; we are able to show that the chain sense reversal progresses spontaneously using the stretched intermediate model. A nonlinear DNA model is used to describe the origins and mechanism of base rotation in the stretched intermediate state of DNA. We also propose an experiment that can verify the existence of a stretched intermediate state during B-Z transition, thus opening up fresh grounds for experimentation in this field.     

A model of the strain-induced B-Z transition

The B-to-Z transition in supercoiled circular DNA is modeled as a strain-induced nonlinear excitation process. Using a model, in which DNA is regarded as a chain of units with a bistable energy function along the twisting coordinate together with a harmonic inter-unit interaction, we show that a Z region and the accompanying two B-Z junctions of finite width appear naturally as a solution of nonlinear equations, when the strain exceeds a critical value. We examine the B-Z transition behaviour as a function of twist under various situations. We also analyse available experimental results on B-Z transition in supercoiled plasmid with G-C insertions by this mechanistic model in order to estimate the magnitude of model parameters. The energy barrier of the B-Z transition is estimated to be of the order of 1 kcal/mole per base pair. The analysis shows that if the length of the insertion is less than a certain value, the entire insertion converts to Z form at a transition point, but if the insertion is much longer, the B-Z transition exhibits a different behavior, in which part of the insertion flips to Z form and the Z region expands linearly upon changing linking number.     

Theoretical model of the B-Z transition in DNA with an arbitrary sequence

A statistical-mechanical model is suggested that makes it possible to describe the B-Z transition in DNA with an arbitrary sequence of nucleotides. The key point consists in allowance for the fact that each base pair can assume one of the two states with different energies. One of these states corresponds to the standard Z-form with purines in the syn conformation and pyrimidines in the anti conformation. However, in natural DNA sequences such standard base-pair conformations should be interrupted by energetically unfavorable conformations (syn for pyrimidines and anti for purines). Open regions and cruciform structures are also allowed for in the model. The probabilities of formation of the Z-form stretches, open regions and cruciform structures have been calculated for different values of parameters for pBR322 and pAO3 DNA.     

The energetics of the B-Z transition in DNA

The paper deals with the energetics of the transition to left-handed Z form in DNA with an arbitrary base sequence. There is a brief outline of the statistical-mechanical model of the B-Z transition allowing for three possible states of each base pair. The parameters of the model can be determined by comparing the theory with experimental data for the B-Z transition in inserts with given sequences in circular DNA. The model contains six energy parameters, most of which have been determined before. In order to find the remaining parameters of the model and test its adequacy, a number of oligonucleotide sequences were synthesized and inserted into the pUC 19 plasmid. Two-dimensional gel electrophoresis was used to determine the superhelical density at which the inserts adopt the Z form. A statistical-mechanical treatment of these data yielded a complete set of six energy parameters for the B-Z transition. The theoretical assumption that the free energy of Z-form pairs does not depend on the type of adjacent pairs proved to be in agreement with the experimental data.     

Composite cylinder models of DNA: application to the electrostatics of the B-Z transition.

We develop and test a Poisson-Boltzmann model of the electrostatics of the B-Z transition of DNA. Starting from the detailed geometries of the two forms, we compute at each radius the fractions of DNA matter, of volume forbidden (for nonpoint-like ions), and of volume accessible to the center of ions. These radial distributions are incorporated in a composite cylinder model; availability to ions (porosity) and the dielectric constant at each radial distance are then obtained. The phosphate charge is distributed with cylindrical symmetry on two layers at the appropriate radial distances. The porous sheath, between the axis and the charge distribution, provides much more room for ions in B-DNA than in Z-DNA. By using previously developed methods, the Poisson-Boltzmann problem of such cylinders is easily solved. The computational load is small, so that results can be obtained for a large set of salt concentrations and for a number of ionic radii. The variation of the electrostatic free energy difference with salt concentration compares favorably with the experimental value (it is half as large). There is also qualitative agreement with experiments on supercoiled DNA, including a maximum of the free energy difference at submolar salt concentrations. The results for this cylinder with porous sheath are in line with those of the earlier simple planar model and of a plain cylinder with sheath, which is also presented here. They are thus insensitive to details of the model. They support the proposition that the main electrostatic feature of the B-Z transition is the better immersion of the B-DNA phosphates into the solution. They also give confidence in the validity of the Poisson-Boltzmann approach, despite the large salt concentrations involved. Prior studies using an approach based on the potential of mean force are discussed.     

Cooperative chiral order in the B-Z transition in random sequences of DNA.

We present a theory for cooperative chiral order in the transition between right-handed B-DNA and left-handed Z-DNA. This theory, based on the random-field Ising model, predicts the characteristic length scale of Z-DNA segments. This length scale depends on whether the DNA is a homopolymer or a random sequence: it is approximately 4000 nucleotides in a homopolymer but only approximately 25 nucleotides in a random sequence. These theoretical results are consistent with experiments on DNA homopolymers and random sequences.     

Salt dependent premelting base pair opening probabilities of B and Z DNA Poly [d(G-C)] and significance for the B-Z transition

We calculate room temperature thermal fluctuational base pair opening probabilities of B and Z DNA Poly[d(G-C)] at various salt concentrations and discuss the significance of thermal fluctuation in facilitating base pair disruption during B to Z transition. Our calculated base pair opening probability of the B DNA at lower salt concentrations and the probability of the Z DNA at high salt concentrations are in agreement with observations. The salt dependence of the probabilities indicates a B to Z transition at a salt concentration close to the observed concentration.     

An epigenetic role for noncoding RNAs and intragenic DNA methylation

A report on the Keystone Symposium 'Epigenetics: Regulation of Chromatin Structure in Development and Disease', Breckenridge, USA, 11-16 April 2007.     

Effects of 5 cytosine methylation on the B-Z transition in DNA restriction fragments and recombinant plasmids

Alternating (dC-dG)n regions in DNA restriction fragments and recombinant plasmids were methylated at the 5 position of the cytosine residues by the HhaI methylase. Methylation lowers the concentration of NaCl or MgCl2 necessary to cause the B-Z conformational transition in these sequences. Ionic strengths higher than physiological conditions are required to form the Z conformation when the methylated (dC-dG)n tract is contiguous with regions that do not form Z structures, in contrast to the results with the DNA polymer poly(m5dC-dG) . poly(m5dC-dG). In supercoiled plasmids containing (dC-dG)n sequences, methylation reduces the number of negative supercoils necessary to stabilize the Z conformation. Calculations of the observed free energy contributions of the B-Z junction and cytosine methylation suggest that two junctions offset the favorable effect of methylation on the Z conformation in (dC-dG)n sequences (about 29 base-pairs in length). Studies with individual methylated topoisomers demonstrate that increasing Na+ concentration up to approximately 0.2 M inhibits the formation of the Z conformation in the (m5dC-dG)n region of supercoiled plasmids. The results suggest that methylation may serve as a triggering mechanism for Z DNA formation in supercoiled DNAs.     

Application of polyelectrolyte theory to the study of the B-Z transition in DNA (1)

We have used the polyelectrolyte theory to study the ionic strength dependence of the B-Z equilibrium in DNA. A DNA molecule is molded as an infinitely long continuously charged cylinder of radius a with reduced linear charge density q. The parameters a and q for the B and Z forms were taken from X-ray data: aB = 1nm, qB = 4.2, aZ = 0.9 nm and qZ = 3.9. A simple theory shows that at low ionic strengths (when Debye screening length rD much greater than a) the electrostatic free energy difference FelBZ = FelZ - FelB increases with increasing ionic strength since qB greater than qZ. At high ionic strengths (when rD much less than a) the FelBZ would go on growing with increasing ionic strength if the inequality qB/aB greater than qZ/aZ were valid. In the converse case when qZ/qB greater than aZ/aB the FelBZ value decreases with increasing salt concentration at high ionic strength. Since X-ray data correspond to the latter case, theory predicts that the FelBZ value reaches a maximum at an intermediate ionic strength of about 0.1 M (where rD approximately a). We also performed rigorous calculations based on the Poisson-Boltzmann equation. These calculations have confirmed the above criterion of nonmonotonous behaviour of the FelBZ value as a function of ionic strength. Different theoretical predictions for the B-Z transition in linear and superhelical molecules are discussed. Theory predicts specifically that at a very low ionic strength the Z form may prove to be more stable than the B form.(ABSTRACT TRUNCATED AT 250 WORDS)     

BAL 31 nuclease as a probe in concentrated salt for the B-Z DNA junction

The BAL 31 nuclease, an extracellular nuclease from A. espejiana, specifically recognizes and cleaves the salt induced conformational junction between B and Z-DNA. Short segments of (dC-dG) left-handed Z-helix, comprising approximately 1% of the total DNA, are specifically detected within two different recombinant plasmids. The BAL 31 enzyme is highly resistant to inactivation by the presence of high concentrations of a variety of electrolytes that stabilize left-handed helices, is active at physiological pH, and can be used to probe both linear and circular DNAs. Additionally, the nuclease cleaves left-handed (dC-dG)n only very poorly, if at all. Thus, the BAL 31 nuclease can be utilized as a probe for helical junctions and consequently for segments of left-handed DNA that might exist within predominantly right-handed naturally occurring genomes.     

A unified theory of the B-Z transition of DNA in high and low concentrations of multivalent ions

We showed recently that the high-salt transition of poly[d(G-C)]. poly[d(G-C)] between B-DNA and Z-DNA (at [NaCl] = 2.25 M or [MgCl(2)] = 0.7 M) can be ascribed to the lesser electrostatic free energy of the B form, due to better immersion of the phosphates in the solution. This property was incorporated in cylindrical DNA models that were analyzed by Poisson-Boltzmann theory. The results are insensitive to details of the models, and in fair agreement with experiment. In contrast, the Z form of the poly[d(G-m5C)] duplex is stabilized by very small concentrations of magnesium. We now show that this striking difference is accommodated quantitatively by the same electrostatic theory, without any adjustable parameter. The different responses to magnesium of the methylated and nonmethylated polymers do not come from stereospecific cation-DNA interactions: they stem from an experimentally derived, modest difference in the nonelectrostatic component of the free energy difference (or NFED) between the Z and B forms. The NFED is derived from circular DNA measurements. The differences between alkaline earth and transition metal ions are explained by weak coordination of the latter. The theory also explains the induction of the transition by micromolar concentrations of cobalt hexammine, again without specific binding or adjustable parameters. Hence, in the case of the B-Z transition as in others (e.g., the folding of tRNA and of ribozymes), the effect of multivalent cations on nucleic acid structure is mediated primarily by nonspecific ion-polyelectrolyte interactions. We propose this as a general rule for which convincing counter-examples are lacking.     

Reaction conditions affect the specificity of bromoacetaldehyde as a probe for DNA cruciforms and B-Z junctions.

The reaction of bromoacetaldehyde (BAA) was investigated further with recombinant plasmids containing tracts of (CG)16, in pRW756, or (CA)32, in pRW777, which adopt left-handed Z-structures under the influence of negative supercoiling. The cruciform structures adopted by the inverted repeat sequences near the replication origins of the pBR322 vectors served as internal controls for the number of unpaired bases. The extent of reaction with the B-Z junctions and the cruciforms was dependent on the reaction and analysis conditions, the method of preparation of BAA, ionic conditions, and the amount of negative supercoiling. In contrast to the previous results of Kang and Wells, B-Z junctions in addition to cruciforms do react with BAA. However, more forcing conditions are required to detect this reaction since B-Z junctions appear to be less reactive than the single stranded loops of cruciforms. The site of reaction with DNA was readily mapped with high precision at the nucleotide level. Also, a simple method is described for determining the concentration of BAA as well as its intrinsic reactivity in a given ionic medium.     

The lumbosacral transition: morphologic variations and their functional significance

The articular facet of a superior articular process of the sacrum is directed backward, inward, and upward with marked variations. 4 angles characterize the orientation of this facet: a) The relative angle of tilt: i.e. the angle between the articular facet and the upper end-plate of the sacrum, measured in a sagittal plane. b) The absolute angle of tilt: i.e. the angle between the articular facet and the horizontal plane, measured in a sagittal plane. c) The tilted part-angle of opening: i.e. the angle between the articular facet and the sagittal plane, measured in a plane parallel to the upper end-plate of the sacrum. d) The horizontal part-angle of opening: i.e. the angle between the articular facet and the sagittal plane, measured in a horizontal plane. These 4 angles are determined by characteristic straights within the articular facet and certain reference planes (upper end-plate of the sacrum, horizontal plane, sagittal plane). Only 2 intersecting straights suffice for an adequate determination of a geometrical plane; therefore, if we know the relative angle of tilt and the tilted part-angle of opening, we are able to construct or to calculate the absolute angle of tilt as well as the horizontal part-angle of opening by using the range of inclination of the sacrum. The shape as well as the orientation of the articular facets at the superior articular processes of the sacrum do not depend on the inclination of the pelvis nor on the inclination of the sacrum nor on the range of the lumbosacral angle. Only the absolute angle of tilt shows a reference to the inclination of the sacrum because the relative angle of tilt shows a certain constancy. The orientation of the articular facets is slightly influenced by static moments, but considerably determined by dynamical requirements. At spines with irregular numbers of praesacral vertebrae, the orientation of the lumbosacral articular facets do not differ from the orientation of these facets at spines with the regular number of 24 praesacral vertebrae. This, however, does not prove right at spines, that have a lumbosacral "transitional vertebra". Such lumbosacral transitional vertebrae detract much from the stability of the lumbosacral region of the spine.     

Theoretical prediction of base sequence effects in DNA. Experimental reactivity of Z-DNA and B-Z transition enthalpies

Molecular modeling is used to study the sequence dependence of conformation and stability within helically regular duplex Z-DNA. The variations of conformation that are found are sufficiently important to be classified as a new type of polymorphism within the Z family. It is also demonstrated that certain sequences can adopt more than one of these polymorphic forms. Comparison with experimental studies of chemical reactivity within a natural DNA fragment, forced into a left-handed conformation, suggests that the results of our modeling may be used to explain the chemical reactivity observed. Comparison of the Z results with similar studies of the B form allow enthalpies of transition to be calculated as a function of base sequence.     

How proteins search for their specific sites on DNA: the role of DNA conformation

It is known since the early days of molecular biology that proteins locate their specific targets on DNA up to two orders-of-magnitude faster than the Smoluchowski three-dimensional diffusion rate. An accepted explanation of this fact is that proteins are nonspecifically adsorbed on DNA, and sliding along DNA provides for the faster one-dimensional search. Surprisingly, the role of DNA conformation was never considered in this context. In this article, we explicitly address the relative role of three-dimensional diffusion and one-dimensional sliding along coiled or globular DNA and the possibility of correlated readsorption of desorbed proteins. We have identified a wealth of new different scaling regimes. We also found the maximal possible acceleration of the reaction due to sliding. We found that the maximum on the rate-versus-ionic strength curve is asymmetric, and that sliding can lead not only to acceleration, but also in some regimes to dramatic deceleration of the reaction.     

Possible role of long-range interactions in the proton transfer of model DNA systems

The possible role of the long-range interactions has been examined within the semiempirical approach for model doubly stranded DNA systems involving the screw symmetry operation. The interaction energy terms seem to be sensitive to the sequence of base pairs. The essential influence of long-range corrections to the proton transfer potential was found resulting in a remarkably more unsymmetrical energy curve. Only in the case of the singlet excited (n, pi*) electronic state of the A...T base pair, more symmetrical potential is predicted. It is concluded that highly polar sugar-phosphate species are of significant importance for the interaction energy components as well as for related proton transfer processes.