Guest Editor: A. S. Sidorenko Beilstein J. Nanotechnol.2016,7, 1434–1446.https://doi.org/10.3762/bjnano.7.135 Received 21 May 2016,
Accepted 23 Sep 2016,
Published 12 Oct 2016
Osmium tetroxide 2,2’-bipyridine (OsBp) is known to react with pyrimidines in ssDNA and preferentially label deoxythymine (T) over deoxycytosine (C). The product, osmylated DNA, was proposed as a surrogate for nanopore-based DNA sequencing due to OsBp’s “perfect” label attributes. Osmylated deoxyoligos translocate unassisted and measurably slow via sub-2 nm SiN solid-state nanopores, as well as via the alpha-hemolysin (α-HL) pore. Both nanopores discriminate clearly between osmylated and intact nucleobase; α-HL was also shown to discriminate between osmylated T and osmylated C. Experiments presented here confirm that the kinetics of osmylation are comparable for short oligos and long ssDNA and show that pyrimidine osmylation is practically complete in two hours at room temperature with less than 15 mM OsBp. Under the proposed labeling conditions: deoxyoligo backbone degradation measures less than 1/1,000,000; false positives such as osmylated deoxyadenine (A) and osmylated deoxyguanine (G) measure less than 1/100,000; false negatives, i.e., unosmylated C measure less than 1/10,000; and unosmylated T must measure substantially lower than 1/10,000 due to the 27-fold higher reactivity of T compared to C. However, osmylated C undergoes degradation that amounts to about 1–2% for the duration of the labeling protocol. This degradation may be further characterized, possibly suppressed, and the properties of the degradation products via nanopore translocation can be evaluated to assure base calling quality in a DNA sequencing effort.
At the time the human genome was completed in 2001 [1,2], it was becoming clear that the relationships between genes, noncoding regions, epigenetic modifications, and the transcriptome were intertwined and led to health vs disease in complex patterns. A large effort in academia, government, and industry pushed DNA sequencing to next-generation technologies [3] that flourished exponentially due to technological advancements in parallelization, miniaturization, and bioinformatics. DNA sequencing was and still is an extension of the Sanger sequencing-by-synthesis (SBS) method [4]. Nevertheless, SBS limitations [5], and the need for single-molecule sequencing at a fast speed, low cost, and accurate base-calling has led to the exploration of nanopore-based solutions and other platforms [6,7].
Physical scientists embraced nanopore technology for a plethora of applications [8-10]. The advantage of protein pores is that they are well defined and highly reproducible, whereas the advantage of solid-state nanopores is that they can be manufactured to desirable dimensions. Protein pores such as α-HL, a modified version of the Mycobacterium smegmatis porin A (MspA), and the Ph29 connector channel have been investigated as single-molecule sensing devices for ssDNA, RNA, dsDNA, and proteins [11-13]. The concept of nanopore-based sequencing, patented in 1998 [14], is based on applying a potential across an open pore embedded within an insulating membrane that separates two compartments filled with electrolyte. Influenced by the electric field, the electrolyte ions traverse the pore producing a constant current. Also led by the applied field, a nucleic acid in one compartment moves through the pore to the other compartment and obstructs the current in a sequence information-rich manner [15-17].
Nanopore-based sequencing of ssDNA relies on current modulation attributed to the single-file passing of the nucleotide bases A, G, C, or T through the pore. Typically, nanopores interact with a short sequence, and therefore the observed current obstruction may not be attributed to a single base [11,18]. From a chemical point of view, the bases are not that different from each other and yield comparable current levels, as shown by experiments where the ssDNA was immobilized within the pore [19-21]. Discrimination is further reduced when the nucleic acid traverses the pore in speeds that correspond microseconds per base − speeds that are too fast for detection by current instrumentation [11]. To remedy the situation, work mostly with α-HL and MspA has yielded engineered versions thereof with improved DNA sensing capabilities [12,22,23], as well as the incorporation of Ph29 polymerase, or other appropriate enzymes at the lip of the nanopores. This acts to slow down translocation by binding to the strand and moves it one base at a time through the pore [24,25].
To resolve the issues of similar chemical structure and unreadable speed of translocation [26], we considered a different approach, namely selective labeling of ssDNA. Historically, selective labeling was evaluated in the 1960s in an effort to image DNA conjugated with contrast agents by electron microscopy in order to elucidate sequence from the images [27-29]. This was tested within this work and it was found that one of the contrast agents, a 1:1 mixture of osmium tetroxide and 2,2’-bipyridine, exhibited extraordinary properties, and was considered to be the “perfect” label (Scheme 1). These attributes are (i) room temperature reactivity at mM concentration, (ii) no loss of label upon standing, (iii) superior selectivity for T over C by a factor of 27, (iv) undetectable degradation of the DNA strand under the labeling conditions, and (v) undetectable reactivity towards the purines (false positives). Most importantly OsBp reactivity towards T or C was found to be independent of composition, sequence, length, and secondary structure of the DNA strand, as shown in [30-32] and also the results herein. For example, oligothymidylates, such as T8, T15 and T20, exhibit fully labeled product formation that is identical to the rate of labeling of an oligoadenylate with one T or, even to osmylation of the monomer, the deoxythymidylic acid triphosphate, 5’dTTP [30]. The across-the-board equal reactivity leads to predictable osmylation of any DNA without prior sequence knowledge. Conjugation of OsBp to a pyrimidine leads to a new chromophore that absorbs in the range of 300 to 320 nm, a range where DNA does not absorb. This last feature enabled the development of a quality control UV–vis assay to confirm the extent of labeling [30-32].
It was hypothesized that a nanopore of suitable size would differentiate between osmylated and bare base due to their major differences in size and structure. Indeed it was recently experimentally confirmed that both 1.6 nm wide SiN nanopores [33] and α-HL [34] allow translocation of an osmylated oligo at strikingly slow, measurable speeds. Current recordings from experiments using 20 nucleotide long deoxyoligos with sequence A10XA9, where X = A, C, T, 5-MeC, U (deoxyuridine), solidly support the conclusion that α-HL discriminates between osmylated and intact base, as well as among the different osmylated pyrimidines tested. This last observation yields a tentative nanopore-based methodology whereby the osmylated target strand may provide sequencing information on T and C, and the osmylated complementary strand may provide information on A and G.
Undoubtedly a labeled nucleic acid, used as a surrogate for sequencing via electron microscopy or nanopore technology, is required to meet stringent qualifications. These qualifications should ascertain a “perfect” match between intact and osmylated, generally speaking “labeled” nucleic acid, so that sequencing of the labeled nucleic acid can accurately reproduce the sequence of the target polymer. These properties, defined elsewhere [26], include: (i) high solubility in water/electrolyte, (ii) less than 1/10,000 in false positives, (iii) less than 1/10,000 in false negatives, and (iv) negligible degradation of the DNA backbone. In addition, any side reaction that occurs during labeling should be identified and quantified. Experiments are described here to show that under the proposed osmylation conditions, these qualifications are met.
Results and Discussion
Earlier studies identified a chromophore, resulting from the formation of the pyrimidine/OsBp bond, which absorbs in the range of 300 to 320 nm where DNA has practically no absorbance [30-32]. These studies also selected two wavelengths and their absorbance ratio (R = A(λ = 312 nm)/A(λ = 272 nm), i.e., the absorbance at 312 nm divided by the absorbance at 272 nm) as optimal for monitoring the osmylation of either T or T + C. Osmylation can be monitored at 312 nm directly, but determining R cancels sampling errors. R is a direct measure of the extent of osmylation and can be obtained spectrophotometrically after purification/removal of the excess label. In the presence of OsBp, the values of R are directly obtained from the DNA peak when the reaction mixture is analyzed and monitored by high performance liquid chromatography (HPLC) or capillary electrophoresis (CE). It turns out that both analytical tools clearly resolve DNA from the OsBp label, which appears as a single peak with the methods used in this report. For the more typical cases where the oligo/DNA has 5 or more pyrimidines, the resolution between intact and osmylated nucleic acid is lost, and the substrate/product appears as a single broad peak. Figure 1 illustrates R for the reaction of a 32 nucleotide long oligo with 12.6 mM OsBp in water at 25 °C, as monitored by HPLC. Upon HPLC analysis, OsBp elutes very early, and the osmylated product elutes late, allowing the visualization of complete osmylation (plateau), as well as of any side reactions (additional peaks) that could be attributed to degradation products. Under these conditions, the reaction is quite fast and all pyrimidines have reacted after about 40 min, as evidenced by the plateau. Moreover, no new peaks accumulate over time, as shown by stability studies for up to two weeks at room temperature.
Table 1:
List of oligos/DNA tested in this study with abbreviations and sequences; all bases deoxy including deoxyuridine (U). R1 and R2 represent the ratio R value upon using protocol A or B, respectively (see Experimental).
Earlier studies unexpectedly showed the kinetic independence of the labeling process on sequence, composition and length among a large number of tested deoxyoligos with length 8-mer to 99-mer, including the monomers. Note that only fully labeled product formation is kinetically independent of the above parameters, but not the disappearance of the intact oligo, as its kinetics are statistically defined by the number of pyrimidines in the molecule [30]. The absorbance of the osmylated product varies proportionally with the fraction of the osmylated pyrimidines over the total number of bases (Ntot), and the corresponding equation R = 2.01·[T(OsBp)+C(OsBp)]/Ntot was obtained using an oligo training set [31]. Osmylation kinetics of M13mp18, a 7249 nucleotide long circular ssDNA, concluded that this DNA osmylates just like any oligo and led to the conjecture that other ssDNAs may follow suit.
Figure 2 illustrates the reactivity of M13mp18, and directly compares it with that of a 20 nucleotide long deoxyoligo, the PCR Primer “16S RNA for” (see Table 1). This oligo was selected with T and C composition similar to M13mp18, so that R has a comparable value for both. The experiment was conducted by reacting each material with the same concentration of OsBp, initially 3.4 mM to practically osmylate Ts. After 2 h, OsBp was added so that the final concentration was 8.5 mM – a concentration that is suitable for labeling Cs relatively slowly so that the process can be monitored by HPLC or CE. Reactions were also conducted with 13.6 mM OsBp and the reaction mixtures were purified ahead of analysis. It is striking that under the above three different osmylation conditions, the kinetic behavior (see ascending curves, and plateaus) of the oligo and the long DNA are practically identical, suggesting that the same protocol may be used predictably for any ssDNA. The fact that all the data points lie somewhat lower for the oligo compared to M13mp18 is due to the somewhat lower fraction of T and C in the oligo. The ratio T/Ntot = 0.30 for the 20-mer and 0.33 for M13mp18; (T + C)/Ntot equals 0.50 for the 20-mer and 0.54 for M13mp18 [31].
Upper limits for purine osmylation and DNA backbone degradation (false positives)
Capillary zone electrophoresis (CZE) methods were extensively used to study the osmylation kinetics with nucleic acids [30,31]. The characteristic of CZE is that the reagent, OsBp, migrates early with a migration time (mt) that corresponds to neutral molecules, whereas the nucleic acid migrates late and the osmylated nucleic acid migrates in between, when resolution between substrate and product is feasible. More extensive osmylation leads to earlier mt values. The advantage of CZE is that typically nucleic acids share comparable mt irrespective of length and composition, which leads to the convenience of using a single method for all tested materials. The disadvantage is that CZE may not show tentative degradation of the backbone or plausible side reactions. Therefore, ion exchange (IE) HPLC with DNA Pac column from Dionex (see Experimental) was exploited. To our knowledge, this is one of the best columns to resolve oligos primarily based on length and also based on composition. Figure 3 illustrates the last point, namely the resolution of four 20 nucleotide long deoxyoligos that differ only in the 11th nucleotide: A10XA9 with X = A, T, G, and C, a separation not achieved by CZE.
Figure 4a illustrates the HPLC profile of A20 as the major peak preceded by a number of impurities of decreasing relative abundance. These impurities are magnified 40-fold in Figure 4b; they are the oligos with n < 20, where the lower the n, the earlier the rt. In Figure 4, the red trace represents the intact A20 and blue trace represents A20 after 12 h osmylation in 14.2 mM OsBp at room temperature. The two traces are superimposable, suggesting that A20 does not react and does not degrade under these conditions. Since the small peak after the main component measures 0.3%, one may conclude that a 0.2% peak would also be detectable and estimate that any reaction/degradation would be detectable at or above 0.2%. A comparable behavior is shown by A10GA9 in Figure 5a,b. Here again osmylation of this oligo does not produce any changes and/or any peaks that would have been detected at levels at or above 0.2%, as per the argument made above for A20. Since incubation is 6-times longer than protocol B requires and A20 has 20As, then the observation of 0.2% product would have corresponded to a 0.002/20/6 = 1.7/100,000 upper limit for A osmylation and 0.002/6 = 3.3/10,000 = 0.033% for G osmylation in A10GA9.
A second major conclusion from Figure 4 and Figure 5 is that oligos of n−1 length are clearly resolved with this method, which will be exploited later to show that an oligo composed of all four bases exhibits no detectable degradation to shorter oligos during osmylation by protocol B. The conclusion for undetectable osmylation of G is in contrast to observations published elsewhere [36], and it is attributed to the different conditions used in the two studies.
In contrast to the nonreactivity observed with A20 and A10GA9 (see Figure 4 and Figure 5), osmylated A10TA9 exhibits a completely different HPLC profile compared to the intact oligo (Figure 6). Figure 6 illustrates that the intact 20-mer (red trace) is replaced by two peaks (blue trace) with an approximate 2:1 ratio upon osmylation. Figure 7 (magnification of Figure 6) illustrates that the first major impurity (red) is also replaced by two peaks (blue) and this pattern continues up to about n = 11; namely, every impurity peak that includes a T reacts with OsBp and forms two isomeric products in an approximate 2:1 ratio. The formation of two isomers for T(OsBp) was observed by CE and reported earlier [30,31]. Specifically, two isomers are formed by addition of OsBp to the double bond from the top or the bottom of the pyrimidine ring. The HPLC analysis of the osmylated A10CA9 (not shown) is comparable to Figure 6 and Figure 7.
The above experiment with oligo A10GA9 (Figure 5) indicated that false positive and backbone degradation on a G base would have been detectable at levels at or above 0.033%. In order to assess a tentative G reactivity at higher sensitivity, we monitored the osmylation of an oligo (A5GA5GA5GA2) with 3 Gs by CE. The reaction was attempted in the presence of 12.6 mM OsBp and instead of 2 h, it was monitored for three days (Figure 8). Assuming a limit of 0.2% in product detectability, and accounting for 3 Gs and a 36-fold longer incubation, the upper limit for G reactivity with OsBp under protocol B conditions is 0.002/3/36 = 0.002% or 2/100,000. From this experiment, the limit for A-osmylation is at or less than 0.002/17/36 = 0.3/100,000, in agreement with the earlier experiment (Figure 4). These calculations of 0.3/100,000 for A reactivity and 2/100,000 for G reactivity are upper limits of detectability, and not actual levels of reactivity. This is due to limitations in instrumentation and the specifics of the experiment. These limits are also consistent with undetectable backbone degradation of the over 70 oligos that we have osmylated so far. It is likely that false positives for purine reaction and backbone degradation of purines under the harshest of our osmylation conditions (2 h with <15 mM OsBp in water at room temperature) are below 1/100,000.
To assess reactivity of OsBp towards the backbone of an oligo composed of all four nucleobases, including osmylated pyrimidines, we experimented with an oligo that happened to be of exceptional purity. This oligo is a 32-mer, 3’-end-adenylated PCR primer pGEX3’, abbreviated as pGEX3’A9 (Table 1). Figure 9 shows the IE HPLC profile of this 32-mer with the main peak eluting close to 18 min and several impurities eluting between the main peak and rt = 16 min. Since this method is the identical method we used for the analysis of the 20-mers shown in Figures 3–6, any short product oligo with n > 10 is expected to elute in the range of 5 min < rt < 16 min. To push the detectability limit, analysis was done on a large sample of intact oligo with an area of 20,468 HPLC units (red trace in Figure 9). A second aliquot from this oligo batch was “super”-osmylated (3 times longer than protocol B) and analyzed immediately after a sample of the intact. The osmylated oligo, as expected, elutes ahead of the intact as a broader peak, due to multiple topoisomers. Even though the tail of the osmylated peak extends towards the 15 min range, any peak produced as the result of backbone breakdown would have been detectable in the range of 5 min < rt < 15 min. Figure 9 shows a magnification of this range and assigns HPLC area to 2 units, attributable to a small visible peak that is present in both materials. This experiment suggests that degradation of the backbone of a 32 nucleotide long oligo that incorporates osmylated pyrimidines and 31 plausible positions for backbone cleavage would be visible at the level of 2/20,135/3/31 = 1/1,000,000, and suggests that OsBp reactivity towards a DNA backbone, osmylated or not, is lower than 1/1,000,000. As previously discussed, this limit is the result of our current instrumentation and other experiments may push this limit even lower.
Instability of C(OsBp)
There is a literature precedent that suggests that osmylated C deamidates and hydrolyzes to osmylated U [37]. In this work, ion channel recordings for pyrimidine(OsBp) were obtained within a sequence of oligodeoxyadenylates (A10XA9, where X is deoxypyrimidine), using α-HL nanopores. The response of the nanopore to the presence of the osmylated pyrimidine is dramatically different compared to an intact base. Still C(OsBp) and U(OsBp) exhibit comparable relative current obstruction (Ir/Io), but different dwell times (τ) at 120 mV in pH 7 buffered 1M KCl [34]. Specifically, Ir/Io = 0.12 for both, and corrected per unit mean translocation time τ = 312.5 and τ = 422.5 μs, for C(OsBp) and U(OsBp), respectively [34,38]. In contrast, Ir/Io = 0.08 and τ = 102.5 μs for T(OsBp) was obtained under the same conditions. In the context of nanopore-based osmylated ssDNA sequencing, the conversion of C(OsBp) to U(OsBp) may be inconsequential, because these two exhibit, relatively speaking, comparable properties, while both are clearly distinguishable from T(OsBp). Hence a mix-up between U(OsBp) and T(OsBp) may be as unlikely as a mix-up between C(OsBp) and T(OsBp). Still the conversion is undesirable in the context of nanopore-based sequencing of osmylated RNA, where riboC and riboU are the most abundant pyrimidines, especially if the ion-channel recordings of the osmylated riboC and riboU mimic the trends observed with their deoxy counterparts (see above).
To test the tentative transformations of C(OsBp) to U(OsBp), and of C(Me)(OsBp) to T(OsBp), as well as any other instability during the labeling process, we set up osmylation experiments in 13.6 mM OsBp using A10CA9, A10C(Me)A9, A10CA4 and evaluated the stability of the osmylated product over time using both IE HPLC and CE with the long 50 μm × 72 cm capillary (see the Experimental section). In addition, A10C(Me)(OsBp)A9, A10T(OsBp)A9 and A10U(OsBp)A9 were prepared and their stability was tested over time as well (see Supporting Information File 1 Figures S1, S2, and S3, respectively). These materials were also used for identification of tentative byproducts from the presumed hydrolysis of C(OsBp) to U(OsBp) and C(Me)(OsBp) to T(OsBp). We confirmed that C(Me)OsBp, T(OsBp), and U(OsBp) are stable over a couple of weeks at room temperature, whereas C(OsBp) both within the 20-mer and the 15-mer degrades over time. It was also established that the degradation of C(OsBp) is in the range of about 1–2% during the protocol B osmylation and that this degradation does not involve the loss of OsBp.
In the first set of experiments, the reaction of A10CA9 with OsBp was monitored by CE using a long capillary. There is a clear decrease in the substrate peak and a clear increase in the product peak. The two isomeric products are not resolved in this case, so only one product peak was observed. Osmylation after 2 h in 13.6 mM OsBp yielded a CE profile shown in Figure 10 where a small broad peak, measuring 1% of the total, appears at a later mt. Stability studies for 20 days at −5 °C and for an additional 5.5 h at 25 °C yielded a steady increase of the peak area. The stability studies under the above conditions with A10C(Me)(OsBp)A9, A10T(OsBp)A9, and A10U(OsBp)A9 showed no new peak and no other detectable changes (see Supporting Information File 1, Figures S1, S2, and S3, respectively).
In an attempt to establish whether or not the observed degradation can be attributed to deamidation and hydrolysis of C(OsBp) to U(OsBp), we explored the stability of A10C(OsBp)A4 by CE using the same long capillary. The starting oligo was completely consumed in 13.6 mM OsBp and formed two product peaks migrating ahead of the intact oligo. A comparable result was obtained by osmylation of A10UA4. The two product peaks are attributed to the two isomers that are now resolved by CE because of the shorter length, 15-mer vs 20-mer. Still CE partially resolves the two isomers of A10C(OsBp)A4 from the two isomers of A10U(OsBp)A4; upon mixing, three peaks were observed instead of the expected four (see Figure 11). Over a period of three days at room temperature degradation was observed that measured less than 3% (see Figure 12).
The above reaction (A10CA4 + OsBp) was reran using 5.5 mM OsBp (in order to slow down the kinetics) and monitored by IE HPLC. The disappearance of the intact oligo and the simultaneous appearance of two peaks eluting earlier was observed. These two peaks are attributed to the expected A10C(OsBp)A4 isomers (see Supporting Information File 1, Figure S4). In addition to the expected products, a small peak eluting after the starting material continued to increase with time. After 2 h, this peak measured about 2% of the total material. It should be noted that IE HPLC did not resolve A10C(OsBp)A4 from A10U(OsBp)A4. Hence, the observed, minor side product is not deamidation. The observed late peak grows over time and becomes the major component in room temperature stability samples. Similar observations were made with A10C(OsBp)A9, but not with any of the other oligos tested; hence only C(OsBp) degrades to a product not likely to be U(OsBp). A several-weeks-aged solution of A10C(OsBp)A9 is shown in the HPLC profiles of Figure 13 and illustrates the presence of a major decomposition product. Based on the literature, this degradation product may have been trans-uridine glycol [37], but that transformation requires that the OsBp moiety is labile. Figure 13 shows the profiles obtained at 272 nm and 312 nm and clearly indicates that each peak/product has the expected R = 0.11, strongly suggesting that OsBp is not lost upon standing. Comparable profiles are shown in Figure 13 for aged solutions of A10T(OsBp)A9 with no detectable changes. The stability of A10U(OsBp)A4 was also tested over a period of days and showed no detectable changes, suggesting that a transformation of C(OsBp) via U(OsBp) to the observed major degradant is not likely. It should be noted that all of the above stability tests were done in pure water, as the osmylation reaction is carried out in pure water. Further investigation into the stability of C(OsBp) containing oligos as well as the stability of C(Me)(OsBp) oligos should be done, and major degradants should be identified by NMR and mass spectrometry, but the corresponding reactions need to be carried out in buffers, because most of this instability is known to be pH-dependent [37]. For the purposes of sequencing osmylated DNA via nanopores, one may conjecture that tentative degradation of osmylated C under current labeling conditions of 2 h incubation time in 12–15 mM OsBp at room temperature yields 1–2% of alternative material that still carries the OsBp group. How this affects the quality of base calling still remains to be determined.
Even though there is no precedent in the literature that OsBp could react in a bifunctional way and yield inter- or intra-strand dimerization and/or polymerization of the reacting nucleic acid, this was investigated possibility as follows: pGEX3’-A9 was employed, which is a 32-mer long deoxyoligo, and both aged and fresh OsBp stock solutions were used for labeling. The reactions were monitored by IE HPLC and the method was adjusted to show the appearance of tentative 64 nucleotide long byproducts. No such long oligos were detected, even when the reaction was conducted with the highest concentration of OsBp for three times longer incubation compared to protocol B.
Estimate for unreacted pyrimidines (false negatives)
To ensure reproducibility, a reaction must be exploited under pseudo-first-order conditions. Under such conditions, i.e., 20–25-fold excess of OsBp over the reacting nucleobase in DNA, reactivity depends on OsBp concentration, but it is independent of the actual concentration of the reacting nucleobase (N). Hence, the fraction of unreacted N (pyrimidine in this case) can be calculated using the first-order rate equation (Equation 1) from the observed rate of reaction, k, for a specific OsBp concentration, and the incubation time, t.
where N is unreacted, N0 is the initial concentration, and t is the incubation time. The half-life of a reaction (t1/2) is by definition when N/N0 = 0.5, or ln(0.5) = 0.693, with t1/2 = 0.693/k, or k = 0.693/t1/2. Therefore, Equation 1 transforms to Equation 2 as:
Equation 2 provides an easy way to calculate multiples of half-lives necessary to reach a certain level of reactivity. For example, for N/No = 1/1000 incubation time t = 10 half-lives, for N/No = 10,000, t = 13.3 half-lives and for N/N0 = 1/100,000, t = 16.6 half-lives, etc.
Equation 2 applies to any first-order or pseudo-first-order reaction. Hence the incubation time t necessary to reach a fraction of unreacted T or C equal to 1/10,000 is 13.3 half-lives of rate for T or C, respectively. We have repeatedly found that T osmylates with a rate of about 0.045 per min in the presence of 3.2 mM OsBp (protocol A) and C osmylates with a rate of about 0.08 per min with 12.6 mM OsBp (protocol B). These rates yield half-lives of about 15 min for T and 9 min for C at the corresponding OsBp concentrations. Hence 13.3 half-lives for C osmylation corresponds to 120 min and this is the incubation time recommended for protocol B. Evidently, after employing protocol B, intact T is orders of magnitude less than 1/10,000, because of its higher reactivity compared to C. On the contrary, protocol A aims to yield the minimum amount of osmylated C, while osmylating over 90% T. This is why protocol A exploits a 60 min incubation time with 3.2 mM OsBp, which yields from Equation 1: ln(N/N0) = −2.7 and N/N0 = 0.067, i.e., 6.7% intact T, or 93.3% osmylated T. It is advisable that new stock solutions of OsBp, as well as aged OsBp solutions, are evaluated for consistent reactivity.
Conclusion
The experiments reported in this study strongly support the proposal that OsBp is, from a chemistry point-of-view, an outstanding agent to label pyrimidines in deoxyoligonucleotides. Whether or not OsBp is equally fit for translocation experiments that yield sequencing information remains to be shown. Our experimental set-up allows the determination of upper levels of purine osmylation and DNA backbone degradation under the proposed conditions. Approximate upper limits of 1/100,000, 1/100,000, and 1/1,000,000 were obtained for G-osmylation, A-osmylation, and DNA backbone cleavage, respectively. These upper limits were obtained using oligodeoxynucleotides and are likely to be valid with oligoribonucleotides and long nucleic acids because the osmylation reaction appears to be confined to the actual pyrimidine ring and independent of strand length and composition. Under the current protocol B, for practically 100% osmylation of T + C, levels of unreacted C are estimated at 1/10,000 and levels of unreacted T are negligible. Whereas osmylated deoxy T, (Me)C, and U are stable at room temperature over weeks, osmylated C exhibits degradation of 1–2% during the labeling protocol B. Osmylated C continues to degrade in water over time to products all carrying the OsBp moiety; the major degradant being inconsistent with C to U deamidation/hydrolysis. Due to C(OsBp) degradation, osmylated DNA needs to be stored at −20 °C. In summary, the results presented here provide improved protocols for the preparation of osmylated DNA, an overview of the upper limits in false positives and false negatives, and define the limits for accurate representation of the parent DNA in sequencing efforts.
Experimental
Materials
HPCE grade solution of 50 mM sodium tetraborate (pH 9.3) was purchased from Agilent Technologies. A 4% aqueous osmium tetroxide solution in 2 mL ampules was purchased from Electron Microscopy Sciences. 2,2’-bipyridine, 99+% (bipyridine) was purchased from Acros Organics. Oligos were ordered from Integrated DNA Technologies (IDT), some desalted and others PAGE purified. Circular ssDNA M13mp18 was purchased from Bayou Labs [35]. Oligos and DNA were diluted with DNase/RNase-free water (from MP Biomedical) to desired μg/μL concentration or to 100 μM stock solutions. Oligo solutions were stored at −20 °C. The purity of purchased oligos was tested by CE analysis and was typically better than 85%. The oligos used in this study, abbreviations, values of R of the corresponding osmylated derivative, and sequences are listed in Table 1.
Preparation of OsBp stock solutions and protocols for osmylation reactions
Typically, 1:1 mixtures of OsO4 and 2,2’-bipyridine were prepared. A stock solution of OsBp at 15.75 mM OsO4 was prepared by mixing a 2 mL 4% OsO4 (ampule) with 18 mL of water that contained the equivalent 15.75 mM of 2,2’-bipyridine. The solubility tests in water were performed visually using weighted amounts of 2,2’-bipyridine, and its solubility was estimated at about 21 mM at room temperature. The solubility of 2,2’-bipyridine was also tested in solutions with increasing concentration of OsO4, and did not measurably increase over 21 mM, confirming that OsBp does not form a very stable complex. Typically, stock solutions of OsBp were prepared as 1:1 mixtures, and the use of saturated 2,2’-bipyridine solutions was abandoned [34] because of the difficulty in assuring that all the suspended 2,2’-bipyridine was removed. Due to OsBp reactivity, paper filters should not be used to filter OsBp, and plastic vials should also not be used to centrifuge OsBp. The stock solutions were dispensed in small, clear glass vials, sealed with parafilm, and kept at −20 °C until use; no detectable change in reactivity was observed with solutions that were stored sealed at 4 °C for four weeks. The experiments were initiated by mixing the OsBp stock solution and the oligo stock solution directly in a CE glass vial or HPLC vial or another stoppered glass vial filled with the appropriate amount of distilled water. The final concentration of [OsO4] = [2,2’-bipyridine], reported as [OsBp] is listed with each experiment. No buffer was added and reaction mixtures were incubated at 25 ± 2 °C. Analyses were carried out automatically. The aliquoting time together with the product distribution, in the form of an electropherogram or chromatogram, were recorded by the CE or HPLC. Spin columns (TC-100 FC, TrimGen) were used to remove excess OsBp according to the manufacturer’s instructions. Practically 100% recovery of labeled oligo and removal of up to 15 mM OsBp down to detectability levels after one or two purifications was achieved. Protocol A refers to 60 min incubation with 3.2–3.4 mM OsBp at room temperature and protocol B refers to 120 min incubation with 12.6–15 mM OsBp at room temperature. An earlier version of protocol B exploited 12 h labeling at room temperature, but recent results have shown it to be unnecessarily long. Protocol A results in over 90% T osmylation and less than 8% C-osmylation and protocol B results in 100% T + C osmylation.
Capillary electrophoresis methods and analyses
Analyses of the reaction mixtures were carried out using an Agilent 1600 CE instrument equipped with a diode array detector (DAD) and Chemstation software, Rev.B.04.02 SP1, for data acquisition and processing. Only glass-type CE vials were used, as OsBp was found to react with polyurethane vials, lowering its effective concentration. Untreated fused silica capillaries (short: 50 μm × 40 cm and long: 50 μm × 72 cm) with an extended light path were purchased from Agilent Technologies.
The analysis of reaction mixtures was conducted using capillary zone electrophoresis (CZE) methods in 50 mM sodium tetraborate (pH 9.3) with 20 kV for the short capillary and 30 kV for the long capillary. The capillary temperature was typically 25 °C. The method associated with the short capillary was 18 min long and included a 3 min automatic buffer conditioning before sample analysis. The method associated with the long capillary was about 30 min long and included a 5 min wash. OsBp migrates early together with any other neutral compound present, the triphosphates of the nucleotides migrate last, and oligos/DNA in between. OsBp-labeled products migrate between the OsBp peak and the corresponding DNA peak. The resolution between the starting material and osmylated product, is among other parameters, a function of the oligo length and the number of reacting bases. With the longer DNA, a shift of the peak to earlier migration is detected as a function of the reaction progress. In CE profiles, the more abundant topoisomer migrates first ([30] and Figure 11), whereas in HPLC profiles, the most abundant topoisomer elutes last (Figure 6).
High performance liquid chromatography methods and analyses
As a second analytical tool, we used an Agilent 1100/1200 LC HPLC equipped with a binary pump, DAD and a 1290 Infinity Autosampler/Thermostat, and Chemstation software Rev.B.04.01 SP1 for data acquisition and processing. Both the autosampler and the column compartment have individual temperature control. The method most suitable for oligo analysis is ion exchange (IE) with a salt gradient at 30 °C in the column compartment. A Dionex HPLC column, DNAPac 200PA, with a 2.0 μm inner diameter and 25 cm path length was used for the resolution of oligos based on both length and composition. Typically, a flow of 0.5 mL/min was used and a linear gradient system with solvent A of 25 mM Tris HCl at pH 7.0 and solvent B of 1.0 M NaCl in solvent A; the initial elution conditions were typically at 75% for A and 25% for B. The actual gradient slope, at about 1–1.5% B per minute, varied depending on the length of the oligos that needed to be resolved. Both CE and HPLC peaks were detected and identified using the DAD in the UV–vis region 200–450 nm and the electropherograms and chromatograms were recorded at several wavelengths, including 260 nm, 272 nm and 312 nm.
Supporting Information
Supporting Information File 1:
Additional figures are included to illustrate the stability of oligos with C(Me)(OsBp), T(OsBp) and U(OsBp) and the osmylation reaction of a 15-mer oligo monitored by IE HPLC.
Anastassia Kanavarioti is the founder and director of Yenos Analytical LLC, a company engaged in the development and manufacturing of osmylated nucleic acids.
Sanger, F.; Nicklen, S.; Coulson, A. R. Proc. Natl. Acad. Sci. U. S. A.1977,74, 5463–5467. doi:10.1073/pnas.74.12.5463
Return to citation in text:
[1]
Fuller, C. W.; Middendorf, L. R.; Benner, S. A.; Church, G. M.; Harris, T.; Huang, X.; Jovanovich, S. B.; Nelson, J. R.; Schloss, J. A.; Schwartz, D. C.; Vezenov, D. V. Nat. Biotechnol.2009,27, 1013–1023. doi:10.1038/nbt.1585
Return to citation in text:
[1]
Bayley, H.; Martin, C. R. Chem. Rev.2000,100, 2575–2594. doi:10.1021/cr980099g
Return to citation in text:
[1]
Venkatesan, B. M.; Bashir, R. Nat. Nanotechnol.2011,6, 615–624. doi:10.1038/nnano.2011.129
Return to citation in text:
[1]
Reiner, J. E.; Balijepalli, A.; Robertson, J. W. F.; Campbell, J.; Suehle, J.; Kasianowicz, J. J. Chem. Rev.2012,112, 6431–6451. doi:10.1021/cr300381m
Return to citation in text:
[1]
Haque, F.; Li, J.; Wu, H.-C.; Liang, X.-J.; Guo, P. Nano Today2013,8, 56–74. doi:10.1016/j.nantod.2012.12.008
Return to citation in text:
[1]
Maglia, G.; Heron, A. J.; Stoddart, D.; Japrung, D.; Bayley, H. Methods Enzymol.2010,475, 591–623. doi:10.1016/S0076-6879(10)75022-9
Return to citation in text:
[1]
[2]
[3]
Butler, T. Z.; Pavlenok, M.; Derrington, I. M.; Niederweis, M.; Gundlach, J. H. Proc. Natl. Acad. Sci. U. S. A.2008,105, 20647–20652. doi:10.1073/pnas.0807514106
Return to citation in text:
[1]
[2]
Jing, P.; Haque, F.; Shu, D.; Montemagno, C.; Guo, P. Nano Lett.2010,10, 3620–3627. doi:10.1021/nl101939e
Return to citation in text:
[1]
Church, G.; Deamer, D. W.; Branton, D.; Baldarelli, R.; Kasianowicz, J. Characterization of individual polymer molecules based on monomer-interface interactions. U.S. Patent 5,795,782, Aug 18, 1998.
Return to citation in text:
[1]
Kasianowicz, J. J.; Brandin, E.; Branton, D.; Deamer, D. W. Proc. Natl. Acad. Sci. U. S. A.1996,93, 13770–13773. doi:10.1073/pnas.93.24.13770
Return to citation in text:
[1]
Akeson, M.; Branton, D.; Kasianowicz, J. J.; Brandin, E.; Deamer, D. W. Biophys. J.1999,77, 3227–3233. doi:10.1016/S0006-3495(99)77153-5
Return to citation in text:
[1]
Meller, A.; Nivon, L.; Brandin, E.; Golovchenko, J.; Branton, D. Proc. Natl. Acad. Sci. U. S. A.2000,97, 1079–1084. doi:10.1073/pnas.97.3.1079
Return to citation in text:
[1]
Laszlo, A. H.; Derrington, I. M.; Ross, B. C.; Brinkerhoff, H.; Adey, A.; Nova, I. C.; Craig, J. M.; Langford, K. W.; Samson, J. M.; Daza, R.; Doering, K.; Shendure, J.; Gundlach, J. H. Nat. Biotechnol.2014,32, 829–833. doi:10.1038/nbt.2950
Return to citation in text:
[1]
Ashkenasy, N.; Sánchez-Quesada, J.; Bayley, H.; Ghadiri, M. R. Angew. Chem., Int. Ed.2005,44, 1401–1404. doi:10.1002/anie.200462114
Return to citation in text:
[1]
Stoddart, D.; Heron, A. J.; Mikhailova, E.; Maglia, G.; Bayley, H. Proc. Natl. Acad. Sci. U. S. A.2009,106, 7702–7707. doi:10.1073/pnas.0901054106
Return to citation in text:
[1]
Wolna, A. H.; Fleming, A. M.; An, N.; He, L.; White, H. S.; Burrows, C. J. Isr. J. Chem.2013,53, 417–430. doi:10.1002/ijch.201300022
Return to citation in text:
[1]
Ervin, E. N.; Barrall, G. A.; Pal, P.; Bean, M. K.; Schibel, A. E. P.; Hibbs, A. D. Bionanosci.2014,4, 78–84. doi:10.1007/s12668-013-0119-0
Return to citation in text:
[1]
Stoddart, D.; Ayub, M.; Höfler, L.; Raychaudhuri, P.; Klingelhoefer, J. W.; Maglia, G.; Heron, A.; Bayley, H. Proc. Natl. Acad. Sci. U. S. A.2014,111, 2425–2430. doi:10.1073/pnas.1312976111
Return to citation in text:
[1]
Wendell, D.; Jing, P.; Geng, J.; Subramaniam, V.; Lee, T. J.; Montemagno, C.; Guo, P. Nat. Nanotechnol.2009,4, 765–772. doi:10.1038/nnano.2009.259
Return to citation in text:
[1]
Lieberman, K. R.; Cherf, G. M.; Doody, M. J.; Olasagasti, F.; Kolodji, Y.; Akeson, M. J. Am. Chem. Soc.2010,132, 17961–17972. doi:10.1021/ja1087612
Return to citation in text:
[1]
Branton, D.; Deamer, D. W.; Marziali, A.; Bayley, H.; Benner, S. A.; Butler, T.; Di Ventra, M.; Garaj, S.; Hibbs, A.; Huang, X.; Jovanovich, S. B.; Krstic, P. S.; Lindsay, S.; Ling, X. S.; Mastrangelo, C. H.; Meller, A.; Oliver, J. S.; Pershin, Y. V.; Ramsey, J. M.; Riehn, R.; Soni, G. V.; Tabard-Cossa, V.; Wanunu, M.; Wiggin, M.; Schloss, J. A. Nat. Biotechnol.2008,26, 1146–1153. doi:10.1038/nbt.1495
Return to citation in text:
[1]
[2]
Beer, M.; Moudrianakis, E. N. Proc. Natl. Acad. Sci. U. S. A.1962,48, 409–416. doi:10.1073/pnas.48.3.409
Return to citation in text:
[1]
Chang, C.-H.; Beer, M.; Marzilli, L. G. Biochemistry1977,16, 33–38. doi:10.1021/bi00620a006
Return to citation in text:
[1]
Kanavarioti, A.; Greenman, K. L.; Hamalainen, M.; Jain, A.; Johns, A. M.; Melville, C. R.; Kemmish, K.; Andregg, W. Electrophoresis2012,33, 3529–3543. doi:10.1002/elps.201200214
Return to citation in text:
[1]
[2]
[3]
[4]
[5]
[6]
[7]
[8]
[9]
[10]
Kanavarioti, A. Labeled Nucleic Acids: A Surrogate for Nanopore-based Nucleic Acid Sequencing. U.S. Patent Application 14/944,888, Nov 18, 2015.
Return to citation in text:
[1]
[2]
[3]
Henley, R. Y.; Vazquez-Pagan, A. G.; Johnson, M.; Kanavarioti, A.; Wanunu, M. PLoS One2015,10, e0142155. doi:10.1371/journal.pone.0142155
Return to citation in text:
[1]
Ding, Y.; Kanavarioti, A. Beilstein J. Nanotechnol.2016,7, 91–101. doi:10.3762/bjnano.7.11
Return to citation in text:
[1]
[2]
[3]
[4]
Branton, D.; Deamer, D. W.; Marziali, A.; Bayley, H.; Benner, S. A.; Butler, T.; Di Ventra, M.; Garaj, S.; Hibbs, A.; Huang, X.; Jovanovich, S. B.; Krstic, P. S.; Lindsay, S.; Ling, X. S.; Mastrangelo, C. H.; Meller, A.; Oliver, J. S.; Pershin, Y. V.; Ramsey, J. M.; Riehn, R.; Soni, G. V.; Tabard-Cossa, V.; Wanunu, M.; Wiggin, M.; Schloss, J. A. Nat. Biotechnol.2008,26, 1146–1153. doi:10.1038/nbt.1495
Fuller, C. W.; Middendorf, L. R.; Benner, S. A.; Church, G. M.; Harris, T.; Huang, X.; Jovanovich, S. B.; Nelson, J. R.; Schloss, J. A.; Schwartz, D. C.; Vezenov, D. V. Nat. Biotechnol.2009,27, 1013–1023. doi:10.1038/nbt.1585
Butler, T. Z.; Pavlenok, M.; Derrington, I. M.; Niederweis, M.; Gundlach, J. H. Proc. Natl. Acad. Sci. U. S. A.2008,105, 20647–20652. doi:10.1073/pnas.0807514106
22.
Ervin, E. N.; Barrall, G. A.; Pal, P.; Bean, M. K.; Schibel, A. E. P.; Hibbs, A. D. Bionanosci.2014,4, 78–84. doi:10.1007/s12668-013-0119-0
23.
Stoddart, D.; Ayub, M.; Höfler, L.; Raychaudhuri, P.; Klingelhoefer, J. W.; Maglia, G.; Heron, A.; Bayley, H. Proc. Natl. Acad. Sci. U. S. A.2014,111, 2425–2430. doi:10.1073/pnas.1312976111
Kanavarioti, A.; Greenman, K. L.; Hamalainen, M.; Jain, A.; Johns, A. M.; Melville, C. R.; Kemmish, K.; Andregg, W. Electrophoresis2012,33, 3529–3543. doi:10.1002/elps.201200214
Church, G.; Deamer, D. W.; Branton, D.; Baldarelli, R.; Kasianowicz, J. Characterization of individual polymer molecules based on monomer-interface interactions. U.S. Patent 5,795,782, Aug 18, 1998.
Maglia, G.; Heron, A. J.; Stoddart, D.; Japrung, D.; Bayley, H. Methods Enzymol.2010,475, 591–623. doi:10.1016/S0076-6879(10)75022-9
12.
Butler, T. Z.; Pavlenok, M.; Derrington, I. M.; Niederweis, M.; Gundlach, J. H. Proc. Natl. Acad. Sci. U. S. A.2008,105, 20647–20652. doi:10.1073/pnas.0807514106
Kanavarioti, A.; Greenman, K. L.; Hamalainen, M.; Jain, A.; Johns, A. M.; Melville, C. R.; Kemmish, K.; Andregg, W. Electrophoresis2012,33, 3529–3543. doi:10.1002/elps.201200214
Maglia, G.; Heron, A. J.; Stoddart, D.; Japrung, D.; Bayley, H. Methods Enzymol.2010,475, 591–623. doi:10.1016/S0076-6879(10)75022-9
18.
Laszlo, A. H.; Derrington, I. M.; Ross, B. C.; Brinkerhoff, H.; Adey, A.; Nova, I. C.; Craig, J. M.; Langford, K. W.; Samson, J. M.; Daza, R.; Doering, K.; Shendure, J.; Gundlach, J. H. Nat. Biotechnol.2014,32, 829–833. doi:10.1038/nbt.2950
Kanavarioti, A.; Greenman, K. L.; Hamalainen, M.; Jain, A.; Johns, A. M.; Melville, C. R.; Kemmish, K.; Andregg, W. Electrophoresis2012,33, 3529–3543. doi:10.1002/elps.201200214
Kanavarioti, A.; Greenman, K. L.; Hamalainen, M.; Jain, A.; Johns, A. M.; Melville, C. R.; Kemmish, K.; Andregg, W. Electrophoresis2012,33, 3529–3543. doi:10.1002/elps.201200214
Kanavarioti, A.; Greenman, K. L.; Hamalainen, M.; Jain, A.; Johns, A. M.; Melville, C. R.; Kemmish, K.; Andregg, W. Electrophoresis2012,33, 3529–3543. doi:10.1002/elps.201200214
Kanavarioti, A.; Greenman, K. L.; Hamalainen, M.; Jain, A.; Johns, A. M.; Melville, C. R.; Kemmish, K.; Andregg, W. Electrophoresis2012,33, 3529–3543. doi:10.1002/elps.201200214
Kanavarioti, A.; Greenman, K. L.; Hamalainen, M.; Jain, A.; Johns, A. M.; Melville, C. R.; Kemmish, K.; Andregg, W. Electrophoresis2012,33, 3529–3543. doi:10.1002/elps.201200214
Kanavarioti, A.; Greenman, K. L.; Hamalainen, M.; Jain, A.; Johns, A. M.; Melville, C. R.; Kemmish, K.; Andregg, W. Electrophoresis2012,33, 3529–3543. doi:10.1002/elps.201200214
Kanavarioti, A.; Greenman, K. L.; Hamalainen, M.; Jain, A.; Johns, A. M.; Melville, C. R.; Kemmish, K.; Andregg, W. Electrophoresis2012,33, 3529–3543. doi:10.1002/elps.201200214
Branton, D.; Deamer, D. W.; Marziali, A.; Bayley, H.; Benner, S. A.; Butler, T.; Di Ventra, M.; Garaj, S.; Hibbs, A.; Huang, X.; Jovanovich, S. B.; Krstic, P. S.; Lindsay, S.; Ling, X. S.; Mastrangelo, C. H.; Meller, A.; Oliver, J. S.; Pershin, Y. V.; Ramsey, J. M.; Riehn, R.; Soni, G. V.; Tabard-Cossa, V.; Wanunu, M.; Wiggin, M.; Schloss, J. A. Nat. Biotechnol.2008,26, 1146–1153. doi:10.1038/nbt.1495
Kanavarioti, A.; Greenman, K. L.; Hamalainen, M.; Jain, A.; Johns, A. M.; Melville, C. R.; Kemmish, K.; Andregg, W. Electrophoresis2012,33, 3529–3543. doi:10.1002/elps.201200214