1. Introduction
Therapeutic oligonucleotides (ONs) provide an opportunity for treating serious, life-threatening diseases with limited options using traditional low molecular weight molecules and antibody drugs [
1]. Although there are many ongoing clinical trials and some ON drugs [
2], ON therapeutics has not quite lived up to the early expectations, especially regarding the time frame for development. There has been a substantial number of promising results [
1,
2] but it is recognized that the most severe limitation in ON therapeutics is their delivery (or rather lack thereof) to the chosen site of action [
3,
4].
Moreover, antisense oligonucleotide (AON) action based on a 1:1 binding stoichiometry may require a relatively high dosage of ON to silence the target mRNA transcript effectively [
5]. Efficacy of therapeutic ONs is more readily achieved if catalytic cleavage of the target RNA is obtained, which can occur if native enzymes (e.g., RNAse H for antisense and RNA-induced silencing complex (RISC) for siRNA) can be recruited. Chemically modified antisense oligonucleotides often do not retain the ability to activate native enzymes. It is also not completely clear if in some circumstances the activity and/or concentration of the host enzymes can become limiting.
ON-based artificial nucleases (OBANs) [
6,
7,
8,
9,
10,
11] have the potential to affect RNA levels and thereby also diseases, while at the same time being independent of native enzymes with the freedom to use many available ON modifications. RNA sequence and secondary structures affect the reactivity of the phosphodiester bonds and studies have been done to evaluate structural parameters to support the development of sequence-specific RNA cleaving agents that mimic natural restriction enzymes [
12,
13]. RNA oligonucleotides are considerably more reactive in the single stranded form (also as bulges or loops) than as a duplex, and the cleavage of the phosphodiester bonds has a strong dependence on the sequence [
14,
15,
16]. Furthermore, the molecular environment within the RNA bulge, including the size, bulge sequence, stacking of the nucleobases, and orientation of the phosphodiester bonds within, has a considerable influence on the overall reactivity of the phosphodiester bonds. It is also preferable to cleave the target RNA in a central part of the RNA sequence so that RNA fragments can readily dissociate [
17,
18,
19].
We have designed artificial nucleases based on the concept that the modified oligonucleotide part of the artificial nuclease recognizes the target RNA upon hybridization and forces the RNA to form a bulge. This concept enabled the development of artificial enzymes capable of cleaving target RNA sequences with turnover [
20,
21,
22]. Previous efforts have also resulted in the development of peptide nucleic acid-based artificial nucleases (PNAzymes) [
23,
24,
25,
26]. PNAzymes using Cu
2+ as a co-factor were shown to act as RNA restriction enzymes with high site selectivity and half-lives of cleavage down to 15 to 30 min, at a 1:1 ratio of PNAzyme to RNA [
27,
28]. It was also found that the rate of cleavage with the Cu
2+ PNAzymes was highly dependent on the sequence of the non-paired bulge sequence that were formed in the target RNA [
27]. PNAzymes with Zn
2+ as a co-factor was shown to cleave target RNA within 3- and 4-nt bulges with a half-life as low as 11 h [
26]. However, these zinc-based systems were not studied with respect to how the bulge sequence may affect the rate and site selectivity of the RNA cleavage. In the present study, we report on how the size and sequence of formed 3- and 4-nt bulges affect the rate and site selectivity of RNA cleavage by a Zn
2+-dependent PNAzyme. This study was motivated by the fact that we found the bulge sequence to have a substantial effect on the rate of cleavage in Cu
2+-based systems [
27]. In our current development towards Zn
2+-based artificial nucleases that are more biocompatible than copper-based systems, it is essential to investigate different RNA sequences in order to determine how this can effect PNAzyme action.
2. Results and Discussion
Metal ions, the stacking interactions within an RNA/PNA duplex, and the formed RNA bulge may have a structural role affecting the reactivity of phosphodiester bonds within the bulge, making the activity of PNAzyme complexes quite sensitive to the sequence of the RNA molecule and also dependent on the metal co-factor used. For the evaluation of structural factors in the bulge, we chose to utilize a PNAzyme [
26] that has been shown to cleave a model of the Leukemia-related bcr/abl mRNA in the presence of zinc ions. This PNAzyme cleaves the target RNA in a catalytic fashion with a half-life of down to 11 h when an AAAA bulge is formed [
26]. Scission of the bcr/abl mRNA takes place in several positions within the bulged out RNA. As a starting point for this study, we selected systems where the RNA/
PNAzyme1 complex will form 4-nt bulges in the RNA part (
Figure 1).
Four-bulge forming target RNAs: We performed a relatively systematic variation of the initial 4-adenosine bulge forming RNA substrate (
RNA1) and investigated how the bulge sequence influenced the rate and selectivity of the RNA cleavage reaction conducted by the Zn
2+-dependent
PNAzyme1. Thus, 21 additional RNA sequences (
RNA2–
22,
Table 1) were separately treated with
PNAzyme1. The reactions were performed using equimolar amounts of
PNAzyme1 and RNA (both at 4 μM) in the presence of 100 μM Zn
2+ (pH 7.4 and 37 °C) following conditions described in previous studies [
26]. Aliquots were withdrawn from the reaction mixtures at 3 and 24 h. Variation of the 4-adenosine bulge sequence by replacing one adenine base with one of the other three natural nucleobases (
RNA2–
4,
RNA5–
7,
RNA11–
13, and
RNA14–
16) revealed a clear dependence of the cleavage rate on the bulge composition (
Table 1). Further variation of the bulge sequence (
RNA17–
22) confirmed the substantial influence of the bulge sequence on the rate of cleavage.
Replacement of the adenosine at position 5, adjacent to the GT wobble base pair (
RNA2–
4), decreased the rate of cleavage regardless of the base substituted. The main cleavage site remained the same and was located between positions 5 and 6, i.e., the nucleotide at position 6 provided the 2′-hydroxyl nucleophile for the phosphodiester cleavage reaction and resulted in the formation of a 2′, 3′-cyclic phosphodiester and release of the 3′-end fragment with a free 5′-hydroxyl. Replacing the adenosine at position 6 (
RNA5) with guanosine was somewhat tolerated, while the cleavage rate dropped substantially if the purine bases were replaced by a pyrimidine in this position (
RNA6–
7). The change to a guanosine at position 5 was also better tolerated than the change to a pyrimidine base (
RNA2–
4). These effects were confirmed by the data obtained with
RNA8–
10 and suggest that either stacking between base 5 and 6 and/or similar interaction with the neocuproine affected the rate of cleavage, and that stacking to the G in the wobble base pair may be important. This is in clear contrast to the results obtained for position 7, where a base change to a pyrimidine (
RNA12–
13) in the target RNA led to an even higher rate of cleavage (change to G (
RNA11) in this position did not affect cleavage efficiency much). The cleavage rate is less sensitive to variation at position 8 in the bulge, although G and A seem to be less favoured.
RNA13 (AUAA),
RNA15 (CAAA),
RNA17 (UUAA), and
RNA20 (AUGA) gave the most cleavage after 24 h, closely followed by a few more sequences. The sequence preferences with respect to rate of cleavage are summarized in
Figure 2.
The best site selectivity in the 4-nt bulge systems was achieved with the bulge sequence -ACAA- (
RNA12), showing more than 70% selectivity for cleavage at one site (
Figure 3), which was an improvement compared to the AAAA bulge system [
26]. The ACAA bulge was also cleaved at a higher rate, though not as efficiently as, e.g., the AUAA bulge. Several similarities in the preferences for bulge sequences can be seen when comparing the rate of cleavage with Cu
2+-based PNAzymes [
27]. Although differences were larger for the latter, the zinc- and copper-based PNAzymes seem to show a similar nucleobase-related rate dependence at positions 5–7 of the RNA bulge and AUAA was also the most favoured sequence for the Cu
2+-based cleavage. In contrast, position 8 suggested a preference for A in copper-based systems, while the zinc-based systems equally accept pyrimidine bases.
Three-bulge forming target RNAs: OBANs can display different preferences for different RNA bulge sizes, and the sequence dependence has not been studied for 3-nt bulges. To explore this further, we investigated the influence of the 3-nt bulge sequence on RNA target cleavage by the Zn
2+-based
PNAzyme1. The starting point was again an all-adenosine bulge where each one of the adenines were systematically replaced by other nucleobases. This revealed an even higher dependence of the rate of cleavage on the bulge composition than that observed for the 4-nt bulges (
Table 2).
The 3A bulge RNA (
RNA23,
Table 2) was cleaved at a lower rate than the corresponding 4A bulge target (
RNA1,
Table 1). The results of this experiment, taken on their own, could indicate that the smaller size of the RNA bulge affects the cleavage efficiency negatively. However, as further experiments with other sequences were continued, this became less evident. For the nucleotide adjacent to the wobble base pair (position 5) adenine is clearly the preferred base (cf
RNA23 and
RNA27,
RNA24 and
RNA30,
RNA26 and
RNA37). For the central position (position 6) the presence of a uracil base in the RNA target appears to give a higher rate of cleavage (
RNA24). This is the position that provides the 2′-hydroxyl nucleophile at the main cleavage reaction site, as evidenced by the main fragments formed. Changing the adenine to guanine (
RNA25) gives a smaller enhancement, while change to a cytosine (
RNA26) led to a slightly decreased rate.
Replacing the adenine in position 7 (
RNA28) led to a notable increase in the rate of cleavage. Out of the
RNA23–
28 target RNAs, the AUA bulge (
RNA24) gave the best site selectivity with more than 80% of the cleavage taking place at the phosphodiester connecting A-5 and U-6. Further investigations were continued by systematically changing the bases within the AUA bulge region of
RNA24. Changing the base at position 5 neighbouring the GT wobble base pair generated a decrease in the cleavage activity (
RNA29–
31). This correlated with the 4-nt bulge systems where adenine was the preferred choice in the corresponding position (
Figure 2). Replacing adenine in position 7 with any of the other three standard bases pushed the rate of cleavage upwards. These 3-nt bulge PNAzyme complexes (
RNA32–
34) displayed rates comparable to the fastest 4-nt bulge (AUAA,
RNA13) system. Three additional RNAs (
RNA35–
37) were also evaluated, but cleavage activity decreased. The two RNA targets that seemed to be cleaved most efficiently were those with either GUA (
RNA32) or UUA (
RNA34) sequences in the bulge, as is schematically shown in
Figure 4.
There is some, but not complete, correlation between the cleavage rate and site selectivity of
PNAzyme1-mediated cleavage reactions. The highest site selectivity of 3-nt bulge systems studied was achieved with the bulge sequences AUA (
RNA24), GUA (
RNA32), and UUA (
RNA34), which resulted in over 80% relative site selectivity (e.g., GUA in
Figure 5). This is a significant improvement compared to other reported 3-nt [
26] systems and to the 4-nt Zn
2+-based systems in
Table 1.
3. Materials and Methods
Peptide nucleic acid monomers, Fmoc-PNA-A(Bhoc)-OH, Fmoc-PNA-G(Bhoc)-OH, Fmoc-PNA-C(Bhoc)-OH and Fmoc-PNA-T-OH, were purchased from Link Technologies Ltd. (Glasgow, UK). Rink Amide resin (ChemMatrix, 0.47 mmol/g) was purchased from Biotage (Uppsala, Sweden). Solvents and reagents for solid-phase synthesis; Iris Biotech GmbH (Marktredwitz, Germany) (2,3-diaminopropionic acid (Fmoc-L-Dap(Mtt)-OH), Fmoc-L-Lys(Boc)-OH, trifluoroacetic acid (TFA), triisopropylsilane (TIS), dichloromethane (DCM) and
N-methylpyrrolidone (NMP)), Acros Organics (NJ, US) (diethyl ether), Sigma-Aldrich (St. Louis, MO, USA) (piperidine, phenyl chloroformate, 4-methylmorpholine), Merck-Millipore (Burlington, MA, USA) (ethyl cyano(hydroxyimino) acetate (Oxyma) and
N,
N’-diisopropylcarbodiimide (DIC)), Alfa Aesar (Haverhill, MA, USA) (2,6-lutidine), VWR (Radnor, PA, USA) (acetonitrile, MeCN), and Applied Biosystems (Foster City, CA, USA) (acetic anhydride). 5-amino-2,9-dimethyl 1,10 phenanthroline was synthesized as previously reported [
30]. All chemicals used in the RNA cleavage assays were of molecular biology grade. All RNA oligomers were purchased from Dharmacon (Lafayette, CO, USA). RNA oligomers 5′-AGAGUUCAAAGCCC-3′ (
RNA23) and 5′-AGAGUUCAUAGCCC-3′ (
RNA24) were bought as IE-HPLC-purified and desalted. All other RNAs were deprotected according to manufacturer′s protocol (2’-ACE protecting groups), purified by the Ion-Exchange High Performance Liquid Chromatography (IE-HPLC) and desalted by Reversed Phase HPLC (RP-HPLC). IE-HPLC purifications were performed using an analytical Thermo Scientific DNAPac PA-100 BioLC (4 × 250 mm) column (Waltham, MA, USA). A linear gradient of 0–35% buffer B over 15 min was used with a flow rate of 1.5 mL/min at 60 °C, (Buffer A) 20 mM NaOAc in 30% MeCN/aq., (Buffer B) 20 mM NaOAc, 0.4 M LiClO
4 in 30% MeCN/aq. UV detection was carried out at 260 nm. Purified RNAs were desalted by reverse-phase HPLC on a Jasco HPLC system using the Supelco Discovery BIO Wide Pore C18-5.5 μm (250 × 4.6 mm) column (Sigma-Aldrich). A linear gradient of 0–37% buffer B over 20 min was used with a flow rate of 1 mL/min at 50 °C, (Buffer A) 50 mM triethylammonium acetate (TEAA) in water (pH 6.5) and (Buffer B) 50 mM TEAA in water (pH 6.5)-acetonitrile (1:1,
v/
v). Purified and desalted RNAs were lyophilized three times before use and stored as frozen solution. Concentrations of RNA sequences were determined by UV absorption at 260 nm on a Varian Cary 300 UV-vis dual beam spectrophotometer (Varian) and calculated from extinction coefficients obtained by the nearest neighbour approximation [
31].
3.1. Synthesis of Peptide Nucleic Acid-Based Artificial Ribonuclease (PNAzyme1)
Chimeric PNA oligomer was prepared on a Biotage Initiator+ Alstra microwave peptide synthesizer. The sequence was assembled automatically on Rink Amide ChemMatrix resin (0.47 mmol/g) on a 10-μmol scale in a 5-mL reactor vial using 9-fluorenylmethoxycarbonyl (Fmoc) chemistry under inert gas (N2). Fmoc deprotection was performed at room temperature (RT) in two stages by treating the resin with piperidine-NMP (1:4, v/v) for 3 min followed by piperidine-NMP (1:4, v/v) for 10 min. The resin was then washed with NMP (×4). PNA couplings were performed using 4 eq. of PNA monomer (or Lys or Dap monomers), 4 eq. Oxyma, and 4 eq. DIC in NMP. A coupling time of 6 min at 75 °C (microwave) was employed and the resin was washed with NMP (×4). This was followed by a capping step using NMP-lutidine-acetic anhydride (89:6:5, v/v/v) for 2 min and then washing with NMP (×4). After the synthesis cycles were completed, the resin was washed with NMP (×5) and DCM (×6), and then dried. Prior to post-conjugations, the methyltrityl (Mtt) protection group of 2,3-diaminopropionic acid was removed by subjecting the solid-supported PNA to 1% trifluoroacetic acid (TFA) in dichloromethane (DCM) for 5 × 1 min, followed by washing with DCM and NMP. 5-amino-2,9-dimethyl-1,10-phenanthroline (10 mg, 45 μmol) was dissolved in DCM (200 μL). 4-methylmorpholine (NMM, 10 μL, 90 μmol) and phenyl chloroformate (5 μL, 40 μmol) was added to the solution, and the reaction mixture was agitated for 1 h. The mixture was partitioned between DCM and H2O (200 μL). The organic phase was added to the solid support bound PNA with 50 μL NMM and reaction was incubated overnight at room temperature. The support was filtered and washed with NMP and DCM. PNA was cleaved from the solid support using a cocktail of TFA-H2O-TIS (95:2.5:2.5, v/v/v/v) for 1.5 h at room temperature. The PNA product was evaporated to dryness by a nitrogen flow, diluted with deionized (MilliQ) water and evaporated to dryness under reduced pressure. The crude product was partitioned between diethyl ether and water (×3), and the water phase was evaporated to dryness.
RP-HPLC purification and analysis was carried out on a Ascentis Express Supelco Peptide ES-C18 column (2.7 μm, 150 × 4.6 mm) with a linear gradient elution of 27% B over 30 min at 60 °C, using a flow rate of 1.0 mL/min and UV detection at 260 nm. The following solvent system was used: solvent A, water containing 0.1% TFA; solvent B, 50% MeCN: water containing 0.1% TFA (1:1; v/v). Collected products were evaporated to dryness and lyophilized from water (×3). Mass spectrometry of the collected products was performed on a ESI-TOF mass spectrometer in positive ion mode using a solution of acetonitrile-water (1:1, v/v) containing 0.1% formic acid. ES-MS m/z: mass calculated for PNAzyme1 C145H180N72O38 [M]+4 calcd 3539, found 3539.
3.2. RNA Cleavage Assay Using Anion Exchange (IE) HPLC
RNA cleavage reactions were carried out in sealed tubes immersed in a thermostated water bath (37 °C). These experiments were performed in 10 mM HEPES buffer at pH 7.4, containing 0.1 M NaCl and 100 μM Zn
2+ (110 μM Zn
2+ , 10 μM ethylenediaminetetraacetic acid (EDTA), resulting in 100 μM accessible zinc ions).
PNAzyme1 is known to cleave phosphodiester bonds efficiently in the presence of copper ions [
27]. EDTA was used to eliminate possible copper contamination. All reactions were performed in equimolar concentrations (4 μM) of substrate RNA and PNAzyme. The samples were prepared at room temperature by adding the appropriate amounts of deionized (MilliQ) water, buffer concentrate, substrate RNA, metal ion stock solution, and PNAzyme to achieve the final concentrations. After the addition of all components, at 3-h and 24-h time intervals, 70-μL aliquots were withdrawn from the reactions and quenched by adding 30 μL of a solution of 1 mM EDTA in 30% aq. MeCN. The samples were analysed by anion exchange HPLC (IE-HPLC). IE-HPLC analysis was carried out on a Dionex NucleoPac PA-100 (4 × 250 mm) column with a linear gradient elution of 0–45% buffer B over 30 min at 60 °C. A flow rate of 1.5 mL/min was used and UV detection was carried out at 260 nm. The following solvent system was used: (A) 20 mM NaOAc in 30% aq. MeCN and (B) 20 mM NaOAc, 0.4 M LiClO
4 in 30% aq. MeCN. Cleavage of RNA substrates were obtained by the quantification of the remaining RNA and the sum of all the formed fragments detected in the IE-HPLC analysis. Further examples of chromatograms can be found in
Supplementary materials (
Figures S1 and S2).
3.3. Determination of RNA Cleavage Sites
RNA cleavage reactions were performed as above. The main RNA fragments were collected from IE-HPLC runs and desalted as described earlier. Mass spectrometry of the collected fractions was performed in negative ion mode using a solution of MeCN–water (1:1, v/v). ES-MS m/z: RNA12 (-ACAA-) Fragment1 C47H61N19O32P4 [M]−2 calcd 1528, found 1527, Fragment2 (2′,3′-cyclic phosphate) C96H118N40O68P10 [M]−2 calcd 3230, found 3229. RNA13 (-AUAA-) Fragment1 C47H61N19O32P4 [M]−2 calcd 1528, found 1527, Fragment2 (2′,3′-cyclic phosphate) C96H117N39O69P10 [M]−2 calcd 3231, found 3230. RNA24 (-AUA-) Fragment1 C47H61N19O32P4 [M]−2 calcd 1528, found 1527, Fragment2 (2′,3′-cyclic phosphate) C86H105N34O63P9 [M]−2 calcd 2902, found 2901. RNA32 (-GUA-) Fragment1 C47H61N19O32P4 [M]−2 calcd 1528, found 1527, Fragment2 (2′,3′-cyclic phosphate) C86H105N34O64P9 [M]−2 calcd 2918, found 2917.