Yujun Zhou1, Ryohei Kajino1, Seiichiro Ishii2, Kenji Yamagishi2, Yoshihito Ueno1,3,4,5. 1. Graduate School of Natural Science and Technology 1-1 Yanagido Gifu 501-1193 Japan uenoy@gifu-u.ac.jp +81-58-293-2919 +81-58-293-2919. 2. Department of Chemical Biology and Applied Chemistry, College of Engineering, Nihon University 1 Nakagawara, Tokusada, Tamuramachi Koriyama Fukushima 963-8642 Japan. 3. Faculty of Applied Biological Sciences 1-1 Yanagido Gifu 501-1193 Japan. 4. United Graduate School of Agricultural Science 1-1 Yanagido Gifu 501-1193 Japan. 5. Center for Highly Advanced Integration of Nano and Life Sciences (G-CHAIN), Gifu University 1-1 Yanagido Gifu 501-1193 Japan.
Oligonucleotide therapeutics are emerging as a novel solution for treating a variety of genetic and intractable diseases, so-called “undruggable” diseases that cannot be targeted by small molecule or antibody-based drugs.[1] Antisense oligonucleotides (ASOs), which are composed of single-stranded DNA-like oligonucleotides, represent one kind of oligonucleotide therapeutics, and constitute the majority of oligonucleotide therapeutics approved by 2020.[2-6] ASOs are designed to target disease-related pre-mRNAs or mature mRNAs. After forming the DNA/RNA duplexes inside the cells, they are expected to alter the splicing process of pre-mRNAs to control the synthesis of the proteins encoded, or to induce a selective cleavage of the RNA strand processed by ribonuclease H (RNase H), resulting in the prevention of translation into the encoded proteins.[7]However, there are still some challenges left with regard to ASOs composed of natural monomers (deoxynucleotides or ribonucleotides). It is known that a long half-life of therapeutics is necessary for eliciting sufficient pharmacological effects. In the case of ASOs, the half-life depends on the resistance of these ASOs against numerous nucleases present in the blood circulation. Moreover, it is also essential to reduce the side effects of ASOs by weakening their cytotoxicity, and by enhancing their binding affinity to the target RNAs to suppress off-target effects. As a measure of improving the properties of ASOs, various chemically modified nucleoside analogs have been developed, some of which have already been granted approval for clinical usage. For example, one kind of modification on phosphodiesters, phosphorothioate (PS), which is used in mipomersen and nusinersen, among others, could not only increase the nuclease resistance, but also enhance the intracellular internalization of ASOs, with the cytotoxicity of PS also determined in some of the cases.[2,3,8-10]As in our previous reports, inserting an aminoalkyl side chain into the 4′- or 5′-site of nucleosides could significantly increase the nuclease resistance of the resulting oligonucleotides.[11-15] Due to the decrease of nuclease resistance when the terminal amino moiety was protected by acetyl groups, it is assumed that the positive charge on the amino moiety of the aminoalkyl side chain acts as an obstacle to nucleases binding to oligonucleotides, and thus inhibits the degradation of oligonucleotides. However, combined with the 2′-O-methyl (2′-OMe) modification, it was confirmed that the 4′-aminoalkyl modification would cause a decrease in thermal stability, which was related to a deterioration in binding affinity of RNA targets.[12] It was also revealed that the thermal instability could be improved by altering the aminoalkyl side chain from the 4′- to the 5′-site, and, compared with the (R)-5′-aminopropyl-2′-OMe-modified nucleoside analog (shown in Fig. 1, A1), the (S)-5′-aminopropyl-2′-OMe-modified nucleoside analog (shown in Fig. 1, A2), which remained appropriate thermal stability with robust nuclease resistance, was considered the most potent potential candidate for siRNA in our studies so far.[15]
Fig. 1
The chemical structure of (R)-5′-aminopropyl-2′-OMe-uridine (A1), (S)-5′-aminopropyl-2′-OMe-uridine (A2), (S)-5′-C-aminopropyl-thymidine (A3) and (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine (A4).
Based on the results obtained by siRNA studies, we addressed whether the (S)-5′-aminopropyl modification was also suitable for ASOs to enhance nuclease resistance. However, it was reported that the existence of the 2′-OMe modification would act as an obstruction during RNase H-mediated cleavage of the RNA strand.[16,17] Therefore, instead of employing the 2′-OMe modification, we focused on the application of the 2′-arabinofluoro (2′-araF or 2′-β-fluoro) modification. The synthesis and properties of 2′-araF-modified DNAs (FANAs) have been reported by Damha and co-workers.[18] According to these reports, the incorporation of a fluorine atom into the 2′-β-site of the sugar moiety showed no inhibition towards RNase H-mediated cleavage, as well as increased thermal stability of the FANA/RNA duplexes compared with the natural ones.[19,20] In addition, 2′-deoxynucleosides, the natural monomers of DNA, are known not to interfere with the action of RNase H.In this study, we designed two novel thymidine analogs (Fig. 1): (S)-5′-C-aminopropyl-thymidine (A3) and (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine (A4), envisioning that the dual modification would improve the thermal stability of the DNA/RNA duplexes, and enhance the nuclease resistance of ASOs without inhibiting the cleavage of the RNA strand mediated by RNase H. The following sections outline the synthesis of the nucleoside analogs we designed, and the functional properties of oligonucleotides containing these analogs were further evaluated.
Results and discussion
Synthesis of nucleoside analogs
The synthesis route of phosphoramidites corresponding to (S)-5′-C-aminopropyl-thymidine (A3) is shown in Scheme 1. First, compound 1 was prepared according to a previously reported method.[21,22] Hydroboration and oxidation of the allyl moiety in compound 1 quantitatively produced compound 2. Next, site-selective tosylation was achieved by controlling the reagent equivalent, which gave compound 3 in 81% yield, and the tosyl group was replaced with the azide group to afford compound 4 in 80% yield. Then, the 5′-hydroxy function, considered as a secondary alcohol, was protected with the DMTr group to give compound 5 in 91% yield, via a reported method utilizing DMTrOTFA reagent.[23] After converting the azide groups into an amino function by the Staudinger reaction, the resulting amino function was protected with TFA group to give compound 6 in 89% yield. Finally, deprotection of the 3′-O-silyl group of compound 6 by treatment with TBAF afforded compound 7 in 92% yield, followed by a standard phosphorylation procedure to give the corresponding phosphoramidite 8 in 86% yield.
The synthesis route of phosphoramidites corresponding to (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine (A4), similar to that of A3, is shown in Scheme 2. Compound 9 was prepared via a reported method, followed by oxidization and Hosomi–Sakurai allylation to give compound 11 and 12 in 71% and 4% yield, respectively.[18] Configurations of the 5′-carbons in 11 and 12 were determined via nuclear Overhauser effect spectroscopy (NOESY) measurements and 1H NMR. After deprotecting the 3′-O-silyl groups of 11 and 12, the TIPDS group was used to lock the conformations of the sugar moieties of 11 and 12, resulting in compound 20 and 21 in 83% and 84% yields, respectively (Scheme 3). As a result of the NOESY measurements, the NOE was observed between H-3′ and H-5′ of 21, but not in 20, which suggests that compound 11, the main product generated by Hosomi–Sakurai allylation, was the desired (S)-isomer, while compound 12 was the other isomer.[15,24] This conclusions from NOESY measurements were consistent with the results of 1H NMR; the spin–spin coupling constant between H-4′ and H-5′ was measured as 1.84 Hz in 20 and 9.16 Hz in 21, which suggested the conformation of H-4′ and H-5′ were of the gauche- and anti-form, respectively. The pure (S)-isomer 11 was then converted into compound 13 by hydroboration and oxidation of the allyl moiety, followed by site-selective tosylation of the 5′-hydroxy function to give compound 14 in 81% yield. Notably, the solvent for tosylation of 13 was DCM instead of MeCN due to its poor solubility, which was assumed to result from the incorporation of the fluorine atom. The tosyl group was replaced by azidation to afford compound 15 in 76% yield, followed by the 5′-hydroxy function being protected by the DMTr group to give compound 16 in 93% yield. After converting the azide group into amine by the Staudinger reaction, the resulting amine was protected by TFA group to give compound 17 in 96% yield. Deprotection of the 3′-O-silyl group of 17 afforded compound 18 in 95% yield, followed by a standard phosphorylation procedure to give the corresponding phosphoramidite 19 in 89% yield.
Synthetic route of the 3′,5′-TIPDS derivatives 20 and 21.Reagents and conditions: (a) (1) TBAF/THF, THF, r.t., 2 h; (2) TIPDSCl2, pyridine, r.t., 24 h; 20: 83%; 21: 87%.
Oligonucleotide synthesis
The modified nucleoside analogs A3 and A4 were incorporated into DNA oligomers utilizing phosphoramidites 8 and 19, respectively, by a DNA/RNA synthesizer. After the synthesis, to prevent the additional reaction of acrylonitrile with the 5′-C-aminopropyl groups, the CPG beads were treated with 10% dimethylamine in MeCN at room temperature for 5 min, followed by rinsing with MeCN to selectively remove cyanoethyl groups. Then, the DNA oligomers were cleaved from the CPG beads and deprotected by treatment with concentrated NH3 solution for 12 h at 55 °C. The oligomers were purified by 20% denaturing polyacrylamide gel electrophoresis (PAGE) in presence of 7 M urea. The corresponding RNA oligomers used in this study were also prepared by a DNA/RNA synthesizer. In contrast to the DNA oligomers, the RNA oligomers were cleaved from CPG beads and deprotected by treatment with concentrated NH3 solution/40% methylamine (1 : 1, v/v) for 10 min at 65 °C after the synthesis was complete. Then, 2′-O-TBDMS groups in RNA oligomers were removed by incubation with Et3N·3HF (125 μL) in DMSO (100 μL) for 1.5 h at 65 °C. The reaction was quenched with 0.1 M TEAA buffer (pH 7.0), and the mixture was desalted using a Sep-Pak C18 cartridge. Both DNA and RNA oligomers were finally purified by 20% PAGE. The sequences of DNA and RNA used in this study are shown in Tables 1–3, 5, 6, S1 and S2 and Fig. S8.†
The DNA oligomers used in enzymatic stability assays
Abbreviation of DNA
Sequencea
DNA 7
5′-F-d(TCTTTCTCTTTCCCTT)-3′
DNA 8
5′-F-d(TCTTTCTCTTTCCCTT)-3′
DNA 9
5′-F-d(TCTTTCTCTTTCCCTT)-3′
F, T (bold) and T (italic) denote fluorescein, (S)-5′-C-aminopropyl-thymidine (A3) and (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine (A4), respectively.
The DNA/RNA duplexes used in RNase H-dependent cleavage assays
Abbreviation of DNA/RNA duplex
Abbreviation of DNA or RNA
Sequencea
Duplex 16
DNA 1
5′-d(TCTTTCTCTTTCCCTT)-3′
RNA 6
5′-F-r(AAGGGAAGAGAAAGA)-3′
Duplex 17
DNA 2
5′-d(TCTTTCTCTTTCCCTT)-3′
RNA 6
5′-F-r(AAGGGAAGAGAAAGA)-3′
Duplex 18
DNA 3
5′-d(TCTTTCTCTTTCCCTT)-3′
RNA 6
5′-F-r(AAGGGAAGAGAAAGA)-3′
F , T (bold) and T (italic) denote fluorescein, (S)-5′-C-aminopropyl-thymidine (A3) and (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine (A4), respectively.
Thermal stability of duplexes
To investigate the impact of nucleoside analogs A3 and A4 on the thermal stability of the DNA/RNA duplexes, temperature-induced melting was measured by ultraviolet (UV) spectroscopy in 10 mM sodium phosphate buffer (pH 7.0) containing 100 mM NaCl. The 50% melting temperature (Tm) values are shown in Table 1. As a result, the Tm value of the natural DNA/RNA duplex 1 was 60.0 °C, while that of modified duplexes 2 and 3 were 57.3 °C and 59.1 °C, respectively. Additionally, the incorporation of A3 decreased the Tm value of the duplex by 0.7 °C per modification, while the incorporation of A4 decreased the Tm value by 0.2 °C per modification, which was considered as stable as the natural duplex 1. On the other hand, although both A3 and A4 showed a slight negative effect on thermal stability, it was revealed that the introduction of the 2′-araF modification improved the thermal stability of the duplexes.
Sequence of DNA/RNA duplexes, DNAs, RNAs, and Tm values of duplexes
Abbreviation of DNA/RNA duplex
Abbreviation of DNA or RNA
Sequencea
Tmb (°C)
ΔTmc (°C)
ΔTmc (°C) per mod.
Duplex 1
DNA 1
5′-d(TCTTTCTCTTTCCCTT)-3′
60.0
—
—
RNA 1
5′-r(AAGGGAAGAGAAAGA)-3′
Duplex 2
DNA 2
5′-d(TCTTTCTCTTTCCCTT)-3′
57.3
−2.7
−0.7
RNA 1
5′-r(AAGGGAAGAGAAAGA)-3′
Duplex 3
DNA 3
5′-d(TCTTTCTCTTTCCCTT)-3′
59.1
−0.9
−0.2
RNA 1
5′-r(AAGGGAAGAGAAAGA)-3′
T (bold) and T (italic) denote (S)-5′-C-aminopropyl-thymidine (A3) and (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine (A4), respectively.
The Tm value were measured in 10 mM sodium phosphate buffer (pH 7.0) containing 100 mM NaCl. The concentrations of the duplexes were 3 μM. All measurements were carried out three times, and the data are shown as an average value in Table 1.
ΔTm represents [Tm(duplex 2 or 3) − Tm(duplex 1)].
T (bold) and T (italic) denote (S)-5′-C-aminopropyl-thymidine (A3) and (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine (A4), respectively.The Tm value were measured in 10 mM sodium phosphate buffer (pH 7.0) containing 100 mM NaCl. The concentrations of the duplexes were 3 μM. All measurements were carried out three times, and the data are shown as an average value in Table 1.ΔTm represents [Tm(duplex 2 or 3) − Tm(duplex 1)].To determine the reason for the difference in Tm in detail, the thermodynamic parameters of duplexes 1–3 were calculated. The thermodynamic parameters from these calculation were based on the slope of a 1/Tmvs. ln(CT/4) plot, where CT (1, 3, 6, 12, 15, 21, 30 and 60 μM) was the total concentration of each single strand, the results of which are shown in Table 2. The ΔG° value of the natural duplex 1 was −23.0 kcal mol−1 at 25 °C and −18.3 kcal mol−1 at 37 °C, whereas that of modified duplexes 2 and 3 were −27.4 and −21.6 kcal mol−1 at 25 °C, −20.6 and −17.3 kcal mol−1 at 37 °C, respectively. The ΔH° value of duplexes 1–3 were −140.0, −196.8 and −129.1 kcal mol−1, while the ΔS° value were −392.6, −568.6 and −360.6 cal mol−1 K−1, respectively. These results indicated that the thermal destabilization of duplex 2 was related to the unfavorable entropy, whereas duplex 3 had a more favorable entropy owing to the incorporation of A4.
Thermodynamic parameters of duplexesa
Abbreviation of DNA/RNA duplex
Abbreviation of DNA or RNA
Sequenceb
ΔH° (kcal mol−1)
ΔS° (cal mol−1 K−1)
ΔG° (kcal mol−1)
Duplex 1
DNA 1
5′-d(TCTTTCTCTTTCCCTT)-3′
−140.0
−392.6
−23.0
(25 °C)
RNA 1
5′-r(AAGGGAAGAGAAAGA)-3′
−18.3
(37 °C)
Duplex 2
DNA 2
5′-d(TCTTTCTCTTTCCCTT)-3′
−196.8
−568.6
−27.4
(25 °C)
RNA 1
5′-r(AAGGGAAGAGAAAGA)-3′
−20.6
(37 °C)
Duplex 3
DNA 3
5′-d(TCTTTCTCTTTCCCTT)-3′
−129.1
−360.6
−21.6
(25 °C)
RNA 1
5′-r(AAGGGAAGAGAAAGA)-3′
−17.3
(37 °C)
Thermodynamic parameters were calculated based on the slop of a 1/Tmvs. ln(CT/4) plot.
T (bold) and T (italic) denote (S)-5′-C-aminopropyl-thymidine (A3) and (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine (A4), respectively.
Thermodynamic parameters were calculated based on the slop of a 1/Tmvs. ln(CT/4) plot.T (bold) and T (italic) denote (S)-5′-C-aminopropyl-thymidine (A3) and (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine (A4), respectively.Additionally, the base-discriminating ability of the nucleoside analog A3 and A4 in the DNA/RNA duplexes was examined utilizing duplexes 4–15 composed of DNAs 4–6 and RNAs 2–5. The results obtained are shown in Table 3 and S1† in detail. The ΔTm values, calculated from [Tm(duplexfully matched) − Tm(duplexmismatched)], of the natural duplexes ranged from −5.3 to −13.8 °C, whereas those of the modified duplexes 8–11 and 12–15 were ranged from −6.0 to −15.8 °C and from −4.0 to −14.6 °C, respectively. These results indicated that the base-discriminating abilities of nucleoside analogs A3 and A4 were the same as or better than those of the natural ones.
Base-discriminating ability of the nucleoside analog A1 and A2 in DNA/RNA duplexa
Tm/ΔTmb (°C)
Abbreviation of DNA/RNA duplex
X
Y
A
U
G
C
Duplex 4–7
Thymidine
47.4 (−)
35.9 (−11.5)
42.1 (−5.3)
33.6 (−13.8)
Duplex 8–11
A3
47.1 (−)
31.3 (−15.8)
41.1 (−6.0)
32.4 (−14.7)
Duplex 12–15
A4
47.3 (−)
34.0 (−13.3)
43.3 (−4.0)
32.7 (−14.6)
The sequences of duplexes are 5′-d(GGCTAX̲AATGTCG)-3′ for the DNA strand and 5′-r(CGACAUUY̲UAGCC)-3′ for the RNA strand, shown in Table S1 (see ESI) in detail. Underlined letters denote mismatched base pairs.
The Tm value were measured in 10 mM sodium phosphate buffer (pH 7.0) containing 100 mM NaCl. The concentrations the duplexes were 3 μM.
The sequences of duplexes are 5′-d(GGCTAX̲AATGTCG)-3′ for the DNA strand and 5′-r(CGACAUUY̲UAGCC)-3′ for the RNA strand, shown in Table S1 (see ESI) in detail. Underlined letters denote mismatched base pairs.The Tm value were measured in 10 mM sodium phosphate buffer (pH 7.0) containing 100 mM NaCl. The concentrations the duplexes were 3 μM.
Circular dichroism (CD) spectra
To determine how the incorporation of the nucleoside analogs A3 and A4 affected the thermal stability in comparison with the natural ones, CD spectra of the global structures of the duplexes 1–3 was measured, and the results are shown in Fig. 2. It is known that DNA/RNA duplexes adopt an A-type structure with negative and positive peaks at 210 and 270 nm, respectively, consistent with the result for the natural duplex 1. Since the CD spectra of the modified duplexes 2 and 3 was similar to that of duplex 1, it is considered that duplexes 2 and 3 form the A-type helix as well. However, the local maxima and minima of the modified duplexes were lower than those of the natural duplex. These variations in CD spectra indicated that there was a moderate change in the secondary structure of the modified duplexes due to the introduction of (S)-5′-C-aminopropyl modification, which was hence considered the reason for the moderate thermal destabilization observed for the modified duplexes. In addition, the CD spectra of duplexes 2 and 3 were highly similar to each other, resulting in the reverse effect with respect to the entropy values, as revealed by the calculation of the thermodynamic parameters, which was not related to the secondary structure of the duplexes.
Fig. 2
CD spectra of the native duplex 1 and the modified duplex 2 and 3 in a buffer containing 10 mM sodium phosphate (pH 7.0) and 100 mM NaCl at 25 °C. The concentration of the duplexes was 4 μM.
Molecular modeling study
In addition to the findings obtained based on the thermodynamic parameters and the CD spectra, the effect of the 2′-araF modification was assumed to play a vital role in determining the thermal stability. It is known that nucleosides adopt the north-type (N-type) C3′-endo pucker in DNA/RNA duplexes, which adopt the A-like helix structure. Therefore, it is assumed that there might be a change in sugar pucker resulting from the incorporation of the fluorine atom. To confirm this assumption, the features of nucleoside analogs A3 and A4, together with the natural thymidine, were calculated based on molecular dynamics (MD) simulation.The root mean square deviations (RMSDs) of heavy atoms were calculated relative to the reference structure attained from the energy minimization calculation (Fig. S1†). As the results shown in Table 4, the wide-range distribution of RMSDs of each nucleoside analog indicates that the fluctuation of the nucleoside analogs was more pronounced than that of the natural thymidine (Fig. S2†), indicating that nucleoside analogs require extra energy while forming duplexes. Next, the time-dependence of sugar puckering was calculated using pseudorotational phase angle (P) (Fig. S3†), based on which the percentages in the N-type sugar conformation of each nucleoside were calculated as well. As shown in Table 4, the percentage of A4 in the N-type sugar conformation was apparently higher than that of the natural one, suggesting that the 2′-araF modification might cause a shift of sugar pucker into the N-type sugar conformation. In addition, the introduction of the (S)-5′-C-aminopropyl modification has almost no effect on sugar pucker because of the inconspicuous disparity in the N-type abundance ratios of A3 and thymidine.
The features of thymidine, A3 and A4 obtained from MD simulation
Thymidine
A3
A4
RMSDa (Å)
0.55 (±0.20)
1.03 (±0.41)
0.83 (±0.33)
N-typeb (%)
12.5
14.1
19.5
The mean values and standard deviations of RMSD (Å).
The percentages of the N-type sugar conformation were calculated based on a criteria that considers the sugar conformation as the N-type if 0° < P < 90°and 270° < P < 360°.
The mean values and standard deviations of RMSD (Å).The percentages of the N-type sugar conformation were calculated based on a criteria that considers the sugar conformation as the N-type if 0° < P < 90°and 270° < P < 360°.Therefore, according to the molecular modeling study, the molecular fluctuation caused by the introduction of the (S)-5′-C-aminopropyl side chain would result in unfavorable conditions for the DNA/RNA duplex formation and thus cause thermal destabilization, whereas the 2′-arabinofluoro modification suppressed thermal destabilization by preferably adopting the N-type sugar pucker. These molecular modeling-based results were consistent with the findings on the thermodynamic parameters and the CD spectra.
Nuclease resistance of DNAs
Since natural nucleotides can easily be disassembled by various nucleases in the blood circulation, it is necessary to enhance the nuclease resistance of oligomers for therapeutic application. To investigate whether the nuclease resistance of DNA oligomers was improved by the incorporation of the nucleoside analogs A3 and A4, an enzymatic stability assay was performed utilizing bovine serum (BS). For this, the fluorescein-labeled single-stranded DNAs 7–9 (Table 5) were dissolved in OPTI-MEM, and subsequently incubated with 3% BS at 37 °C for the required time. Reaction mixtures were analyzed after various incubation times (0, 0.5, 1, 5, and 30 min, and 1, 3, 6 h) by 20% PAGE in the presence of 7 M urea, and then quantified by Luminescent Image analyzer LAS-4000 (Fujifilm). As shown in Fig. 3, the natural DNA 7 was completely degraded in BS within 3 h of incubation, whereas the modified DNAs 8 and 9 remained >80% and >90% of the complete strand, respectively, even after 6 h of incubation. Therefore, we could show that the incorporation of nucleoside analogs A3 and A4 could enhance the nuclease resistance of DNA oligomers. The experimental results were also consistent with previous studies.[15] It is presumed that the positive charges carried by the amino group of the aminopropyl side chain would mask the phosphodiester backbone of the oligomers and electrostatically impede approaching nucleases, thereby improving the nuclease resistance.
Fig. 3
PAGE analysis of single-stranded DNAs (8 μM) treated in OPTI-MEM containing 3% bovine serum. Fluorescein-labeled DNA 7–9 (300 pmol) were incubated in OPTI-MEM containing 3% bovine serum, and the reaction mixtures at various incubation time (0, 0.5, 1, 5, 30, 60, 180, 360 min) were analyzed by 20% PAGE containing 7 M urea.
F, T (bold) and T (italic) denote fluorescein, (S)-5′-C-aminopropyl-thymidine (A3) and (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine (A4), respectively.
The ability for RNase H activation
To investigate whether the modified duplexes containing the nucleoside analogs A3 and A4 have the same ability for RNase H activation as the natural ones, RNase H-mediated cleavage assays were performed utilizing RNase H purified from Escherichia coli. The natural duplex 16 was composed of DNA 1 and fluorescein-labeled RNA 6, and the modified duplexes 17 and 18 were composed of DNAs 2 and 3, respectively, as well as RNA 6 (Table 6). These duplexes were dissolved in a buffer containing 50 mM Tris–HCl (pH 8.0), 75 mM KCl, 3 mM MgCl2 and 10 mM dithiothreitol. Then, diluted RNase H solution (60 unit per L in H2O) was added, and subsequently the mixture was incubated at 37 °C for the required times (0, 1, 5, 15 and 30 min, and 1, 2, 4 h). The reaction mixtures were analyzed after various incubation times by 20% denaturing PAGE, and then quantified by Luminescent Image analyzer LAS-4000 (Fujifilm). As shown in Fig. 4, after 4 h of incubation, almost no complete strand of RNA 6 remained in the natural duplex 16, while in duplexes 17 and 18, the residual prevalence of complete RNA 6 was 88% and 46%, respectively. In other words, compared to the natural oligomer, the incorporation of nucleoside analog A3 significantly impeded RNase H-mediated cleavage of the RNA strand in the DNA/RNA duplex, while the impact of A4 was comparatively inconspicuous. According to these results, although the introduction of the (S)-5′-C-aminopropyl side chain could enhance the nuclease resistance of oligomers considerably, as desired, it might have inhibited the interaction between RNase H and the DNA/RNA duplexes as well. In addition, it was revealed that the combination with the 2′-araF modification could improve the ability to activate RNase H, which was in line with our expectations.
Fig. 4
PAGE analysis of DNA/RNA duplexes treated in a buffer containing 50 mM Tris–HCl (pH 8.0), 75 mM KCl, 3 mM MgCl2, 10 mM dithiothreitol, and RNase H from E. coli. The duplex 16–18 were incubated in a buffer containing 50 mM Tris–HCl (pH 8.0), 75 mM KCl, 3 mM MgCl2 and 10 mM dithiothreitol, and then diluted RNase H solution (60 unit per L in H2O) was added, and subsequently the mixture was incubated at 37 °C for the required time. The reaction mixtures at various incubation time (0, 1, 5, 15, 30, 60, 120, 240 min) were analyzed by 20% PAGE.
F , T (bold) and T (italic) denote fluorescein, (S)-5′-C-aminopropyl-thymidine (A3) and (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine (A4), respectively.However, it still remains the possibility that the functional properties of the modified duplexes might depend on the number and position of nucleoside analogs. For this reason, DNA 10–13, which contain only one analog A4 in different positions, were synthesized. Thermal stability of the DNA/RNA duplexes composed of DNA 10–13 and RNA 1 were evaluated and RNase H-mediated cleavage assays were also performed using the duplexes composed of DNA 10–13 and RNA 6. The Tm values and images of PAGE analysis were shown in Table S2 and Fig. S8,† respectively. The Tm values indicated that duplexes 19–22 would be thermally stable enough in physiological conditions, and duplexes 23–26 showed the same behavior as the natural duplex in RNase H-mediated cleavage assays. Meanwhile, it was revealed that the introduction positions of nucleoside analogs should be optimized, and the minimum length of gap, the distance between each nucleoside analog, for RNase H binding should be investigated in further studies.
Conclusions
In summary, this is the first time that the synthesis of (S)-5′-C-aminopropyl- and (S)-5′-C-aminopropyl-2′-arabinofluoro-modified thymidines have been accomplished and reported. The properties of the DNA oligomers containing these thymidine analogs were also evaluated. It was revealed that the (S)-5′-C-aminopropyl-2′-arabinofluoro-thymidine was thermally stable enough compared with the natural one, while the incorporation of (S)-5′-C-aminopropyl-thymidine caused thermal destabilization. The difference in thermal stability might have resulted from a moderate change in the secondary structure of the DNA/RNA duplexes, as well as a molecular fluctuation in nucleoside structures derived from the introduction of the (S)-5′-C-aminopropyl side chain, and the variation in sugar pucker derived from the 2′-arabinofluoro modification. Meanwhile, the incorporation of these nucleoside analogs could significantly enhance the nuclease resistance of DNA oligomers exposed to BS. Moreover, the DNA/RNA duplexes containing the (S)-5′-C-aminopropyl-2′-arabinofluoro-thymidine showed a superior ability to activate RNase H-mediated cleavage of the RNA strand. Hence, the (S)-5′-C-aminopropyl-2′-arabinofluoro-modified nucleoside analogs might be potential candidates for the application of RNase H-dependent antisense oligonucleotides as therapeutics, although the introduction position should be optimized, and their cytotoxicity and cell membrane permeability need to be evaluated in further studies.
Experimental section
General remark
All chemicals and dry solvents (DCE, DCM, DMF, MeCN, THF, and pyridine) were obtained from commercial sources and used without any further purification. Thin layer chromatography (TLC) was performed on silica gel plates precoated with fluorescent indicator with visualization by UV light or by dipping into a solution of 5% (v/v) concentrated H2SO4 in a mixture of p-anisaldehyde and methanol, and then heating. Silica gel (63–210 mesh) was used for column chromatography. 1H NMR (400 or 500 MHz), 13C {1H} NMR (101 MHz), 19F NMR (376 MHz), 31P NMR (162 MHz) were recorded on 400 or 500 MHz NMR equipment. Acetonitrile-d3, CDCl3, or DMSO-d6 was used as a solvent for obtaining NMR spectra. Chemical shifts (δ) are given in parts per million (ppm) from CDCl3 (7.26 ppm) for 1H NMR spectra and CDCl3 (77.2 ppm) for 13C NMR spectra. The abbreviations s, d, t, q, and m signify singlet, doublet, triplet, quadruplet, and multiplet, respectively. High resolution mass spectra (HRMS) were obtained in positive ion electrospray ionization (ESI-TOF) mode. The purity of each oligomers which were synthesized in this study was confirmed by HPLC and MALDI-TOF/MS.
9-Borabicyclo[3.3.1]nonane (9-BBN, 0.5 M in THF, 62 mL) was added dropwise to a solution of compound 1 (2.69 g, 5.16 mmol) in THF (27 mL) under argon atmosphere at 0 °C, and the reactant was stirred for 17 h at room temperature. Then, 2 N NaOH solution (46 mL) and 30% aqueous hydrogen peroxide solution (11 mL) were added dropwise to the reactant in turn. After stirring for 15 min, the reactant was quenched with 30% Na2S2O3 solution, and extracted with ethyl acetate and water; the organic layer was washed with brine, dried by Na2SO4, filtered and concentrated. The crude material was purified by column chromatography (80% ethyl acetate in hexane) to afford desired product 2 as a white solid (3.78 g, quant.). 1H NMR (400 MHz, CDCl3) δ: 1.08 (s, 9H), 1.42–1.62 (m, 4H), 1.85 (d, J = 0.92 Hz, 3H), 2.23 (m, 2H), 3.08–3.10 (m, 1H), 3.52–3.65 (m, 2H), 3.80 (t, J = 2.3 Hz, 1H), 4.46 (d, J = 1.8 Hz, 1H), 6.24 (t, J = 6.9 Hz, 1H), 7.37–7.47 (m, 6H), 7.49 (d, J = 0.92 Hz, 1H), 7.62–7.66 (m, 4H), 8.54 (s, 1H); 13C{1H} NMR (101 MHz, CDCl3) δ: 12.6, 19.2, 27.1, 29.2, 31.9, 39.9, 63.0, 71.2, 74.9, 87.5, 90.3, 111.0, 128.0, 130.2, 133.3, 133.6, 135.9, 136.0, 137.7, 150.5, 163.7; HRMS (ESI-TOF) m/z calcd for C29H38N2NaO6Si [M + Na]+, 561.2397, found 561.2371.
N,N-Diisopropylethylamine (DIPEA, 0.64 mL, 3.66 mmol) and 2-cyanoethyl N,N-diisopropylchlorophosphoramidite (CEPCl, 0.32 mL, 1.44 mmol) were added in a solution of compound 7 (0.500 g, 0.72 mmol) in THF (5.00 mL) in turn under argon atmosphere. The reactant was stirred for 1.5 h at room temperature, and then extracted with ethyl acetate and saturated NaHCO3; the organic layer was washed with brine, dried by Na2SO4, filtered and concentrated. The crude material was purified by column chromatography (60% ethyl acetate in hexane) to afford desired product 8 as a white solid (0.556 g, 86%). 31P NMR (162 MHz, CDCl3) δ: 148.9, 150.4; HRMS (ESI-TOF) m/z calcd for C45H55F3N5NaO9P [M + Na]+, 920.3587, found 920.2581.
(S)-5′-C-Allyl-3′-O-[(1,1-dimethylethyl)diphenylsilyl]-2′-β-fluoro-thymidine (11) and (R)-5′-C-allyl-3′-O-[(1,1-dimethylethyl)diphenylsilyl]-2′-β-fluoro-thymidine (12)
50% DMSO in DCM (44 mL) was added to a mixture of compound 9 (4.38 g, 8.79 mmol) and 1-(3-(dimethylamino)propyl)-3-ethylcarbodiimide hydrochloride (EDC·HCl, 6.75 g, 35.2 mmol) under argon atmosphere, and the mixture was stirred at 0 °C. Then, CHCl2CO2H (0.36 mL, 4.40 mmol) was added dropwise to the mixture and stirred for 3 h at 0 °C. The reactant was extracted with ethyl acetate and saturated NaHCO3; the organic layer was washed with brine, dried by Na2SO4, filtered and concentrated. The residue was dissolved in DCM (87 mL) under argon atmosphere, and the mixture was stirred at 0 °C. Then, allyltrimethylsilane (7.00 mL, 44.11 mmol) and boron trifluoride-ethyl ether complex (5.00 mL, 39.81 mmol) were added to the mixture. The reactant was stirred for 4 h at 0 °C, then extracted with ethyl acetate and saturated NaHCO3; the organic layer was washed with brine, dried by Na2SO4, filtered and concentrated. The crude material was purified by column chromatography (28.5% ethyl acetate in hexane) to afford desired product 11, of which the Rf value on TLC was little lower than that of another product 12, as a white solid (3.36 g, 71%). 1H NMR (400 MHz, CDCl3) δ: 1.10 (s, 9H), 1.41 (d, J = 5.5 Hz, 1H), 1.88 (s, 3H), 2.01–2.11 (m, 1H), 2.31–2.35 (m, 1H), 3.56 (s, 1H), 3.89 (dd, J = 2.1 Hz, 7.1 Hz, 1H), 4.54 (d, J = 15.1 Hz, 1H), 4.83 (dd, J = 2.3 Hz, 51.3 Hz, 1H), 5.10–5.22 (m, 2H), 5.67–5.78 (m, 1H), 6.30 (q, J = 2.5 Hz, 22.7 Hz, 1H), 7.18 (s, 1H), 7.40–7.49 (m, 6H), 7.65–7.69 (m, 4H), 8.29 (s, 1H).Another product 12 as a white solid (0.188 g, 4%).1H NMR (400 MHz, CDCl3) δ:1.10 (s, 9H), 1.56 (broad s, 1H), 1.86 (d, J = 1.2 Hz, 3H), 2.05–2.10 (m, 2H), 3.21–3.26 (m, 1H), 3.80 (t, J = 4.0 Hz, 1H), 4.39–4.46 (m, 1H), 4.96–5.11 (m, 3H), 5.58–5.68 (m, 1H), 6.30 (dd, J = 3.8 Hz, 18.6 Hz, 1H), 7.32 (s, 1H), 7.40–7.70 (m, 10H), 8.07 (broad s, 1H).
N,N-Diisopropylethylamine (DIPEA, 0.85 mL, 4.87 mmol) and 2-cyanoethyl N,N-diisopropylchlorophosphoramidite (CEPCl, 0.45 mL, 2.00 mmol) were added in a solution of compound 18 (0.713 g, 1.00 mmol) in THF (7.00 mL) in turn under argon atmosphere. The reactant was stirred for 1.5 h at room temperature, and then extracted with ethyl acetate and saturated NaHCO3; the organic layer was washed with brine, dried by Na2SO4, filtered and concentrated. The crude material was purified by column chromatography (60% ethyl acetate in hexane) to afford desired product 19 as a white solid (0.813 g, 89%). 19F NMR (376 MHz, CDCl3) δ: −119.58 (broad s); 31P NMR (162 MHz, CDCl3) δ: 152.1, 151.7; HRMS (ESI-TOF) m/z calcd for C45H54F4N5NaO9P [M + Na]+, 938.3493, found 938.3500.
Tetrabutylammonium fluoride (TBAF, 0.39 mL of a 1 M solution in THF) was added in a solution of compound 11 (0.095 g, 0.18 mmol) in THF (3.9 mL) under argon atmosphere, and the reactant was stirred for 2 h at room temperature. The reactant was extracted with ethyl acetate and saturated NaHCO3; the organic layer was washed with brine, dried by Na2SO4, filtered and concentrated. The residue was dissolved with pyridine (1.10 mL) and 1,3-dichloro-1,1,3,3-tetraisopropyldisiloxane (TIPDSCl2, 0.11 mL, 0.34 mmol) was added at 0 °C. The reactant was stirred for 24 h at room temperature and then extracted with chloroform and saturated NaHCO3; the organic layer was washed with brine, dried by Na2SO4, filtered and concentrated. The crude material was purified by column chromatography (100% chloroform) to afford desired product 20 as a white solid (0.082 g, 83%). 1H NMR (400 MHz, CDCl3) δ: 0.96–1.12 (m, 28H), 1.92 (d, J = 0.92 Hz, 3H), 2.38–2.56 (m, 2H), 3.61–3.69 (m, 1H), 4.05 (m, 1H), 4.45 (m, 1H), 5.07–5.24 (m, 3H), 5.75–5.85 (m, 1H), 6.23 (dd, J = 6.0 Hz), 7.27 (s, 1H), 8.92 (s, 1H); 13C{1H} NMR (101 MHz, CDCl3) δ: 12.5, 12.7, 12.8, 13.2, 13.6, 16.9, 17.0, 17.1, 17.49, 17.53, 37.4, 68.6, 73.3, 73.5, 79.7, 79.8, 81.0, 81.2, 93.5, 95.5, 110.8, 118.5, 133.7, 135.9, 150.4, 163.7; 19F NMR (376 MHz, CDCl3) δ: −120.62 (d, J = 92.42 Hz); HRMS (ESI-TOF) m/z calcd for C25H43FN2NaO6Si2 [M + Na]+, 565.2541, found 565.2532.
Tetrabutylammonium fluoride (TBAF, 0.80 mL of a 1 M solution in THF) was added in a solution of compound 12 (0.1817 g, 0.34 mmol) in THF (7.3 mL) under argon atmosphere, and the reactant was stirred for 2 h at room temperature. The reactant was extracted with ethyl acetate and saturated NaHCO3; the organic layer was washed with brine, dried by Na2SO4, filtered and concentrated. The half of residue was dissolved with pyridine (1.00 mL) and 1,3-dichloro-1,1,3,3-tetraisopropyldisiloxane (TIPDSCl2, 0.11 mL, 0.34 mmol) was added at 0 °C. The reactant was stirred for 24 h at room temperature and then extracted with chloroform and saturated NaHCO3; the organic layer was washed with brine, dried by Na2SO4, filtered and concentrated. The crude material was purified by column chromatography (100% chloroform) to afford desired product 21 as a white solid (0.067 g, 87%). 1H NMR (400 MHz, CDCl3) δ: 0.90–1.26 (m, 28H), 1.95 (d, J = 0.92 Hz, 3H), 2.27–2.34 (m, 1H), 2.50–2.55 (m, 1H), 3.51 (dd, J = 5.3 Hz, 9.2 Hz, 1H), 4.00–4.05 (m, 1H), 4.41 (dd, J = 5.3 Hz, 26.4 Hz, 1H), 4.98–5.16 (m, 3H), 5.88–5.99 (m, 1H), 6.05 (dd, J = 3.2 Hz, 21.1 Hz, 1H), 7.21 (dd, J = 1.4 Hz, 2.3 Hz, 1H), 8.74 (s, 1H); 13C{1H} NMR (101 MHz, CDCl3) δ: 12.7, 13.2, 13.3, 13.6, 16.9, 17.1, 17.6, 17.8, 17.9, 39.7, 74.8, 79.6, 79.9, 85.1, 110.4, 118.1, 134.0, 136.7, 150.2, 163.5; 19F NMR (376 MHz, CDCl3) δ: −117.58 (q, J = 23.11 Hz, 46.23 Hz); HRMS (ESI-TOF) m/z calcd for C25H43FN2NaO8Si2 [M + Na]+, 565.2541, found 565.2554.
Solid-phase oligonucleotide synthesis
The synthesis was carried out with a DNA/RNA synthesizer by the phosphoramidite method. After the synthesis, the RNA oligomers were cleaved from CPG beads and deprotected by treatment with concentrated NH3 solution/40% methylamine (1 : 1, v/v) for 10 min at 65 °C, while the DNA oligomers were treated with concentrated NH3 solution for 12 h at 55 °C. Notably, before the treatment of the NH3 solution, the CPG beads were treated with 10% dimethylamine in MeCN for 5 min followed by a rinse with MeCN to selectively remove cyanoethyl groups, if there are analogs introduced in oligomers. Then, 2′-O-TBDMS groups in RNA oligomers were removed by Et3N·3HF (125 μL) in DMSO (100 μL) for 1.5 h at 65 °C. The reaction was quenched with 0.1 M TEAA buffer (pH 7.0) and the mixture was desalted using a Sep-Pak C18 cartridge. The oligomers were purified by 20% PAGE containing 7 M urea to give highly purified oligonucleotides.
MALDI-TOF/MS analysis of ONs
The spectra were obtained with a time-of-flight mass spectrometer equipped with a nitrogen laser (337 nm, 3 ns pulse). A solution of 3-hydroxypicolinic acid (3-HPA) and diammonium hydrogen citrate in H2O was used as a matrix. Data of synthetic ONs: DNA 1: m/z = 4712.78 (calcd for C154H203N38O104P15 [M − H]−, 4713.60); DNA 2: m/z = 4941.01 (calcd for C166H231N42O104P15 [M − H]−, 4943.23); DNA 3: m/z = 5012.98 (calcd for C166H227F4N42O104P15 [M − H]−, 5014.16); DNA 5: m/z = 4044.76 (calcd for C131H168N50O77P12 [M − H]−, 4043.31); DNA 6: m/z = 4062.75 (calcd for C131H167FN50O77P12 [M − H]−, 4061.76); DNA 7: m/z = 5250.91 (calcd for C181H228N39O113P16 [M − H]−, 5252.09); DNA 8: m/z = 5479.15 (calcd for C193H256N43O113P16 [M − H]−, 5481.17); DNA 9: m/z = 5551.11 (calcd for C193H252F4N43O113P16 [M − H]−, 5552.23); DNA 10: m/z = 4787.83 (calcd for C157H209FN39O104P15 [M − H]−, 4789.59); DNA 11: m/z = 4787.83 (calcd for C157H209FN39O104P15 [M − H]−, 4790.08); DNA 12: m/z = 4787.83 (calcd for C157H209FN39O104P15 [M − H]−, 4789.44); DNA 13: m/z = 4787.83 (calcd for C157H209FN39O104P15 [M − H]−, 4790.53); RNA 1: m/z = 5298.85 (calcd for C160H193N80O100P15 [M − H]−, 5300.40); RNA 6: m/z = 5836.99 (calcd for C187H218N81O109P16 [M − H]−, 5838.61). The other ONs, which have been used in this study but not mentioned above, were synthesized in the previous study.[14]
Reverse-phase HPLC analysis of ONs
The purity of each synthesized ONs in this study was evaluated by reverse-phase HPLC (C18G 5 μm, 150 × 4.6 mm SS). Elution started from 100% buffer A, followed by a linear gradient to 35% buffer B in 20 min at a flow rate of 1.0 mL min−1 (buffer A: 5% MeCN in 0.1 M TEAA (pH 7.0); buffer B: 50% MeCN in 0.1 M TEAA (pH 7.0).)
Thermal denaturation study
The solution containing 3.0 μM DNA/RNA duplex in a buffer of 10 mM sodium phosphate (pH 7.0) containing 100 mM NaCl was heated at 90–100 °C and then cooled gradually to room temperature, and used for the thermal denaturation study. Thermally induced transitions were monitored at 260 nm with a UV/vis spectrometer fitted with a temperature controller in quartz cuvettes with a path length of 1.0 cm. The sample temperature was increased by 0.5 °C min−1. The thermodynamic parameters of the duplexes on duplex formation were determined by calculations based on the slope of a 1/Tmvs. ln(CT/4) plot, where CT (1, 3, 6, 12, 15, 21, 30 and 60 μM) is the total concentration of single strands.
CD spectroscopy
All CD spectra were recorded at 25 °C. The following instrument settings were used: resolution, 1 nm; response 1.0 s; speed, 50 nm min−1; accumulation, 10.The initial geometry of the natural thymidine, (S)-5′-C-aminopropyl-thymidine and (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine were modeled by a molecular modeling software, SYBYL.[25] The geometry of each nucleoside was optimized and the electrostatic potential (ESP) of each nucleoside was calculated by Gaussain16 at the HF/6-31G* level of theory.[26] The restrained ESP (RESP) charges were derived using the Antechamber RESP fitting procedure.[27] The generalized AMBER force field (GAFF) was used to natural and modified nucleosides.[28] Each of the nucleoside was solvated in a cubic box with TIP3P water model and periodic boundary conditions were applied.[29] The size of cubic box was determined to be 43.0 Å lengths based on the molecular size of (S)-5′-C-aminopropyl-2′-β-fluoro-thymidine. Chloride ions were added into the solution to neutralize the charge of the systems. The energy minimization calculation was performed using the SANDER module of AMBER18. The calculation of the solute molecule was carried out by 2500 circles of the steepest descent method followed by 2500 circles of the conjugate gradient method, applying a restraint force constant of 2.0 kcal (mol Å2) −1 to the nucleoside. Then, the calculation of the entire system was carried out by 5000 circles of the conjugate gradient method without any restraints. After that, the temperature of the system was heated gradually from 0 K to 300 K over a period of 100 ps of NVT dynamics with a restraint force constant of 2.0 kcal (mol Å2) −1 to the nucleoside, followed by 50 ps of NPT equilibration at 300 K and 1 bar pressure with the same restraint force. After the equilibration phase, 500 ns MD simulations were performed in an NPT ensemble at 300 K and 1 bar pressure without any restraints. During all MD simulations, a time step of 2 fs was used. The cutoff distance towards the van der Waals force was 10 Å, the SHAKE algorithm was used to constrain all bonds involving hydrogen atoms, and the Particle Mesh Ewald (PME) method was used for long-range electrostatic interaction.[30,31] The temperature was regulated by Langevin dynamics with a collision frequency of 2.0 ps−1.[32] All the MD simulations were performed by the PMEMD module of AMBER18.[33] Trajectory analysis was done for 50 000 snapshots, every 10 ps, sampled during the full 500 ns simulation.
Nuclease resistance of single-stranded DNA
Fluorescein labeled DNAs (300 pmol) were dissolved in OPTI-MEM (37.5 μL) and used for the serum stability test. 1.2 μL bovine serum was added, and subsequently the mixture was incubated at 37 °C for the required time. Aliquots of 2.4 μL were diluted with a stop solution (10 mM EDTA in formamide 10 μL). Samples were subjected to electrophoresis in 20% PAGE containing 7 M urea and quantified by Luminescent Image analyzer LAS-4000 (Fujifilm).
RNase H assay
The DNA/RNA duplexes used for RNase H assay were prepared by mixing the DNAs (600 pmol) with fluorescein labeled complementary target RNA (3000 pmol) in 75 μL of 50 mM Tris–HCl (pH 8.0) containing 75 mM KCl, 3 mM MgCl2 and 10 mM dithiothreitol, followed by heating at 90–100 °C for 5 min and cooling gradually to room temperature. Then, 70 μL diluted RNase H solution (60 unit per L in H2O) was added, and subsequently the mixture was incubated at 37 °C for the required time. Aliquots of 5 μL were diluted with 100% formamide (10 μL). Samples were subjected to electrophoresis in 20% PAGE containing 7 M urea and quantified by Luminescent Image analyzer LAS-4000 (Fujifilm).
Authors: Merrill D Benson; Márcia Waddington-Cruz; John L Berk; Michael Polydefkis; Peter J Dyck; Annabel K Wang; Violaine Planté-Bordeneuve; Fabio A Barroso; Giampaolo Merlini; Laura Obici; Morton Scheinberg; Thomas H Brannagan; William J Litchy; Carol Whelan; Brian M Drachman; David Adams; Stephen B Heitner; Isabel Conceição; Hartmut H Schmidt; Giuseppe Vita; Josep M Campistol; Josep Gamez; Peter D Gorevic; Edward Gane; Amil M Shah; Scott D Solomon; Brett P Monia; Steven G Hughes; T Jesse Kwoh; Bradley W McEvoy; Shiangtung W Jung; Brenda F Baker; Elizabeth J Ackermann; Morie A Gertz; Teresa Coelho Journal: N Engl J Med Date: 2018-07-05 Impact factor: 91.245
Authors: Wen Shen; Cheryl L De Hoyos; Hong Sun; Timothy A Vickers; Xue-Hai Liang; Stanley T Crooke Journal: Nucleic Acids Res Date: 2018-03-16 Impact factor: 16.971