Salt Disproportionation in the Solid State: Role of Solubility and

Oct 21, 2016 - Vyazovkin , S.; Wight , C. A. Model-free and model-fitting approaches to kinetic analysis of isothermal and nonisothermal data Thermoch...
0 downloads 0 Views 3MB Size
Article pubs.acs.org/molecularpharmaceutics

Salt Disproportionation in the Solid State: Role of Solubility and Counterion Volatility Naveen K. Thakral,†,‡,§ Robert J. Behme,†,∥ Aktham Aburub,† Jeffrey A. Peterson,† Timothy A. Woods,† Benjamin A. Diseroad,† Raj Suryanarayanan,‡ and Gregory A. Stephenson*,† †

Lilly Corporate Center, Eli Lilly and Company, Indianapolis, Indiana 46285, United States Department of Pharmaceutics, College of Pharmacy, University of Minnesota, Minneapolis, Minnesota 55455, United States



S Supporting Information *

ABSTRACT: Disproportionation propensity of salts (HCl, HBr, heminapadisylate) and adipic acid cocrystal of corticotropin releasing hormone receptor-1 antagonist was studied using model free kinetics. Using thermogravimetic weight loss profile or heat flow curves from differential scanning calorimetry, an activation energy plot for salts and cocrystal was generated based on model free kinetics. This activation energy of disproportionation provided qualitative information about the solid state salt stability. To ensure the stability throughout the shelf life, “prototype” formulations of salts and cocrystal in tablet form were stored at 40 °C and several water vapor pressures. Disproportionation kinetics were studied in these prototype tablet formulations using two-dimensional X-ray diffractometry. Formulations containing the adipic acid cocrystal or heminapadisylate salt did not show disproportionation of API when stored at 40 °C/75% RH for 300 days. On the other hand, formulations containing HCl or HBr salt disproportionated. Though isostructural, the disproportionation propensity of HBr and HCl salts was quite different. The HCl salt highlighted the important role that volatility of the counterion plays in the physical stability of the formulations. Solution state stability (i.e., in dissolution medium) of salts and cocrystal was also assessed and compared with solid state stability, by determining their solubility at different pH’s, and intrinsic dissolution rate. KEYWORDS: salt, cocrystal, disproportionation, model free kinetics, counterion, volatilization, solubility, intrinsic dissolution, X-ray diffractometry



back to the less soluble, un-ionized form.6 This reaction can be brought about when the pH of the aqueous microenvironment is more basic than the pHmax of the salt (of the weakly basic drug). Trace levels of basic impurities in the excipients, for example in magnesium stearate, can cause such a pH shift. Salt disproportionation may adversely affect both bioavailability and chemical stability of an active pharmaceutical ingredient (API).7,8 Early formulations for “first-in-human” use are quite simple: commonly drug in capsule. While selecting the form of the API (this term encompasses both the physical and chemical form) for commercial development, physical and chemical stability typically over a shelf life of two years must be ensured.9 Disproportionation issues detected in late stages of development can dramatically increase the cost of development and may result in significant delays.

INTRODUCTION With the advent of high throughput screening (HTS) and combinatorial chemistry in sourcing the drug leads, empirical screening methods have become less important.1 HTS methods are efficient in selection of new drug candidates, but often result in high molecular weight lipophilic molecules with low aqueous solubility.2 An analysis of approved drugs by the FDA in the past 30 years indicated a trend of continuous rise in lipophilicity (log P), and decrease in aqueous solubility.3 Following oral administration, absorption of these compounds is expected to be dissolution rate and/or solubility limited. These active pharmaceutical ingredients are usually weak electrolytes, and their corresponding acidic or basic salts exhibit improved aqueous solubility, and hence potentially improved oral absorption.4,5 Also, during the “lead optimization stage” of drug discovery, the drug candidate must have adequate aqueous solubility to achieve sufficiently high exposure to observe a toxic effect in animals to establish its margin of safety. This can be achieved by making salts of the weak electrolytes. However, due to processing conditions or reaction with formulation ingredients, salts disproportionate, that is, dissociate and revert © 2016 American Chemical Society

Received: Revised: Accepted: Published: 4141

August 11, 2016 October 3, 2016 October 20, 2016 October 21, 2016 DOI: 10.1021/acs.molpharmaceut.6b00745 Mol. Pharmaceutics 2016, 13, 4141−4151

Molecular Pharmaceutics



During lead optimization, when only a limited amount of material is available, simple slurry of the candidate salt form in water, measurement of the solution pH, and identification of the form in equilibrium with water can be useful. Similarly, pH−solubility profile of a candidate salt form, where the residual crystalline form in equilibrium with water is identified, provides an early indication of the potential for salt disproportionation. We further advocate the use of “prototype” formulations to provide more meaningful information, as disproportionation in our experience is observed within the first 2 weeks and almost certainly within one month of equilibration.9 Thermal analysis using model free kinetics10−15 may provide qualitative information about the solid state salt stability, and can be used during formulation development to ensure the stability of salt throughout the shelf life of the formulation. Similarly, our current simplistic model based primarily upon solution state equilibria within the water contained in the formulation and the solubility of the constituents of the formulation, has been useful in predicting disproportionation in the solid formulations.16 There are limitations in our understanding of the chemistry of the system that need to be incorporated into the model if it is to be relied upon. The appropriate modifications to the model will be the subject of a future publication. In this investigation, we will specifically address the role of counterion volatility on salt disproportionation reaction. In the present study, a poorly water-soluble, weakly basic compound (corticotropin releasing hormone receptor-1 antagonist; CRH-1 (Figure 1), three of its salt forms (HCl, HBr, heminapadisylate) and a cocrystal (adipic acid cocrystal) were used to study disproportionation/dissociation occurring in the solid state.

Article

EXPERIMENTAL SECTION

Materials. CRH-1 free base was synthesized by Eli Lilly and Company’s Research Laboratories and provided for the study. The HBr salt was prepared by suspending 545 mg of free base in 3 mL of isopropyl alcohol and stirring at 60 °C. To this mixture was added 150 μL of 48% HBr acid (pKa −9), and it clarified quickly before solid precipitation. The precipitation continued with the addition of excess diethyl ether. The mixture was slurried and evaporated at ambient temperature. The heminapadisylate salt was prepared by suspending 528.8 mg of free base in 2 mL of acetone to which 340 μL of 3.8 M naphthalene-1,5-disulfonic acid (pKa −3) in methanol was added. Within 30 s, there was precipitation from solution, and 2 mL of methanol was added with further precipitation noted. The liquid was decanted and solids were dried under nitrogen. The solids were suspended in 10 mL of ethanol and slurried overnight at ambient temperature before vacuum filtration and drying. The HCl salt was prepared by addition of 47 μL of HCl acid (pKa −7) to an acetone solution of the free base (420 mg in 2 mL) at a low temperature (0−10 °C). The solution was then warmed to 20 °C and filtered. The adipic acid cocrystal was prepared by suspending 5.8 g of free base in 11.5 mL of ethyl acetate. To this slurry was added 2.02 g of adipic acid in 1:4 acetone:methanol (v/v) solution. The clear solution was then filtered into a vial containing seeds of the adipic acid cocrystal. While some seeds persisted, the bulk dissolved so the liquor was concentrated at 55 °C by ∼50%. The clarified solution was filtered again, seeded again, and cooled to rt overnight. Approximately 1.5 g of adipic acid cocrystal formed along the bottom of the vial and was isolated and dried under vacuum at 65 °C. Microcrystalline cellulose (MCC; Avicel PH-102; FMC BioPolymer), hydroxypropyl methyl cellulose (HPMC; Methocel K4M; Dow Chemical Company), polyvinylpyrrolidone (PVP; Kollidon 30; BASF), and magnesium stearate (Fisher Scientific Company; lot # 740042), were used as received. Thermogravimetric Analysis with Simultaneous Differential Scanning Calorimetry (TGA/sDSC) or StandAlone DSC. Samples were heated usually in open 40 μL aluminum pans without lids (Mettler Toledo TGA/SDTA851e) except in one case where a lid was hermetically sealed. TGA results were acquired at a variety of heating rates from as low as 2 °C/min to as high as 20 °C/min. For each TGA weight loss profile at a desired heating rate, an analysis of an empty pan under the same conditions was subtracted from the sample TGA thermal curve. Performance qualification of the TGA balance was conducted with standard weights. The experimental weight loss observed using sodium tartrate dihydrate and barium chloride dihydrate as standards was in excellent agreement with the theoretical weight loss. The TGA temperature axis was calibrated using the melting points of indium, tin, and zinc melting point standards assessed by simultaneous DSC analysis. Samples were heated at 2, 5, 10, 20, and 100 °C/min. Performance qualification of the heat flow axis and the temperature axis utilized indium, tin, and zinc at heating rates of 2, 5, 10, 20, 50, 100, and 200 °C/min. The onset melting point and apparent heat of fusion data from multiple melting point standards at all heating rates were entered into the FlexCal software. The elaborate calibration procedure ensured that regardless of heating rate over wide temperature ranges any instrumental

Figure 1. Corticotropin releasing hormone receptor-1 antagonist, 4-(4-chloro-5-(2,6-dimethyl-8-(pentan-3-yl)imidazo[1,2-b]pyridazin3-yl)thiazol-2-yl)morpholine; (C20H26ClN5OS), molecular weight 419.98.

One particularly interesting aspect of the study is that two of the salts, the hydrobromide and the hydrochloride, are isostructural. This affords the rather unique opportunity to compare the disproportionation tendency of two salts from structures of nearly identical packing and highlight the important role that volatility of the counter substance plays in the physical stability in formulations. The objectives of the present research were (i) to assess/ predict the disproportionation in dry neat salts by model free kinetics using thermogravimetric/differential scanning calorimetry data, (ii) study the disproportionation kinetics in the solid state in “prototype” formulations containing different salts and the cocrystal of CRH-1, at elevated temperature and humidity conditions, (iii) generate single crystal data of HCl and HBr salts of CRH-1 to verify their isostructural nature, and (iv) study the effect of counterion volatility on disproportionation kinetics. 4142

DOI: 10.1021/acs.molpharmaceut.6b00745 Mol. Pharmaceutics 2016, 13, 4141−4151

Article

Molecular Pharmaceutics

The structure was solved and refined using the Bruker SHELXTL Software Package, using the space group P2(1)c, with Z = 4 for the formula unit, C20H27BrClN5OS. The final anisotropic full-matrix least-squares refinement on F2 with 271 variables converged at R1 = 7.08% for the observed data and wR2 = 22.16% for all data. The goodness-of-fit was 1.025. The largest peak in the final difference electron density synthesis was 0.687 e−/Å3, and the largest hole was −0.495 e−/Å3 with an RMS deviation of 0.088 e−/Å3. On the basis of the final model, the calculated density was 1.400 g/cm3 and F(000), 1032 e−. The pentyl group is disordered within the crystal structure, giving rise to a high peak in the final difference map and relatively high wR2. The remainder of the molecule and the counterion refined well. The structure was sufficient for reasonable calculation of its theoretical X-ray powder diffraction pattern and calculation of its true density, and to confirm the isostructural nature of the HCl and HBr salts. The structure parameters at room temperature are reported. Crystal Structure Determination: HCl Salt. A clear, colorless, prismatic-like specimen of C20H27Cl2N5OS, approximate dimensions 0.020 mm × 0.070 mm × 0.080 mm, was used for the X-ray crystallographic analysis. The data set was collected at 296 K using a copper source and wavelength of 1.5418 Å. The X-ray intensity data were measured. A total of 2832 frames were collected. The total exposure time was 2.36 h. The frames were integrated with the Bruker SAINT software package using a narrow-frame algorithm. The integration of the data using a monoclinic unit cell yielded a total of 10341 reflections to a maximum θ angle of 64.22° (0.86 Å resolution), of which 3704 were independent (average redundancy 2.792, completeness = 95.9%, Rint = 6.56%, Rsig = 8.39%) and 2525 (68.17%) were greater than 2σ(F2). The final cell constants of a = 13.1687(4) Å, b = 16.6379(5) Å, c = 11.5914(3) Å, β = 114.243(2)°, and volume = 2315.70(12) Å3 are based upon the refinement of the XYZ centroids of 2321 reflections above 20σ(I) with 9.914° < 2θ < 120.4°. Data were corrected for absorption effects using the multiscan method (SADABS). The ratio of minimum to maximum apparent transmission was 0.653. The calculated minimum and maximum transmission coefficients (based on crystal size) are 0.7777 and 0.9490. The structure was solved and refined using the Bruker SHELXTL Software Package, using the space group P21/c, with Z = 4 for the formula unit, C20H27Cl2N5OS. The final anisotropic full-matrix least-squares refinement on F2 with 271 variables converged at R1 = 9.18% for the observed data and wR2 = 29.72% for all data. The goodness-of-fit was 1.087. The largest peak in the final difference electron density synthesis was 0.959 e−/Å3 and the largest hole was −0.482 e−/Å3 with an rms deviation of 0.159 e−/Å3. On the basis of the final model, the calculated density was 1.309 g/cm3 and F(000), 960 e−. The pentyl group is disordered within the crystal structure, giving rise to a high peak in the final difference map and relatively high wR2. There was no improvement by collecting data for a longer time or at a lower temperature. The rest of the molecule and counterion refined well, and the site of protonation was found in the difference map. The structure was not deemed of sufficient quality for deposition in the Cambridge Structural Database. The structure was sufficient for reasonable calculation of its theoretical X-ray powder diffraction pattern and calculation of its true density, and to confirm the isostructural nature of the HCl and HBr salts. The room temperature structure’s parameters were reported.

thermal lag was compensated, and differences detected for materials could be attributed to kinetics. Model Free Kinetics (MFK) and Software. Model free kinetics, sometimes called isoconversion kinetics, was applied to the conversion curves of the TGA weight loss profiles (or DSC profiles) obtained at different rates. The MFK module calculates apparent activation energy for a family of curves obtained at different heating rates at each isoconversion point. The MFK application is a module within the Mettler Toledo software suite (STARe software; version 13.00a) that assumes that the kinetic process is the same across the family of curves at each isoconversion point. However, it does not assume a single process across the entire conversion profile. Weight loss profiles at different heating rates were plotted together. In each thermal curve, the weight loss profiles were truncated to focus on the specific weight loss region of interest. The software enabled normalization of the data by converting the weight loss to percent conversion profiles. Model free kinetics software calculated activation energies at each isoconversion point and plotted the isoconversion activation energy values versus percent converted. Model free kinetics in the advanced MFK software option requires at least three curves based on different temperature programs. The curves can be isothermal or dynamic or include a combination of isothermal and dynamic segments. The curves are evaluated by numerically calculating the minimum of the following integral. n

I(Eα) =

n

∑∑ i=1 j≠1

J1{Eα , Ti(tα)} J2 {Eα , Tj(tα)}

where J1{Eα , Ti(tα)} =



∫t −Δα e E /(RT(t)) dt α

i

α

J2 {Eα , Tj(tα)} =



∫t −Δα e E /(RT (t)) dt α

j

α

where T is temperature [K], t is time [sec], α is the fraction converted, J1 and J2 are subintegrals, Eα is activation energy as a function of conversion, and R is the universal gas constant. Crystal Structure Determination. A clear, colorless, prismatic-like specimen of C20H27BrClN5OS, approximate dimensions 0.060 mm × 0.070 mm × 0.080 mm, was used for the X-ray crystallographic analysis. The data set was collected at 296 K using a copper source and wavelength of 1.5418 Å. The total exposure time was 2.82 h. A total of 3383 frames were collected and integrated with the Bruker SAINT software package using a narrow-frame algorithm. The integration of the data using a monoclinic unit cell yielded a total of 14134 reflections to a maximum θ value of 66.71° (0.84 Å resolution), of which 4023 were independent (average redundancy 3.513, completeness = 95.8%, Rint = 6.28%, Rsig = 5.98%) and 3238 (80.49%) were greater than 2σ(F2). The final cell constants of a = 13.3301(5) Å, b = 16.8060(6) Å, c = 11.6385(4) Å, β = 114.292(2)°, and volume = 2376.47(15) Å3 are based upon the refinement of the XYZ centroids of 6111 reflections above 20σ(I) with 7.275° < 2θ < 130.3°. Data were corrected for absorption effects using the multiscan method (SADABS). The ratio of minimum to maximum apparent transmission was 0.696. The calculated minimum and maximum transmission coefficients (based on crystal size) are 0.7210 and 0.7792. 4143

DOI: 10.1021/acs.molpharmaceut.6b00745 Mol. Pharmaceutics 2016, 13, 4141−4151

Article

Molecular Pharmaceutics

Two Dimensional X-ray Diffractometry (2D-XRD). Intact tablets were exposed, at room temperature, to Co Kα radiation (1.78899 Å; 35 kV × 40 mA) in a two-dimensional X-ray diffractometer (D8 Discover 2D, Bruker with a 140 mm diameter window VÅNTEC-500 detector). XRD patterns were collected for 100 s, using a 0.8 mm collimator set at 4° angle of incidence and an area detector (angular range 36°) set at an angle of diffraction at 8° 2θ. The irradiated area can be described by an ellipse with a major axis of 13.91 mm and minor axis of 0.97 mm. To keep the X-ray irradiated area of tablet constant, five random points on the centerline of the upper surface of each tablet were analyzed at each time point (Figure S1). Since no internal standard was used, it was necessary to detect and correct for any short- and long-term instrumental drift. On each day before and at the end of the experiment, the instrument was calibrated using corundum (SRM 1976a, NIST). By setting the Ω at 20°, and 2θ at 40° (scan time 50 s), the integrated intensity of the 41.02° 2θ (2.55 Å) peak was determined. The coefficient of variation of such measurements was found to be 0.92%. There was no measurable short-term instrumental drift during the time of analysis of each sample. Quantification of Free Base in Intact Tablets: Generation of Standard Curves. Tablets containing known compositions of free base and salt (0−40% w/w; total concentration 40% w/w) and placebo (60% w/w) were prepared (11 compositions; 3 tablets of each composition). Two dimensional XRD was used to collect patterns from 5 locations (along the center line of the top surface) in each of the 3 tablets. The intensity (integrated) of the 7.6° 2θ peak of the free base was determined and averaged for the 15 determinations and plotted as a function of free base concentration. This experiment was performed separately for HCl and HBr salts. The relationship between integrated peak intensity (y) and free base fraction (x; % w/w) was expressed as y = 2.1024x + 0.9086 in the case of HBr salt and y = 2.0834x + 5.1613 in the case of HCl salt. These standard curves were used to determine the free base concentrations, formed due to disproportionation of HBr or HCl salt in tablets. Determination of Microenvironmental pH by Slurry Method. To determine the microenvironmental pH, 500 mg of HBr or HCl salt was suspended in 5 mL of distilled water. On similar lines, slurries of the formulations of both salts were prepared by suspending 500 mg of each formulation in 5 mL of distilled water. The pH of the slurry was measured (pH meter, Mettler Toledo) in an instrument calibrated with buffers of pH 2, 4, 7, and 10. After measuring the pH at time “0”, each slurry vial was stored in an oven maintained at 40 °C. The vials were withdrawn at specific time intervals, and the pH was measured. pKa Determination. Neat compound was weighed into a sample vial and titrated in triplicate from pH 2 to 10 using Sirius Automated Titrator Instrument (GLpKa with DPAS) in the presence of cosolvent (methanol). The approximate wt % of methanol for the three titrations was 27, 36, and 47. Yasuda− Shedlovsky extrapolation was used to extrapolate the pKa values obtained using different concentrations of cosolvent to 0% cosolvent.19 Solubility Determination. Equilibrium solubility of CRH-1 salts and cocrystals was determined at different pH values. For solubility determination at lower pH values (≤3), different molar concentrations of corresponding acids forming counterions were used. One exception to this was in the case of adipic acid cocrystal, where the solubility study below the pHmax was

Synchrotron X-ray Diffractometry (SXRD). The powder samples, ∼20 mg, were filled in aluminum pans (used for differential scanning calorimetry, TA Instruments), crimped hermetically at room temperature, and placed in a specially fabricated sample holder. Experiments were performed in the transmission mode in the 17-BM-B beamline at Argonne National Laboratory (Argonne, IL, USA). A monochromatic X-ray beam [wavelength 0.72808 Å; beam size 250 μm (horizontal) × 160 μm (vertical)] and a two-dimensional area detector (XRD-1621, PerkinElmer) were used. A triplebounce channel-cut Si single crystal monochromator with [111] faces polished was used, which limited the line broadening to its theoretical low limit, i.e., the Darwin width. The flux of the incident X-ray was 8 × 1011 photons/s at 17 keV. Calibration was performed using an Al2O3 standard (SRM 674a, NIST). Using a stepper motor, the sample was oscillated (±1 mm from the center along the horizontal axis) during data collection. Each sample was scanned 30 times, with an exposure time of 1 s for each scan, and the results were averaged. The raw images were integrated to yield one-dimensional d-spacing (Å) or 2θ (degree) scans using the FIT2D software developed by A. P. Hammersley of the European Synchrotron Radiation Facility.17,18 Commercially available software (JADE 2010, Material Data, Inc.) was used for determining the integrated peak intensities. Preparation of Tablets. APIs and excipients were passed through sieve # 170 (D90 = 90 μm). API (40% w/w) and placebo (60% w/w) were mixed manually. The placebo composition was MCC (70.0% w/w), HPMC (23.5% w/w), PVP (5.0% w/w), and magnesium stearate (1.5% w/w). The blend (100 mg) was filled in a circular stainless steel die (8 mm). Tablets were prepared using a universal material testing machine (Zwick/Roell, Zwick GmbH & Co., KG, Ulm, Germany), equipped with 8 mm diameter flat-faced punches. During compression, the lower punch was stationary while the upper punch moved at a constant speed of 10 mm/min. Tablets were compressed under ambient conditions (25 °C; 45% RH) to 100 MPa. The upper punch displacement was measured using a linear variable differential transformer (LVDT) position gauge, attached to the compression platens. The upper punch was allowed to retract, as soon as the desired compression force was achieved. For ejection of the compressed tablets from the die, the lower punch was removed and a force was applied to the upper punch using a hydraulic press (Carver model C laboratory press, Menomonee Falls, WI, USA). Three different tablets of each formulation composition were prepared. Solid State Stressed Stability Testing. Three tablets of each formulation were kept in open pans in glass chambers maintained at different relative humidity (RH) conditions, using saturated salt solutions (lithium chloride, 11% RH; potassium acetate, 22% RH; magnesium chloride hexahydrate, 32% RH; sodium chloride, 75% RH) at 40 °C. Tablets containing the HCl salt were stored at 11, 32, and 75% RH conditions, tablets containing HBr salt were stored at 22 and 75% RH conditions, while tablets containing heminapadisylate salts and adipic acid cocrystal were stored only at 75% RH. Tablets were withdrawn at specific time intervals and were analyzed using two-dimensional X-ray diffractometry. After analyses, the tablets were placed back in the chamber. The top surface of each tablet was marked so as to analyze the exposed surface at each time point. 4144

DOI: 10.1021/acs.molpharmaceut.6b00745 Mol. Pharmaceutics 2016, 13, 4141−4151

Article

Molecular Pharmaceutics determined in HCl solution as it was not possible to achieve low pH using adipic acid. Solubility was also determined in phosphate buffer (pH 4, 6, and 8) and deionized water. Excess solids were weighed in HPLC vials, and solvent was added. Samples were rotated using a laboratory shaker (Barnstead Thermolyne Labquake Rotisserie Shaker) for 24 h at RT. After 24 h, samples were filtered using Ultrafree-MC centrifugal filter devices (0.22 μm) using a centrifuge (5415C Eppendorf centrifuge) for 5 min. The pH of the filtered samples was determined using a calibrated pH meter (Orion 7110 BH), dilutions were made using 50/50 (acetonitrile/water) + 0.1% trifluoroacetic acid, and concentrations were determined using HPLC. Intrinsic Dissolution Rate (IDR). To assess the effect of salt and cocrystal formation on dissolution rate, CRH-1, its different salts and cocrystal were subjected to IDR analysis. An accurately weighed quantity of free base, salts, or cocrystal (n = 3, each) was transferred into the die cavity of a Wood’s apparatus. The powder was compressed in a hydraulic press (Carver press) at 2,000 psi for 1 min. The pellets thus formed were mounted on the shafts of a dissolution apparatus (Varian, VK7000). The dissolution medium used was 500 mL of 0.01 N HCl (37 °C; vacuum degassed before using), and the experiment was conducted for 4 h at 100 rpm. Data was collected every minute up to 4 h using an inline UV probe (head size 5 mm; wavelength 285 nm).

Figure 2. pH−solubility profiles (HCl, HBr, heminapadisylate salts, and adipic acid cocrystal). The curves were generated from the experimental data.

when pH < pHmax, and can be described mathematically for the HCl and HBr salts as shown below. Region 2, solubility description when pH < pHmax for a diprotic weak base:



RESULTS AND DISCUSSION Solution Stability of the Solid State Forms of CRH-1. The solubility of a weak base is influenced by its solubility product (Ksp), pKa(s), intrinsic solubility (free base solubility), and pH. To calculate the pH−solubility profile of a weak base, only the solubility product, ionization constant(s), and intrinsic solubility are needed to be known. Region 1 of the pH−solubility profile of a diprotic weak base is where solubility of the free base (intrinsic solubility) is limiting, that is, when pH > pHmax, and can be described mathematically as shown below (details in Supporting Information). Region 1, solubility description when pH > pHmax for a diprotic weak base: S = [BaseH+2 2] + [BaseH+] + [Base]

(1)

⎛ H+2 ⎞ H+ S = [Base]⎜⎜ + + 1⎟⎟ K a2 ⎝ K a1K a2 ⎠

(2)

(4)

⎛ ⎞ K a1H+ ⎟[Counterion−] K sp = S⎜⎜ 2 ⎟ + + H K H K K + + ⎝ a1 a1 a2 ⎠

(5)

⎛ H+2 + K H+ + K K ⎞ 1 a1 a1 a2 ⎟ S = K sp⎜⎜ + ⎟ [Counterion−] K H a1 ⎝ ⎠

(6)



[Counterion ] = S + [HC]

(7)

⎛ 10−2pH + 10−(pKa1+ pH) + 10−(pKa1+ pKa2) ⎞ ⎟⎟ S = K sp⎜⎜ ⎝ ⎠ 10−(pKa1+ pH) 1 × [Counterion−]

(8)

Here [HC] is the concentration of acid (i.e., HCl, HBr) in which the solubility experiment was conducted. Since S was determined experimentally, and the concentration of acid in which the solubility experiment was conducted is known, the concentration of counterion as a function of pH was plotted (Figure S2) and an equation that describes the concentration of counterion as a function of pH was developed. This equation (concentration of counterion as a function of pH) was then used to replace [Counterion−] in eq 8. Since the napadisylate salt is a hemisalt, two drug molecules are associated with one napadisylate molecule. This is mathematically described below:

As [Base] = S0, S = S0(10(pKa1+ pKa2 − 2pH) + 10(pKa2 − pH) + 1)

K sp = [BaseH+][Counterion−]

(3)

where S is solubility, S0 is intrinsic solubility, and Ka1 and Ka2 are the acid dissociation constants of the diprotic weak base. The pH−solubility profile is identical, irrespective of the salt (HCl, HBr, heminapadisylate) or cocrystal (adipic acid), when pH > pHmax (Figure 2). The data was fitted (Excel Solver function, MS Office 2003) using eq 3 and pKa1, pKa2, and intrinsic solubility were estimated to be 0.6, 4.8, and 0.0016 mg/mL; respectively. The pKa values are consistent with those determined using potentiometric titration and in silico (0.4 and 4.8 for pKa1 and pKa2, respectively). Region 2 of the pH−solubility profile of a diprotic weak base is where solubility of the protonated form is limiting, that is,

K sp = [BaseH+]2 [Counterion−2] ⎛

⎞2 ⎟ [Counterion−2] K sp = S ⎜ 2 ⎟ + + ⎝ H + K a1H + K a1K a2 ⎠ 2⎜

4145

(9)

K a1H+

(10)

DOI: 10.1021/acs.molpharmaceut.6b00745 Mol. Pharmaceutics 2016, 13, 4141−4151

Article

Molecular Pharmaceutics 2 ⎛ ⎞0.5⎛ H+ + K a1H+ + K a1K a2 ⎞ K sp ⎟ ⎜ S=⎜ ⎟ ⎟ K a1H+ ⎝ [Counterion−2] ⎠ ⎜⎝ ⎠

(11)

[Counterion−] = (0.5)S + [H2C]

(12)

of the HCl salt is likely due to common ion effect (0.01 N HCl dissolution medium). The IDR of the HCl salt could be described by two linear profiles (Figure 3a). From 0 to 10 min, the drug dissolution rate was ∼101 μg/min, but during the time period of 20− 240 min, its dissolution rate dropped to ∼20 μg/min. The dissolution rate during 20−240 min was similar to that of the free base, which was ∼16 μg/min. This strongly suggested that HCl salt disproportionated to the free base after ∼10−20 min in the dissolution medium. A similar trend was observed with the HBr salt (Figure 3b). In Figure 3b, IDR data is provided only up to 90 min to enable the visualization of the disproportionation/dissociation. Complete IDR profiles are presented in Figure S3. The IDR profile of heminapadisylate salt does not indicate disproportionation in dissolution medium, as seen in solid state stability studies, whereas the HCl and HBr salt forms and the adipic acid cocrystal forms do disproportionate/dissociate shortly after exposure to the dissolution medium. In the case of HCl and HBr salts, there is a substantial difference in their solubilities, as compared to free base. The pHmax of the HCl and HBr salts (1.0 and 1.3 respectively) is lower than the pH of the dissolution medium (2.0), hence the free base form is the thermodynamically favored form. On the other hand, pHmax of heminapadisylate salt is 2.5, which is higher than that of the dissolution medium. Consequently, the heminapadisylate salt form is thermodynamically favored to exist in its ionized salt form in the dissolution medium. Interestingly, the adipic acid cocrystal has the coformer present in an un-ionized form, and as such it does not have a true pHmax per se. Under these conditions, adipic acid is more soluble than the free base and results in its transformation to the less soluble free base form. Studies of CRH-1 Physical Stability in the Solid State. X-ray Diffraction. The two-dimensional SXRD patterns and corresponding one-dimensional patterns of CRH-1 free base, HCl, HBr, heminapadisylate salts, and adipic acid cocrystal are depicted in Figures S4A and S4B, respectively. The calculated XRD patterns of HCl and HBr salts are presented in Figure S4C. The HCl and the HBr salts of CRH-1 are isostructural, wherein the unit cell dimensions and the molecular arrangement within the unit cell are virtually identical, with the exception that in one structure the counterion is a chloride ion and in the other it is a bromide ion (Figure 4). Due to the lack of hydrogen bond donating groups in the molecule itself, the only formal hydrogen bond found is between the counterion and the

⎛ ⎞0.5 K sp S=⎜ ⎟ ⎝ [Counterion−2] ⎠ ⎛ 10−2pH + 10−(pKa1+ pH) + 10−(pKa1+ pKa2) ⎞ ⎟⎟ ×⎜⎜ ⎝ ⎠ 10−(pKa1+ pH)

(13)

Here [H2C] is the concentration of diprotic napadisylic acid in which solubility experiment for heminapadisylate salt was conducted. Below pHmax, the pH−solubility profiles of the HCl and HBr salts were generated (Excel Solver function, MS Office 2003) using eq 8. The pH−solubility profile of the heminapadisylate salt was generated using eq 13. A summary of the solubility product, pHmax, and acid dissociation constants can be seen in Table 1. Table 1. Solubility Product, pHmax, and pKa(s) for the HCl, HBr, and Heminapadisylate Saltsa salt

pHmax

Ksp

pKa1

pKa2

HCl HBr heminapadisylate

1.0 1.3 2.5

2.8 × 10−3 [M]2 7.7 × 10−4 [M]2 8.0 × 10−10 [M]3

0.4 0.6 2.1

4.6 4.8 4.8

Ksp and pKa(s) values were estimated from fitting pH−solubility data using eqs 8 (HCl and HBr salts) and 13 (heminapadisylate salt). a

The pH−solubility profile of the adipic acid cocrystal overlaps with that of the HCl (Figure 2). This is expected above pHmax since the observed solubility is intrinsic solubility limited. Below pHmax, the adipic acid cocrystal solubility study was conducted in HCl solutions since it was not possible to achieve low pH using adipic acid (adipic acid has low solubility). Hence, the cocrystal was converted to the HCl salt. The intrinsic dissolution rates (IDR) of CRH-1 free base, its salts, and cocrystal were determined. This provided a measure of their propensity for disproportionation during dissolution. The IDR (μg/min·cm2; 1−10 min) of CRH-1 was found to be 40.8 (±1.63); for HCl salt 202.4 (±15.31); for HBr salt 232.8 (±13.46); for heminapadisyalte salt 63.9 (±0.48); and for adipic acid cocrystal 164.3 (±7.71). Though the solubility of the HCl salt was more than that of the HBr salt, the lower IDR

Figure 3. (a) IDR of CRH-1-HCl and CRH-1 free base. Disproportionation was observed in HCl salt after ∼10−20 min during IDR analysis. (b) Comparison of IDR of CRH-1 salts and adipic acid cocrystal up to 90 min showing disproportionation. 4146

DOI: 10.1021/acs.molpharmaceut.6b00745 Mol. Pharmaceutics 2016, 13, 4141−4151

Article

Molecular Pharmaceutics

So in this case, the crystalline phase is incongruently melting at, during, and above its melting point. If the liquid phase cools and spontaneously crystallizes, the crystals are the free form of the drug if the counterion volatilized or may be a solution of molten drug and counterion if the counterion did not completely volatilize. Another possible scenario is disproportionation as was observed in the HCl salt. Spontaneous disproportionation begins with complete loss of the volatile HCl gas, resulting in collapse of the crystalline phase, forming a molten liquid of the free base. For such a salt, incongruent melting (disproportionation) and liquefaction occurs below the theoretical congruent melting point of the crystalline material. This can be practically determined as explained in Figure 5. The upper panel of the

Figure 4. Microscopic image (20×; Olympus BX51) and unit cell diagram of the HBr salt (a and c) and HCl salt (b and d) of CRH-1. The unit cells of the two salt forms viewed along their lattice c-axis demonstrates the similarity (isostructural).

protonated imidazo[1,2-b]pyridazin-3-yl nitrogen atom. In the HCl salt, the Cl to N(19) distance was 2.97 Å, having a hydrogen bond motif designation of D1,1(2)a, whereas for the larger bromide ion distance to the adjacent nitrogen atom was 3.15 Å [Br(1)···N(19)]. This is a discrete hydrogen-bonding pattern, where there is a single acceptor and a single donor atom, with each of the two atoms participating in the hydrogenbonding motif.20 Also this is considered to be the simplest form of hydrogen bonding that consists merely of a protonated drug molecule associated with anion, due to the scarcity of hydrogen bond donating groups. Isostructural pharmaceutical salts are known to exhibit substantial differences in physicochemical properties.21−23 It is therefore of interest to see if isostructural salts exhibit similar disproportionation behavior. Thermal Analysis of Neat API. In the early stages of development, thermal analysis is particularly useful in the characterization of crystalline salts. Thermally stable crystalline salts melt congruently with no change in composition across the melting temperature range, and its molten liquid retains the counterion. On cooling below the melting point, the salt recrystallizes. In other words, the salt is thermally stable both in its crystalline phase and as a molten liquid. This process is reversible and repeatable assuming that crystallization kinetics are favorable and crystallization occurs at a reasonable rate. The fragility of the supercooled melt and its proclivity to crystallize has been the subject of much interest in the literature. Congruently melting crystalline salts melt consistently at a single temperature with a constant heat of fusion. This type of behavior is well documented in the pharmaceutical literature, which is often assumed to be what occurs during conventional thermal analysis. We observed the congruent melting of CRH-1 free base at different heating rates (Figure S5). On the other hand, some salts may be stable in the solid state, but quite unstable when melted. If the counterion is volatile, it will begin to vaporize with the onset of melting and continue until the crystalline solid has melted completely. While this unstable liquid phase may not negatively impact the measurement of the melting point, it seriously affects the heat of fusion measurement because it is the sum of the heat of fusion and the heat of disproportionation. In addition, there will be contribution from the heat of vaporization of the counterion.

Figure 5. TGA/DSC demonstrating spontaneous volatile loss of the HCl counterion. In the upper half of the figure are the TGA curves while the lower half of the figure shows simultaneously collected DSC curves. The red curves are results for material heated in an open pan while black curves are for material in a hermetically sealed pan.

figure shows TGA weight loss profiles while the lower panel contains the simultaneously collected DSC heat flow signals. In an open pan, spontaneous disproportionation results with a weight loss equivalent to the monohydrochloride (8%) beginning at about 150 °C and continuing up to 240 °C. The simultaneously collected DSC curve shows an endothermic event in this same temperature range. When the HCl salt in a hermetically sealed pan was analyzed, volatile loss (endothermic event) was inhibited, that is, these thermal events did not occur until above 250 °C. The initial volatile loss of HCl fills the head space of the hermetic pan and prevents further reaction. As the temperature increases, eventually the pressure limit is exceeded with sudden loss of about 8% weight, attributed to complete loss of HCl, over a narrow temperature range. This experiment suggests that the back pressure in the hermetic pan inhibits disproportionation as has occurred in the open pan and perhaps allows the crystalline HCl salt to approach its “true” congruent melting point. When the HCl salt is subjected to unencumbered heating in an open pan, spontaneous disproportionation ensues and incongruent melting occurs at about 100 °C or well below the inhibited melting point in the hermetic pan. Volatile loss attributed to the disproportionation of HCl gas in Figure 5 was confirmed by simultaneous TGA/mass spectrometry (Figure S6). DSC heat flow curves of the HBr salt heated at different rates show differences in onset temperature or melting point depression 4147

DOI: 10.1021/acs.molpharmaceut.6b00745 Mol. Pharmaceutics 2016, 13, 4141−4151

Article

Molecular Pharmaceutics

Figure 6. DSC heat flow curves of (A) CRH-1 HBr salt and (B) CRH-1 HCl salt, heated in open pans at five different heating rates.

at slower heating rates (Figure 6A). This is due to incongruent melting (disproportionation) of the HBr salt. The HBr and HCl salts (Figure 6B) share this common feature. Model free kinetics (MFK) was applied to the conversion curves of the TGA weight loss profiles of HCl and HBr salts. An initial evaluation using MFK is a rather small investment in time to gain some important empirical information. The MFK module calculates apparent activation energy for a family of curves at each isoconversion point. MFK is especially useful in the evaluation of solid state reactions in materials that may undergo more than one type of process or when the mechanism of the reaction changes throughout the heating range. Materials can be analyzed more generally without an in depth knowledge of the underlying processes. For example, it would be plausible to hypothesize that disproportionation and volatile loss of either HCl or HBr from the parent molecule would be more facile from the surfaces of the drug particles and that the rate of loss might change as diffusion of the gases from inside particles lowers the rate. It might also be expected that the rate of volatile loss would be different initially from crystalline particles than from a collapsed liquid that could form throughout the complete process. In the bottom panel of Figure 7, the activation energy is plotted as a function of the extent of conversion. There was a gradual increase in the apparent activation energy through the first 25% disproportionation of the HCl salt, and then it leveled out. Between 10 and 90% conversion, the activation energy was 122 ± 23 kJ/mol. The weight loss profiles for HBr salt occurred at higher temperatures, indicating that it was more thermally stable than HCl salt. The activation energy for HBr was 210 ± 10 kJ/mol, roughly twice that for the HCl salt. MFK of the HCl salt by TGA weight loss profiles was corroborated using the simultaneous DSC heat flow curves. The DSC profiles produced an activation energy plot similar to that produced by the TGA data (Figure S7). It is both interesting and important that similar kinetic information was obtained for a volatile counterion to provide some confidence that MFK disproportionation kinetics could be obtained by DSC even for a nonvolatile counterion. TGA/sDSC of heminapadisylate salt showed that the melting endotherm and the weight loss occurred approximately over the same temperature range. This was independent of the heating rate (Figure S8). More importantly, the temperature at which the weight loss was initiated was lower than the

Figure 7. TGA weight loss profile (top panel; acquired at several heating rates) was used to create conversion curve (mid panel), showing % converted as a function of time. Activation energy plots (bottom panel), created using conversion curve, reveal that the activation energy for the HBr salt disproportionation is about twice that of the HCl salt.

temperature of the melting endotherm onset. These results strongly suggest that decomposition is initiated before the onset of melting. TGA weight loss profiles for the heminapadisylate material occurred at progressively lower temperatures as heating rates decreased. MFK applied to these time and temperature dependent weight loss profiles for heminapadisylate salt revealed the conversion rate was not constant across the entire range. It was higher at the outset and through about 50% of the weight loss and then decreased. The activation energy for heminapadisylate was 360 ± 44 kJ/mol, much higher than for HCl or the HBr salt (Figures S9 and S10). TGA/sDSC of the adipic acid cocrystal (Figure S11) indicates that endothermic melting is complete before the sample decomposes (weight loss in TGA). Therefore, there is no effect of heating rate on its melting transition. Based on thermal analyses, the adipic acid cocrystal, though with the lowest melting point, was the most thermally stable physical form. Counterion Volatility. The Litmus Paper Test. To confirm the volatility of HCl counterion, as observed during DSC, we conducted a simple litmus paper experiment. A saturated 4148

DOI: 10.1021/acs.molpharmaceut.6b00745 Mol. Pharmaceutics 2016, 13, 4141−4151

Article

Molecular Pharmaceutics solution of sodium chloride was filled in a small glass vial, and this open vial was placed in a bigger glass vial, containing 200 mg of powdered tablet of HCl or HBr salt formulation. The open small glass vial containing saturated salt solution was used to create an environment of ∼75% RH inside the capped large glass vial. A blue litmus paper was suspended in the glass vial, and the vial was capped and kept at 80 °C. After 24 h of storage, the blue litmus paper had turned red in the vial containing the HCl salt, due to volatilization of HCl gas. On the other hand, in the vial containing the HBr salt, no color change of the litmus paper was observed (Figure 8).

Figure 9. Free base (%, w/w) formed in tablets containing HBr salt (n = 3) stored at 40 °C and 22 or 75% RH. Inset shows early time points for onset of disproportionation.

It is interesting to note that, irrespective of the RH conditions, the disproportionation reached a plateau after few days of storage. Slurry pH measurements were done as a surrogate to microenvironmental pH measurement. At time “0”, the slurry pH of the formulation containing the HBr salt was found to be 2.28, which is higher than the slurry pH of the HBr salt alone, which was found to be 1.63. A drop in the slurry pH of the formulation started and continued until it reached a constant pH of ∼1.65 within 26 h (Figure S12). Dissolution of the dissociated HBr counterion from the salt within the mobile water may have brought down the microenvironmental pH and stopped further disproportionation reaction, resulting in the observed plateau in the disproportionation kinetics of HBr salt. The formulation containing the HCl salt, when subjected to accelerated stability testing at 40 °C and 11%, 32%, or 75% RH conditions, disproportionated to the free base. The disproportionation kinetics are depicted in Figure 10. In the case

Figure 8. (a) No change in litmus paper color in vial containing HBr salt, and (b) blue litmus paper changed to red in the vial containing the HCl salt.

Physical Stability in “Prototype” Formulations by Quantitative XRPD. The tablets containing prototype formulations were subjected to accelarated stability testing, in an open pan. As predicted by thermal analysis, formulation containing the adipic acid cocrystal did not show any evidence of free base formation, even after storage for 300 days at 40 °C/75% RH. This was a bit surprising, given the solution state data, but consistent with thermal data obtained under a dry environment without the presence of excipients. Hence, the cocrystal appears to be suitable as a crystalline API. However, there would be concern for its dissociation in the solution state and potential in vivo influence. In the case of the heminapadisylate salt, thermal analysis showed degradation even before completion of its melting, and predicted dispropotionation propensity. But with no free base formation, the formulation containing heminapadisylate salt was also found to be stable after storage for 300 days at 40 °C/75% RH. This may be due to the fact that the heminapadisylate salt may be unstable at very high temperatures, conditions not relevant to storage of pharmaceutical dosage forms. This physical stability data is also consistent with the solution stability data, which indicates that the heminapadisylate salt form is stable even when in equilibrium with water at a pH of 2. The formulation containing the HBr salt, when subjected to accelerated stability testing at 40 °C and 22% or 75% RH conditions, was unstable and disproportionated to the free base. The disproportionation kinetics are depicted in Figure 9. In tablets stored at 75% RH, the disproportionation rate was the highest during the first 2 days of storage. On the other hand, in tablets stored at 22% RH, disproportionation could be discerned after 2 days. This is likely due to the much lower water content in the sample stored at the lower RH. Following storage at 75 and 22% RH for 2 days, the water content in the tablets containing HBr salt was ∼10% and ∼0.5% w/w respectively (determined by the water sorption experiment). Disproportionation is reported to be a solution mediated transformation.16,24,25

Figure 10. Free base (%, w/w) formed in tablets containing HCl salt (n = 3) stored at 40 °C and 11, 32, or 75% RH. Inset shows early time points for onset of disproportionation.

of samples stored at 32 and 75% RH, disproportionation was observed within 4 h of storage. But in the case of samples stored at 11% RH, no disproportionation was observed within 24 h of storage. Table 2 contains the slope values of the trend lines of initial and later time points in disproportionation kinetics data. In contrast to the HBr salt, where a plateau was observed in the disproportionation, the HCl salt disproportionation 4149

DOI: 10.1021/acs.molpharmaceut.6b00745 Mol. Pharmaceutics 2016, 13, 4141−4151

Article

Molecular Pharmaceutics

was self-limiting and occurred only up to ∼9%. We hypothesize that the HBr formed following disproportionation and dissolved in the sorbed water, and its ionization caused a dramatic lowering of the microenvironmental pH and stopped further disproportionation reaction. Thus, the reaction byproduct makes disproportionation self-limiting. But in the case of the HCl salt, the disproportionation reaction did not exhibit a plateau. In this case, HCl left the “system” by volatilization. As a result, the disproportionation reaction continued. Interestingly, no disproportionation was observed in the case of the heminapadisylate salt and adipic acid cocrystal in solid state. However, the cocrystal dissociated readily in the dissolution medium during IDR analysis. If the salt to free base solubility ratio is very high, the pHmax is reduced to well below the pKa of the compound. When the salt form is taken to a pH where the free base is the thermodynamically stable form, disproportionation occurs spontaneously. This was observed in the case of the HCl and HBr salts. Therefore, the objective of salt selection should be to identify counterions which yield the desired rather than very high aqueous solubility. Very high salt solubility (with respect to the free base) can translate to an increase in physical instability and lead to salt → free base conversion.

Table 2. Slope of Disproportionation Kinetics Plots of Formulation, Containing HCl Salt slope RH (%)

initial time points

later time points

11 32 75

0.31 1.98 3.39

0.02 0.05 0.05

reaction does not plateau. We shall discuss this in detail later. The initial time points in the curve (Figure 10) indicate rapid disproportionation. A similar trend was also observed in the slurry pH measurement that was done as a surrogate to microenvironmental pH measurement (Figure S13). In later time points during stability study, HCl vaporization may be the driving force for the disproportionation (Figure 10) and may be the rate limiting step. Volatilization of HCl, following disproportionation of HCl salt, has been reported in the literature.26−30 Hydrochloride and methanesulfonate salts of triazole antifungal agent of the α-styrylcarbinol class (pKa = 1.7−2.2) were prepared, and stored at 60 °C. A loss of 12.7% of chloride ion was observed after 6 h of storage, in the case of its HCl salt. On the other hand, the methanesulfonate salt with 5% w/w water was stable after 3 weeks of storage.26 Other examples also indicate that the primary driving force for disproportionation reactions in HCl salts is liberation of HCl gas. The liberated hydrogen chloride gas either escapes from the system or reacts with other components of the formulation.27,28 A proprietary hydrochloride salt, REV-6000-A(SS), was found to liberate hydrogen chloride gas and thereby cause the rusting of tablet punches and dies. The thermogravimetric analysis of the salt indicated the loss of a mole of hydrogen chloride between 147 and 175 °C.29 Volatilization of HCl may also influence the chemical stability of an API. Lyophilization of quinapril HCl, in the absence of pH control, produced an amorphous mixture of salt and free base due to volatilization of a small amount of HCl. As amorphous free base is chemically less stable, the overall degradation rate increased relative to quinapril HCl salt alone.30 HCl and HBr are salts of mineral acids, and moisture associated with these salts can be very acidic. When compressed as a solid dosage form, these acidic salts share the available moisture with other basic formulation excipients. As compared to sulfonate or carboxylate salts, salts of these mineral acids can produce a highly unfavorable stability environment in dosage forms.31 Though HCl and HBr salts of CRH-1 are isostructural, their disproportionation kinetics were quite different. The released HBr remaining in the system caused a pronounced lowering of the microenvironmental pH below pHmax and stopped further disproportionation. But in the case of the HCl salt, the released HCl volatilized and escaped the system, keeping the microenvironmental pH well above its pHmax. HCl vaporization became the driving force for continuous (or larger degree) of disproportionation as HCl was quantitatively removed from the system.



CONCLUSION



ASSOCIATED CONTENT

1. Disproportionation propensity of salts (HCl, HBr, heminapadisylate) and adipic acid cocrystal of corticotropin releasing hormone receptor-1 antagonist was studied using MFK. The DSC heat flow curves of HCl, HBr, and heminapadisylate salts showed incongruent melting (disproportionation). Based on thermal analysis, though the adipic acid cocrystal had the lowest melting point, it was the most stable crystal form. 2. To ensure stability in solid state, tablets of “prototype” formulations containing salts and cocrystals were subjected to accelerated stability studies. As predicted by MFK, formulations containing the adipic acid cocrystal did not show disproportionation of API even after 300 days of storage. Similarly, formulations containing the HCl or HBr salts disproportionate as predicted by MFK. Though isostructural, the disproportionation kinetics of HBr and HCl salts were quite different, as the HCl counterion is volatile. 3. Formulation containing the napadisyalte salt also did not disproportionate after 300 days of storage, but this contradicted the findings of model free kinetics. Finally, solution state stability of salts and cocrystal was also assessed by determining their solubility at different pH values and intrinsic dissolution rate using 0.01 N HCl. In the present study, all the salts and cocrystal, except the heminapadisylate salt, disproportionated in dissolution medium.



S Supporting Information *

SIGNIFICANCE Disproportionation of salts of weakly basic drugs can be modulated by microenvironmental pH. While the HCl and HBr salts of CRH-1 are isostructural, they exhibit a profound difference in their disproportionation behavior. When stored at 40 °C/75% RH, in the case of the HBr salt, disproportionation

The Supporting Information is available free of charge on the ACS Publications website at DOI: 10.1021/acs.molpharmaceut.6b00745. IDR profiles and solubility, stability, SXRD, DSC, TGA/MS, TGA/sDSC, and pH studies (PDF) 4150

DOI: 10.1021/acs.molpharmaceut.6b00745 Mol. Pharmaceutics 2016, 13, 4141−4151

Article

Molecular Pharmaceutics



(15) Zhou, D.; Schmitt, E. A.; Zhang, G. G.; Law, D.; Vyazovkin, S.; Wight, C. A.; Grant, D. J. Crystallization kinetics of amorphous nifedipine studied by model-fitting and model-free approaches. J. Pharm. Sci. 2003, 92 (9), 1779−1792. (16) Merritt, J. M.; Viswanath, S. K.; Stephenson, G. A. Implementing quality by design in pharmaceutical salt selection: a modeling approach to understanding disproportionation. Pharm. Res. 2013, 30 (1), 203−217. (17) Hammersley, A. P. Fit2D: an introduction and overview; ESRF internal report, ESRF97HA02T; 1997. (18) Hammersley, A. P.; Svensson, S. O.; Hanfland, M.; Fitch, A. N.; Hausermann, D. Two-dimensional detector software: from real detector to idealised image or two-theta scan. High Pressure Res. 1996, 14 (4−6), 235−248. (19) Avdeef, A.; Comer, J. E.; Thomson, S. J. pH-Metric log P. 3. Glass electrode calibration in methanol-water, applied to pKa determination of water-insoluble substances. Anal. Chem. 1993, 65 (1), 42−49. (20) Motherwell, W. S.; Shields, G. P.; Allen, F. H. Graph-set and packing analysis of hydrogen-bonded networks in polyamide structures in the Cambridge Structural Database. Acta Crystallogr., Sect. B: Struct. Sci. 2000, 56 (5), 857−871. (21) Wood, P. A.; Oliveira, M. A.; Zink, A.; Hickey, M. B. Isostructurality in pharmaceutical salts: How often and how similar? CrystEngComm 2012, 14 (7), 2413−2421. (22) Ong, W.; Cheung, E. Y.; Schultz, K. A.; Smith, C.; Bourassa, J.; Hickey, M. B. Sodium and potassium salts of bumetanide trihydrate: Impact of counterion on structure, aqueous solubility and dehydration kinetics. CrystEngComm 2012, 14 (7), 2428−2434. (23) Galcera, J.; Molins, E. Effect of the counterion on the solubility of isostructural pharmaceutical lamotrigine salts. Cryst. Growth Des. 2009, 9 (1), 327−334. (24) Guerrieri, P.; Jarring, K.; Taylor, L. S. Impact of counterion on the chemical stability of crystalline salts of procaine. J. Pharm. Sci. 2010, 99 (9), 3719−3730. (25) Hsieh, Y. L.; Taylor, L. S. Salt stability-effect of particle size, relative humidity, temperature and composition on salt to free base conversion. Pharm. Res. 2015, 32 (2), 549−561. (26) Maurin, M. B.; Addicks, W. J.; Rowe, S. M.; Hogan, R. Physical chemical properties of alpha styryl carbinol antifungal agents. Pharm. Res. 1993, 10 (2), 309−312. (27) Elder, D. P.; Delaney, E. D.; Teasdale, A.; Eyley, S.; Reif, V. D.; Jacq, K.; Facchine, K. L.; Oestrich, R. S.; Sandra, P.; David, F. The utility of sulfonate salts in drug development. J. Pharm. Sci. 2010, 99 (7), 2948−2961. (28) Lin, S. L.; Lachman, L.; Swartz, C. J.; Huebner, C. F. Preformulation investigation I: relation of salt forms and biological activity of an experimental antihypertensive. J. Pharm. Sci. 1972, 61 (9), 1418−1422. (29) Narurkar, A. N.; Purkaystha, A. R.; Sheen, P. C. Effect of various factors on the corrosion and rusting of tooling material used for tablet manufacturing. Drug Dev. Ind. Pharm. 1985, 11 (8), 1487−1495. (30) Guo, Y.; Byrn, S. R.; Zografi, G. Effects of lyophilization on the physical characteristics and chemical stability of amorphous quinapril hydrochloride. Pharm. Res. 2000, 17 (8), 930−936. (31) Gould, P. L. Salt selection for basic drugs. Int. J. Pharm. 1986, 33 (1), 201−217.

AUTHOR INFORMATION

Corresponding Author

*E-mail: [email protected]. Present Addresses

§ N.K.T.: Amgen Inc., One Amgen Center Drive, Thousand Oaks, CA 91320, United States. ∥ R.J.B.: Eurofins Lancaster Laboratories PSS, 2425 New Holland Pike, Lancaster, PA 17601, United States.

Notes

The authors declare no competing financial interest.



ACKNOWLEDGMENTS N.K.T. was sponsored by the Lilly Innovation Fellowship Award. The work was partially supported by the William and Mildred Peters Endowment Fund. The authors thank Jeremy Merritt for the consultations and enlightening discussion. Use of the Advanced Photon Source was supported by the U.S. Department of Energy, Office of Science, Office of Basic Energy Science, under Contract No. DE-AC02-06CH11357. Parts of this work were carried out in the Characterization Facility, University of Minnesota, which receives partial support from NSF through the MRSEC program.



REFERENCES

(1) Lipinski, C. A.; Lombardo, F.; Dominy, B. W.; Feeney, P. J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Delivery Rev. 2012, 64, 4−17. (2) Thayer, A. M. Finding solutions. Chem. Eng. News 2010, 88 (22), 13−18. (3) Crew, M. A new year for solubility enhancement. Drug Dev. Delivery 2014, 14 (1), 22. (4) Serajuddin, A. T. Salt formation to improve drug solubility. Adv. Drug Delivery Rev. 2007, 59 (7), 603−616. (5) Engel, G. L.; Farid, N. A.; Faul, M. M.; Richardson, L. A.; Winneroski, L. L. Salt form selection and characterization of LY333531 mesylate monohydrate. Int. J. Pharm. 2000, 198 (2), 239−247. (6) Stephenson, G. A.; Aburub, A.; Woods, T. A. Physical stability of salts of weak bases in the solid-state. J. Pharm. Sci. 2011, 100 (5), 1607−1617. (7) Rohrs, B. R.; Thamann, T. J.; Gao, P.; Stelzer, D. J.; Bergren, M. S.; Chao, R. S. Tablet dissolution affected by a moisture mediated solid-state interaction between drug and disintegrant. Pharm. Res. 1999, 16 (12), 1850−1856. (8) Hsieh, Y. L.; Yu, W.; Xiang, Y.; Pan, W.; Waterman, K. C.; Shalaev, E. Y.; Shamblin, S. L.; Taylor, L. S. Impact of sertraline salt form on the oxidative stability in powder blends. Int. J. Pharm. 2014, 461 (1), 322−330. (9) Baertschi, S. W.; Alsante, K. M.; Reed, R. A. Pharmaceutical Stress Testing: predicting drug degradation; CRC Press: 2011. (10) Vyazovkin, S.; Wight, C. A. Kinetics in solids. Annu. Rev. Phys. Chem. 1997, 48 (1), 125−149. (11) Vyazovkin, S.; Wight, C. A. Model-free and model-fitting approaches to kinetic analysis of isothermal and nonisothermal data. Thermochim. Acta 1999, 340, 53−68. (12) Vyazovkin, S. Reply to “What is meant by the term ‘variable activation energy’when applied in the kinetics analyses of solid state decompositions (crystolysis reactions)? Thermochim. Acta 2003, 397 (1), 269−271. (13) Vyazovkin, S. Thermal analysis. Anal. Chem. 2006, 78 (12), 3875−3886. (14) Zhou, D.; Schmitt, E. A.; Zhang, G. G.; Law, D.; Wight, C. A.; Vyazovkin, S.; Grant, D. J. Model-free treatment of the dehydration kinetics of nedocromil sodium trihydrate. J. Pharm. Sci. 2003, 92 (7), 1367−1376. 4151

DOI: 10.1021/acs.molpharmaceut.6b00745 Mol. Pharmaceutics 2016, 13, 4141−4151