Ab initio investigation of the molecular structure of methyl methoxymethyl phosphonate, a promising nuclease-resistant alternative of the phosphodiester linkage

Novel isosteric, isopolar phosphonate analogs of oligonucleotides: preparation and properties

A synthetic approach leading to novel-type modified oligothymidylates containing an isosteric, isopolar, enzyme-stable C3'-O-P-CH(2)-O-C4'' phosphonate alternative to phosphodiester internucleotide bond was elaborated. The suitable monomers were prepared from 4'-phosphonomethoxy derivatives of alpha-L-threo and beta-D-erythro-2',5'-dideoxythymidine, which were considered interesting as structurally related to nucleoside 5'-monophosphates. The phosphotriester method was applied to the automated synthesis of both homooligomeric phosphonate 15-mer chains and alternating phosphonate-phosphate constructs. The fully modified homooligomers did not hybridize while homooligomers with alternating sequences containing alpha-L-threo-configured units (but not beta-D-erythro-) showed a significant decrease in T(m) values in comparison with natural dT(15). For a comparative study, phosphodiester 4'-CH(3)-substituted oligothymidylate was synthesized and physical studies (NMR, CD, MDS modeling) were undertaken to shed more light on the changes in conformational behavior arising from the chosen structural alterations.

Altered structural fluctuations in duplex RNA versus DNA: a conformational switch involving base pair opening

DNA and RNA are known to have different structural properties. In the present study, molecular dynamics (MD) simulations on a series of RNA and DNA duplexes indicate differential structural flexibility for the two classes of oligonucleotides. In duplex RNA, multiple base pairs experienced local opening events into the major groove on the nanosecond time scale, while such events were not observed in the DNA simulations. Three factors are indicated to be responsible for the base opening events in RNA: solvent-base interactions, 2'OH(n)-O4'(n+1) intra-strand hydrogen bonding, and enhanced rigid body motion of RNA at the nucleoside level. Water molecules in the major groove of RNA contribute to initiation of base pair opening. Stabilization of the base pair open state is due to a 'conformational switch' comprised of 2'OH(n)-O4'(n+1) hydrogen bonding and a rigid body motion of the nucleoside moiety in RNA. This rigid body motion is associated with decreased flexibility of the glycosyl linkage and sugar moieties in A-form structures. The observed opening rates in RNA are consistent with the imino proton exchange experiments for AU base pairs, although not for GC base pairs, while structural and flexibility changes associated with the proposed conformational switch are consistent with survey data of RNA and DNA crystal structures. The possible relevance of base pair opening events in RNA to its many biological functions is discussed.

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.

-CH2- lengthening of the internucleotide linkage in the ApA dimer can improve its conformational compatibility with its natural polynucleotide counterpart

The complete family of ApA phosphonate analogues with the internucleotide linkage elongated by insertion of a -CH2- group was prepared and the hybridisation and structural properties of its members in interaction with polyuridylic acid were investigated using an original 2D Raman approach. Except for the conformationally restricted A(CH)pA(2'3'endo-5') modification, all of the isopolar, non-isosteric analogues form triplex-like complexes with poly(rU) at room temperature, in which two polymer strands are bound by Watson-Crick and Hoogsteen bonds to a central pseudostrand consisting of a 'chain' of A-dimers. For all of these dimers, the overall conformation of the triplexes was found to be similar according to their extracted Raman spectra. A simple semi-empirical model was introduced to explain the observed dependency of the efficiency of triplex formation on the adenine concentration. Apparently, for most of the modifications studied, the creation of a stable complex at room temperature requires the formation of a central pseudostrand, consisting of several adenine dimers. Molecular dynamics calculations were finally performed to interpret the differences in 'cooperative' behaviour between the different dimers studied. The results indicate that the exceptional properties of the Ap(CH2)A(3'-5') dimer could be caused by the 3D conformational compatibility of this modified linkage with the second (Hoogsteen) poly(rU) strand.

Vibrational analysis of phosphorothioate DNA: II. The POS group in the model compound dimethyl phosphorothioate [(CH3O)2(POS)]-

The results of Raman and Infrared (IR) spectroscopic investigations on the vibrational modes of dimethyl phosphorothioate (DMPS) anion, [(CH3O)2(POS)]-, are reported. Ab initio calculations of the vibrational modes, the IR and Raman spectra and the interatomic force constants of DMPS were performed. A normal mode calculation was performed and the results were used to calculate the potential energy distribution for the vibrational modes. This analysis shows that in DMPS the P-S stretching mode has a frequency of about 630 cm-1 and an angle bending mode involving the sulfur atom has a frequency of about 440 cm-1. The proposed vibrational mode assignments will serve as marker bands in the conformational studies of phosphorothioate oligonucleotides which play a central role in the novel antisense therapeutic paradigm.