- Journal List
- HHS Author Manuscripts
- PMC4312256
ROLE OF BENZYL ALCOHOL IN THE UNFOLDING AND AGGREGATION OF INTERFERON α-2A†
Abstract
Benzyl alcohol (BA) is the most widely used antimicrobial preservative in multi-dose protein formulations, and has been shown to cause protein aggregation. Our previous work on a model protein cytochrome c demonstrated that this phenomenon occurs via partial unfolding. Here, we examine the validity of these results by investigating the effect of BA on a pharmaceutically relevant protein, interferon α-2a (IFNA2). IFNA2 therapeutic formulations available on the pharmaceutical market contain BA as a preservative. Isothermal aggregation kinetics and temperature scanning demonstrated that BA induced IFNA2 aggregation in a concentration dependent manner. With increasing concentration of BA, the apparent aggregation temperature of IFNA2 linearly decreased. Denaturant melts measured using intrinsic protein fluorescence and that of the ANS dye indicated that IFNA2 stability decreased with increasing BA concentration populating a partially unfolded intermediate. Changes in NMR chemical shifts and hydrogen exchange rates identified the structural nature of this intermediate, which correlated with an aggregation “hot-spot” predicted by computational methods. These results indicate that BA induces IFNA2 aggregation by partial unfolding rather than global unfolding of the entire protein, and is consistent with our earlier conclusions from model protein studies.
INTRODUCTION
Protein-based pharmaceuticals comprise nearly half of all parenteral products available worldwide1,2. One third of these products are available as multi-dose formulations3, which are advantageous for patient compliance as well as economics. In addition to a shelf-life stability requirement of 18–24 months4, multi-dose formulations also necessitate the inclusion of antimicrobial preservatives (APs)3,5 to combat the growth of bacteria and other microbes during repeated contact between the solution and syringe needle.
In recent years, the study of APs in protein pharmaceutics has become increasingly important due to the fact that these compounds have been shown to cause protein aggregation. One of the first reports demonstrated that the addition of various aromatic and aliphatic alcohols caused aggregation of human growth hormone6. In recent studies, APs have been shown to induce the aggregation of numerous protein systems, including interleukin-1 receptor antagonist7,8, interferon-γ9, and an antibody10. Such protein aggregates in formulations can affect the amount and efficacy of the delivered drug as well as promote undesirable immunologic responses in patients11–17. Therefore, an understanding of the molecular mechanisms underlying AP-induced protein aggregation and the development of strategies to minimize such aggregation are critical in developing stable multi-dose drugs.
Previously published work from our laboratory on the model protein cytochrome c (Cyt c) demonstrated that local unfolding of the protein, rather than the global unfolding, leads to protein aggregation18–20. This local unfolding identified an aggregation “hot-spot” that could be modified to increase protein stability and decrease aggregation. It is important to test these principles learned from model protein studies on a pharmaceutically relevant protein to determine whether these conclusions are general to AP-induced aggregation of proteins. For this purpose, we chose to study interferon α-2a (IFNA2) (Figure 1). IFNA2 has been shown in the literature to aggregate under formulation conditions21–24. Benzyl alcohol (BA) has been used as a preservative in IFNA2 formulations. Whether BA causes IFNA2 aggregation has not been explored. We show that IFNA2 is susceptible to increased aggregation in the presence of BA, and that the aggregation mechanism proceeds through a partially unfolded intermediate. Structural nature of this intermediate has been obtained using 2D NMR, which gives information at the resolution of individual amino acids. Using such high resolution structural techniques is relatively rare in characterizing pharmaceutical proteins.
Molecular structure of interferon α-2a (IFNA2) (PDB ID: 1ITF). The protein is α-helical in nature with helices organized according to color: helix A, residues 11–21 (orange); helix B, residues 52–68 (green); helix B’, residues 70–75 (yellow); helix C, residues 78–100 (purple); helix D, residues 110–132 (cyan); helix E, residues 137–157 (blue). Residues 22–51 comprise the AB-loop, and residues 40–43 form a 310 helix (red).
MATERIALS AND METHODS
Materials
Synthetic cDNA corresponding to IFNA2 was obtained from Operon (Huntsville, Alabama), and was cloned into the pET-SUMO expression vector25,26 (a generous gift from Christopher Lima, Sloan-Kettering Institute). Protein was expressed in Escherichia coli BL21(DE3) cells, and the soluble protein was purified using a Nickel Sepharose 6 Fast Flow column (GE Healthcare Life Sciences, Pittsburgh, Pennsylvania). The SUMO tag was cleaved using the Ulp1 protease, leaving no additional amino acids. The final protein sequence was identical to that of the pharmaceutical protein, which was confirmed by mass spectrometry, circular dichroism, and NMR. Details of IFNA2 expression, purification, and characterization have been described in our earlier publication 27.
Size Exclusion Chromatography
To monitor the effects of BA on protein aggregation, IFNA2 in formulation buffer (10 µM in 0.01 M ammonium acetate, 0.12 M sodium chloride, pH 5) was incubated at 50°C and samples were taken at desired intervals. Samples were centrifuged prior to injection on the HPLC column, and only soluble monomer was detected. Monomer concentration was estimated by injecting 70µL onto a TSKgel 5 µM G3000SWxl column (Tosoh Bioscience LLC, San Francisco, California) on an Agilent 1100 HPLC (Santa Clara, California). The mobile phase used was 0.01 M ammonium acetate, 0.12 M sodium chloride, pH 5, at a flow rate of 1 mL min−1. Absorbance at 280 nm was used to determine the monomer content.
Isothermal Incubation Experiments
IFNA2 (10 µM) in formulation buffer was incubated at the desired temperature, and the changes in optical density at 350 nm were measured as a function of the incubation time28. Buffer and protein do not absorb at this wavelength. The aggregation kinetics were monitored until the signal reached a plateau. At longer incubation times, the aggregates started to settle down to the bottom of the cuvette, resulting in decreased optical density. At that point, the experiment was stopped.
Thermal Scanning Method
The midpoint aggregation temperature (TmAgg) of the protein was measured on an UV-Visible spectrophotometer (Agilent Technologies, Santa Clara, California). The temperature was increased at a rate of 1°C/step followed by 90 sec equilibration, and changes in the optical density at 450 nm were monitored29. TmAgg was determined as the temperature at which half the maximum optical density was reached. For these experiments, 10 µm IFNA2 in formulation buffer was used with varying BA concentration.
Denaturant Melts
Guanidinium chloride (GdmCl) was used as the denaturant. Protein solutions at varying GdmCl concentrations were prepared and equilibrated for one hour at room temperature before measuring changes in fluorescence signals as a function of the denaturant concentration. Concentration of the denaturant was determined using refractive index measurements30. For measuring changes in the 280 nm (aromatic) fluorescence, 3 µM IFNA2 in formulation buffer was used with varying concentrations of BA. Fluorescence spectra were obtained using a Quantamaster spectrofluorometer (Photon Technology International, Birmingham, New Jersey) with 4 nm slits. The ΔG values were determined by fitting the changes in optical signals at different denaturant concentrations to a two-state unfolding model31,32. The m-value, which is the slope of linear variation of ΔG with denaturant concentration, was obtained by fitting the data to a Santero-Bolen two-state equation31,32.
ANS Fluorescence
1-Anilinonaphthalene-8-Sulfonic Acid (ANS) (Sigma, St. Louis, MO) at a concentration of 100 µM was used to monitor partial unfolding of IFNA2. The same experimental setup described above for the denaturant melts was employed for ANS fluorescence experiments. Fluorescence intensity was measured using a Quantamaster spectrofluorometer. Excitation and emission wavelengths used were 360 nm and 480 nm, respectively, and slit widths were 4 nm. Samples were corrected for ANS and BA fluorescence in the absence of protein. All samples were normalized using a rhodamine cuvette (Starna Cells, Atascadero, CA).
NMR
BA Titration
2D 15N-1H heteronuclear single quantum coherence (HSQC) experiments were run to monitor changes in IFNA2 amide crosspeaks as a function of BA concentration using a Varian Inova 600 MHz NMR instrument equipped with a cryoprobe (Palo Alto, California). IFNA2 was singly labeled with 15N using M9 minimal media and expressed using cold induction for 16 hours. For these NMR experiments, 100 µM IFNA2 and a gradient of BA concentrations were used. Samples included 0%, 0.3%, 0.6%, and 0.9% v/v BA. HSQC spectra were collected over eight hours. Changes in the crosspeak positions (chemical shifts) were calculated using the nmrDraw package (Frank Delaglio, National Institutes of Health). IFNA2 residue assignments available in the literature33 were used for this purpose.
Hydrogen Exchange (HX)
2D 15N-1H heteronuclear single quantum coherence (HSQC) experiments were run to monitor changes in IFNA2 amide crosspeaks as a function of BA concentration using a Varian Inova 900 MHz NMR instrument equipped with a cryoprobe (Palo Alto, California). For these experiments, 200 µM IFNA2 (0.05 M deuterated acetic acid, 0.02% sodium azide, pDr 3.5) was used. A time-zero HSQC of IFNA2 was recorded in aqueous buffer (0.05 M acetic acid, 0.02% sodium azide, pH 3.5) prior to exchange. Buffer exchange was carried out on a Q-Sepharose size-exclusion column (bed volume 4.5 ml), and HSQC spectra were recorded every hour. Changes in peak positions and volumes were calculated using the nmrDraw package.
RESULTS
BA Induces IFNA2 Aggregation
Previous work has demonstrated that BA causes aggregation of pharmaceutical proteins over a course of several days and months when stored at 4°C or at room temperature6–10,34–36. Because of these long incubation times, screening the effect of various BA concentrations on protein aggregation is not feasible on a convenient laboratory time scale. Therefore, we accelerated aggregation kinetics by performing isothermal incubation studies at a higher temperature, following the protocols used in earlier studies6–10,34–36. Using elevated temperatures to accelerate protein aggregation is becoming a commonly used method for scanning the effect of a large number of solution conditions on protein stability and aggregation37–39. In addition, our earlier work on the aggregation of the model protein Cyt c showed that the mechanism by which BA and other APs induce protein aggregation is identical at low temperature as well as at high temperature18–20. Using size-exclusion chromatography, we show that BA caused the IFNA2 monomer loss as a function of the incubation time at 50°C (Figure 2A). After one day of incubation in the absence of BA, approximately 90% of IFNA2 monomer remained in solution. A similar result was seen in the presence of 0.3% v/v BA. However, no monomer was detected after eight hours in the presence of 0.9% v/v BA. We then examined the effect of BA concentration on the real-time isothermal kinetics of BA-induced protein aggregation by monitoring the increase in optical density due to increases in solution turbidity28. With the increase of BA concentration, the IFNA2 aggregation kinetics was accelerated (Figure 2B).
Effects of BA on IFNA2 aggregation under isothermal conditions. (A) Soluble monomer remaining in solution during 24 hours of incubation at 50°C in the presence of no BA (closed circles), 0.3% v/v BA (open circles), and 0.9% v/v BA (closed triangles). Error bars indicate triplicate data. Individual time points are normalized with respect to no BA at hour zero. (B) Aggregation kinetics at 42°C with varying concentrations of BA: 0% v/v (closed circles), 0.9% v/v (open circles), 1.5% v/v (closed triangles), and 2% v/v (open triangles). Inset shows absorbance traces for buffer and IFNA2, indicating that neither IFNA2 nor the buffer have absorbance at 350 nm.
BA Decreases the Temperature at Which IFNA2 Aggregates
An alternative method used in the literature to probe the effects of various solution conditions is to monitor their influence on the temperature at which the protein aggregates during thermal scanning6–10,18–20,34–39. For this purpose, we measured the optical density at 450 nm29. At this wavelength, the protein and the buffer do not absorb (Figure 2B Inset), and hence the observed changes in the optical density can be attributed solely to protein aggregation. The signal initially increased with increasing temperature and reached a plateau. At higher temperatures, the optical density decreased as protein aggregates began to settle to the bottom of the cuvette. We performed this thermal scanning experiment in the presence of varying concentrations of BA to measure its effect on the aggregation temperature, TmAgg, of IFNA2 (Figure 3A). With the addition of BA, IFNA2 melted at lower temperatures, resulting in a decrease of the TmAgg. In the absence of BA, the TmAgg was 63.9±0.9°C, whereas the addition of 2% v/v BA decreased the TmAgg to 42.1±0.3°C, indicating that the presence of BA enhanced IFNA2 aggregation.
Aggregation Temperature Decreases Linearly with Increase in BA Concentration
The aggregation temperature, TmAgg, obtained from thermal scanning experiments was plotted as a function of BA concentration (Figure 3B). A linear correlation was observed between the increase in BA concentration and the decrease in TmAgg of IFNA2. This linear correlation is very similar to the commonly observed general protein unfolding behavior where thermodynamic stability (ΔG) or the folding rate (kf) decreases linearly with the addition of denaturant31,32,40. This correlation is significant in that it is a qualitative measure of how effective BA is in inducing aggregation of IFNA2, reflected in the value of the slope of TmAgg decrease with increase in BA concentration. The corresponding slope for BA indicates that for every percent increase of BA, the TmAgg of IFNA2 decreases by 10.5°C (Figure 3B).
BA Causes Destabilization of IFNA2
In order to monitor the effect of BA on IFNA2 thermodynamic stability, we followed changes in protein intrinsic fluorescence as a function of BA and denaturant concentration (Figure 4A). With the addition of increasing concentrations of BA, the calculated ΔG decreased significantly (Figure 4B). For example, in the absence of BA, the calculated ΔG for IFNA2 was 10.2±0.6 kcal/mol. With the addition of 2% v/v BA, the ΔG value dropped to 3.8±0.3 kcal/mol. These denaturant melts indicate that BA causes protein destabilization. Changes in the corresponding m-values were also observed (Table 1). The m-value is a measure of the accessible surface area that is exposed upon protein unfolding41. One possible explanation for the decrease in m-value while maintaining the same midpoint (Cm) is the presence of an intermediate42.
Destabilization of IFNA2 as a function of BA concentration at room temperature. (A)IFNA2 denaturant melts in the presence of varying concentrations of BA: 0% v/v (closed circles), 0.9% v/v (open circles), 1.5% v/v (closed triangles), and 2% v/v (open triangles). (B) The free energy of unfolding (ΔG) of IFNA2 as a function of BA.
Table I
Parameters obtained from fitting denaturant melts in Figure 4 to a 2-state unfolding model. The decrease in m-value with no significant change in Cm implies the presence of an intermediate.
| BA (% v/v) | ΔG (kcal/mol) | m (kcal/mol/M [GdnHCl]) | Cm (M[GdnHCl]) |
|---|---|---|---|
| 0.0% | 10.2 ± 0.6 | −2.6 ± 0.2 | 3.9 ± 0.02 |
| 0.9% | 8.4 ± 0.4 | −2.1 ± 0.1 | 4.0 ± 0.09 |
| 1.5% | 5.5 ± 0.9 | −1.4 ± 0.2 | 3.9 ± 0.02 |
| 2.0% | 3.8 ± 0.3 | −0.9 ± 0.1 | 4.2 ± 0.16 |
BA Preferentially Populates an Intermediate
In order to further test the increased population of an intermediate in the presence of BA, we used the hydrophobic fluorescent dye ANS. ANS fluorescence has been used to detect partially unfolded intermediate states in the presence of GdmCl43. The denaturant concentration at which maximum ANS fluorescence was observed was used for comparative purposes. With increasing concentrations of BA, the denaturant concentration at which ANS intensity was highest shifted to lower values, indicating an earlier onset of partial protein unfolding and increased intermediate population (Figure 5). For example, in the presence of 0% v/v BA, denaturant concentration at which ANS fluorescence was maximum was observed at 3.7 M GdmCl. In contrast, 1.5% v/v BA shifted the maximum ANS fluorescence to 3.0 M GdmCl.
Intermediate detection using ANS as a function of BA concentration in the presence of the denaturant GdmCl at room temperature. (A) Increase in fluorescence emission intensity of ANS as a function of BA concentration, normalized using a rhodamine cuvette. Inset shows representative raw ANS intensity curves as a function of GdmCl concentration at 0% v/v BA (closed circles) and 0.9% v/v BA (open circles). (B) Comparison of midpoint denaturant concentration at maximum ANS intensity with the addition of BA.
BA Perturbs IFNA2 Locally
To better understand the specific effects of BA on IFNA2 structure, 2D 1H-15N HSQC NMR spectra were recorded at various BA concentrations. With an increase in BA concentration, numerous 1H-15N crosspeaks were perturbed, which were monitored through chemical shift changes. These changes were used to identify the local regions in IFNA2 that were most susceptible to influences by BA (Figure 6A). For all concentrations of BA, 1H chemical shift changes were less than 0.1 ppm. However, multiple residues demonstrated significant chemical shift alterations (i.e., greater than 0.1 ppm) in the 15N dimension. The most substantial changes were observed in residues within the AB-loop and helix A (Figures 1 and and6A).6A). For example, between 0% and 0.3% v/v BA, the 15N chemical shift of F27 increased by nearly 0.4 ppm. As BA concentration was increased, changes in chemical shifts were observed across the entire protein, indicating that the unfolding propagates to other protein regions.
2D 15N-1H HSQC NMR spectra of IFNA2 in the presence of BA. (A) Chemical shift changes as a function of BA concentration in the 15N dimension. The changes in crosspeak positions represent the absolute changes from 0% v/v BA. Top panel: 0.3% v/v BA; middle panel: 0.6% v/v BA; bottom panel: 0.9% v/v BA. The black rectangle shows the region with the highest chemical shift changes compared to 0% v/v BA. (B) HX rates in IFNA2 plotted as a function of residue position. The rate values represent the relative increase of individual amide exchange rates in the presence of 0.9 % v/v BA when compared to 0% v/v BA. Colored bars indicate helices: orange (helix A), red (310 helix), green (helix B), yellow (helix B’), purple (helix C), cyan (helix D), dark blue (helix E). Black rectangles show maximum HX rate changes compared to 0% v/v BA.
HX experiments in the presence and absence of BA were conducted in order to monitor the effects of BA on structural stabilities around individual amides (Figure 6B). Nearly all the residues in the AB-loop exchanged within the deadtime of the experiment (2 min), including residue F27, demonstrating that the AB-loop is the most dynamic region of IFNA2. Residues L17, A19, M21, and L30 demonstrated significantly increased exchange rates in the presence of 0.9% v/v BA. These residues are located within the C-terminus of helix A and the AB-loop. Additionally, W76, D82, and K83 in helix C showed increased exchange with the addition of BA. In the three-dimensional structure of IFNA2, L17 docks against K83 and has a contact area of 45.3Å2 (contact maps were determined using the online server http://ligin.weizmann.ac.il/cma/). Residue M21 and W76 have a similar relationship, with a contact area of 33.4Å2 (Figure 7).
Amino acid contacts within IFNA2 three-dimensional structure. (A) Critical residues in the AB-loop are not protected through intra-protein contacts and are highly labile in the native state. (B) Residue L17 in helix A shares a contact area of 43.5 Å2 with K83 in helix C. (C) Residue M21 in helix A shares a contact area of 33.4 Å2 with W76 in helix C.
In light of the HX results that show partial protein unfolding, we used computational programs available in the literature to identify the protein regions that may act as “hot-spots” for aggregation44–46. Protein aggregation is a multi-molecular reaction that is governed by hydrophobic interactions, secondary structural elements, electrostatics, and other physical properties44,47,48. Aggregation prediction programs, such as TANGO44 and AGGRESCAN45, take these parameters into account and determine the specific regions that are aggregation-prone. The 3D analysis program ZipperDB46 uses the available structure of a protein to determine its propensity for forming cross-β structure, the most common structure found in protein aggregates. These programs were used in conjunction to identify three common aggregation hot-spots in IFNA2 that were predicted by all the programs (Figure 8). One of these predicted hot-spots is within the highly dynamic AB-loop, which has very little protection in the three-dimensional structure of IFNA2. Another identified hot-spot is located within helix A, a highly ordered α-helical region. BA appears to affect multiple areas of IFNA2 to varying degrees, but imparts large influence on the AB-loop and within helix A, as evidenced by significant changes observed in chemical shifts and HX rates for many amino acids in this region (Figure 6). The third predicted hot-spot is within helix D, which did not show significant changes in either the titration experiment or HX, demonstrating the importance of performing experimental studies using high-resolution techniques such as NMR rather than relying upon computational predictions.
DISCUSSION
Multi-dose formulations require the inclusion of APs to combat the growth of bacteria and other microbes during administration and storage3–5. However, recent studies have suggested a link between the presence of APs and protein aggregation3–6,9,10,34–36. The ideal method of understanding aggregation mechanisms is to study how various APs induce the aggregation of a pharmaceutical protein of interest using biophysical methods. In our previous work, we used a model protein, Cyt c, which has been well-characterized in the literature, to monitor preservative-induced protein aggregation mechanisms in response to various stresses18–20. We demonstrated that BA, the most commonly used preservative in liquid formulation, caused protein aggregation, and the aggregation was preceded by a local unfolding event rather than global unfolding of the entire protein. This result was significant in that it was the first study to identify the interaction of BA with proteins as mechanistically local in nature, rather than global. Using NMR and other high-resolution techniques provided a clear advantage in structurally identifying the partial unfolding, when compared with other AP-induced aggregation studies performed in the literature.
In the present study, we tested the conclusions obtained from our model system on a pharmaceutically relevant protein interferon α-2a. IFNA2 is one of 13 interferon-α subtypes found in humans and belongs to the cytokine family of proteins, which plays critical roles in the innate immune response. Interferons are naturally produced as a result of viral or microbial infections and activate a non-specific, self-perpetuating response by binding to cell surface receptors to initiate a downstream cascade of anti-proliferative and anti-viral factors49–53. As a result, interferon-based pharmaceuticals are used to treat numerous debilitating diseases including various leukemias and hepatitis C54.
IFN2A therapeutic formulations have been known to aggregate for nearly two decades21–24. Numerous studies have been performed to demonstrate the occurrence of IFNA2 aggregation in therapeutic formulations. For example, in the late 1990s, Hochuli23 tested IFNA2 formulations for aggregates and determined that the inclusion of human serum albumin (HSA) increased the aggregation propensity of IFNA2. As a result, HSA was removed from the formulation design and replaced with polysorbate 80 as a stabilizer. However, IFNA2 formulations still contain a measurable population of aggregates. In this work, we have shown that BA, the AP used in IFNA2 formulations, also causes protein aggregation, adding a new dimension to IFNA2 aggregation phenomenon under formulation conditions. The 2D-NMR solution structure for IFNA2 has already been solved33, which provided a clear advantage in expediting the structural studies to monitor changes on an amino acid level induced by BA.
The influence of BA on IFNA2 stability and aggregation was monitored with numerous biophysical techniques. BA was shown to induce IFNA2 aggregation in a concentration-dependent manner (Figures 2 & 3). Thermodynamic stability (ΔG) of IFNA2 decreased in the presence of increasing concentrations of BA, with a constant Cm but a changing m-value (Figure 4). The m-value is related to the differences in the accessible surface areas of the protein’s folded and unfolded states. Changes in this value suggest a deviation from two-state behavior, and may indicate the presence of a partially unfolded intermediate42. Using the hydrophobic fluorescent dye ANS, we demonstrated that IFNA2 contains such an intermediate that is preferentially populated in the presence of BA (Figure 5). The identity of this intermediate was determined using high resolution 2D-NMR and HX techniques, demonstrating that the AB-loop and helix A were specifically influenced by the presence of BA in solution (Figure 6).
Multiple aggregation propensity programs were utilized to determine potential aggregation “hot-spots” in IFNA2 (Figure 8). Two of the predicted “hot-spots” were significantly perturbed in the presence of BA (Figure 6), and these regions share significant surface contacts in the folded protein structure (Figure 7). It is possible that unfolding of one or both of these regions leads to the formation of the partially unfolded intermediate. BA may be influencing this region by simply disrupting the contacts between these critical residues and initiating a local unfolding event that populates an aggregation-prone species. These data indicate that the AB-loop and the terminus of helix A may play a major role in initiating IFNA2 aggregation.
Another possible interpretation is the disruption of these contacts causes a domain swap, leading to new inter-helical connections in a native-like configuration. Significant chemical shift changes were observed in helix A and the AB-loop, indicating that the amino acids in this region move with the addition of BA. A number of other protein examples demonstrate the relationship between domain swapping and the formation of protein aggregates55,56. Additionally, Cyt c has been shown to assemble into higher order species in the presence of ethanol through the successive swapping of helices from one monomer to another monomer57.
Similar to Cyt c18–20, stabilizing the aggregation “hot-spots” of IFNA2 via specific mutations may reduce BA-induced aggregation. Previous work using alanine scanning demonstrated that a number of residues located in the AB-loop are critical for receptor binding, including F27, L30, and K3158–60. Accordingly, introducing a stabilizing mutation at this location could alter protein activity significantly. On the other hand, IFNA2 has been pegylated in recent formulations in order to increase the drug lifetime by protecting the protein from proteolytic degradation, and has been shown to inhibit aggregation61,62. Pegylation occurs at six locations in IFNA2 including at K31 in the AB-loop. This positional isomer has very low affinity to the receptor compared to unmodified IFNA2 (520nM and 5nM, respectively). However, the opposite effects of pegylation on protein stability and activity seem to cancel each other due to increased net protein concentration available for function. It is possible that mutating residues within the AB-loop may confer stability at the cost of potency. Mutating a pharmaceutical protein to minimize aggregation must be done carefully to ensure that the modified drug product maintains efficacy and safety.
Understanding the specific effects of preservatives on protein systems is critical to the development of safer and more potent protein pharmaceuticals. Our results suggest that the BA influences model proteins and pharmaceutical proteins in a similar manner. In-depth studies on other protein formulations in the presence of APs are required in order to further test the principles we have demonstrated here for IFNA2 and for Cyt c. Probing the specific mechanism of interaction between APs and proteins will be useful in the development of stable pharmaceuticals.
ACKNOWLEDGMENTS
The authors thank John Carpenter, Theodore Randolph, David Bain, and LaToya Jones Braun for many helpful discussions and critical comments. This work was funded by the University of Colorado Skaggs School of Pharmacy and Pharmaceutical Sciences. Regina Bis was partially supported by a NIH Leadership training grant in Pharmaceutical Biotechnology (T32GM008732) and a predoctoral fellowship from the PhRMA Foundation (AWD-120487). Rocky Mountain Regional 900 MHz NMR Facility was funded by the National Institutes of Health instrumentation Grant P41GM068928.
ABBREVIATIONS
| ANS | 1-Anilinonaphthalene-8-Sulfonic Acid |
| AP | antimicrobial preservative |
| BA | benzyl alcohol |
| Cm | midpoint denaturant concentration in protein denaturant melt |
| Cyt c | cytochrome c |
| ΔG | Gibbs free energy of protein unfolding |
| GdmCl | guanidinium chloride |
| HSQC | heteronuclear single quantum coherence |
| HX | hydrogen exchange |
| IFNA2 | interferon alpha-2a |
| m-value | slope of the linear variation of ΔG with denaturant concentration |
| NMR | nuclear magnetic resonance |
| pDr | measured pH of a deuterated solution using a standard hydrogen electrode |
| TmAgg | midpoint aggregation temperature |
Footnotes
This manuscript is an invited contribution to the Special issue of the Journal of Pharmaceutical Sciences on “Two Decades of Publishing Excellence in Pharmaceutical Biotechnology”.
Facebook
Twitter
Google+