Gelation Mechanism of Erythromycin Ethylsuccinate During Crystallization

In this paper, the gelation mechanism of erythromycin ethylsuccinate (EES) during crystallization is investigated for the first time. The generated semisolid gel-like phase exhibited a 3D fibrillar network morphology and the typical rheological properties of gels. The fibers inside the gel-like phase were confirmed to be new types of EES solvates using powder X-ray diffraction, thermogravimetric analysis/differential scanning calorimetry, and gas chromatography. The gelation and crystallization regions in EES-1-propanol solid–liquid phase diagram were determined. Analyses of solvent parameters showed that moderate solvent polarity may promote EES gelation. Fourier transform infrared spectra, nuclear magnetic resonance spectra, and scanning electron microscopy analyses indicated that through intermolecular hydrogen bonds, EES and solvent molecules assembled into fibers via crystallographic mismatch branching growth. The fibers intertwined into a 3D network microstructure and formed a gel-like phase, completely immobilizing the solution.


Introduction
Solution crystallization is widely used in the production of pharmaceutical crystals. During the crystallization of active pharmaceutical ingredients (APIs) under various conditions, several studies have reported the occurrence of gelation phenomenon [1,2]. Gelation of some APIs, such as clopidogrel sulfate [3] and valsartan [4], is rather likely to form a gel-like phase in their supersaturated state, which is an unwanted process in industrial crystallization. However, gelation may have some advantages. New perspectives propose that gelation allows APIs to be used as molecular gel dosage forms and allows the control of crystal habits of APIs. Hence, the study of APIs gelation and the gelation mechanism is of great importance.
Gelation refers to a phenomenon where, under certain conditions, the solute molecules in a solution self-assemble into a 3D network structure that immobilizes the whole liquid phase [5]. At the macroscopic scale, the gelled solution loses its flowability and turns into a semisolid soft matter. Gelation is driven by the assembly of solute molecules through multiple reversible noncovalent interactions [6][7][8], including hydrogen bonds, π-π stacking, van der Waals forces, electrostatic, dipole-dipole, and hydrophobic interactions.
Similar to crystallization, gelation is also a common phase separation phenomenon. During crystallization, macroscopic phase separation occurs, leading to a solid-liquid phase equilibrium, whereas, during gelation, only a semisolid phase exists because of a microscopic phase separation. Both crystallization and gelation are driven by supersaturation. Methods that result in supersaturation or add external energy to induce the formation of new phases can result in gelation or crystallization [9][10][11][12][13][14], such as cooling a hot solution, pH switch, ultrasonic treatment, and optical radiation.
Dynamic arrest is regarded as the main mechanism of physical gelation [15], which means that a new phase induced by solute molecules is trapped in a kinetically favorable gel-like state rather than a thermodynamic stable state. The metastable state is the outcome of the competition between solute-solute interaction and solvent-solute interaction [16] and occurs in the metastable region in its phase diagram. For APIs, most studies focused on the metastable region of crystallization. Only a few reports involved the metastable equilibrium phase diagrams that contain the gel-like phase; meanwhile, gelation mechanism remains poorly understood.
In this paper, the gelation of a widely used erythromycin analog antibiotic, erythromycin ethylsuccinate (EES, CAS No. 1264-62-6, Fig. 1) is studied for the first time. First, the EES gel-like phase was characterized and the solvate nature of the fibers inside the phase was confirmed. Then, the EES gelation region was developed in the solvent system of 1-propanol. Finally, Fourier transform infrared (FTIR) spectra, nuclear magnetic resonance (NMR) spectra, and solvent parameters analyses were conducted to unravel the gelation mechanism at a molecular level. The results of this study provide useful guidelines to understand the gelation mechanism of other APIs.

Materials
A white crystalline powder of EES (C 43 H 75 NO 16 , MW 862.07) with a mass fraction purity of over 98% was purchased from Wuhan Yuancheng Gongchuang Science and Technology Ltd, China. The organic solvents (acetone, 1-propanol, etc., analytical reagent grade) were purchased from Tianjin Jiangtian Chemical Technology Co., Ltd, China, and were used without further purification. Deionized water was purchased from Tianjin QingYuanquan Purified Water Co., Ltd, China. Deuterated reagents (acetone-d6 and D 2 O, 99.9% deuterium content) were purchased from Tianjin Heowns Biochem Technology Co., Ltd, China.

Preparation of Gel-Like Phase
Method 1: Solventing-out. To achieve a gel-like phase, 1.5 g of EES powder was dissolved completely in 10 mL acetone and a transparent solution was formed at 20 °C. Under moderate mechanical stirring (300 r/min), 6 mL water was added as an anti-solvent at a constant flowrate (1.00 mL/min) using a peristaltic pump, and then the whole system turned into a gel-like state. Gelation also occurred in a 1-propanol-water mixture at different volume ratios.
Method 2: Cooling. 1.5 g EES powder was added into a 4 g water-acetone mixture (1 mol/mol) and heated to 51 °C to make the powder dissolve completely. Then, the hot solution was kept at room temperature to form a gel-like phase. Gelation also occurred in 1-propanol using a similar method.
Method 3: A gel-like phase was also formed when a moderate amount of EES powder was added to certain solvents. In a typical gelation process, 0.37 g EES powder was added into a 4 g water-acetone mixture (1 mol/mol) at 20 °C under 300 r/min stirring. After 5 min, the solution became turbid and then a gel-like phase emerged. Gelation also occurred in 1-propanol when the amount of EES powder was changed.

Characterization
Micrographs were acquired by scanning electron microscopy (Hitachi TM3000 tabletop SEM, Japan) and polarized optical microscopy (Olympus BX51, Japan). Powder X-ray Diffraction (PXRD) data were collected by a PXRD instrument (Rigaku D/MAX-2500, Japan) under the following conditions: 8°/min scan speed, 0.02° step size, 2°-50° scan range, and Cu Kα radiation (λ = 1.5406 Å). Rheological experiments were conducted on a rheometer (Anton Paar MCR 302, Austria) at 20 °C. Dynamic frequency sweep (DFS) was conducted in an angular frequency range of 100-0.1 rad/s with a strain amplitude of 0.05%, and dynamic strain sweep (DSS) was conducted in a strain range of 0.01-100% at 10 rad/s frequency. Thermogravimetric analysis (TGA) and differential scanning calorimetry (DSC) were conducted on a TGA instrument (Mettler Toledo TGA/DSC 1/SF, Switzerland) and DSC instrument (Mettler Toledo DSC 1/500, Switzerland), respectively. 5-10 mg of samples was loaded in a crucible and heated from 25 to 250 °C at 10 °C/min under N 2 atmosphere. FTIR spectra were determined on a spectrometer (Bruker Tensor 27, Switzerland) with a resolution of 1 cm −1 , and the solid powder was tableted with KBr and scanned in the range of 4000-400 cm −1 . NMR spectra were determined on an NMR spectrometer (Bruker Avance III 400 MHz, Switzerland) in a D 2 O and acetone-d6 mixture (1 mol/mol), which included 13 C-NMR spectra and concentration variation 1 H-NMR spectra, with a resolution of 0.45 Hz and sensitivity of 250:1.

Solubility Measurement
Dynamic Method 4 g of 1-propanol (designated as m 2 ) was poured into a 50-mL jacketed glass vessel, and a thermostat bath was used to regulate the experimental temperature. At a stirring rate of 300 r/min, a certain small amount of EES powder was continuously added to the solvent to approach saturation. The last added amount was not more than 0.01 g. The total amount of EES powder consumed when the slightly turbid solution did not turn transparent any more was recorded as m 1 . Each measuring point was repeated thrice, and the solubility of EES was calculated by Eq. (1): The temperature ranged from 10 to 50 °C, and 5 equidistant points were considered.
Static Method The cloudiness of the slightly turbid solution mentioned above increased with time, suggesting a phase transition after the EES dissolution. Thus, the static method was used to acquire the solubility data and the equilibrium solid phase was tested. The experiments were carried out in Erlenmeyer flasks, which were shaken in a thermostat shaking bath at 300 r/min. Excessive EES powder and 15 g 1-propanol were added into the Erlenmeyer flasks, and samples were taken every 2 h. After 8 h, a liquid-solid equilibrium was reached, and then the shaking was stopped for 0.5 h for the suspension to separate into two phases. A total sample of 5 mL was extracted with a syringe equipped with a 0.22-μm polyvinylidene fluoride filter. The sample was placed in an evaporating dish, and the weight of the dish (X 0 ) and the total weight of the dish and solution (X 2 ) were recorded. The liquid samples were evaporated in a vacuum oven at 50 °C and weighed every 12 h until a constant weight (X 1 ) was achieved. The PXRD pattern of the wet solid phase was recorded. The solubility of the equilibrium solid phase was calculated by Eq. (2): The samples were measured thrice for each Erlenmeyer flask, and the average value was used to calculate the solute solubility. The temperature ranged from 10 to 50 °C over 5 equidistant points.

Gelation Points Measurement
Cooling Method An appropriate amount of EES powder (m 1 ) and 6 g of 1-propanol (m 2 ) were added into a 50-mL jacketed glass vessel within a thermostat bath. Under a 300 r/min stirring rate, the vessel was heated at 1 °C/min to dissolve EES completely, and then, it was cooled at the same rate. The temperature was recorded when the solution became discernibly turbid. This was assigned as the gelation point of EES, after which the sample would completely turn into a gel-like phase. The amount of EES powder was then changed, and the experiments were repeated to obtain a series of gelation points at different concentrations. The concentration of the solution was calculated by Eq. (1).

Gelation Phenomenon and Characterization of the EES Gel-like Phase
The EES gel-like phase obtained through three methods showed a white opaque appearance and exhibited no gravitational flow under tube inversion test. When heated, the gel-like phase was separated into solid and liquid phases and eventually turned into a transparent solution. When the solution was cooled, the gel-like state reappeared, as shown in Fig. 2. The transition temperatures depended on the solution concentration and solvent composition. Long, thin, and flexible fibers were observed in the microscope images of the gel-like phase, and their sizes varied with the preparation methods ( Fig. S1-S3). The variation of the microscale morphology (Fig. 3) of the gel-like phase was tracked throughout the gelation process by a microscope. On the initial stage, the emergence of a few single fibers caused The mechanical properties of the gel-like phase were measured by dynamic rheology characterization. The DSS curves of the EES gel-like phase are shown in Fig. 4a. Both storage modulus (G′) and loss modulus (G″) remained constant with the increase in strain from 0.01 to 0.1%, which is the linear viscoelastic region (LVR) of the gel-like phase. The low strain limit (0.1%) demonstrated the crystalline nature of the fibers [17]. The DFS was performed within the LVR, as shown in Fig. 4b. In the whole frequency range, G′ was larger than G″, and both were almost independent of frequency, meaning that the sample was unable to flow and remained in a solid-like phase regardless of the bearing force. These rheological properties are similar to those of typical gel substances [18].

Analysis of PXRD Patterns
The morphology of the fibers was maintained when the EES gel-like phase was dried in air for 1 h, while the immobilized solution was completely dried out. The PXRD patterns of EES fibers are shown in Fig. 5, where the series of diffraction peaks in the PXRD patterns indicate the crystalline nature of the fibers. The fibers formed in different solvents showed the same PXRD patterns, indicating identical packing modes of the EES molecules in the crystal structure. The obtained Bragg lattice spacing values (d) of 1.99, 0.97, 0.63, 0.55 nm closely followed the ratio of 1:1/2:1/3:1/4, suggesting the existence of an 1D lamellar structure in the nanofibers of the gel-like phase with a layer distance of 1.99 nm [19,20].
Compared with the PXRD pattern of EES crystals, the PXRD patterns of the fibers showed two new peaks at 4.44° and 13.96°, and the peaks at 8.15°, 11.05°, 15.22°, 17.00°, and 20.68° disappeared, suggesting that the crystalline form of the fibers is different from that of the EES crystals and has not been reported. After drying the samples in a vacuum oven at 50 °C for 24 h, the PXRD patterns of the dried samples (Fig. S4) were different from those of the fibers but   similar to that of the EES crystals. The SEM image of a dried sample (Fig. S5) shows a similar morphology to that of the EES crystals. The changes indicate that the fibers transformed to EES crystals. Figure 6 shows the thermal analysis of EES fibers. The TGA curve of the fibers formed in 1-propanol showed a weight loss of 2.55 wt%, and the DSC curve had an endothermic peak at the corresponding temperature range, indicating that the sample was 1-propanol solvate. For the fibers formed in the water-acetone mixture, a weight loss of 0.90 wt% and an endothermic peak appeared in TGA and DSC curves in the same temperature range, respectively, proving that the sample was also a kind of EES solvate.

Thermal Analysis
The headspace sampling gas chromatography (GC) analysis verified the existence of acetone and 1-propanol in the fibers formed in water-acetone and 1-propanol, respectively (Figs. S6 and S7). The existence of water in the fibers formed in the water-acetone mixture was proved by volumetric Karl Fischer titration (see Supporting Information). The results obtained from the above two sections indicate that the fibers formed in different solvents were EES solvatomorphs, which were discovered for the first time. The gelation of the EES in its solution was due to the formation and cross-linking of the fibrous solvates.

Phase Diagram
The solubility data in 1-propanol were measured by both static and dynamic methods, corresponding to the blue and red curves in Fig. 7. The equilibrium solid phase in the static method was proved to be EES 1-propanol solvate (Fig. S8), demonstrating the solubility data (blue curve) corresponding to the solubility of EES 1-propanol solvate. At the same temperature, the solubility of the EES solvate is lower than that of the EES crystals. Thus, with the addition of the EES crystals into 1-propanol, the EES crystals transformed into a solvate when the concentration of the EES in solution exceeded the solubility of the solvate. In a series of typical experiments, solutions represented by points from L 1 to J 1 were prepared. It showed that because few fibrous crystalline solvates crystallized out under low supersaturation, the solution could not gel in the region between the blue and red curves, which can be regarded as the crystallization region. Gelation would inevitably occur when the concentration is above the red curve, which can be regarded as the gelation region. The lower the supersaturation in the gelation region, the more difficult and time-consuming it is for gelation to occur. Thus, it is indicated that a high degree of supersaturation promotes gelation. The gelation method 3 described in Experimental section is hence termed as "trans-crystal gelation," and the solubility of the EES crystals can be regarded as the critical gelation concentration. The black curve in Fig. 7 corresponds to the gelation points of EES measured by cooling. Compared with the trans-crystal gelation region, the solution gels fast in the cooling gelation region, and the critical concentration of cooling gelation is higher than that of trans-crystal gelation.

Analysis of Solvent Parameters
The influence of solvent properties on the EES gelation was investigated by solvent screening. The solutions of EES in acetone, methanol, ethanol, 1-propanol, 2-propanol, DMF, and NMP gelled entirely with added water, while the solution in acetonitrile gelled only partially. As shown in Table 1, the solvents that facilitated gelation have something in common: a dielectric constant range of 18.3-37.5, δ p (dipole interactions) of 6.1-18 MPa 1/2 in Hansen solubility parameters, and π * (polarizability parameter) of 0.48-0.88 in Kamlet-Taft solubility parameters [21]. Outside these ranges, it is impossible for gelation to occur in the corresponding solvents.
The matching degree of solvent-solute intermiscibility can be described by the residual of solubility Δδ. A smaller value of Δδ means a better solubility of a solute in a solvent. Δδ can be calculated from Eq. (3): where δ d and δ h represent the dispersion interactions and hydrogen bond interactions, respectively, while subscripts 1 and 2 correspond to EES and solvent, respectively. Table 1 ( shows the Δδ of the studied gelation solvents, and it ranges from 4.3 to 14 MPa 1/2 , which is smaller than those of other solvents, indicating that the gelation solvents have a better dissolving capacity for EES. Combined with the analysis of solubility parameters, it can be concluded that a high solubility of EES in the solvents is necessary for gelation. A low solubility results in the separation of a small amount of fibers from the solution, which makes it hard to trap quantities of solvent molecules and ultimately leads to crystallization.

FTIR
FTIR characterization was performed on the EES crystals and fibers. The peaks at 1741 and 1698 cm −1 in Curve a of Fig. 8 represent the C=O stretching vibration of ester groups and carbonyl of EES, respectively. For fibers, the former peak shifted to 1736 cm −1 , while the latter disappeared, which could have been masked by the former. It was speculated that ester groups participated in the intermolecular hydrogen bonds as H acceptors, accepting H from the solvent molecules to form the EES solvate fibers. Furthermore, the red shift of C-O-C vibration peaks in the ester groups and glucosidic bonds also imply the participation of these O atoms in the supramolecular self-assembly of EES solvate fibers as H acceptors (Table S1).

3
The wide peak of -OH stretching vibration at 3453 cm −1 of EES crystals indicates that the EES molecules interacted with each other through hydrogen bonds between hydroxyls. For fibers, the -OH stretching vibration peak split into two peaks, indicating the variation of intermolecular hydrogen bonds between the hydroxyl groups of the EES molecules after fibers formation. Some of the EES-EES intermolecular hydrogen bonds between the hydroxyls were substituted by the new EES-solvent intermolecular hydrogen bonds in the supramolecular fibers [23,24], as the small size makes the solvent molecules bond easily with the EES molecules. Moreover, both the out-of-plane bending vibration peak of -OH and the bending vibration peak of -OH on the cladinose ring exhibited a red shift in the fibers, demonstrating the participation of these hydroxyls in intermolecular hydrogen bonds.
Because of the presence of many intermolecular hydrogen bonds between the EES and solvent molecules, they assembled together to form 1D lamellar structural units, which assembled into 2D fibers quickly, and many branches were generated around the fibers, as shown in Fig. 9. Finally, the fibers intertwined and cross-linked with each other to construct a 3D self-assembled fibrillar network structure. Thus, the formation and growth of the EES-solvent supramolecular fibers is in accordance with the mechanism of crystallographic mismatch branching [25].

NMR
Concentration variation 1 H NMR spectra are shown in Fig. 10. With the increase in EES concentration, the system changed from a solution to gel-like state, which led to the weakening of peak intensity and peak splitting, suggesting that the mobility of molecules decreased, and protons could not produce individual sharp signals as they were in the aggregated state in the gel-like phase [26]. The H atoms of H 2 O (δ of ~ 4.30 ppm) show downfield chemical shifts with the increasing EES concentration, indicating the participation of H 2 O molecules in the EES-solvent hydrogen bonds as H donors.
The 13 C NMR spectra (see Fig. S10) were also recorded for the EES solution (1.875 mg/mL) and gel-like state (3.75 mg/mL). The chemical shift at nearly δ =211.5 ppm (Fig. 11a) corresponds to the C atoms of carbonyl groups. As the amount of solvents was much larger, we assigned the high peak to acetone-d6 and the low peak to the carbonyl group on the EES ring. The downfield shifts of these peaks indicate the participation of carbonyl groups as H acceptors in the EES-solvent intermolecular hydrogen bonds. In addition, the peaks of -COOC-in ester groups around 175 ppm (Fig. 11b) showed minute downfield shifts (0.06-0.1 ppm), suggesting the involvement of ester groups in the EES-solvent intermolecular hydrogen bonds.

Conclusions
Gelation phenomenon of EES in single or mixed solvents was reported for the first time in this paper. Analyses of PXRD, TGA/DSC, and GC indicate that the fibers formed during the gelation process were solvates of EES, which were discovered for the first time. The phase diagram of the gelation region and crystallization region was obtained for EES in 1-propanol. High supersaturation was found to be the reason for the formation of numerous fibrous solvate and gelation. Solvent parameters analysis indicated that solvent polarity was also a key factor for gelation. The results showed that gelation can be avoided by controlling the supersaturation and selecting appropriate solvents. FTIR and NMR analyses suggested that the dominant driving force for the formation of EES solvate fibers was the intermolecular hydrogen bonds between EES and solvent molecules. The crystalline fibers intertwined with each other to construct a 3D fibrillar network, immobilizing solvent molecules and leading to the formation of a gel-like phase.