Base sequence, local helix structure, and macroscopic curvature of A-DNA and B-DNA

Given a specified DNA sequence and starting with an idealized conformation for the double helix (A-DNA or B-DNA), the dependence of conformational energy on variations in the local geometry of the double helix can be examined by computer modeling. By averaging over all thermally accessible states, it is possible to determine 1) how the optimum local structure differs from the initial idealized conformation and 2) the energetic costs of small structural deformations. This paper describes such a study. Tables are presented for the prediction of helix twist angles and base pair roll angles for both A-DNA and B-DNA when the sequence has been specified. Local deviations of helix parameters from their average values can accumulate to produce a net curvature of the molecule, a curvature that can be sharp enough to be experimentally detectable. As an independent check on the method, the calculations provide predictions for the longitudinal compressibility (Young's modulus) and the average torsional stiffness, both of which are in good agreement with experimental values. In examining the role of sequence-dependent variations in helix structure for the recognition of specific sequences by proteins, we have calculated the energy needed to deform the self-complementary hexanucleotide d(CAATTG) to match the local geometry of d(GAATTC), which is the sequence recognized by the EcoRI restriction endonuclease. That energy would be sufficient to reduce the binding of the incorrect sequence to the protein by over 2 orders of magnitude relative to the correct sequence.     

Principal component analysis of DNA oligonucleotide structural data

The microstructure of a DNA helix is characterized by several base pair and base step parameters such as twist, rise, roll, propeller twist, etc., in addition to conformational parameters such as the backbone and the glycosidic torsion angles. Among these only a few, which are independent of all others and of each other, may be used to precisely characterize the helix. The problem however is to identify these independent parameters. We have used principal component analysis to identify a relatively small set of independent parameters, with which to characterize each DNA helix. We show that these principal components clearly discriminate between A and B DNA helical types. The calculations further suggest that the microstructure of a DNA helix is better characterized using dinucleotides.     

Statistical Mechanical Approach for Predicting the Transition to Non-B DNA Structures in Supercoiled DNA

Abstract

Supercoiling causes global twist of DNA structure and the supercoiled state has wide influence on conformational transition. A statistical mechanical approach was made for prediction of the transition probability to non-B DNA structures under torsional stress. A conditional partition function was defined as the sum over all possible states of the DNA sequence with basepair 1 and basepair n being in B-form helix and a recurrence formula was developed which expressed the partition function for basepair n with those for less number of pairs. This new definition permits a quick enumeration of every configuration of secondary structures. Energetic parameters of all conformations concerned, involving B-form, interior loop, cruciform and Z-form, were included in the equation. The probability of transition to each non-B conformation could be derived from these conditional partition functions. For treatment of effects of superhelicity, supercoiling energy was considered, and a twist of each conformation was determined to minimize the supercoiling energy. As the twist itself affects the transition probability, the whole scheme of equations was solved by renormalization technique. The present method permits a simultaneous treatment of serveral types of conformations under a common torsional stress.

A set of energetic parameters of DNA secondary structures has been chosen for calculation. Some DNA sequences were submitted to the calculation, and all the sequences that we submitted gave stable convergence. Some of them have been investigated the critical supercoil density for the transition to non-B DNA structures. Even though the reliability of the set of parameters was not enough, the prediction of secondary structure transition showed good agreement with reported observation. Hence, the present algorithm can estimate the probability of local conformational change of DNA under a given supercoil density, and also be employed to predict some specific sequences in which conformational change is sensitive to superhelicity.     

100ps molecular dynamic simulation of d(TATCACC)2.

A heptanucleotide sequence d(TATCACC)2 from OR3 region of bacteriophage lambda is considered sufficient for the recognition of Cro protein. We present here results on molecular dynamic simulations on this sequence for 100 ps in 0.02 ps interval. The simulations are done using computer program GROMOS. The conformational results are averaged over each ps. The IUPAC torsional parameters for 100 conformations are illustrated using a wheal and a dial systems. Several other stereochemical parameters such as H-bonding lengths and angles, sugar puckers, helix twist and roll angles as also distances between opposite strand phosphorus are depicted graphically. We find that there is rupture of terminal H-bonds. The bases are tilted and shifted away from the helix axis giving rise to bifurcated H-bonds. H-bonds are seen even in between different base pairs. The role of these dynamic structural changes in the recognition of OR3 operator by Cro protein is discussed in the paper.     

Directed binding

We propose a novel physical mechanism to describe the mode of processive propagation of twoheaded kinesin motor proteins along microtubule (MT) filaments. Binding and unbinding of the kinesin heads to and from the MT filament play a crucial role in producing movement. The chemical energy of adenosine triphosphate hydrolysis is used in large part for the unbinding process of kinesin from the MT filament. Importantly, in our model, the binding of each head is to be directionally oriented to the MT filament. Therefore, we treat the two motor domains (heads) as extended objects that are connected with each other by a neck region that contains the kinesin dimerization domain. The head domains recognize tubulin binding sites by feeling the two-dimensional periodic potential from the MT surface and are also subjected to thermal noise. Using experimentally determined results regarding physical parameters of the walk, we develop a simple mathematical and mechanical model in which directed binding of the heads to tubulin results in a directed twist of the molecule, probably in the neck linker region, away from its relaxed state. Unbinding of the head from the filament relaxes the twist and defines the propagation direction. We showed that there must be at least two torsional springs (one for every head) involved that can store elastic energy. Consequently, in our model, it is the internal structure both of the relaxed and tensed-up state and the transition mode between them that define the walking direction of kinesin. We present calculations based on the model that are in good quantitative agreement with experimental observations for kinesin.     

100ps Molecular Dynamic Simulation of d(TATCACC)2

Abstract

A heptanucleotide sequence d(TATCACC)2 from OR3 region of bacteriophage X is considered sufficient for the recognition of Cro protein. We present here results on molecular dynamic simulations on this sequence for 100 ps in 0.02 ps interval. The simulations are done using computer program GROMOS. The conformational results are averaged over each ps. The IUPAC torsional parameters for 100 conformations are illustrated using a wheal and a dial systems. Several other stereochemical parameters such as H-bonding lengths and angles, sugar puckers, helix twist and roll angles as also distances between opposite strand phosphorus are depicted graphically. We find that there is rupture of terminal H-bonds. The bases are tilted and shifted away from the helix axis giving rise to bifurcated H-bonds. H- bonds are seen even in between different base pairs. The role of these dynamic structural changes in the recognition of OR3 operator by Cro protein is discussed in the paper.     

Base interaction and conformation manifestations of repetitive nucleotide sequences

To understand why different nucleotide sequences prefer different double helical conformations and to predict conformational behaviour of definite sequences the base-base interaction energy in regular helices consisting of A:U, A:T, G:C and I:C (hypoxanthine-cytosine) base pairs was calculated. Interaction energy was assumed to be a function of eight conformational parameters: H, the distance between adjacent pairs along helix axes; tau, turn angle of one pair relative to the neighbouring one; angles between base planes in a pair (TW, propeller twist and BL, buckle) and position of pairs with respect to helix axes (D and SL, displacements in the plane normal to helix axes, and TL and RL, inclinations to this plane, tilt and roll, respectively). For H and tau characteristic of A- and B-families of nucleic acid conformations (2.5 A less than H less than or equal to 3.5 A, 30 degrees less than or equal to tau less than or equal to 45 degrees) the ranges of conformational parameters corresponding to energy values close to minimal ones (valleys) and correlations between conformational parameters were revealed. Valleys for different sequences largely coincide but have distinctive characteristics for each sequence. Reasons for base pair planarity distortion in double-stranded helices were considered. The calculations permit to account for A-phility of G:C sequences and B-phility for A:T sequences. The valley for I:C sequence branches. This corresponds to A:T-like behaviour in some cases and G:C-like in the others.     

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.     

Conformational search of antisense nucleotides

A preliminary MMFF implementation of selenium atom parameters necessary to model the nucleoside 1 is reported. X-ray structures of two compounds 1 and 2 have been used as references. Ab initio methods have been adopted for checking torsional energy profile and charge distribution. Monte Carlo calculations and energy minimization in solvation complete the conformational search.     

The helical repeat of DNA at high temperature.

The increasing number of studies on thermophilic organisms addressed the question of DNA double helix parameters at high temperature. The present study shows that the helix rotation angle per base pair omega of an unconstrained DNA decreases linearly upon temperature increase, up to the premelting range. In the ionic conditions tested, this rule extends to temperatures up to 85 degrees C, which is a common growth temperature for many hyperthermophilic organisms. In addition, the torsional constant K of DNA decreases with temperature, indicating that the energy required to modify the DNA twist is lower at high temperature. These findings have several implications for people working on the structure and enzymology of DNA at high temperature.     

Periodic structurally similar oligomers are found on one side of the axes of symmetry in the lac, trp, and gal operators

Three well-defined E. coli operator regions were examined for recurring conformational deviation from a regular B-DNA helix. All three, the lac, trp, and gal, show repeats of the same set of neighboring helical twist angles. These angles recur with a periodicity equal to the helix periodicity on one side of the operator's axes of symmetry. The probability that their occurrence is random was found to be extremely small. Therefore, we propose that in addition to specific bases, repeating twist angles patterns are likely to be among the local parameters involved in repressor-operator recognition.     

Conformational analysis and molecular modelling of the branching point of amylopectin

The conformational properties of 6(2) alpha-D-glucosylmaltotriose have been studied using energy calculations that include van der Waals interactions, hydrogen bond stabilization, exo-anometric effect and torsional potential contributions. The calculations focused mainly on the conformational properties displayed at the alpha(1----6) linkage within the tetrasaccharide for which the conformational space is reported. The tetrasaccharide molecule was then considered as a model compound of the branching point in amylopectin. From molecular modelling, some basic structural features associated with branching were clearly established. It was found that, among the low energy arrangements, the side chain would fold back onto the main backbone, thereby producing dense three-dimensional structures in which a 'parallel' arrangement is achieved. The branching between two strands of the double helix, as found in the crystalline moiety of A and B starches, was further investigated. It was found that one particular set of conformations about the glycosidic linkages in the two different strands, could result in an arrangement such that strands could be connected through an alpha(1----6) glycosidic linkage, with a minimum of distortion. The three-dimensional features derived from the molecular modelling agree with the physical properties and mode of biogenesis within the starch granule; they are in accord with a 'cluster' type of structure.     

Periodic Structurally Similar Oligomers are Found on One Side of the Axes of Symmetry in the lac,trp, and gal Operators

Abstract

Three well-defined E. coli operator regions were examined for recurring conformational deviation from a regular B-DNA helix. All three, the lac, trp, and gal, show repeats of the same set of neighboring helical twist angles. These angles recur with a periodicity equal to the helix periodicity on one side of the operator's axes of symmetry. The probability that their occurrence is random was found to be extremely small. Therefore, we propose that in addition to specific bases, repeating twist angle patterns are likely to be among the local parameters involved in repressor-operator recognition.     

The influence of complementary base pair interactions on secondary structure formation and conformational transitions in polynucleotides]

The calculations have been carried out of interaction energy between complementary base pairs of nucleic acids in the function of conformational parametres of double helix (Arnott's parameters) by the method of atom-atom potential functions. Interaction energy as a function of conformational parametres is valley-like and varies little along the bottom of the valley. The regions of interaction energy minima are compared with experimentally determined conformational parametres of nucleic acid double helices. On the basis of calculation results the pathways of conformational transitions between different forms of double-helical polynucleotides are discussed.     

Assignments of 31P NMR resonances in oligodeoxyribonucleotides: origin of sequence-specific variations in the deoxyribose phosphate backbone conformation and the 31P chemical shifts of double-helical nucleic acids

It is now possible to unambiguously assign all 31P resonances in the 31P NMR spectra of oligonucleotides by either two-dimensional NMR techniques or site-specific 17O labeling of the phosphoryl groups. Assignment of 31P signals in tetradecamer duplexes, (dTGTGAGCGCTCACA)2, (dTAT-GAGCGCTCATA)2, (dTCTGAGCGCTCAGA)2, and (dTGTGTGCGCACACA)2, and the dodecamer duplex d(CGTGAATTCGCG)2 containing one base-pair mismatch, combined with additional assignments in the literature, has allowed an analysis of the origin of the sequence-specific variation in 31P chemical shifts of DNA. The 31P chemical shifts of duplex B-DNA phosphates correlate reasonably well with some aspects of the Dickerson/Calladine sum function for variation in the helical twist of the oligonucleotides. Correlations between experimentally measured P-O and C-O torsional angles and results from molecular mechanics energy minimization calculations show that these results are consistent with the hypothesis that sequence-specific variations in 31P chemical shifts are attributable to sequence-specific changes in the deoxyribose phosphate backbone. The major structural variation responsible for these 31P shift perturbations appears to be P-O and C-O backbone torsional angles which respond to changes in the local helical structure. Furthermore, 31P chemical shifts and JH3'-P coupling constants both indicate that these backbone torsional angle variations are more permissive at the ends of the double helix than in the middle. Thus 31P NMR spectroscopy and molecular mechanics energy minimization calculations appear to be able to support sequence-specific structural variations along the backbone of the DNA in solution.     

Conformational studies of (2''-5'') polynucleotides: theoretical computations of energy, base morphology, helical structure, and duplex formation.

A detailed theoretical analysis has been carried out to probe the conformational characteristics of (2'-5') polynucleotide chains. Semi-empirical energy calculations are used to estimate the preferred torsional combinations of the monomeric repeating unit. The resulting morphology of adjacent bases and the tendency to form regular single-stranded structures are determined by standard computational procedures. The torsional preferences are in agreement with available nmr measurements on model compounds. The tendencies to adopt base stacked and intercalative geometries are markedly depressed compared to those in (3'-5') chains. Very limited families of regular monomerically repeating single-stranded (2'-5') helices are found. Base stacking, however, can be enhanced (but helix formation is at the same time depressed) in mixed puckered chains. Constrained (2'-5') duplex structures have been constructed from a search of all intervening glycosyl and sugar conformations that form geometrically feasible phosphodiester linkages. Both A- and B-type base stacking are found to generate non-standard backbone torsions and mixed glycosyl/sugar combinations. The 2'- and 5'-residues are locked in totally different arrangements and are thereby prevented from generating long helical structures.     

The torsional state of DNA within the chromosome

Virtually all processes of the genome biology affect or are affected by the torsional state of DNA. Torsional energy associated with an altered twist facilitates or hinders the melting of the double helix, its molecular interactions, and its spatial folding in the form of supercoils. Yet, understanding how the torsional state of DNA is modulated remains a challenging task due to the multiplicity of cellular factors involved in the generation, transmission, and dissipation of DNA twisting forces. Here, an overview of the implication of DNA topoisomerases, DNA revolving motors, and other DNA interactions that determine local levels of torsional stress in bacterial and eukaryotic chromosomes is provided. Particular emphasis is made on the experimental approaches being developed to assess the torsional state of intracellular DNA and its organization into topological domains.     

The mammalian postorbital bar as a torsion-resisting helical strut

The mammalian skull is asymmetrically loaded during mastication because most of these animals chew on only one side at a time. This loading regime tends to twist the braincase relative to the rostral, tooth bearing part of the skull at the zone of potential weakness between the orbits. This torsional effect is exaggerated, and a postorbital bar is present, in those animals with very large masseter and pterygoid muscles. The lines of action of these muscles are oriented at large angles to the long axis of the skull in lateral view, providing large components of force that twist the skull segments relative to one another. When the temporalis is the dominant muscle, the torsional effect is usually less important, and the bar is absent, because this muscle acts at a smaller angle to the skull axis. The postorbital bar exhibits the predicted three dimensional spatial orientation required to resist torsional forces: it is a segment of an imaginary 45° helix that is wound around the skull, if the skull is idealized as a cylinder. This orientation is significant because, in general, maximum compressive and tensile shear stresses lie along 45° helices on a cylinder loaded in torsion; to resist torsion, material should be placed far from the axis of torsion and along a helix oriented at 45° to the deforming forces. Each half of a supraorbital ridge is also a segment of a 45° helix that is perpendicular to the helix passing through the postorbital bar. This model suggests that the postorbital bar is loaded in compression on the chewing side and in tension on the non-chewing side; the supraorbital ridge is loaded in tension on the chewing side and in compression on the non-chewing side.