Scientists & Organizations
Contributors and affiliated organizations for this study.
| Kristina Flavier | FreeThink Technologies, Inc. |
| Chetan Shende | FreeThink Technologies, Inc. |
| Justin McCabe | FreeThink Technologies, Inc. |
| Kenneth C. Waterman | FreeThink Technologies, Inc. |
| Christian Reidy | Waters Corporation (Milford, MA) |
| Makda Araya | Waters Corporation (Milford, MA) |
| Balasubrahmanyam Addepalli | Waters Corporation (Milford, MA) |
| Koley Hall | Waters Corporation (Golden, CO) |
Abstract
A case study was carried out that confirms the applicability of ASAPprime® for accurate determination of shelf life of lyophilized oligonucleotides using a lipid-conjugated antisense oligonucleotide reference standard from Waters Corporation. This oligonucleotide was exposed to a range of statistically designed conditions, with the resulting purity data found to be well fit using the ASAPprime® process of combining isoconversion times (times to hit the specification limit of 10% loss of purity) with a moisture-modified Arrhenius equation which accounts for both temperature and relative humidity impact.
Introduction
Oligonucleotides are a growing class of drug molecules that are based on nucleic acids1. As with most biological molecules, their stability is sensitive to environmental factors such as temperature, humidity, and oxygen. Lyophilization (freeze drying) is a widely used technique to improve the stability and shelf life of drug products by removing the majority of the water, thereby reducing molecular mobility and corresponding reactivity.
The Accelerated Stability Assessment Program (ASAP) has been successfully applied to numerous small molecule, peptide, and biological drug products2–6. In the ASAP process, samples are stressed at a range of temperatures and relative humidities outside of packaging for up to one month, and the isoconversion time (i.e., time to reach the failure point) is determined at each condition. The isoconversion times are used to calculate the temperature and moisture sensitivities of degradation using a moisture-modified Arrhenius equation (Eq. 1).
(1)
Here, k is the isoconversion rate constant (i.e., the inverse of the time to fail), A is the collision frequency, Ea is the activation energy or temperature sensitivity, R is the gas constant, T is the absolute temperature, B is the moisture sensitivity, and RH is the equilibrium relative humidity to which the sample is exposed. Based on these parameters, the ASAP model is used to predict shelf life at the intended storage conditions and in different packaging configurations.
To verify that ASAP can be used to model the long-term stability behavior of oligonucleotides, FreeThink conducted an ASAP case study using Waters’ antisense oligonucleotide (ASO) lipid-conjugated LC-MS reference standard. Purity loss was monitored using a hydrophilic interaction liquid chromatography (HILIC) UPLC method.
Methods
Materials and Accelerated Stressing
A lipid-conjugated ASO LC-MS reference standard (Part Number 186010747, Table 1) was provided by Waters Corporation as 1 nmol per vial. For each ASAP condition, one vial was placed open into a canning jar with a second vial containing the appropriate saturated salt for humidity control (Table 2). The canning jars were sealed and placed into the appropriate ovens in a staggered fashion such that all samples were removed at the end of the study. Control samples were stored refrigerated and unopened for the duration of the study before being prepared for UPLC analysis along with stressed samples. All samples were analyzed in a pre-established random sequence to minimize any impact of chromatographic drift (systemic bias) during the analysis. Sample purity was calculated as a percentage of the total peak area.
| Component | Elemental Composition | Average Mass (Da) |
|---|---|---|
| 16mer | 5′ d 5-Pal-*-MOE-MeC-*-MOE-G-*-MOE-MeC-*C*G*A*T*A*A*G*G*T*A*-MOE-MeC-*-MOE-A-*-MOE-MeC 3′ | 6050 |
* denotes PS bonds, MOE denotes methoxy ethyl 2′ modifications, and MeC denotes a 5-methyl cytidine residue.
| Temp (°C) | Saturated Salt | % RH | Time (d) |
|---|---|---|---|
| Control | NA | NA | 0 (3 repeats) |
| 25 | LiCl | 11 | 12 |
| 30 | CsF | 3 | 13, 21 |
| 35 | KF | 25 | 6, 21 |
| 40 | LiCl | 11 | 4, 9, 21 |
| 45 | KF | 21 | 2, 5, 16 |
| 50 | LiCl | 11 | 2, 5, 15 |
| 50 | NaI | 29 | 1 |
Chromatography
Stressed and control samples were reconstituted by adding 50 µL of 50:50 acetone:acetonitrile (ACN) and gently vortexing to ensure complete dissolution. A 15 µL aliquot of the resulting solution was transferred to a vial containing 985 µL of 50:50 acetone:ACN and mixed thoroughly. All samples were placed in an autosampler for analysis using an ACQUITY UPLC H-Class System with a Sample Manager with Flow Through Needle (SM-FTN) and a GTxResolve Premier BEH Amide, 300 Å, 1.7 µm, 2.1 × 100 mm column under HILIC conditions7. The sample temperature was 5 °C, the flow rate was 0.3 mL/min, and the needle wash was 50:50 ACN:H2O. The detection wavelength was 260 nm. All data were acquired using Empower™ 3 software. The column temperature was set at 75 °C and the injection volume was 10 µL. The Mobile Phases A and B were 90:10 and 10:90 ACN:100 mM ammonium acetate (v:v), respectively. The gradient was 0, 90:10; 0.67, 90:10; 4.67, 22:78; 5, 22:78; 5.33, 90:10; and 13.33 min, 90:10 A:B.
ASAPprime® Modeling
Chemical purity data for the oligonucleotide were modeled using ASAPprime® software (version 7.0). Calculations used 50,000 Monte Carlo simulations to determine error bars.
Results and Discussion
Loss of purity in the lipid-conjugated ASO was modeled using a default fit (linear or bilinear if a single line has an R2 < 0.90) with a standard deviation of all samples assumed to be 1% based on the reproducibility determined for the controls. Control samples had an average initial purity of 88%, so with an allowed decrease of 10% purity, a specification limit of 78% was used for the ASAP modeling. The fitting of these data at each condition is shown in Figure 1.

The resulting fit parameters to the moisture-modified Arrhenius equation (Eq. 1) are given in Table 3. As can be seen, loss of purity is well fit by this equation, as evidenced by R2 and Q2 values of 0.982 and 0.945, respectively. Q2 is the predicted residual error sum of squares (or PRESS R2) and is a sensitive indicator of how well the fitted data predict points not used in the fitting8. The residuals plots of the model fit (Figure 2) indicate that the experimental isoconversion times (expressed as the logarithm of the inverse isoconversion times, ln k) are close to those derived from the model.
| ln A | Ea (kcal/mol) | B | R2 | Q2 |
|---|---|---|---|---|
| 27.5 ± 3.7 | 17.2 ± 2.2 | 0.066 ± 0.011 | 0.982 | 0.945 |

There is limited published information about the temperature and humidity dependence of oligonucleotide degradation. The International Conference on Harmonisation (ICH) uses an assumption that six months of storage at 25 °C (accelerated conditions) is predictive of a two-year shelf life at 2–8 °C. This equates to an Ea of 11.4 kcal/mol. Loss of purity in the lipid-conjugated ASO exhibits an Ea value of 17 kcal/mol. While this is lower than that of the average small molecule stability (27 kcal/mol), it is still higher than the ICH standard.
The high B term (0.07) versus the average 0.04 for small molecules indicates that degradation is significantly impacted by relative humidity. With lyophilized oligonucleotides often stored in impermeable (glass) packaging, the relative humidity is set by any residual moisture from the lyophilization process with no further moisture ingress. In this case, that residual moisture level will have a large impact on the overall shelf life of the oligonucleotide.
Figure 3 shows the predicted long-term loss of purity for the lipid-conjugated ASO calculated at 5 °C/1% RH and −20 °C. It should be noted that for temperatures below water’s freezing point (0 °C), relative humidity (i.e., the percent of saturated vapor pressure in the air) does not have a definition. As can be seen, the mean shelf life is six months at 5 °C (refrigerated), while at −20 °C (frozen), this oligonucleotide is predicted to be stable for more than three years.

Conclusion
A three-week ASAP case study was conducted on an oligonucleotide lyophile to demonstrate that the chemical stability of RNA-based materials could be effectively modeled using the combination of isoconversion (time to hit the specification limit) and a moisture-modified Arrhenius equation. The mathematical modeling shows with good statistics that loss of oligonucleotide purity is well-fit by a moisture-modified Arrhenius equation, supporting the use of ASAP for this class of molecules. Based on this case study, the lipid-conjugated ASO is predicted to exceed 10% purity loss after approximately six months at 5 °C, even under low humidity conditions. This oligonucleotide therefore requires frozen (−20 °C) storage to maintain acceptable purity. At that condition, the oligonucleotide is stable for more than three years.
References
- Vinjamuri PV, Pan J, Peng P. “A review on commercial oligonucleotide drug products.” J. Pharm. Sci. 2024, 113 (7), 1749–1768. DOI: 10.1016/j.xphs.2024.04.021.
- Qiu F, Scrivens G. “Accelerated predictive stability (APS): fundamentals and pharmaceutical industry practices, 1st ed.” Academic Press. 2018. DOI: 10.1016/C2014-0-02298-8.
- McMahon ME, Abbott A, Babayan Y, Carhart J, Chen CW, Debie E, Fu M, Hoaglund-Hyzer C, Lennard A, Li H, Mazzeo T, McCaig L, Pischel S, Qiu F, Stephens D, Timpano R, Webb D, Wolfe C, Woodlief K, Wu Y. “Considerations for updates to ICH Q1 and Q5C stability guidelines: embracing current technology and risk assessment strategies.” The AAPS Journal. 2021, 23, 1–9. DOI: 10.1208/s12248-021-00641-6.
- Waterman KC, Swanson JT, Lippold BL. “A scientific and statistical analysis of accelerated aging for pharmaceuticals. Part 1: accuracy of fitting methods.” J Pharm Sci. 2014, 103 (10), 3000–3006. DOI: 10.1002/jps.24075.
- González-González O, Ramirez IO, Ramirez BI, O’Connell P, Ballesteros MP, Torrado JJ, Serrano DR. “Drug stability: ICH versus accelerated predictive stability studies.” Pharmaceutics. 2022, 14 (11), 2324. DOI: 10.3390/pharmaceutics14112324.
- Waterman R, Lewis J, Waterman KC. “Accelerated stability modeling for peptides: a case study with bacitracin.” AAPS PharmSciTech. 2017, 18 (5), 1692–1698. DOI: 10.1208/s12249-016-0635-7.
- Finny AS, Johnston T, Sauner B, Addepalli B, Lauber MA. “Reproducible HILIC for Denaturing and Non-Denaturing Oligonucleotide Analysis Using GTxResolve Premier BEH Amide Columns.” Waters Application Note. 2024.
- Nguyen T. Quan, “The Prediction Sum of Squares as a General Measure for Regression Diagnostics.” Journal of Business & Economic Statistics 6, no. 4. 1988. 501–4. DOI: 10.1080/07350015.1988.10509698.

