Explicit solvent molecular dynamics simulation of duplex formed by the modified oligonucleotide with alternating phosphate/phosphonate internucleoside linkages and its natural counterpart

Impact of the internucleoside linkage modification by inserting a methylene group on the ability of the modified oligonucleotide to hybridize with a natural DNA strand was studied by fully solvated molecular dynamics (MD) simulations. Three undecamer complexes were analyzed: natural dT(11).dA(11) duplex as a reference and two its analogs with alternating modified and natural linkages in the deoxyadenosine chain. The isopolar, non-isosteric modified linkages were of 5'-O-PO(2)-CH(2)-O-3' (5'PC3') or 5'-O-CH(2)-PO(2)-O-3' (5'CP3') type. Simulations were performed by using the AMBER 5.0 software package with the force field completed by a set of parameters needed to model the modified segments. Both modifications were found to lead to double helical complexes, in which the thymidine strand as well as deoxyriboses and unmodified linkages in the adenosine strand adopted conformations typical for the B-type structure. For each of the two conformational richer modified linkages two stable conformations were found at 300 K: the -ggt and ggt for the 5'PC3' and ggg, tgg for the 5'CP3', respectively. Both modified chains adopted helical conformations with heightened values of the inclination parameter but without affecting the Watson-Crick hydrogen bonds.     

Adenylyl-(5′, 2′)-adenosine 5′-Phosphate ((2′-5′) pA-A) in Crude Crystals of 5′-Adenylic Acid Isolated from Nuclease P1 (Penicillium Nuclease) Digest of Technical Grade Yeast RNA

Adenylyl (5′,2′)-adenosine 5′-phosphate ((2′-5′)pA-A) was detected in crude crystals of 5′-AMP prepared from Penicillium nuclease (nuclease P₁) digest of a technical grade yeast RNA. While (3′–5′)A-A was split by nuclease P₁, spleen phosphodiesterase, snake venom phosphodiesterase or alkali, (2′–5′)A-A was not split by a usual level of nuclease P₁ or spleen phosphodiesterase. Nuclease P₁ digests of 12 preparations of technical grade yeast RNA tested were confirmed to contain (2′–5′)pA-A. Its content was about 1 to 2% of the AMP component of each RNA preparation. As poly(A) was degraded completely by the Penicillium enzyme into 5′-AMP without formation of any appreciable amount of (2′–5′)pA-A, the technical grade RNA is supposed to contain 2–5′ phosphodiester linkages in addition to 3′–5′ major linkages.     

Studies on Alternate Strand Triple Helix Formation by Oligodeoxyribonucleotides Containing A 3′-3′ Phosphodiester Bond

Abstract

ODNs containing a 3′-3′ phosphodiester linkage as inversion of polarity motif have been shown to cooperatively bind to 5′-(purine)_m(pyrimidine)_n-3′ type duplexes by specific alternate strand recognition of the adjacent oligopurine domains. An NMR study has been undertaken to investigate the role of the 3′-3′ linked nucleosides and their nearest neighbours in the stabilization of the triple helical complexes.     

Interaction of molecules with nucleic acids. II. Two pairs of families of intercalation sites, unwinding angles, and the neighbor-exclusion principle

Intercalation-site geometries are generated for a tetramer duplex extracted from B-DNA. Glycosidic angles and puckers of the deoxyribose sugar groups bonded to base pairs BP₁ and BP₄, namely, those at either end of the tetramer duplex, are assumed to be those of B-DNA to insure continuity. All possible geometrical conformations for combinations of C(2′)-endo, C(3′)-endo, C(2′)-exo, and C(3′)-exo sugar puckers are determined for the tetranucleotide backbone. Those with minimum energy are selected as candidates for intercalation sites. Calculations reveal two pairs of physically meaningful families of intercalation sites which occur in two distinct regions, I and II, of helical angles which orient BP₂ relative to BP₃ and with the helical axis disjointed between these base pairs. For each site I and II within BP₂ and BP₃, there are two distinct backbone conformations, A and B, connecting BP₃ to BP₄ or BP₁ to BP₂ which do not disrupt backbone conformations connecting BP₂ to BP₃. Hence two pairs, IA and IB, and IIA and IIB, of intercalation sites exist in which the sugar puckers along the backbone of the tetramer alternate from C(2′)-endo to C(3′)-endo on the backbone (5′p3′) connecting BP₂ to BP₃. The glycosidic angles of the C(3′)-endo sugar χ₃^γ are, coincidentally, 80° ± 2° for both conformations γ = A and B connecting BP₃ to BP₄ along the phosphate backbone (5′p3′). Consistent with the theoretical results, the experimental unwinding angles can be grouped into two categories with absolute values of 18° and 26°. The theoretical unwinding angles for sites IA and IB of 16° and for sites IIA and IIB of 20° occur for a displacement of -0.8 Å in the helical axes of BP₂ and BP₃ and for a 100% G·C composition, with a decrease depending on the amount of A·T base pairs present. Ratios of theoretical unwinding angles of sites I and II, which range from 0.75 to 0.84 for the two principal sites, compare well with the experimental value of 0.71. The theoretical results, in agreement with experimental observation, provide a new interpretation of the nature and conformation of the possible binding sites. Conformations obtained from these studies of intercalation sites in a tetramer duplex are used to rationalize the well-known neighbor-exclusion principle. The possibility of violation of this principle is demonstrated by the existence of two families of physically meaningful conformations. Conformations of unconstrained dimer duplexes are also obtained, one of which corresponds to the experimental crystal structure of ethidium–dinucleoside complexes, but these cannot be joined to the B-DNA structure. Backbone conformations of the tetramer duplex can be constructed until the base-pair separation reaches 8.25 Å, which may limit the molecules that can intercalate.     

Structural Properties of Hybrid Triplex of Polycation Deoxyribonucleic S-methylthiourea (DNmt) Strands with a Complementary DNA Strand,Probed by Nanosecond Molecular Dynamics

Abstract

The replacement of phosphodiester linkages of the polyanion DNA with S-methylthiourea linkers provides the polycation deoxyribonucleic S-methylthiourea (DNmt). Molecular dynamics studies to 1,220 ps of the hybrid triplex formed from octameric DNmt strands d(Tmt)₈ with a complementary DNA oligomer strand d(Ap)₈ have been carried out with explicit water solvent and Na counterions under periodic boundary conditions using the CHARMM force field and the Ewald summation method. The Watson-Crick and Hoogsteen hydrogen-bonding patterns of the A/T tracts remained intact without any structural restraints for triplex structures throughout the simulation. The duplex portion of the triplex structure equilibrated at a B-DNA conformation in terms of the helical rise and other helical parameters. The dynamic structures of the DNmt·DNA·DNmt triplex were determined by examining histograms from the last 800 ps of the dynamics run. These included the hydrogen-bonding pattern (sequence recognition), three-centered bifurcating occurrences, minor groove width variations, and bending of tracts for the hybrid triplex structures. Together with the Watson-Crick hydrogen-bondings, the strong Hoogsteen hydrogen-bondings, the partially maintained three-centered bifurcatings in the Watson-Crick pair, and the medium-strength three-centered bifurcatings in the Hoogsteen pair suggest that the hybrid triplex is energetically favorable as compared to a duplex with similar base stacking, van der Waals interactions, and helical parameters. This is in agreement with our previously reported thermody- namic study, in which only triplex structures were observed in solution. The bending angle measured between the local axis vectors of the first and last helical axis segments is about 20° for the Watson-Crick portion of the averaged structure. Propeller twist (associated with three-centered hydrogen-bonding) up to ?30°, native to DNA AT base pairing, was also observed for the triplex structure. The sugar pseudorotation phase angles and the ring rotation angles for the DNA strand are within the C3′-endo domain and C2′-endo domain for the DNmt strand. Water spines are observed in both minor and major grooves throughout the dynamics run. The molecular dynamics simulations of the structural properties of DNmt·DNA·DNmt hybrid triplex is compared to the DNG·DNA·DNG hybrid triplex (In DNG the -O-(PO₂-)-O- linkers in DNA is replaced by -NH-C(=N₂)-NH-).     

Left-handed Helical Polynucleotides with D-Sugar Phosphodiester Backbones

Naturally occurring polynucleotides have right-handed helical confrontations in the solid state¹ and in solution². Poly(dI-dC)poly(dI-dC) was found to form a left-handed helix in spite of the D-sugar backbone. Also, L-adenylyl-(3′–5′)-L-adenosine synthesized by Tazawa et al⁴. takes up the left-handed stacked conformation. We had synthesized a dinucleoside monophosphate, 8,2′-anhydro-8-mercapto-9-β-D-arabinofuranosyladenine phosphoryl-(3′–5′)-8,2′-anhydro-8-mercapto-9-β-D-arabinofuranosyladenine (A^spA^s) (molecular structure Ia; see also ref. 5) and this compound has a left-handed stacked conformation. The two bases in Ia, having the D-sugar backbone, stacked along the left-handed helical axis; these bases are fixed at ?_CN = ?108° (syn-anti region) by the anhydro linkages.     

Molecular Conformation of 2′-Deoxy-2′-methylidene-cytidine: A Potent Antineoplastic Nucleoside

Abstract

2′-Deoxy-2′-methylidenecytidine (DMDC), a potent inhibitor of the growth of tumor cells, was crystallized with two different forms. One is dihydrated (DMDC·2H₂O) and the other is its hydrochloride salt (DMDC·HCLl). Both crystal and molecular structures have been determined by the X-ray diffraction method. In both forms the glycosidic and sugar conformations are anti and C(4′)-exo, respectively, whereas the conformation about the exocyclic bond is trans for DMDC·2H₂O and gauche ⁺ for DMDC·HCl. Proton nuclear magnetic resonance data of DMDC indicate a preference for the anti C(4′)-exo conformation found in the solid state. These molecular conformations were compared with the related pyrimidine nucleosides. When the cytosine bases are brought into coincidence, DMDC displays the exocyclic C(4′)-C(5′) bond located on the very close position to those of pyrimidine nucleosides with typical overall conformations. On the other hand, the hydroxyl O(3′)-H groups are separated by ca. 3 Å in the cases of DMDC and other pyrimidine nucleosides which have the C(2′)-endo sugar conformation. This result may be useful for the implication about the mechanism of the biological activity of DMDC.     

E. COLI RNASE HI AND THE PHOSPHONATE-DNA/RNA HYBRID: MOLECULAR DYNAMICS SIMULATIONS

A model for the complex between E. coli RNase HI and the DNA/RNA hybrid (previously refined by molecular dynamics simulations) was used to determine the impact of the internucleotide linkage modifications (either 3′–O–CH₂–P–O–5′ or 3′–O–P–CH₂–O–5′) on the ability of the modified-DNA/RNA hybrid to create a complex with the protein. Modified internucleotide linkages were incorporated systematically at different positions close to the 3′-end of the DNA strand to interfere with the DNA binding site of RNase H. Altogether, six trajectories were produced (length 1.5). Mutual hydrogen bonds connecting both strands of the nucleic acids hybrid, DNA with RNase H, RNA with RNase H, and the scissile bond with the Mg⁺⁺ · 4H₂O chelate complex (bound in the active site) were analyzed in detail. Many residues were involved in binding of the DNA (Arg88, Asn84, Trp85, Trp104, Tyr73, Lys99, Asn100, Thr43, and Asn16) and RNA (Gln76, Gln72, Tyr73, Lys122, Glu48, Asn44, and Cys13) strand to the substrate-binding site of the RNase H enzyme. The most remarkable disturbance of the hydrogen bonding net was observed for structures with modified internucleotide linkages positioned in a way to interact with the Trp104, Tyr73, Lys99, and Asn100 residues (situated in the middle of the DNA binding site, where a cluster of Trp residues forms a rigid core of the protein structure).     

Nuclear magnetic resonance study of the impact of ribose 2′-O-methylation on the aqueous solution conformation of cytidylyl-(3′ → 5′)-cytidine

In order to obtain information about the conformational characteristics at the nearestneighbor level in the 2′-O-methylated region of t-RNA, as well as in the bizarre 5′-terminus of eucaryotic mRNA, a detailed nuclear magnetic resonance study of 2′-O-methyl-cytidylyl-(3′ → 5′)-cytidine (CmpC) was conducted. Proton spectra were recorded at 270 MHz in the Fourier mode in D₂O solutions, 0.01M, pD 7.3 in the temperature range 5–80°C. Complete accurate sets of nmr parameters were derived for each of the nucleotidyl units by a combination of homo-nuclear decouplings and simulation iteration methods. The data were translated into conformational parameters using procedures developed in earlier studies from these laboratories. It is shown that the ribofuranose ring exists at a ²E ? ³E equilibrium with clear preference [(75–80)%] for the ³E mode. The C(4′)-C(5′) and C(5′)-O(5′) bonds form a stable conformational network with outspoken preference for conformers in which Ψ₁, Ψ₂ ? 60° and ?₂ ? 180°. The orientation of the 3′-phosphate and 2′-O-methyl groups is such that ?₁′ ? 210° and ?″ ? 60°. The phosphodiester bonds are flexible and shift trends for base, H(1′), and H(5″) suggest the existence of a conformational blend of right-handed stack (g^?g^?), left-handed stack (g⁺g⁺), and unstacked arrays (tg^? and tg⁺). Elevation of temperature perturbs the ²E ? ³E equilibrium accompanied with modest depopulation of ψ₁, ψ₂ ? 60° and ?₂ ? 180° conformers. The major effect of elevation of temperature is in the increase of unstacked arrays at the expense of g^?g^? and g⁺g⁺ conformers. The shift trend of Cmp-H(3′) with temperature shows that torsional variation about O(3′)-P is facilitated by increase in temperature and the preferred rotamer about O(3′)-P in the unstacked form is t (ω₁′ = 180°). A detailed comparison of the aqueous solution conformations of CpC and CmpC reveals that 2′-O-methylation causes: (i) a reduction in the magnitude of _χ1; (ii) an increase in the population of ³E pucker at the 3′-nucleotidyl unit; and (iii) modest perturbations in the O(3′)-P and P-O(5′) bond conformations. Comparison of the aqueous solution conformations of AmpA and CmpC makes clear that the conformational properties of pyrimidine-pyrimidine and purine-purine dimers which carry a 2′-O-methylated 3′-nucleotidyl unit are significantly different.     

Crystal Structure and Molecular Conformation of 4-Thiouracil-2′-Trifluorothioacetamide -3′, 5′-Diacetyl-β-D-Riboside

Abstract

4-thiouracil-2′-trifluorothioacetamide-3′, 5′-diacetyl-β-D-riboside is one of the modified thiouracil analogs synthesized in our institute. The determination of the crystal and molecular structure of this compound was carried out with a view to study the conformation of the molecule in the solid state as well as to investigate the conformations of the trifluoroacetamide and the acetyl substituents of the ribose and their effects on the conformation of the ribose ring. Crystals of 4-thiouracil-2′-trifluorothioacetamide-3′,5′- diacetyl-β-D-riboside are orthorhombic, space group P2₁ 2₁ 2₁, with cell dimensions a= 15.351 (2), b= 15.535 (1), c= 8.307 (1) Å, V=1981.0 (7) ų, Z=4, D_m= 1.53, D_c=1.527 g/c.c. and μ=30.1cm -¹. The structure was determined using CuKα (λ, =1.5418 Å) at a temperature T of 297K, with 2333 reflections, which were collected on a Enraf-Nonius CAD-4 diffactometer, out of which 2249 (I ≥20) were considered observed. The structure was determined by direct methods using MULTAN and refined by full matrix least squares method to a final reliability factor of 0.054 and a weighted R factor of 0.079. The nucleoside is in the anti conformation [X_CN =51.4 (5)°], the ribose has the unusual C (2′) endo -C (1′) exo (²T₁), and a g+ conformation [ψ=47.5 (4)] across C(4′)-C(5′) bond. The pseudorotation angle P is 152.8 (4) ° and the amplitude of pucker τ_m of 42.7 (3)°. The average C-F bond distance is 1.308 Å. There is no base pairing and the typical base-base hydrogen bonded interactions are not present in this structure. On the other hand, a hydrogen bonded dimer is formed involving C(3′) - H(3′)… O (2) and N(3) -H (N3) … O (Al) hydrogen bonds joining the base, ribose ring and the acetyl group. The trend towards longer exocyclic bonds at the acetyl centers in compounds with strongly electronegative aglycones, is also exhibited in this compound, with C(3′)-O(3′) and C(5′)-0(5′) being much longer than C(1′)-O(4′). The acetyl groups also take part in C-H…O hydrogen bonding with the acetyl oxygen atom OA2.     

The Synthesis of Cyclonucleotides with Fixed Glycosidic Bond Linkages as Putative Agonists for P2-Purinergic Receptors

Abstract

Cyclonucleotides with fixed glycosidic bond linkages were investigated as possible ligands for purinoceptors in PC12 cells. P₂Y₂-purinoceptors were not activated by the ATP analogue, 8,2′-thioanhydroadenosine-5′-triphosphate (4) and only weakly by the UTP analogue, 2,2′-anhydrouridine-5′-triphosphate (6). However, both analogues were agonists for P2X₂-purinoceptors although the potencies were approximately 30-fold less than that of the parent nucleotides.     

Hydrolysis of some mRNA 5′-Cap Analogs Catalyzed by the Human Fhit Protein - and Lupin ApppA Hydrolases

Abstract

Hydrolysis of the following four cap analogs: m⁷G(5′)ppp(5′)A, m⁷G(5′)ppp(5′)m⁶A, m⁷G(5′)ppp(5′)m^2′OG and m⁷G(5′)ppp(5′)2′dG catalyzed by homogeneous human Fhit protein and yellow lupin Ap₃A hydrolase has been investigated. The hydrolysis products were identified by HPLC analysis and the K_m and V_max values calculated based on the data obtained by the fluorimetric method.     

Substrate/Inhibitor Properties of Human Deoxycytidine Kinase (dCK) and Thymidine Kinases (Tk1 And Tk2) Towards the Sugar Moiety of Nucleosides,Including O′-Alkyl Analogues

Abstract

Nucleoside analogues with modified sugar moieties have been examined for their substrate/inhibitor specificities towards highly purified deoxycytidine kinase (dCK) and thymidine kinases (tetrameric high-affinity form of TK1, and TK2) from human leukemic spleen. In particular, the analogues included the mono-and di-O′-methyl derivatives of dC, dU and dA, syntheses of which are described. In general, purine nucleosides with modified sugar rings were feebler substrates than the corresponding cytosine analogues. Sugar-modified analogues of dU were also relatively poor substrates of TK1 and TK2, but were reasonably good inhibitors, with generally lower K_i values vs TK2 than TK1. An excellent discriminator between TK1 and TK2 was 3′-hexanoylamino-2′,3′-dideoxythymidine, with a K_i of ～600 μM for TK1 and ～0.1 μM for TK2. 3′-OMe-dC was a superior inhibitor of dCK to its 5′-O-methyl congener, consistent with possible participation of the oxygen of the (3′)-OH or (3′)-OMe as proton acceptor in hydrogen bonding with the enzyme. Surprisingly α-dT was a good substrate of both TK1 and TK2, with K_i values of 120 and 30 μM for TK1 and TK2, respectively; and a 3′-branched α-L-deoxycytidine analogue proved to be as good a substrate as its α-D-counterpart. Several 5 ′-substituted analogues of dC were     

Disaccharide Pyrimidine Nucleosides and Their Derivatives: A Novel Group of Cell-Penetrating Inhibitors of Poly(ADP-Ribose) Polymerase 1

Nearly 30 synthetic nucleosides were tested with human recombinant poly(ADP-ribose) polymerase 1 as potential inhibitors of this enzyme. The most active compounds were some disaccharide analogues of thymidine: 3′-O-β-D-ribofuranosyl-5-iodo-dUrd (2d; IC₅₀ = 45 μM), 3′-O-β-D-ribofuranosyl-2′-deoxythymidine (2e; IC₅₀ = 38 μM), and 3′-O-β-D-ribofuranosyl-2′-deoxythymidine oxidized (4; IC₅₀ = 25 μM). These compounds also reduced H₂O₂-induced synthesis of poly(ADP-ribose) in cultured human ovarian carcinoma (SKOV-3) cells in a dose-dependent manner. Furthermore, compounds 2d or 2e until a concentration of 1 mM did not affect growth of SKOV-3 cells, whereas dialdehyde compound 4, as well as thymidine, exhibited a significant cytotoxicity.     

Dynamic Bis-Intercalation of a Homodimeric Thiazole Orange Dye in DNA: Evidence of Intercalator Creeping

Abstract

We have used one and two dimensional exchange ¹H NMR spectroscopy to characterize the dynamics of the binding of a homodimeric thiazole orange dye, 1,1′-(4,4,8,8-tetramethyl-4,8-diaza-undecamethylene)-bis-4-(3-methyl-2,3-dihydro-(benzo-1,3-thiazole)-2-methylidene)-quinolinium tetraiodide (TOTO), to double stranded DNA (dsDNA). The double stranded oligonucleotides used were d-(CGCTAGCG)₂ ( 1 ) and d(CGCTAGCTAGCG)₂ ( 2 ). TOTO binds preferentially to the (5′-CTAG-3′)₂ sites and forms mixtures of 1:1 and 1:2 dsDNA-TOTO complexes with 2 in ratios dependent on the relative amount of TOTO and the oligonucleotide in the sample. The dynamic exchange between preferential binding sites in the case of a 2:1 1 -TOTO mixture is an intermolecular exchange process between two binding sites on different oligonucleotides. In the case of the 1:1 2 -TOTO complex an intramolecular exchange process occur between two different binding sites on the same strand. Both processes were studied. The results demonstrate the ability of TOTO to migrate along a dsDNA strand in an intramolecular exchange process. The migration process (“creeping”) along the DNA strand is 6 times faster than the rate of intermolecular exchange between sites in two different oligonucleotides.     

Conformational features of benzo‐homologated yDNA duplexes by molecular dynamics simulation

yDNA is a base‐modified nucleic acid duplex containing size‐expanded nucleobases. Base‐modified nucleic acids could expand the genetic alphabet and thereby enhance the functional potential of DNA. Unrestrained 100 ns MD simulations were performed in explicit solvent on the yDNA NMR sequence [5′(yA T yA yA T yA T T yA T)₂] and two modeled yDNA duplexes, [5′(yC yC G yC yC G G yC G G)₂] and [(yT5′ G yT A yC yG C yA yG T3′)?(yA5′ C T C yG C G yT A yC A3′)]. The force field parameters for the yDNA bases were derived in consistent with the well‐established AMBER force field. Our results show that DNA backbone can withstand the stretched size of the bases retaining the Watson‐Crick base pairing in the duplexes. The duplexes retained their double helical structure throughout the simulations accommodating the strain due to expanded bases in the backbone torsion angles, sugar pucker and helical parameters. The effect of the benzo‐expansion is clearly reflected in the extended C1′‐C1′ distances and enlarged groove widths. The size expanded base modification leads to reduction in base pair twist resulting in larger overlapping area between the stacked bases, enhancing inter and intra strand stacking interactions in yDNA in comparison with BDNA. This geometry could favour enhanced interactions with the groove binders and DNA binding proteins., 2016. © 2015 Wiley Periodicals, Inc. Biopolymers 105: 55–64, 2016     

Probing the helical forms of Ca2+‐guluronan junction zones in alginate gels by molecular dynamics 1: Duplexes

A molecular dynamics investigation of the helical forms adopted by (1→4)‐α‐L ‐guluronan in explicit water environment was carried out. Single chains and duplexes were modeled at 300 K starting both from 2₁ or 3₂ helical conformations and in the presence of a neutralizing amount of Ca²⁺ ions. All systems were allowed full conformational freedom. The initial perfect helices with integral screw symmetries were lost at the very beginning of simulations and two distinct behaviors were observed: At equilibrium the 2₁ models mostly retained the 2₁ local helical conformations while exploring the 3₂ ones the rest of the time. In duplexes the two chains, which behaved similarly, were well extended and slightly twisted. By contrast, the chains in 3₂ duplex models were dissimilar and explored a much broader conformational space in which 2₁ and 3₂ local helical conformations were dominant and equally represented but the 3₁ and other conformations were also present. The wide variety of conformations revealed in this study is consistent with the general difficulty in obtaining crystals of Ca²⁺‐guluronate with suitable lateral dimensions for crystallographic studies. © 2013 Wiley Periodicals, Inc. Biopolymers 99: 562–571, 2013.     

A Modified Version of the Cornell et al. Force Field with Improved Sugar Pucker Phases and Helical Repeat

Abstract

We have examined some subtle parameter modifications to the Cornell et al. force field, which has proven quite successful in reproducing nucleic acid properties, but whose C2′-endo sugar pucker phase and helical repeat for B DNA appear to be somewhat underestimated. Encouragingly, the addition of a single V₂ term involving the atoms C(sp³)-O-(sp³)-C(sp³)- N(sp²), which can be nicely rationalized because of the anomeric effect (lone pairs on oxygen are preferentially oriented relative to the electron withdrawing N), brings the sugar pucker phase of C2′-endo sugars to near perfect agreement with ab initio calculations (W near 162°). Secondly, the use of high level ab initio calculations on entire nucleosides (in contrast to smaller model systems necessitated in 1994–95 by computer limitations) lets one improve the % torsional potential for nucleic acids. Finally, the O(sp³)-C(sp³)- C(sp³)-O(sp³) V₂ torsional potential has been empirically adjusted to reproduce the ab initio calculated relative energy of C2′- endo and C3′-endo nucleosides. These modifications are tested in molecular dynamics simulations of mononucleosides (to assess sugar pucker percentages) and double helices of DNA and RNA (to assess helical and sequence specific structural properties). In both areas, the modified force field leads to improved agreement with experimental data.     

The occurrence of the syn-C3′ endo conformation and the distorted backbone conformations for C4′-C5′ and P-O5′ in oligo and polynucleotides

Abstract

The nucleoside constituents of nucleic acids prefer the anti conformation (1). When the sugar pucker is taken into account the nucleosides prefer the C2′endo-anti conformation. Of the nearly 300 nucleosides known, about 250 are in the anti conformation and 50 are in the syn-conformation, i.e., anti to syn conformation is 5:1. The nucleotide building blocks of nucleic acids show the same trend as nucleosides. Both the deoxy-guanosine and ribo- guanosine residues in nucleosides and nucleotides prefer the syn-C2′endo conformation with an intra-molecular hydrogen bond (for nucleosides) between the O5′- H and the N3 of the base and, a few syn-C3′endo conformations are also observed. Evidence is presented for the occurrence of the C3′endo-syn conformation for guanines in mis-paired double helical right-handed structures with the distorted sugar phosphate C4′-C5′ and P-O5′ bonds respectively, from g+ (gg) and g- to trans. Evidence is also provided for guanosine nucleotides in left-handed double-helical (Z-DNA) oligo and polynucleotides which has the same syn-C3′endo conformation and the distorted backbone sugar-phosphate bonds (C4′-C5′ and P- O5′) as in the earlier right-handed case.