Kinetics of phosphorothioate oligonucleotide metabolism in biological fluids
Kinetics of phosphorothioate oligonucleotide metabolism in biological fluidsM. Gilar, A. Belenky, D. L. Smisek, A. Bourque and A. S. Cohen*
Hybridon, Inc., 620 Memorial Drive, Cambridge, MA 02139, USA
Received June 3, 1997;Revised and Accepted August 4, 1997
ABSTRACT
The in vitro stability and metabolism of GEM[91, a 25mer phosphorothioate antisense oligonucleotide complementary to the gag mRNA region of HIV-1, was investigated using capillary electrophoresis (CE). The in vitro degradation of the parent compound at 37oC was followed over the course of 120 h in human plasma. A CE method using laser-induced fluorescence detection was able to detect 5'-end intact metabolites including the parent compound extracted from biological fluids. Because the primary metabolic pathway is believed to be via 3'-exonuclease activity, the results of this study were compared with the stability of the compound in a solution containing 3'-exonuclease. The numerical solution of sequential first-order reactions was used to obtain kinetic parameters. Exonuclease digestion of the parent compound, as measured using an automated CE-UV instrument, yielded striking similarities between the two in vitro systems as well as between in vitro and in vivo systems.
INTRODUCTION
A new era of genetic medicine promises to combat disease at the genetic level. One antiviral approach uses chemically modified oligodeoxyribonucleotides as antisense agents that hybridize specifically to a target mRNA sequence resulting in the suppression of mRNA translation to the corresponding protein and thus virus reproduction (1 ). The mechanism of antisense action was first described by Zamecnik and Stephenson (2 ). Zamecnik et al. (3 ) were also first to report an antisense approach to inhibit de novo infection by HIV-1 using an unmodified oligodeoxyribonucleotide (ODN). For successful therapeutic use, the in vivo stability as well as the bioavailability of antisense compounds is crucial. Sulfurization of ODNs is a frequently used approach to increase in vivo stability. The resulting phosphorothioate analogs (SODNs) are more resistant to exonucleases than the unmodified phosphodiester ODNs (4 ). Even though this modification is one of the most conservative, the chemical properties of a phosphorothioate molecule are different from its phosphodiester counterpart. For example, because of increased charge delocalization, phosphorothioates have a lower pKa than phosphodiesters (5 ); phosphorothioates are also more hydrophobic and therefore exhibit more complex secondary structure. Also, the distribution of phosphorothioates in the living organism can be different due to higher binding to proteins and accumulation in organs (6 ).
As more antisense compounds enter human clinical trials, it becomes increasingly important to obtain a fundamental understanding of the kinetics of degradation and metabolism both in vitro and in vivo. To establish a baseline for further study, we have investigated the in vitro stability and metabolism of GEM®91, a 25mer phosphorothioate antisense oligonucleotide complementary to the gag mRNA region of HIV-1, using capillary electrophoresis (CE).
Several papers have been published on in vivo and in vitro stability and pharmacokinetics of phosphorothioate antisense oligonucleotides. Phosphorothioates were found to be stable with little or no metabolism (7 ). Most of these studies used radioisotope labeled 32P, 35S or 14C phosphorothioate oligonucleotides. Slab gel electrophoresis on polyacrylamide gels was used for separation of SODN metabolites with radioisotope detection to estimate the metabolites; however, quantification was very poor. Measurements of SODN stability by radioisotope methods sum not only the parent drug radioactivity but also the radioactivity of shorter metabolites and even free isotopes. Precise measurement of parent SODN stability depends on the reliability and sensitivity of the analytical method used for drug level monitoring. To avoid the disadvantages associated with radioisotope methods, we developed a non-radioisotope method for monitoring SODNs in biological fluids using CE with laser-induced fluorescence detection (CE-LIF) with fluorescently-labeled target compounds (8 ).
Because the primary metabolic pathway in plasma is believed to be by 3'-end exonuclease activity, experiments were conducted both in human plasma and in a solution containing 3'-exonuclease. Our CE-LIF method detects 5'-end intact metabolites including the parent compound from biological fluids. We also developed a CE method with UV detection (CE-UV) that detects both 3'- and 5'-end digested metabolites. For comparison, an unmodified ODN with the same sequence as GEM®91 was run under similar conditions. The kinetics of enzymatic degradation was modeled using a numerical solution of sequential first-order reactions and calculated results compared to measured in vitro digestion data.
MATERIALS AND METHODS
Chemical and reagents
Water and methanol were HPLC grade (J. T. Baker, Inc., Phillipsburg, NJ, USA). Human plasma/serum (male) was obtained from Sigma Chemical (St Louis, MO, USA). Uniflo® 13 mm diameter 0.2 [mu]m cellulose triacetate disposable syringe filters were obtained from Schleicher and Schuell (Keene, NH, USA). Solid-phase extraction (SPE) was performed on silica based weak anion exchange media in a 2 ml SPE tube (Nucleobond AX-5 cartridges, Macherey Nagel Corp., GmbH & Co., Düren, Germany) and silica based reverse phase C18 media in a 1 ml SPE tube (LC-18, Supelco Corp., Bellefonte, PA, USA). Oligomer HPLC purification was performed on Nucleopak PA-100, 250 * 9 mm column (Dionex, Sunnyvale, CA, USA) and polystyrene-divinylbenzene PRP-3, 150 * 4.1 mm column (Hamilton, Reno, NV, USA). Lithium bromide and other chemicals were obtained from Fluka Chemical Corp. (Ronkonkoma, NY, USA). Deionized formamide was purchased from American Bioanalytical (Natick, MA, USA). Phosphodiesterase I, a 3' to 5'-exonuclease enzyme isolated from bovine intestinal mucosa, was obtained from Sigma (St Louis, MO, USA).
The 25mer ODN and SODN were synthesized in our lab with an Expeditetm 8909 NASS synthesizer (PerSeptive Biosystems, Framingham, MA, USA), deprotected, purified and reconstituted in deionized water. Fluorescently labeled primer was obtained from Applied Biosystems (Foster City, CA, USA); T4 DNA ligase and ligase buffer (C#202S) and T4 Kinase (C#201S) were obtained from New England Biolabs (Beverly, MA, USA).
Sample preparation and reaction conditions
Human plasma. Plasma samples spiked with ODN were incubated at 37oC. Aliquots were diluted 1:1 with 0.2 M Tris-PO4, pH 6.3 (EQ), spiked with internal standard, and loaded onto an equilibrated anion exchange (AE) SPE column. The column was washed with 1 * 0.5 ml EQ, 3 * 0.5 ml of 0.2 M NaBr, 0.1 M Tris-PO4, pH 6.3 in 50% formamide (W1), 1 * 0.5 ml of 0.2 M NaBr, 0.1 M Tris-Cl, pH 7.0 (W2) and 1 * 0.2 ml of 2 M NaBr, 0.1M Tris-Cl (pH 8.5) in 10% isopropanol (E1). A 1.5 ml Eppendorf vial was placed under the SPE tube and 1 * 0.5 m of E1 was added. Because the eluate contained a high ionic strength solution which can suppress electrokinetic CE injection, AE-SPE extracts were desalted using reversed phase SPE as follows. The samples were diluted 1:1 with a 2 M triethylammonium acetate (TEAA) solution (pH 7.0) and loaded onto a conditioned C18 cartridge. The cartridge was washed with 2 * 0.5 ml of 0.1 M TEAA, 1 * 0.5 ml 5% ACN, and the DNA was eluted with 1 * 0.5 ml 40% ACN. The extract was lyophilized, reconstituted in 20 [mu]l of 10 mM NaOH and dialyzed by drop dialysis for 60 min over water using Millipore 0.025 [mu]m membranes. Dialyzed samples were subjected to CE.
Plasma samples containing SODNs were left to react at 37oC. The reaction was stopped by freezing, and samples were stored at -20oC. Samples were diluted 10-fold with water and spiked with 30mer internal standard. Samples were then purified from plasma by SPE. A fluorescently-labeled 17mer primer was ligated to the intact 5'-end of the fragments present in the sample including the parent compound. The labeling has two consecutive steps: phosphorylation of the 5'-end of the oligonucleotide sample followed by ligation of the phosphorylated sample with the fluorescently-labeled primer. A more detailed protocol was presented earlier (8 ). Ligation occurs only when the bridge DNA hybridizes to a target oligonucleotide whose 5'-end is intact (Fig. 1 ). Any incompatibility at the 5'-end of the target oligonucleotide with the complementary bridge will prevent ligation to the fluorescent primer. As a result, fluorescently-labeled SODN fragments which differ by one base (due to 3'-end digestion only) can be detected and analyzed at the parts-per-billion level by CE-LIF.
Capillary electrophoresis
CE-LIF. Fluorescently-labeled samples were run on a modular CE system with LIF detection (8 ). An Argon-ion laser (Model 543 100 BS, Omnichrom, Chino, CA, USA) was employed. A 30 kV, 500 [mu]A direct-current high-voltage power supply (Model ER/DM; Glassman, Whitehouse Station, NJ, USA) was used to generate the potential across the capillary. Samples were injected electrokinetically, 10 kV for 2-5 s and then run at 10 kV. Linear polyacrylamide-filled capillaries were prepared in-house by the following manner. Fused-silica capillary tubing (Polymicro Technologies, Phoenix, AZ, USA) with inner diameter of 75 [mu]m, outer diameter of 375 [mu]m, effective length of 10-15 cm, and total length of 27 cm was bifunctionalized with (methylacryloxypropyl) trimethoxysilane (Petrarch Systems, Bristol, PA, USA) and then filled with a degassed solution of 15% polymerizing linear acrylamide in 15% (v/v) formamide media (0.2 M Tris-borate, 2.5 mM EDTA-2Na+ buffer, pH 8.3, containing 7 M urea). Polymerization was achieved using ammonium persulfate/TEMED chemistry. An electric field of 400-500 V/cm was applied resulting in a current of 10-15 [mu]A. CE-UV. A BioFocusr 2000 Capillary Electrophoresis System (BioRad, Hercules, CA, USA) with ultra-violet (UV) detection was used. After each run the polymer solution was replaced. The composition of replaceable polymer solution was as follows. To 2.2 g of BioRad Dynamic Sieving Polymer was added deionized formamide (4.8 ml), 1.2 ml of 10* Tris-boric acid buffer pH 9 (1 M Tris, 0.3 M boric acid), and 1.2 ml of 8.3 M urea solution in 2* Tris-boric acid buffer (0.2 M Tris, 0.06 M boric acid). The Dynamic Sieving Polymer was dissolved, and the solution was degassed under vacuum. The polymer matrix was stable for a week at 4oC. For analysis a 25 cm (20 cm to detection window) coated capillary, 75 [mu]m I.D. * 375 [mu]m O.D. (BioCap oligonucleotide from BioRad) was used. Capillary filling was automated by applying a high pressure N2 pulse for 220 s. After each run the capillary was washed with HPLC-grade water. Injection of samples was 8 s electrokinetically at 13 kV, and samples were run at 15 kV. Separation was performed at 35oC with a liquid thermostated column.
Numerical solution of sequential first-order reactions
The metabolism of oligomers from exonuclease activity can be modeled as a series of first-order reactions:
k25'-1mer
k24'-1mer
k23'-1mer
25mer
-->
24mer
->
23mer
->
22mer
...
1
where ki is the rate constant. If we start with a 25mer and assume first-order kinetics, the rate of parent compound disappearance during the enzymatic reaction is
r = - dc25 / dt = c25k25
2
where c25 is the concentration of 25mer, k25 is the rate constant and t is time. The concentration profile of 25mer as a function of time was solved numerically as:
ci25 = ci-125 + (dci25) - dci25 = k25c25dt
3
The actual concentration ci25 in the first step is equal to c025, and due to a gradually decreasing ci25 value, the concentration step dc25 decreases exponentially with time. For simplicity the integration step dt was chosen as 1 min and the index number i = 1, 2, 3.... Similar calculations were performed to compute the concentration of shorter products noting that the rate of 25mer disappearance dc25 is equal to the rate of 24mer appearance dc24 etc. The calculations were performed using a Microsoft Excel spreadsheet assuming that the kinetic constant k25 is equal to the kinetic constants of subsequent reactions of shorter metabolites (k25=k24 = k23...).
Data quantitation
Quantitation was performed using an internal standard and using a method of internal normalization. Internal normalization was based on the assumption that the sum of all SODN fragments peak height is proportional to the spiked molar concentration:
csp = [sum] ci [approx] Hsp = [sum] Hi
4
where Hi is the peak height of SODN fragments and csp is the spiked concentration of SODN 25mer. Summation of peak heights was possible using LIF detection because of a homoextinction response of SODN fragments from the detector where each fragment at the same molar concentration gives the same signal. The concentrations of the SODN parent compound and metabolism products were calculated as:
ci = Hi / [sum] Hi × csp
5
We assume that the efficiency of fluorescent labeling (ligation with fluorescent primer) is equal for eachSODN fragment and the digestion of SODN from the 5'-end can be neglected. In our experience the shape and efficiency of closely migrating peaks do not significantly differ, and the internal normalization method gave better results than the results obtained by internal standard only.
RESULTS AND DISCUSSION
CONCLUSIONS
The half-lives of SODN and ODN 25mer oligonucleotides (GEMr 91 sequence) were measured both in plasma and in a 3'-exonuclease solution. The in vitro enzymatic degradation in plasma and in a 3'-exonuclease solution are similar indicating that the 3'-exonuclease solution can model the plasma medium.
Phosphodiester metabolism in plasma can be described by first-order kinetics while the phosphorothioate followed first-order kinetics in plasma for only 2 h. During the latter reaction, the reaction rate gradually decreased with time and virtually stopped after the formation of a characteristic right triangle shaped metabolite profile. The same pattern of SODN metabolism was found in plasma and in the 3'-exonuclease solution. Enzyme activity was investigated in both media, and good activity was found for >90 h. Enzyme decomposition during the reaction did not explain the observed decreasing reaction rate. The possibility of SODN-plasma protein interaction was also ruled out as a reason for the decreasing reaction rate. The right triangle shaped metabolite profile is very similar to those typically observed in vivo in plasma. This fact leads us to speculate that this profile is generated by the nature of the metabolic process and does not necessarily depend on the slow efflux of parent drug and metabolites from tissues.
Using the first 2 h to calculate the first-order kinetic constant shows that the 25mer SODN has more than 20 times greater stability than the corresponding 25mer ODN of the same sequence both in plasma and in the 3'-exonuclease solution. Further study of the mechanism of phosphorothioate resistance to enzymatic hydrolysis is ongoing.
REFERENCES
1 Crooke, S.T. and Lebleu, B. (eds) (1993) Antisense Research and Applications. CRC Press, Boca Raton, FL.
2 Zamecnik, P.C. and Stephenson, M.L. (1978) Proc. Natl. Acad. Sci. USA, 75, 280-284.MEDLINE Abstract
3 Zamecnik, P.C., Goodchild, J., Taguchi, Y. and Sarin, P.S. (1986) Proc. Natl. Acad. Sci. USA, 83, 4143-4146.
4 Srinivasan, S.K. and Iversen, P. (1995) J. Clin. Lab. Anal., 9, 129-137.MEDLINE Abstract