Conformational Studies of Nucleic Acids: III. Empirical Multiple Correlation Functions for Nucleic Acid Torsion Angles

Abstract

There are seven significantly variable torsion angles in each monomer unit of a polynucleotide. Because of this, it is computationally infeasible to consider the energetics of all conformations available to a nucleic acid without the use of simplifications. In this paper, we develop functions suggested by and regression fit to crystallographic data which allow three of these torsion angles, α (03′-P-05′-C5′), δ (C5′-C4′-C3′-03′) and ε (C4′-C3′-03′-P), to be calculated as dependent variables of those remaining. Using these functions, the seven independent torsions are reduced to four, a reduction in complexity sufficient to allow an examination of the global conformational energetics of a nucleic acid for the remaining independent torsion angles. These functions are the first to quantitatively relate a dependent nucleic acid torsion angle to several different independent angles. In all three cases the data are fit reasonably well, and in one case, α, the fit is exceptionally good, lending support for the suitability of the functions in conformational searches. In addition, an examination of the most significant terms in each of the correlation functions allows insight into the physical basis for the correlations.     

Conformational studies of nucleic acids: IV. The conformational energetics of oligonucleotides: d(ApApApA) and ApApApA

Utilizing a new method for modeling furanose pseudorotation (D. A. Pearlman and S.-H. Kim, J. Biomol. Struct. Dyn. 3, 85 (1985)) and the empirical multiple correlations between nucleic acid torsion angles we derived in the previous report (D. A. Pearlman and S.-H. Kim, previous paper in this issue), we have made an energetic examination of the entire conformational spaces available to two nucleic acid oligonucleotides: d(ApApApA) and ApApApA. The energies are calculated using a semi-empirical potential function. From the resulting body of data, energy contour map pairs (one for the DNA molecule, one for the RNA structure) have been created for each of the 21 possible torsion angle pairs in a nucleotide repeating unit. Of the 21 pairs, 15 have not been reported previously. The contour plots are different from those made earlier in that for each point in a particular angle-angle plot, the remaining five variable torsion angles are rotated to the values which give a minimum energy at this point. The contour maps are overall quite consistent with the experimental distribution of oligonucleotide data. A number of these maps are of particular interest: delta (C5'-C4'-C3'-O3')-chi (O4'-C1'-N9-C4), where the energetic basis for an approximately linear delta-chi correlation can be seen: zeta (C3'-O3'-P-O5')-delta, in which the experimentally observed linear correlation between zeta and delta in DNA(220 degrees less than zeta less than 280 degrees) is clearly predicted; zeta-epsilon (C4'-C3'-O3'-P), which shows that epsilon increases with decreasing zeta less than 260 degrees; alpha (O3'-P-O5'-C5')-gamma (O5'-C5'-C4'-C3') where a clear linear correlation between these angles is also apparent, consistent with experiment; and several others. For the DNA molecule studied here, the sugar torsion delta is predicted to be the most flexible, while for the RNA molecule, the greatest amount of flexibility is expected to reside in alpha and gamma. Both the DNA and RNA molecules are predicted to be highly polymorphic. Complete energy minimization has been performed on each of the minima found in the energy searches and the results further support this prediction. Possible pathways for B-form to A-form DNA interconversion suggested by the results of this study are discussed. The results of these calculations support use of the new sugar modeling technique and torsion angle correlations in future conformational studies of nucleic acids.     

Nucleic acid model building: the multiple backbone solutions associated with a given base morphology

A constrained model building procedure is used to generate nucleic acid structures of the familiar A-, B-, and Z-DNA duplexes. Attention is focused upon the multiple structural solutions associated with the arrangements of nucleic acid base pairs rather than the optimum sugar-phosphate structure. The glycosyl (chi) and sugar torsions (both the ring puckering and the exocyclic C5'-C4' (psi) torsion) are treated as independent variables and the resulting O3'...O5' distances are used as closure determinants. When such distances conform to the known geometry of phosphate chemical bonding, an intervening phosphorus atom with correct C-O-P valence angles can be located. Four sequential torsion angles--phi', omega', omega and phi--about the C3'-O3'-P-O5'-C5' bonds are then obtained as dependent variables. The resulting structures are categorized in terms of conformation, ranked in potential energy, and analyzed for torsional correlations. The numerical results are quite interesting with implications regarding nucleic acid models constructed to fit less than ideal experimental data. The multiple solutions to the problem are useful for comprehending the conformational complexities of the local sugar-phosphate backbone and for understanding the transitions between different helical forms. According to these studies, unique characterization of a nucleic acid duplex involves more than the determination of its base pair morphology, its sugar puckering preferences, or its groove binding features.     

Vibrational analysis of nucleic acids. V. Force field and conformation-dependent modes of the phosphodiester backbone modeled by diethyl phosphate.

A generalized valence force field is derived for the diethyl phosphate anion [(CH3CH2O)2PO2-] and its deuterium [(CH3CD2O)2PO2-, (CD3CH2O)2PO2- and (CD3CD2O)2PO2-] and carbon-13 [(CH3 13CH2O)2PO2-] derivatives in the stable trans-gauche-gauche-trans conformation. Normal coordinate analysis of the trans-gauche-gauche-trans conformer, which serves as a structural analog of the nucleic acid phosphodiester group, is based on comprehensive infrared and Raman spectroscopic data and vibrational assignments obtained for the diethyl phosphate anion. The generalized valence force field is in good agreement with the scaled ab initio force field of diethyl phosphate and represents significant improvement over earlier modeling of the phosphodiester moiety with dimethyl phosphate. The conformational dependence of skeletal C-C-O-P(O2-)-O-C-C stretching vibrations is also explored. Starting with the trans-gauche-gauche-trans conformation, the frequency dependence of skeletal stretching modes has been obtained by stepwise rotation of the torsion angles of the P-O and C-O bonds corresponding to nucleic acid torsions alpha (P-O5'), beta (O5'-C5'), epsilon (C3'-O3'), and zeta (O3'-P). Both symmetric and antisymmetric phosphoester stretching modes are highly sensitive to P-O and C-O torsions, whereas symmetric and antisymmetric phosphodioxy (PO2-) stretching modes are less sensitive. The present results provide an improved structural basis for understanding previously developed empirical correlations between vibrational marker bands and nucleic acid backbone conformation.     

Determination of nucleic acid backbone conformation by 1H NMR.

In DNA or RNA duplexes, the six-bond C3'-O3'-P-O5'-C5'-C4'-C3' backbone linkage connecting adjacent residues contains six torsion angles (epsilon, zeta, alpha, beta, gamma, delta) but only four protons. This seriously limits the ability to define the backbone conformation by NMR using purely 1H-1H distance geometry (DG) methods. The problem is further compounded by the inability to assign two of the four backbone protons, namely the poorly resolved H5' and H5' protons, and invariably leads to DG structures with poorly defined backbone conformations. We have developed and tested a reliable method to constrain the beta, gamma, and epsilon (and indirectly alpha and zeta) backbone torsion angles by lower-bound NOE distances to unassigned H5'/H5' resonances combined with either 1H line widths or the conservative use of sigma J measurements; the method relies only on 1H 2-D NMR data, does not involve any structural assumptions, and leads to much improved backbone convergence among DG structures. The C4'-C5' torsion angle gamma is constrained by lower-bound NOE distances from H2' and from H6/H8 to any H5'/H5', as well as by sigma JH4, coupling measurements in the 3.9-4.4 ppm region; delta is constrained by H1'-H4' NOE distances and by H3'-H4' and H3'-H2' J couplings in COSY data; epsilon is partially constrained by H3' line width and/or further constrained by subtracting the minimum possible sigma JH3'-H from the observed sigma JH3' (COSY) to arrive at the maximum possible JH3'-P, which is then converted to H3'-P distance bounds. The angle beta is partially constrained via H5'-P and H5'-P distance bounds consistent with the maximum H5'-P and H5'-P J couplings derived from the observed H5' and H5' line widths, while alpha and zeta are indirectly constrained by lower distance bounds on the observed (n)H1' to (n + 1)H5'/H5' NOEs combined with the prior partial constraints on beta, gamma, delta, and epsilon. The combined effects of these additional constraints in determining distance geometry structures have been demonstrated using a 12-base duplex, [d(GCCGTTAACGGC)]2. Coordinate RMSDs per atom between structures refined with these constraints from random-embedded DG structures, from ideal A-DNA, and from B-DNA starting structures were less than 0.4 A for the central 8 base pairs indicating good convergence. All backbone angles for the central 8 base pairs are very well constrained with less than 10 degrees variation in any of the 48 torsion angles.     

DFT studies of the disaccharide, alpha-maltose: relaxed isopotential maps

The disaccharide, alpha-maltose, forms the molecular basis for the analysis of the structure of starch, and determining the conformational energy landscape as the molecule oscillates around the glycosidic bonds is of importance. Thus, it is of interest to determine, using density functionals and a medium size basis set, a relaxed isopotential contour map plotted as a function of the phi(H) and psi(H) dihedral angles. The technical aspects include the method of choosing the starting conformations, the choice of scanning step size, the method of constraining the specific dihedral angles, and the fitting of data to obtain well defined contour maps. Maps were calculated at the B3LYP/6-31+G( *) level of theory in 5 degrees intervals around the (phi(H),psi(H))=(0 degrees ,0 degrees ) position, out to approximately +/-30 degrees or greater, for gg-gg'-c, gg-gg'-r, gt-gt'-c, gt-gt'-r, tg-tg'-c, and tg-tg'-r conformers, as well as one-split gg(c)-gg'(r) conformer. The results show that the preferred conformation of alpha-maltose in vacuo depends strongly upon the hydroxyl group orientations ('c'/'r'), but the energy landscape moving away from the minimum-energy position is generally shallow and transitions between conformational positions can occur without the addition of significant energy. Mapped deviations of selected parameters such as the dipole moment; the C1-O1-C4', H1-C1-O1, and H4'-C4'-O1 bond angles; and deviations in hydroxymethyl rotamers, O5-C5-C6-O6, O5'-C5'-C6'-O6', C5-C6-O6-H, and C5'-C6'-O6'-H', are presented. These allow visualization of the structural and energetic changes that occur upon rotation about the glycosidic bonds. Interactions across the bridge are visualized by deviations in H(O2)...O3', H(O3')...O2, and H1...H4' distances and the H(O2)-O2-C2-C1 and H'(O3')-O3'-C3'-C4' hydroxyl dihedral angles.     

Carbon-13 NMR in conformational analysis of nucleic acid fragments. 3. The magnitude of torsional angle epsilon in d(TpA) from CCOP and HCOP NMR coupling constants.

Carbon-13 NMR spectra of the deoxyribonucleotide d(TpA), 3',5'-cyclic AMP and 3',5'-cyclic dAMP were measured. It is shown that the different substitution of C2' in deoxyribonucleotides versus ribonucleotides does not affect the vicinal C2'-C3'-O3'-P coupling to a measurable extent. Therefore, the same set of Karplus parameters may be used for the C2'-C3-O3'-P couplings in ribonucleotides and in deoxyribonucleotides. Vicinal carbon-phosphorus and proton-phosphorus coupling constants are used to calculate the magnitude of the torsion angle epsilon (C4'-C3'-O3'-P), which amounts to 195(0) in the trans conformer and to 261(0) in the gauche(-) conformer.     

Structural studies on LNA quadruplexes

LNAs (locked nucleic acids) are new DNA analogues with higher binding affinities toward nucleic acids than the canonical counterparts mainly due to the characteristic conformational restriction arising from the 2'-O, 4'-C methylene bridge. In light of the promising therapeutic applications and considering the advantageous characteristics of LNAs, such as their high water solubility, easy handling, and synthetic accessibility through the conventional phosphoramidite chemistry, we undertook a study concerning the capability of these nucleic acid analogues to form quadruplex structures. Particularly, we have been investigating the LNA/DNA chimeras corresponding to the well-known DNA sequences 5-GGTTGGTGTGGTTGG-3', capable of forming an unimolecular quadruplex. This article deals with the study of the sequence 5'-ggTTggTGTggTTgg-3' (upper and lower case letters represent DNA and LNA residues, respectively), which, according to CD spectroscopy, is able to fold into a quadruplex structure.     

Visualisation of a 2'-5' parallel stranded double helix at atomic resolution: crystal structure of cytidylyl-2',5'-adenosine

X-ray crystallographic studies on 3'-5' oligomers have provided a great deal of information on the stereochemistry and conformational flexibility of nucleic acids and polynucleotides. In contrast, there is very little information available on 2'-5' polynucleotides. We have now obtained the crystal structure of Cytidylyl-2',5'-Adenosine (C2'p5'A) at atomic resolution to establish the conformational differences between these two classes of polymers. The dinucleoside phosphate crystallises in the monoclinic space group C2, with a = 33.912(4)A, b = 16.824(4)A, c = 12.898(2)A and beta = 112.35(1) with two molecules in the asymmetric unit. Spectacularly, the two independent C2'p5'A molecules in the asymmetric unit form right handed miniature parallel stranded double helices with their respective crystallographic two fold (b axis) symmetry mates. Remarkably, the two mini duplexes are almost indistinguishable. The cytosines and adenines form self-pairs with three and two hydrogen bonds respectively. The conformation of the C and A residues about the glycosyl bond is anti same as in the 3'-5' analog but contrasts the anti and syn geometry of C and A residues in A2'p5'C. The furanose ring conformation is C3' endo, C2' endo mixed puckering as in the C3'p5'A-proflavine complex. A comparison of the backbone torsion angles with other 2'-5' dinucleoside structures reveals that the major deviations occur in the torsion angles about the C3'-C2' and C4'-C3' bonds. A right-handed 2'-5' parallel stranded double helix having eight base pairs per turn and 45 degrees turn angle between them has been constructed using this dinucleoside phosphate as repeat unit. A discussion on 2'-5' parallel stranded double helix and its relevance to biological systems is presented.     

Wobble dG X dT pairing in right-handed DNA: solution conformation of the d(C-G-T-G-A-A-T-T-C-G-C-G) duplex deduced from distance geometry analysis of nuclear Overhauser effect spectra

We report below on features of the three-dimensional structure of the d(C-G-T-G-A-A-T-T-C-G-C-G) self-complementary duplex (designated 12-mer GT) containing symmetrical G X T mismatches in the interior of the helix. The majority of the base and sugar protons in the 12-mer GT duplex were assigned by two-dimensional nuclear Overhauser effect (NOESY) spectra in H2O and D2O solution. A set of 92 short (less than 4.5-A) proton-proton distances defined by lower and upper bounds for one symmetrical half of the 12-mer GT duplex were estimated from NOESY data sets recorded as a function of mixing time. These experimental distances combined with nucleotide bond length parameters were embedded into Cartesian space; several trial structures were refined to minimize bond geometry and van der Waals and chirality error. Confidence in this approach is based on the similarity of the refined structures for the solution conformation of the 12-mer GT duplex. The G and T bases pair through two imino-carbonyl hydrogen bonds, and stacking is maintained between the G X T wobble pair and adjacent Watson-Crick G X C pairs. The experimental distance information is restricted to base and sugar protons, and hence structural features such as base pair overlap, glycosidic torsion angles, and sugar pucker are well-defined by this combination of NMR and distance geometry methods. By contrast, we are unable to define the torsion angles about the bonds C3'-O3'-P-O5'-C5'-C4' in the backbone of the nucleic acid.     

Flexibility of 3',5' deoxyribonucleoside diphosphates

In 3',5' deoxyribonucleoside diphosphates, in addition to the nature of the base and the sugar puckering, there are six single bond rotations. However, from the analysis of crystal structure data on the constituents of nucleic acids, only three rotational angles, that are about glycosyl bond, about C4'-C5' and about C3'-O3' bonds, are flexible. For a given sugar puckering and a base, potential energy calculations using non-bonded, electrostatic and torsional functions were carried out by varying the three torsion angles. The energies are represented as isopotential energy surfaces. Since the availability of the real-time color graphics, it is possible to analyse these isopotential energy surfaces. The calculations were carried out for C3' exo and C3' endo puckerings for deoxyribose and also for four bases. These calculations throw more light not only on the allowed regions for the three rotational angles but also on the relationships among them. The dependence of base and the puckering of the sugar on these rotational angles and thereby the flexibility of the 3',5' deoxyribonucleoside diphosphates is discussed. From our calculations, it is now possible to follow minimum energy path for interconversion among various conformers.     

Advanced nuclear magnetic resonance lanthanide probe analyses of short-range conformational interrelations controlling ribonucleic acid structures

An advanced method was developed for lanthanide-probe analyses of the conformations of flexible biomolecules such as nucleotides. The new method is to determine structure parameters (such as internal-rotation angles) and population parameters for local conformational equilibria of flexible sites, together with standard deviations of these parameters. As the prominent advantage of this method, the interrelations among local conformations of flexible sites may be quantitatively elucidated from the experimental data of lanthanide-induced shifts and relaxations and vicinal coupling constants. As a structural unit of ribonucleic acids, the molecular conformations and conformational equilibria of uridine 3'-monophosphate in aqueous solution were analyzed. The stable local conformers about the C3'-O3' bond are the G+ (phi' = 281 +/- 11 degrees) and G- (phi' = 211 +/- 8 degrees) forms. The internal rotation about the C3'-O3' bond and the ribose-ring puckering are interrelated; 97 +/- 5% of the C3'-endo ribose ring is associated with the G- form while 70 +/- 22% o the C2'-endo ribose ring is associated with the G+ form. An interdependency also exists between the internal rotation about the C4'-C5' bond and the ribose-ring puckering. These short-range conformational interrelations are probably important in controlling the dynamic aspects of ribonucleic acid structures.     

Conformational properties of branched RNA fragments in aqueous solution

The conformational properties of branched trinucleoside diphosphates ACC, ACG, AGC, AGG, AUU, AGU, AUG, ATT, GUU, and aAUU [XYZ = X(2'p5'Y)3'p5'Z] have been studied in aqueous solution by nuclear magnetic resonance (1H, 13C), ultraviolet absorption, and circular dichroism. It is concluded from these studies that the purine ring of the central residue (X; e.g., adenosine) forms a base-base stack exclusively with the purine or pyrimidine ring of the 2'-nucleotidyl unit (Y; 2'-residue). The residue attached to the central nucleoside via the 3'-5'-linkage (Z; 3'-residue) is "free" from the influence of the other two heterocyclic rings. The ribose rings of the central nucleoside and the 2'- and 3'-residues exist as equilibrium mixtures of C2'-endo (2E)-C3'-endo (3E) conformers. The furanose ring of the central nucleoside (e.g., A) when linked to a pyrimidine nucleoside via the 2'-5'-linkage shows a higher preference for the 2E pucker conformation (e.g., AUG, AUU, ACG, ca. 80%) than those linked to a guanosine nucleoside through the same type of bond (AGU, AGG, AGC, ca. 70%). This indicates some correlation between nucleotide sequence and ribose conformational equilibrium. The 2E-3E equilibrium of 2'-pyrimidines (Y) shows significant, sometimes exclusive, preference (70-100%) for the 3E conformation; 3'-pyrimidines and 2'-guanosines have nearly equal 2E and 3E rotamer populations; and the ribose conformational equilibrium of 3'-guanosines shows a preference (60-65%) for the 2E pucker. Conformational properties were quantitatively evaluated for most of the bonds (C4'-C5', C5'-O5', C2'-O2', and C3'-O3') in the branched "trinucleotides" AUU and AGG by analysis of 1H-1H, 1H-31P, and 13C-31P coupling constants. The C4'-C5' bond of the adenosine units shows a significant preference for the gamma + conformation. The dominant conformation about C4'-C5' and C5'-O5' for the 2'-and 3'-nucleotidyl units is gamma + and beta t, respectively, with larger gamma + and beta t rotamer populations for the 2'-unit. The increased conformational purity in the 2'-residue, compared to the 3'-residue, is ascribed to the presence of an ordered (adenine----2'-residue) stacked state. The favored rotamers about C3'-O3' and C2'-O2' are epsilon- and epsilon'-, respectively. The conformational features of AUU and AGG were compared to those of their constitutive dimers A3'p5'G, A2'p5'G, A3'p5'U, and A2'p5'U and monomers 5'pG and 5'pU.     

Efficient search on energy minima for structure prediction of nucleic acid motifs

Structure prediction of non-canonical motifs such as mismatches, extra unmatched nucleotides or internal and hairpin loop structures in nucleic acids is of great importance for understanding the function and design of nucleic acid structures. Systematic conformational analysis of such motifs typically involves the generation of many possible combinations of backbone dihedral torsion angles for a given motif and subsequent energy minimization (EM) and evaluation. Such approach is limited due to the number of dihedral angle combinations that grows very rapidly with the size of the motif. Two conformational search approaches have been developed that allow both an effective crossing of barriers during conformational searches and the computational demand grows much less with system size then search methods that explore all combinations of backbone dihedral torsion angles. In the first search protocol single torsion angles are flipped into favorable states using constraint EM and subsequent relaxation without constraints. The approach is repeated in an iterative manner along the backbone of the structural motif until no further energy improvement is obtained. In case of two test systems, a DNA-trinucleotide loop (sequence: GCA) and a RNA tetraloop (sequence: UUCG), the approach successfully identified low energy states close to experiment for two out of five start structures. In the second method randomly selected combinations of up to six backbone torsion angles are simultaneously flipped into preset ranges by a short constraint EM followed by unconstraint EM and acceptance according to a Metropolis acceptance criterion. This combined stochastic/EM search was even more effective than the single torsion flip approach and selected low energy states for the two test cases in between two and four cases out of five start structures.     

Molecular conformation of thyroid hormones: structure and binding interactions of thyroxine

In an accurately determined X-ray diffraction study of the thyroid hormone thyroxine (T4), the two independent conformations in the crystal lattice show significant differences in the outer phenyl ring geometry when compared with that of 3,5,3'-triiodothyronine (T3). The major differences between the T4 and T3 structures are a shortened C4'-O4' bond, contraction of the C3'-C4'-C5' angle and an increase in the C3' and C5' angles of T4. These changes can be correlated with the difference in acidity of the 4'-OH of T4 and T3 and help to explain binding affinity differences among thyroactive compounds. The hydrogen bond directionality observed in T4 and other thyroid structures offers an insight into the molecular details of the hormone-receptor site. The conformation of one T4 molecule is cisoid, that of the other transoid, the first such instance of different overall conformations to be found in the same crystal lattice. One T4 molecule has the side chain nearly coplanar with the inner ring, an unusual conformation among thyroid structures.     

Backbone conformation in nucleic acids: an analysis of local helicity through heminucleotide scheme and a proposal for a unified conformational plot

A relationship has been established to express the local helicity of a polynucleotide backbone directly in terms of the virtual bonds spanning the conformationally equivalent heminucleotide repeats, with a view to provide a better understanding of the cumulative effects of all the chemical bond rotational variations on local helicity. Using this, an analysis made with a few oligodeoxynucleotide crystal structures clearly brings forth that it is the concerted movements manifested in the near neighbour correlations between the pair of chemical bonds C4'-C5' and P-O5' and C4'-C3' and P-O3' of the 5' and 3' heminucleotides respectively that are primarily responsible for the observed non-uniform helical twists both in A and B type helical backbones. That these need not be restricted to oligodeoxynucleotides but may be a feature of oligoribonucleotides backbone also is shown from an analysis of helical segments of yeast tRNA(Phe). A proposal of a unified or a grand two dimensional conformational plot which would help visualise succinctly the overall effect of the variations in all the repeating six chemical bonds of a polynucleotide backbone is made. Apart from considerable simplification, the plot affords identification on it regions characteristic of helical, and loop and bend conformations of nucleic acid backbone chain.     

Index to Subjects,Vol. 27

Abstract

A constrained model building procedure is used to generate nucleic acid structures of the familiar A-, B-, and Z-DNA duplexes. Attention is focused upon the multiple structural solutions associated with the arrangements of nucleic acid base pairs rather than the optimum sugar-phosphate structure. The glycosyl (χ) and sugar torsions (both the ring puckering and the exocyclic C5′-C4′ (ψ) torsion) are treated as independent variables and the resulting O3′…O5′ distances are used as closure determinants. When such distances conform to the known geometry of phosphate chemical bonding, an intervening phosphorus atom with correct C-O-P valence angles can be located. Four sequential torsion angles- φ,ω,ω,ω and φ about the C3′-O3′-P-O5′-C5′ bonds are then obtained as dependent variables. The resulting structures are categorized in terms of conformation, ranked in potential energy, and analyzed for torsional correlations. The numerical results are quite interesting with implications regarding nucleic acid models constructed to fit less than ideal experimental data. The multiple solutions to the problem are useful for comprehending the conformational complexities of thelocal sugar-phosphate backbone and for understanding the transitions between different helical forms. According to these studies, unique characterization of a nucleic acid duplex involves more than the determination of its base pair morphology, its sugar puckering preferences, or its groove binding features.     

Analysis of conformational parameters in nucleic acid fragments. II. Co-crystal complexes of nucleic acid bases.

Studies have been made of conformational parameters in co-crystal complexes and compounds of nucleic acid bases in which there is the possibility of formation of hetero-base-pairs. Using published data extracted from the Cambridge structural database, a total of 37 base-pairs were found, of which 25 were hetero-pairs and 12 homo-pairs. These base-pairs were subject to analysis to reveal hydrogen bond parameters, propeller twist, buckle and C1'-C1' separation (or a similar parameter if C1' atoms were not present). Hetero-pairs were found to show larger twists than homo-pairs, the magnitude of twist being unrelated to hydrogen bond parameters or buckle value. The propeller twisting is less pronounced in these nucleic acid bases than in nucleosides, but still has a significant magnitude. Propeller twisting in hetero-pairs is found to be larger than in homo-pairs. Hetero-pairs appear to be formed preferentially in competitive situations.     

Synthesis and properties of 2'-O,4'-C-ethylene-bridged nucleic acids (ENA) as effective antisense oligonucleotides

Novel bicyclo nucleosides, 2'-O,4'-C-ethylene nucleosides and 2'-O,4'-C-propylene nucleosides, were synthesized as building blocks for antisense oligonucleotides to further optimize the 2'-O,4'-C-methylene-linkage of bridged nucleic acids (2',4'-BNA) or locked nucleic acids (LNA). Both the 2'-O,4'-C-ethylene- and propylene-linkage within these nucleosides restrict the sugar puckering to the N-conformation of RNA as do 2',4'-BNA/LNA. Furthermore, ethylene-bridged nucleic acids (ENA) having 2'-O,4'-C-ethylene nucleosides had considerably increased the affinity to complementary RNA, and were as high as that of 2',4'-BNA/LNA (DeltaT(m)=+3 approximately 5 degrees C per modification). On the other hand, addition of 2'-O,4'-C-propylene modifications in oligonucleotides led to a decrease in the affinity to complementary RNA. As for the stability against nucleases, incorporation of one 2'-O,4'-C-ethylene or one 2'-O,4'-C-propylene nucleoside into oligonucleotides considerably increased their resistance against exonucleases to an extent greater than 2',4'-BNA/LNA. These results indicate that ENA is more suitable as an antisense oligonucleotide and is expected to have better antisense activity than 2',4'-BNA/LNA.