Sequence effects on energetic and structural properties of phosphorothioate DNA: a molecular modelling study

Molecular Dynamics Study of the Hybridization between RNA and Modified Oligonucleotides

MicroRNAs (miRNAs) are attractive drug candidates for many diseases as they can modulate the expression of gene networks. Recently, we discovered that DNAs targeting microRNA-22-3p (miR-22-3p) hold the potential for treating obesity and related metabolic disorders (type 2 diabetes mellitus, hyperlipidemia, and nonalcoholic fatty liver disease (NAFLD)) by turning fat-storing white adipocytes into fat-burning adipocytes. In this work, we explored the effects of chemical modifications, including phosphorothioate (PS), locked nucleic acid (LNA), and peptide nucleic acid (PNA), on the structure and energy of DNA analogs by using molecular dynamics (MD) simulations. To achieve a reliable prediction of the hybridization free energy, the AMOEBA polarizable force field and the free energy perturbation technique were employed. The calculated hybridization free energies are generally compatible with previous experiments. For LNA and PNA, the enhanced duplex stability can be explained by the preorganization mechanism, i.e., the single strands adopt stable helical structures similar to those in the duplex. For PS, the S and R isomers (Sp and Rp) have preferences for C2'-endo and C3'-endo sugar puckering conformations, respectively, and therefore Sp is less stable than Rp in DNA/RNA hybrids. In addition, the solvation penalty of Rp accounts for its destabilization effect. PS-LNA is similar to LNA as the sugar puckering is dominated by the locked sugar ring. This work demonstrated that MD simulations with polarizable force fields are useful for the understanding and design of modified nucleic acids.

Iron Oxide Nanoparticles Coated with a Phosphorothioate Oligonucleotide and a Cationic Peptide: Exploring Four Different Ways of Surface Functionalization

The superparamagnetic iron oxide nanoparticles (SPIONs) have great potential in therapeutic and diagnostic applications. Due to their superparamagnetic behavior, they are used clinically as a Magnetic Resonance Imaging (MRI) contrast agent. Iron oxide nanoparticles are also recognized todays as smart drug-delivery systems. However, to increase their specificity, it is essential to functionalize them with a molecule that effectively targets a specific area of the body. Among the molecules that can fulfill this role, peptides are excellent candidates. Oligonucleotides are recognized as potential drugs for various diseases but suffer from poor uptake and intracellular degradation. In this work, we explore four different strategies, based on the electrostatic interactions between the different partners, to functionalize the surface of SPIONs with a phosphorothioate oligonucleotide (ODN) and a cationic peptide labeled with a fluorophore. The internalization of the nanoparticles has been evaluated in vitro on RAW 264.7 cells. Among these strategies, the "«one-step assembly»", i.e., the direct complexation of oligonucleotides and peptides on iron oxide nanoparticles, provides the best way of coating for the internalization of the nanocomplexes.

Theoretical study on the relationship between Rp-phosphorothioation and base-step in S-DNA: based on energetic and structural analysis

Phosphorothioation (PT), previously used in synthetic antisense drugs to arrest the transcription or translation process, is also a novel physiological modification in bacteria DNAs. In the previous study, we reported that Rp-phosphorothioation (Rp-PT) destabilizes B-type helix significantly, using a quantum-mechanics-based energy scoring function developed with a dinucleotide model ( Zhang et al. J. Phys. Chem. B , 2012 , 116 , 10639 - 10648 ). A consequent question surfaces in the field of the phosphorothioated DNA (S-DNA) research: does the endogenous chemical modification interact with the base sequence in the bacterial genomes, e.g., in terms of the most common structure of the B-type helix? In this work, we carried out further energetic analysis on the backbone relative energies calculated with the scoring function according to 16 groups of base-step classifications. Moreover, we conducted molecular dynamics simulations of the B-helical structure with the different base-pair steps, to investigate the detailed structural changes upon the O-/S-substitution. As a result, the Rp-PT modification definitively enhances the stiffness of the backbone and differentiates backbone stability as an interaction with base-steps. Furthermore, certain exceptional sequences such as GT and CC were highlighted in the structural analysis of the sulfur local contacts and relative orientation of double strands, indicating that Rp-PT can cross-talk with particular base-steps. The special effects between the phosphorothioation and base-step may be related to the conservative consensus observed highly frequently in bacterial genomes.

Theoretical study on steric effects of DNA phosphorothioation: B-helical destabilization in Rp-phosphorothioated DNA

Phosphorothioation, with sulfur replacing a nonbridging oxygen of phosphate, has surfaced in bacterial DNA electrophoresis. To understand structural characteristics of the thio-substituted DNA, we have investigated the correlation between the relative energy of phosphate/phosphorothioate linkage and the backbone torsions. The relative energies (R.E.) computed by the quantum mechanical method, the PBE1PBE(CPCM, solvent=water)//PBE1PBE/6-31+G(2df) level of theory, were used to construct energy-scoring functions against backbone torsion variables, resulting in the squared correlation coefficients r(2) of 0.90-0.95. Then, the DNA energy alteration by phosphorothioation is estimated with the relative energy difference (ΔR.E.) between phosphate and phosphorothioate of the phosphate linkages in the DNA crystallographic database (NDB). As a result, Rp-phosphorothioation shifts the relative energy of B-helical structures by 2.7 ± 3.4 kcal/mol, destabilizing about 95% linkages, while Sp-phosphorothioation by -1.4 ± 2.4 kcal/mol, stabilizing over 84% linkages in the data sets. The B-helical destabilization is likely caused by the steric effect between the sulfur atom of Rp-phosphorothioate and the neighboring C-H groups of deoxyribose on the groove wall in B-helix. The unfavorable interaction may be magnified by the increasing rigidness of P-O-involving backbone torsions α and ζ upon the nonbridging phosphorothioations. Since B-helix is the most prevalent DNA double-helical structure and Rp-phosphorothioation is the exclusive configuration in bacteria thio-DNA found to date, the observed stereospecificity-destabilization correlation may reflect a structure-function relationship of biological DNA-phosphorothiation.

The immunostimulatory activity of phosphorothioate CpG oligonucleotides is affected by distal sequence changes

CpG motifs in bacterial DNA activate innate immune cells via toll like receptor 9 (TLR9). Short synthetic oligodeoxynucleotides (ODN) containing a six base CpG motif can mimic the immunostimulatory activity of bacterial DNA. Phosphorothioate (PS) modification of the backbone of ODN makes them more resistant to nuclease degradation and consequently preferable for therapeutic use. Previous studies have shown that the sequence requirements for PS-ODN to have maximal stimulatory activity are more stringent than for normal phosphodiester (PO) ODN. Here we show small sequence changes distal to the CpG motif can affect the activity of PS-ODN whilst having no effect on the activity of PO-ODN. The addition of terminal dG residues and other minor changes to the potently immunostimulatory PS-ODN 1668S caused delayed signalling. The reduction in immunostimulatory activity of PS-ODN was associated with a delay in the activation of MAP kinases.

Geometrical and vibrational features of phosphate, phosphorothioate and phosphorodithioate linkages interacting with hydrated cations: a DFT study

The effect of hexahydrated monovalent and divalent cations on the geometrical and vibrational features of dimethyl phosphate, dimethyl phosphorothioate and dimethyl phosphorodithioate anions (simple suitable model compounds representing the anionic moieties of natural and some modified nucleic acids) was studied. For this purpose, density functional theory (DFT) calculations were carried out at the B3LYP/6-31++G* level. Our results indicate that only K(+) and Mg(2+) prefer to be located in the bisector plane of the PO(2)(-) angle, whereas Li(+) and Na(+) deviate from this plane. Monovalent and divalent cations are slightly deviated from the OPS(-) bisector plane and are found closer to the free oxygen atom. Moreover, the present calculations have shown that in contrast to the general belief, the g(-)g(-) conformer (with respect to the torsion angles defined around the P-O ester bonds) is not always the energetically most favorable. For instance, the g(-)t conformer presents the lowest energy in the case of dimethyl phosphorothioate. The calculated vibrational wavenumbers obtained for dimethyl phosphate and dimethyl phosphorothioate interacting with hydrated sodium counterion, were compared with those previously recorded by Raman scattering and infrared absorption (IR) in aqueous solutions. It has been evidenced that the use of explicit solvent versus dielectric continuum, considerably improves the agreement between the theoretical and observed characteristic wavenumbers.

Synthesis of phosphorothioate oligonucleotides with stereodefined phosphorothioate linkages

A method for solid-phase synthesis of stereodefined PS-oligos via an oxathiaphospholane approach using pure P-diastereomers of nucleoside oxathiaphospholane monomers is described. The oxathiaphospholane monomers are synthesized by phosphitylation of 5'-O-DMTr-N-protected deoxyribonucleosides with 2-chloro-spiro-4,4-pentamethylene-1,3,2-oxathiaphospholane followed by sulfurization. The procedure is general and may be applied to other analogs, depending on the aldehyde (or mercaptoalcohol) used. Starting from an 18O-labeled mercaptoalcohol, the corresponding 18O-labeled phosphitylating reagent and nucleoside monomers can be obtained and used for synthesis of labeled stereodefined PS-oligos, which are useful for studying mechanisms of enzymatic reactions. Details are provided for chromatographic separation of the 5'-O-DMTr-N-protected-deoxyribonucleoside-3'-O-(2-thio-spiro-4,4-pentamethylene-1,3,2-oxathiaphospholane)s into their P-diastereomers, and for manual solid-phase synthesis of PS-oligos. Oxidation of 5'-O-DMTr-N-protected-deoxyribonucleoside-3'-O-(2-thio-spiro-4,4-pentamethylene-1,3,2-oxathiaphospholane)s with selenium dioxide yields their 2-oxo-analogs, which are suitable either for elongation of stereodefined PS-oligos with segments consisting of unmodified nucleotide units possessing phosphate internucleotide linkages, or for generating isotopomeric 18O-labeled PO-oligos of predetermined P-chirality.

Pre-organized structure of viral DNA at the binding-processing site of HIV-1 integrase

The integration of the human immunodeficiency virus type 1 DNA into the host cell genome is catalysed by the viral integrase (IN). The reaction consists of a 3'-processing [dinucleotide released from each 3' end of the viral long terminal repeat (LTR)] followed by a strand transfer (insertion of the viral genome into the human chromosome). A 17 base pair oligonucleotide d(GGAAAATCTCTAGCAGT), d(ACTGCTAGAGATTTTCC) reproducing the U5-LTR extremity of viral DNA that contains the IN attachment site was analysed by NMR using the classical NOEs and scalar coupling constants in conjunction with a small set of residual dipolar coupling constants (RDCs) measured at the 13C/15N natural abundance. The combination of these two types of parameters in calculations significantly improved the DNA structure determination. The well-known features of A-tracts were clearly identified by RDCs in the first part of the molecule. The binding/cleavage site at the viral DNA end is distinguishable by a loss of regular base stacking and a distorted minor groove that can aid its specific recognition by IN.

Effect of phosphorothioate chirality on the grooves of DNA double helices: a molecular dynamics study

Phosphorothioate oligonucleotides (PS-ODNs) have gained considerable attention in drug therapy, primarily as potent antisense or antigene oligomers, which bind to specific DNA or mRNA sequences and lead to transcriptional or translational arrest. These are obtained by substituting one of the anionic oxygen of the phosphate group by a sulfur atom, which introduces chirality to the phosphorus atom of the DNA backbone. In this molecular dynamics simulation study, structural parameters like groove widths, environmental parameters like hydration or cation binding, and electrostatic energy surfaces of both the chiral forms of DNA/PS-DNA duplexes were assessed and compared with that of a normal DNA. Results indicate that, PS-S form with its sulfur atoms facing the minor groove has a widened minor groove, while the scenario is reverse for the PS-R form. Further analysis reveals the existence of several factors like large van der Waals radius of sulfur and the effect it has on its neighboring hydration pattern along with the net electrostatic environment, influencing such structural alterations. This also indicates, for the first time, the effect of absolute phosphorothioate chirality on the global structure of a DNA/PS-DNA hybrid that otherwise resembles a regular B-DNA structure.

NMR and UV studies of 3'-S-phosphorothiolate modified DNA in a DNA : RNA hybrid dodecamer duplex; implications for antisense drug design

High-resolution NMR spectroscopy has been used to establish the conformational consequences of the introduction of a single 3[prime or minute]-S-phosphorothiolate link in the DNA strand of a DNA : RNA hybrid. These systems are of interest as potential antisense therapeutic agents. Previous studies on similarly modified dinucleotides have shown that the conformation of the sugar to which the sulfur is attached shifts to the north (C(3[prime or minute])-endo/C(2[prime or minute])-exo). Comparisons made between NOESY cross-peak intensities, and coupling constants from PE-COSY spectra, for both non-modified and modified duplexes confirm that this conformational shift is also present in the double helical oligonucleotide system. In addition it is noted that in both the dinucleotides and the modified duplex, the conformation of the sugar ring 3[prime or minute] to the site of modification is also shifted to the north. That this pattern is observed in the small monomeric system as well as the larger double helix is suggestive of some pre-ordering of the sequences. The conclusion is supported by consideration of the (1)H chemical shifts of the heterocyclic bases near the site of the modification. The enhanced stability that these conformational changes should bring was confirmed by UV thermal melting studies. Subsequently a series of singly and doubly 3[prime or minute]-S-phosphorothiolate-modified duplexes were investigated by UV. The results are indicative of an additive effect of the modification with thermodynamic benefit being derived from alternate spacing of two modified linkers.

Stereo-enriched phosphorothioate oligodeoxynucleotides: synthesis, biophysical and biological properties

Stereo-enriched [Rp] and [Sp]-phosphorothioate oligodeoxynucleotides are synthesized using oxazaphospholidine derivatized monomers. Three different designs of phosphorothioate oligodeoxynucleotides (PS-oligos), (i) stereo-enriched all-[Rp] or all-[Sp] PS-linkages, (ii) stereo-random mixture of PS-linkages, and (iii) segments containing certain number of stereo-enriched [Rp] and [Sp] PS-linkages ([Sp-Rp-Sp] or [Rp-Sp-Rp]), have been studied. Thermal melting studies of these PS-oligos with RNA complementary strands showed that the binding affinities are in the order [Rp] > [Sp-Rp-Sp]-[Rp-Sp-Rp] > stereo-random > [Sp]. Circular dichroism (CD) studies suggest that the stereochemistry of the PS-oligo does not affect the global conformation of the duplex. The in vitro nuclease stability of these PS-oligos is in the order [Sp] > [Sp-Rp-Sp] > stereo-random > [Rp]. The RNase H activation is in the order [Rp] > stereo-random > [Rp-Sp-Rp] > [Sp] > [Sp-Rp-Sp]. Studies in a cancer cell line of PS-oligos targeted to MDM2 mRNA showed that all oligos had similar biological activity under the experimental conditions employed. Protein- and enzyme-binding studies showed insignificant stereo-dependent binding to proteins. The [Sp] and [Sp-Rp-Sp] chimeric and stereo-random PS-oligos that contained a CpG motif showed higher cell proliferation than [Rp] PS-oligo of the same sequence.