Conformational Analysis of the 3′-5′-Cyclic Dinucleotide d<pApA> by Means of Molecular Mechanics

Abstract

The conformational properties of the cyclic dinucleotide d<(pApA)> were studied by means of molecular mechanics calculations in which a multiconformation analysis was combined with minimum energy calculations. In this approach models of possible conformers are built by varying the torsion angles of the molecule systematically. These models are then subjected to energy minimization; in the present investigation use was made of the AMBER Force field. It followed that the lowest energy conformer has a pseudo-two-fold axis of symmetry. In this conformer the deoxyribose sugars adopt a N-type conformation. The conformation of the sugar-phosphate backbone is determined by the following torsion angles: α⁺, β¹, γ⁺, ?¹ and ζ⁺. The conformation of this ringsystem corresponds to the structure derived earlier by means of NMR spectroscopy and X-ray diffraction. The observation of a preference for N-type sugar conformations in this molecule can be explained by the steric hindrance induced between opposite H3′ atoms when one sugar is switched from N- to S-type puckers. The sugars can in principle switch from N- to S-type conformations, but this requires at least the transition of γ⁺ to γ^?. In this process the molecule obtains an extended shape in which the bases switch from a pseudo-axial to a pseudo-equatorial position. The calculations demonstrate that, apart from the results obtained for the lowest energy conformation, the 180° change in the propagation direction of the phosphate backbone can be achieved by several different combinations of the backbone torsion angles. It appeared that in the low energy conformers five higher order correlations are found. The combination of torsion angles which are involved in changes in the propagation direction of the sugar-phosphate backbone in DNA-hairpin loops and in tRNA are found in the dataset obtained for cyclic d<(pApA)>. It turns out that in the available examples, 180° changes in the backbone direction are localized between two adjacent nucleotides.     

Conformational analysis of a dinucleotide photodimer with the aid of the genetic algorithm.

The solution structure of the photodimer cis,syn-dUp[]dT is derived with the aid of the genetic algorithm. The conformational space available for the molecule is sampled efficiently using the computer program DENISE and tested against a set of constraints available from nmr experiments. The dominant conformation in solution found with this approach can be described by the following combinations of sugar-phosphate backbone torsion angles: epsilon(t), zeta(t), alpha(+), beta(-ac), and gamma(t). The conformation of the sugars and glycosidic torsion angles are S type and syn, respectively. The cyclobutane ring and pyrimidines are puckered. In addition, other conformations that exist in equilibrium with the first are found. It is concluded that the cyclobutane-pyrimidine system is rigid, whereas the sugar-phosphate backbone is flexible. The solution structures are compared with the crystal structure of the strongly related cyano-ethyl ester of cis,syn-dTp[]dT.     

The conformations of locked nucleic acids (LNA)

We have used 2D NMR spectroscopy to study the sugar conformations of oligonucleotides containing a conformationally restricted nucleotide (LNA) with a 2'-O, 4'-C-methylene bridge. We have investigated a modified 9-mer single stranded oligonucleotide as well as three 9- and 10-mer modified oligonucleotides hybridized to unmodified DNA. The single-stranded LNA contained three modifications whereas the duplexes contained one, three and four modifications, respectively. The LNA:DNA duplexes have normal Watson-Crick base-pairing with all the nucleotides in anti-conformation. By use of selective DQF-COSY spectra we determined the ratio between the N-type (C3'-endo) and S-type (C2'-endo) sugar conformations of the nucleotides. In contrast to the corresponding single-stranded DNA (ssDNA), we found that the sugar conformations of the single-stranded LNA oligonucleotide (ssLNA) cannot be described by a major S-type conformer of all the nucleotides. The nucleotides flanking an LNA nucleotide have sugar conformations with a significant population of the N-type conformer. Similarly, the sugar conformations of the nucleotides in the LNA:DNA duplexes flanking a modification were also shown to have significant contributions from the N-type conformation. In all cases, the sugar conformations of the nucleotides in the complementary DNA strand in the duplex remain in the S-type conformation. We found that the locked conformation of the LNA nucleotides both in ssLNA and in the duplexes organize the phosphate backbone in such a way as to introduce higher population of the N-type conformation. These conformational changes are associated with an improved stacking of the nucleobases. Based on the results reported herein, we propose that the exceptional stability of the LNA modified duplexes is caused by a quenching of concerted local backbone motions (preorganization) by the LNA nucleotides in ssLNA so as to decrease the entropy loss on duplex formation combined with a more efficient stacking of the nucleobases.     

Conformational analysis of m4(2)C-m4(2)C-m6(2)A: a chemically modified 3'-acceptor end of tRNA, studied by NMR and CD methods.

A study on the conformation of the title compound, C-C-A, and on its constituent dinucleotides is presented. 1H-NMR spectra at 360 and 500 MHz were completely assigned by decoupling experiments. Computer simulation of the spectra yielded precise proton-proton and proton-phosphorus coupling constant values. The coupling constants are analyzed in terms of torsion angles and of N- and S-type sugar pucker. 31P-NMR spectra gave some information about P-O backbone torsion angles alpha and zeta. CD spectroscopy was used to obtain insight in the base-base interaction. The C(1) and C(2) unit in C-C-A show normal preference for N-type conformation of the sugar ring, whereas the A(3) residue appears rather biased towards the S-conformation. The zeta and alpha backbone torsion angles in the C-C phosphodiester linkage in C-C-A appear to assume normal g-, g- conformation, the zeta, alpha combination in the C-A linkage is proposed to have a g+, t conformation. In the C-C fragment in C-C-A a regular stack is indicated; it is suggested that the C-A part adopts an unusual antiparallel base stack.     

The solution structure of the circular trinucleotide cr(GpGpGp) determined by NMR and molecular mechanics calculation.

The 3'-5' circular trinucleotide cr(GpGpGp) was studied by means of 1D and 2D high resolution NMR techniques and molecular mechanics calculations. Analysis of the J-couplings, obtained from the 1H and 13C-NMR spectra, allowed the determination of the conformation of the sugar rings and of the 'circular' phosphate backbone. In the course of the investigations it was found that the Karplus-equation most recently parametrized for the CCOP J-coupling constants could not account for the measured J(C4'P) of 11.1 Hz and a new parametrization for both HCOP and CCOP coupling constants is therefore presented. Subsequent analysis of the coupling constants yielded 'fixed' values for the torsion angles beta and delta (with beta = 178 degrees and delta = 139 degrees). The value of the latter angle corresponds to an S-type sugar conformation. The torsion angles gamma and epsilon are involved in a rapid equilibrium in which they are converted between the gauche(+) and trans and between the trans and gauche(-) domain respectively. We show that the occurrence of epsilon in the gauche(-) domain necessitates S-type sugar conformations. Given the aforementioned values for beta, gamma, delta and epsilon the ring closure constraints for the ring, formed by the phosphate backbone can only be fulfilled if alpha and zeta adopt some special values. After energy minimization with the CHARMm force field only two combinations of alpha and zeta result in energetically favourable structures, i.e. the combination alpha (t)/zeta(g-) in case gamma is in a gauche(+) and epsilon is in a trans conformation, and the combination alpha (t)/zeta (g+) for the combination gamma (t)/epsilon (g-). The results are discussed in relation to earlier findings obtained for cd(ApAp) and cr(GpGp), the latter molecule being a regulator of the synthesis of cellulose in Acetobacter xylinum.     

1H-, 13C-, 31P-NMR studies and conformational analysis of NADP+, NADPH coenzymes and of dimers from electrochemical reduction of NADP+.

All H,H, H,P and several C,P coupling constants, including those between C-4' and the vicinal phosphorus atom, have been determined for NADP+, NADPH coenzymes and for a 4,4-dimer obtained from one-electron electrochemical reduction of NADP+. From these data the preferred conformation of the ribose, that of the 1,4-dihydronicotinamide rings, and the conformation about bonds C(4')-C(5') and C(5')-O(5') were deduced. The preferred form of the 1,4- and 1,6-dihydropyridine rings and the conformation about the ring-ring junction were also obtained for all the other 4,4- and 4,6-dimers formed in the same reduction. All the dimers show a puckered structure, i.e., a boat form for the 1,4- and a twist-boat for the 1,6-dihydronicotinamide ring; both protons at the ring-ring junctions are equatorial and have preferred gauche orientation. On the contrary, the reduced coenzyme NADPH displays a planar or highly flexible conformation, rapidly flipping between two limiting boat structures. The conformation of the ribose rings, already suggested for the NADP coenzymes to be an equilibrium mixture of C(2')-endo (S-type) and C(3')-endo (N-type) puckering modes, has been reexamined by using the Altona procedure and the relative proportion of the two modes has been obtained. The S and N families of conformers have almost equal population for the adenine-ribose, whereas for the nicotinamide-ribose rings the S-type reaches the 90%. The rotation about the ester bond C(5')-O(5') and about C(4')-C(5'), defined by torsion angles beta and gamma respectively, displays a constant high preference for the trans conformer beta t (75-80%), whereas the rotamers gamma are spread out in a range of different populations. The values are distributed between the gauche gamma + (48-69%) and the trans gamma t forms (28-73%). The gamma + conformer reaches a 90% value in the case of NADP+ and NMN+. The conformations of the mononucleotides 5'-AMP, NMN+ and NMNH were also calculated from the experimental coupling constant values of the literature.     

A conformational analysis of the (+) anti BPDE adduct to the guanine amino group of dCpdG

The (+) anti isomer of benz[a]pyrene diol epoxide (BPDE), 7 beta, 8 alpha-dihydroxy-9 alpha, 10 alpha-epoxy-7,8,9,10-tetrahydrobenz [a] pyrene has been identified as the probable tumorigenic lesion in mammalian systems. It forms a predominant adduct with DNA at N2 of guanine. In order to elucidate its conformation in atomic resolution detail, minimized conformational potential energy calculations were performed for the adduct with dCpdG. A global conformation search involving about 1000 trials was made. The lowest energy conformation had stacking between the hydrocarbon and the adjacent cytidine, in agreement with CD studies on modified GpU and UpG. This conformer differed from the B form most notably in the guanine glycosidic torsion, which is high anti. The next lowest energy form had torsion angles like the B form, with guanine-cytidine stacking. These two conformers differ in energy by only 2.1 kcal./mole, suggesting that their relative stability could easily be reversed in larger polymers, or under specific environmental conditions. Other conformations, with base-hydrocarbon or base-base stacking are also found, at somewhat higher energies. The Z form is at 7.8 kcal./mole. Thus, this adduct stabilizes the B form, in contrast with the N2 linked AAF adduct, which stabilizes the Z conformation.     

Statistical and energetic analysis of side-chain conformations in oligopeptides

The distributions of side-chain conformations in 258 crystal structures of oligopeptides have been analyzed. The sample contains 321 residues having side chains that extend beyond the C beta atom. Statistically observed preferences of side-chain dihedral angles are summarized and correlated with stereochemical and energetic constraints. The distributions are compared with observed distributions in proteins of known X-ray structures and with computed minimum-energy conformations of amino acid derivatives. The distributions are similar in all three sets of data, and they appear to be governed primarily by intraresidue interactions. In side chains with no beta-branching, the most important interactions that determine chi 1 are those between the C gamma H2 group and atoms of the neighboring peptide groups. As a result, the g- conformation (chi 1 congruent to -60 degrees) occurs most frequently for rotation around the C alpha-C beta bond in oligopeptides, followed by the t conformation (chi 1 congruent to 180 degrees), while the g+ conformation (chi 1 congruent to 60 degrees) is least favored. In residues with beta-branching, steric repulsions between the C gamma H2 or C gamma H3 groups and backbone atoms govern the distribution of chi 1. The extended (t) conformation is highly favored for rotation around the C beta-C gamma and C gamma-C delta bonds in unbranched side chains, because the t conformer has a lower energy than the g+ and g- conformers in hydrocarbon chains. This study of the observed side-chain conformations has led to a refinement of one of the energy parameters used in empirical conformational energy computations.     

Effects of initial settings on computational protein–ligand docking accuracies for several docking programs

In this study, the influences of initial settings, i.e. initial conformations, configurations and docking parameters, on docking results were investigated. The conformations used in the study were generated by the CAMDAS program. After the conformational search calculations, five structures were selected from the conformer groups according to their conformation energies and root mean square deviations against crystal structures; for example, the lowest energy conformer, as well as the closest and farthest conformers to the crystal structure, was retrieved. Several docking parameter settings were used (default, high speed, generating 50 poses). In this study, docking calculations were conducted using the GOLD, eHiTS, AutoDock, AutoDock vina, FRED and DOCK programs. The success rates of GOLD, eHiTS and FRED were better than those of AutoDock, AutoDock vina and DOCK. The docking results using the farthest conformations were worse than those obtained using other conformations, indicating that some conformation search for the ligand molecule should be performed before the docking calculations.     

A Conformational Analysis of the (+) Anti BPDE Adduct to the Guanine Amino Group of dCpdG

Abstract

The (+) anti isomer of benz[a]pyrene diol epoxide (BPDE), 7β, 8a-dihydroxy-9α,10α-epoxy- 7,8,9,10-tetrahydrobenz[a]pyrene has been identified as the probable tumorigenic lesion in mammalian systems. It forms a predominant adduct with DNA at N² of guanine. In order to elucidate its conformation in atomic resolution detail, minimized conformational potential energy calculations were performed for the adduct with dCpdG. A global conformation search involving about 1000 trials was made. The lowest energy conformation had stacking between the hydrocarbon and the adjacent cytidine, in agreement with CD studies on modified GpU and UpG. This conformer differed from the B form most notably in the guanine glycosidic torsion, which is high anti. The next lowest energy form had torsion angles like the B form, with guanine-cytidine stacking. These two conformers differ in energy by only 2.1 kcal./mole, suggesting that their relative stability could easily be reversed in larger polymers, or under specific environmental conditions. Other conformations, with base-hydrocarbon or base-base stacking are also found, at somewhat higher energies. The Z form is at 7.8 kcal./mole. Thus, this adduct stabilizes the B form, in contrast with the N²linked AAF adduct, which stabilizes the Z conformation.     

A stereochemical model for phospholipids. II: Conformational analysis and molecular packing of phosphatidylethanolamine (PE)

A partition energy method procedure was applied to select the energetically favoured conformations of phosphatidylethanolamine (PE) as polar constituents of phospholipid molecules. The result indicated a large degree of freedom for the two torsion angles of the ester bond of the phosphate and a gauche, gauche star conformation for the ethane bond.A packing process of the molecule was carried out through a potential energy calculation by considering the conformers selected above, using previously published procedure and conventions. All the arrangements which possess the best packing energy values were characterised by an orientation of the PN dipolar segment parallel to the lattice plain. Rotation of the internal torsion angles and rotation in the eulerian space of the molecule produced differences in the charged groups that interact. An additional minimum was present in the energy packing process of those conformers which have the first torsion angle of the phosphate in a trans conformation. This minimum, which corresponds to an orientation of the molecule orthogonal to the lattice plane, requires a complete neutralisation of the point charges on the system.The results of the calculation underline the importance of changes in the behaviour of the polar group of the phospholipids in the packing process.     

Polypurine/polypyrimidine hairpins form a triple helix structure at low pH.

1D and 2D NMR investigations of the 15 residue deoxynucleotide sequence d(TCTCTC-TTT-GAGAGA) show that above pH = 6.5 the molecule adopts a B-form hairpin conformation. As the pH is lowered below 6.5 molecules progressively associate in pairs to form a partially triple helical, partially single stranded structure in which the bases of the oligopyrimidine d(TC)3 tract from one molecule form Hoogsteen pairs with the d(G-A)3 tract of the other. Imino protons of protonated cytosines can be observed at very low field (approximately 15 ppm). The enthalpy of triplex formation was estimated by NMR techniques to be -16 kcal mol-1. Intense H6 to H3' cross peaks from residues in all three strands suggest the presence of N-type sugars at some but not at all possible sites. Surprisingly strong cross peaks between H5' or H5" and non-exchangeable base protons are also observed. These suggest that certain of the O5'-C5'-C4'-C3' phosphate backbone torsion angles (gamma) are unusual.     

Study of furanose ring flexibility in polynucleotide chains using Raman spectra analysis

The bands within the range of 800–850 cm^?1 of Raman spectra of polynucleotides sensitive to the change in conformation of sugar–phosphate backbone are analyzed theoretically. The bands are interpreted as the appearance of a quasi-local deoxyribose vibrational mode whose frequency is dependent on the ring puckering. The localization region of the vibrational mode is pointed out. The theory establishes a relationship between the observed spectral intensity and the population of deoxyribose conformational states described in the framework of the pseudorotation concept. The approach developed allowed one to describe the band shapes and their temperature behavior, and to determine the pseudorotation potential of deoxyribose in the helix B -form of A · T containing polynucleotides. Using the analysis of Raman spectra of DNA fibers in water–ethanol mixture the deoxyribose flexibility during the B -A transition is investigated in terms of the population of conformers and effective potential. It is shown that N- and S-type deoxyribose conformers are populated in the DNA B -form (those of the S-type are preferable), whereas N-type conformers are primarily populated in the DNA A -form.     

Computation of DNA backbone conformations

We present an algorithm for the computation of 2'-deoxyribose-phosphodiester backbone conformations that are stereochemically compatible with a given arrangement of nucleic acid bases in a DNA structure. The algorithm involves the sequential computation of 2'-deoxyribose and phosphodiester conformers (collectively referred to as a backbone 'segment'), beginning at the 5'-end of a DNA strand. Computation of the possible segment conformations is achieved by the initial creation of a fragment library, with each fragment representing a set of bond lengths, bond angles and torsion angles. Following exhaustive searching of sugar conformations, each segment conformation is reduced to a single vector, defined by a specific distance, angle and torsion angle, that allows calculation of the O(1)' position. A given 'allowed' conformation of a backbone segment is determined based on its compatibility with the base positions and with the position of the preceding backbone segment. Initial computation of allowable segment conformations of a strand is followed by the determination of continuous backbone solutions for the strand, beginning at the 3'-end. The algorithm is also able to detect repeating segment conformations that arise in structures containing geometrically repeating dinucleotide steps. To illustrate the utility and properties of the algorithm, we have applied it to a series of experimental DNA structures. Regardless of the conformational complexity of these structures, we are able to compute backbone conformations for each structure. Hence, the algorithm, which is currently implemented within a new computer program NASDAC (Nucleic Acids: Structure, Dynamics and Conformation), should have generally applicability to the computation of DNA structures.     

A Monte Carlo method for generating structures of short single-stranded DNA sequences

A Monte Carlo method has been developed for generating the conformations of short single-stranded DNAs from arbitrary starting states. The chain conformers are constructed from energetically favorable arrangements of the constituent mononucleotides. Minimum energy states of individual dinucleotide monophosphate molecules are identified using a torsion angle minimizer. The glycosyl and acyclic backbone torsions of the dimers are allowed to vary, while the sugar rings are held fixed in one of the two preferred puckered forms. A total of 108 conformationally distinct states per dimer are considered in this first stage of minimization. The torsion angles within 5 kcal/mole of the global minimum in the resulting optimized states are then allowed to vary by ±10° in an effort to estimate the breadth of the different local minima. The energies of a total of 2187 (3⁷) angle combinations are examined per local conformational minimum. Finally, the energies of all dinucleotide conformers are scaled so that the populations of differently puckered sugar rings in the theoretical sample match those found in nmr solution studies. This last step is necessitated by limitations in the theoretical methods to predict DNA sugar puckering accurately. The conformer populations of the individual acyclic torsion angles in the composite dimer ensembles are found to be in good agreement with the distributions of backbone conformations deduced from nmr coupling constants and the frequencies of glycosyl conformations in x-ray crystal structures, suggesting that the low energy states are reasonable. The low energy dimer forms (consisting of 150–325 conformational states per dimer step) are next used as variables in a Monte Carlo algorithm, which generates the conformations of single-stranded d(CX_nG) chains, where X = A, T and n = 3, 4, 5. The oligonucleotides are built sequentially from the 5′ end of the chain using random numbers to select the conformations of overlapping dimer units. The simulations are very fast, involving a total of 10⁶ conformations per chain sequence. The potential errors in the buildup procedure are minimized by taking advantage of known rotational interdependences in the sugar–phosphate backbone. The distributions of oligonucleotide conformations are examined in terms of the magnitudes, positions, and orientations of the end-to-end vectors of the chains. The differences in overall flexibility and extension of the oligomers are discussed in terms of the conformations of the constituent dinucleotide steps, while the general methodology is discussed and compared with other nucleic acid model building techniques. © 1993 John Wiley & Sons, Inc.     

The influence of helix morphology on co-operative polyamide backbone conformational flexibility in peptide nucleic acid complexes.

A systematic analysis of peptide nucleic acid (PNA) complexes deposited in the Protein Data Bank has been carried out using a set of contiguous atom torsion angle definitions. The analysis is complemented by molecular mechanics adiabatic potential energy calculations on hybrid PNA-nucleic acid model systems. Hitherto unobserved correlations in the values of the (alpha and epsilon) dihedral angles flanking the backbone secondary amide bond are found. This dihedral coupling forms the basis of a PNA backbone conformation classification scheme. Six conformations are thus characterised in experimental structures. Helix morphology is found to exert a significant influence on backbone conformation and flexibility: Watson-Crick PNA strands in complexes with DNA and RNA, that possess A-like base-pair stacking, adopt backbone conformations distinct from those in PNA.DNA-PNA triplex and PNA-PNA duplex P-helix forms. Solvation effects on Watson-Crick PNA backbone conformation in heterotriplexes are discussed and the possible involvement of inter-conformational transitions and dihedral angle uncoupling in asymmetric heteroduplex base-pair breathing is suggested.     

Dinucleoside monophosphate having a high anti conformation: the crystal structure of 8,2'-S-cycloadenylyl-(3'-5')-8,2'-S-cycloadenosine hydrochloride

The crystal and molecular structure of 8,2'-S-cycloadenylyl-(3'-5')-8,2'-S-cycloadenosine (AspAs) hydrochloride has been determined by X-ray method. The conformation of two independent AspAs molecules found in an asymmetric unit are almost identical to each other. The torsion angles concerning the sugar-phosphate backbone are different from those in crystalline dinucleoside monophosphates so far determined by X-rays. Both AspAs molecules are in the sharp bend conformations, i.e. each rotation around P-O bond (omega', omega) is (g-, t) rather than the preferred (g-, g-) or (g+, g+) conformation. There is no intramolecular base stacking or base-pairing but the intermolecular base stacking was found.     

Combined use of spectroscopic and energy calculation methods for the determination of peptide conformation in solution

This paper describes the combined use of energy calculations and spectroscopic data for the determination of peptide conformations in solution. The approach involves (i) experimental measurements of spectroscopic parameters for a molecule, (ii) calculation of these parameters for low-energy conformers previously determined with regard to local fluctuations in conformation and (iii) a random search for statistically weighted combinations of conformers which provide a good agreement between the calculated and experimental data. The above approach was used to study the conformation of a spin-labelled angiotensin molecule (SL-AT). It appears that the C-terminal hexapeptide of SL-AT possesses two geometrically different spatial forms of the backbone in aqueous solution, with mean values for the statistical weight of 0.78 and 0.22, respectively. In contrast, the N-terminal part of the molecule is conformationally labile. The data obtained can be used to describe the conformation of angiotensin in solution.     

NMR solution structure and conformational analysis of the calcium release agent cyclic adenosine 5'-diphosphate ribose (cADPR)

A high-resolution NMR study of the solution structure of the calcium release agent cADPR has been performed. Pseudorotationals analysis reveals that in solution both sugar rings in cADPR adopt predominantly (approximately 75%) South conformations, with the A and N rings adopting approximately 2T3 (C2'-endo(major)-C3'-exo(minor) and 4(3)T (C3'-exo-C4'-endo) conformations, respectively. The backbone torsion angles beta and gamma have also been determined. While the minor North conformers were not observed in the crystal structure of cADPR, the solution values of the major South conformers compare well to those found in crystal structure.     

Conformational characteristics and correlations in crystal structures of nucleic acid oligonucleotides: evidence for sub-states

Sugar phosphate backbone conformations are a structural element inextricably involved in a complete understanding of specific recognition nucleic acid ligand interactions, from early stage discrimination of the correct target to complexation per se, including any structural adaptation on binding. The collective results of high resolution DNA, RNA and protein/DNA crystal structures provide an opportunity for an improved and enhanced statistical analysis of standard and unusual sugar-phosphate backbone conformations together with corresponding dinucleotide sequence effects as a basis for further exploration of conformational effects on binding. In this study, we have analyzed the conformations of all relevant crystal structures in the nucleic acids data base, determined the frequency distribution of all possible epsilon, zeta, alpha, beta and gamma backbone angle arrangements within four nucleic acid categories (A-RNA and A-DNA, free and bound B-DNA) and explored the relationships between backbone angles, sugar puckers and selected helical parameters. The trends in the correlations are found to be similar regardless of the nucleic acid category. It is interesting that specific structural effects exhibited by the different unusual backbone sub-states are in some cases contravariant. Certain alpha/gamma changes are accompanied by C3' endo (north) sugars, small twist angles and positive values of base pair roll, and favor a displacement of nucleotide bases towards the minor groove compared to that of canonical B form structures. Unusual epsilon/zeta combinations occur with C2' (south) sugars, high twist angles, negative values of base pair roll, and base displacements towards the major groove. Furthermore, any unusual backbone correlates with a reduced dispersion of equilibrium structural parameters of the whole double helix, as evidenced by the reduced standard deviations of almost all conformational parameters. Finally, a strong sequence effect is displayed in the free oligomers, but reduced somewhat in the ligand bound forms. The most variable steps are GpA and CpA, and, to a lesser extent, their partners TpC and TpG. The results provide a basis for considering if the variable and non-variable steps within a biological active sequence precisely determine morphological structural features as the curvature direction, the groove depth, and the accessibility of base pair for non covalent associations.