Bayan Alkhawajaa,
Faisal Al-Akayleh*a,
Jehad Nasereddinb,
Muhammad Kamranc,
Tim Woodmand,
Zaid Al-Rubayea,
Nidal Qinnaa,
Mayyas Al-Remawia and
Ali R. Olaimate
aFaculty of Pharmacy and Medical Sciences, University of Petra, Amman, Jordan. E-mail: falakayleh@uop.edu.jo
bDepartment of Pharmaceutical Sciences, Faculty of Pharmacy, Zarqa University, Zarqa, 13110, Jordan
cDepartment of Chemistry, University of Bath, Claverton Down, BA2 7AY, Bath, UK
dDepartment of Life Sciences, University of Bath, Claverton Down, BA2 7AY, Bath, UK
eDrug Directorate, Jordan Food and Drug Administration, Shafa Badran, Amman, 11181, Jordan
First published on 7th May 2024
Therapeutic deep eutectic solvents (THEDSs) are the best exemplification of green alternative formulations of active pharmaceutical ingredients (APIs) that offer superlative properties of APIs. Previously, THEDESs of risperidone, fentanyl and levofloxacin with capric acid (CA) were developed by our group. These APIs share cyclic tertiary amine nuclei. Herein, DESs of two drugs bearing cyclic tertiary amine nucleus, namely, droperidol and aripiprazole, in the presence of CA, were investigated as model drugs. Comprehensive analyses were conducted using liquid-state 1D and 2D NMR and differential scanning calorimetry (DSC) to elucidate the regiochemistry and thermodynamic mechanisms bringing about those THEDESs. Everted gut sac technique was used to study the flux of the developed THEDESs. 1D and 2D NMR techniques analyses revealed the importance of cyclic tertiary amine nuclei in forming interactions with CA. This was confirmed by the downfield shift of the protons proximal to the tertiary amine groups compared to the individual drugs. Diffusion NMR analysis (DOSY) showed a significant reduction in the diffusion coefficient of CA in the mixed system compared with CA in isolation. Thermal analysis of the two drugs revealed that the drugs have a low tendency to recrystallise upon melting but rather vitrify from a melt to form an amorphous solid. Interestingly, the superior absorption and flux of the THEDES formulation of droperidol was demonstrated using the ERIS. Collectively, this work provides a green method to attain liquid formulations of APIs with enhanced pharmacokinetic features.
Owing to their versatility, biocompatibility, and sustainability, DESs are increasingly involved in a multitude of pharmaceutical applications.11,12 Ranging from being solvents for active pharmaceutical ingredients (APIs) to being utilised in drug delivery systems, DESs offered widespread advantages in the pharmaceutical realm. To this end, many recent reports displayed the utility of DESs in enhancing drug solubility of poorly soluble drugs,10,13–15 permeability,9 stability,16 drug delivery, as well as being involved in the synthesis and extraction of active pharmaceutical ingredients.17
More interestingly, DESs of APIs, where the API acts as a component of the DES or IL, or what is known as therapeutic DESs (THEDESs), is another direct application of enhancing the pharmaceutical characteristics in a safer and greener approach.18 Forming the so-called hydrophobic DES, fatty acid-based DESs have been applied in various potential non-aqueous administrations.18–21 Previously, DESs based on fatty acids have been formed with a multitude of APIs, such as menthol,22 fentanyl,23 propranolol,24 risperidone18 and levofloxacin.25 The most significant advantage of DES based on biocompatible fatty acids is the realisation of these properties under mild and solventless conditions.26
The common consensus within the literature is that hydrogen bond interactions are the main interactions bringing about the formation of DES.22 In-depth evaluation of the literature in the field of DESs and THEDESs has highlighted the limited understanding of the structural aspects of such systems, including predicting the propensity to form hydrogen bond interactions based on the chemical structure of APIs, thermodynamic behaviour or physical characteristics related to the lattice properties of the crystals of DES formers.
In our endeavours to generate deep understanding of DESs, we present in this work two novel THEDESs. Studying their chemical structures, some THEDEs with CA, such as fentanyl, risperidone and levofloxacin, share a cyclic tertiary amine moiety (Table 1).18,25 Therefore, we envisaged evaluating the interactions between APIs possessing similar nuclei and CA, primarily using 1D and 2D NMR and DSC techniques. The selection criteria of the drugs in this work, namely droperidol (Drp) and aripiprazole (Arp), was based on their chemical structure; hence, we anticipated similar interactions. Building upon our previous findings, the primary aim of this work is to underline the regiochemistry and thermodynamic interactions taking place in two novel THEDESs with CA. Furthermore, we sought to study the direct application of this system through the evaluation of the absorption and permeation of the THEDSs of Drp using the everted rat intestinal sac (ERIS) technique.
API | Structure | pKa | Solubilitya | Tm (°C) | Medical uses | Ref. |
---|---|---|---|---|---|---|
a Aqueous solubility at 25 °C, all records were obtained from PubChem.27b Solubility at 20 °C. Tm (melting point). | ||||||
Risperidone | 8.76 | 2.16 | 171 (ref. 28) | Treatment of schizophrenia, bipolar mania and psychosis | 29 | |
Fentanyl | 8.99 (ref. 30) | 0.2 | 83–84 | Opioid analgesic for pain management | 30 and 31 | |
Levofloxacin | Basic pKa 8.12 and acidic pKa 6.1 | Sparingly soluble | 225–227 | Broad-spectrum antibacterial drug | 25, 32 and 33 | |
Droperidol (Drp) | 7.46 | 4.21 | 145–146.5 | Management of the nausea and vomiting in surgical and diagnostic procedures | 34 and 35 | |
Aripiprazole (Arp) | 7.6 | 0.045 | 137–140 | Treatment of many forms of mood and psychotic disorders | 36 | |
Capric acid (CA) | 4.9 | 0.062 | 31.5 | — | 37 | |
Methyl caprate (MC) | — | 0.004b | −18.0 | — | 38 |
The plausible interaction behind the formed DES is the hydrogen bond interaction between CA and Drp, as seen previously with DES between CA and menthol.43 To investigate this further, Drp and CA binary mixtures were prepared by increasing the molar ratio of Drp to CA by physical mixing. Amongst the prepared mixtures, D1 (0.9:0.1), D2 (0.8:0.2) and D3 (0.7:0.3) transformed into liquid form, whereas the rest of the binary mixtures transformed into a pasty form (Table 3). Initially, the binary mixtures were investigated by 1D and 2D NMR experiments.
It is well established that hydrogen bonding interactions could cause a downfield shift (higher frequency) of the carboxylic acid protons due to the deshielding effect.44 However, our system couldn't assign the carboxylic acid protons due to their exchangeability (Fig. S1†). Hence, we focused on the change of the chemical shifts of Drp protons. To this end, comparing the chemical shifts between the different mixtures of CA and Drp, most of the protons of Drp exhibited no change in their chemical shifts with the increasing ratio of CA (Fig. S1†).
Nevertheless, a significant downfield shift was observed for the protons located at 3.35 ppm (1) and 2.90 ppm (2), assigned to the protons adjacent to the piperidine ring (Fig. 1A). The positive Δδ (chemical shift deviation) was calculated for the protons at 3.35 ppm (1) and 2.90 ppm (2) of Drp with CA and found to be 0.08 ppm and 0.09 ppm, respectively. These shifts indicate that these protons are either deshielded or the electron density is altered due to the interactions between the piperidine amine and the carboxylic acid of the CA. The observed downfield shift was also reported previously, where the proximal protons of the cyclic tertiary amine of levofloxacin.25
The tertiary amine group of the piperidine ring is expected to interact with CA by forming a hydrogen bond. Therefore, the proximal protons to the piperidine amine are anticipated to be in a new chemical environment due to this newly formed interaction. It is known that the formation of hydrogen bonding with the amine group leads to the increase in the acidity of the proximal protons, and hence, this could explain the observed downfield shift of the protons (1) and (2) (Fig. 1A).45 Their chemical shifts reflect the chemical environment of the protons, the more the electron-withdrawing effect surrounding the protons (more acidic protons), the higher the resonance frequency (higher chemical shift) will be observed. The presence of the neighbouring hydrogen has caused the downfield shift of the proximal protons, as outlined above. In contrast, distal protons to the piperidine ring, such as proton at 3.08 ppm (3), did not exhibit such change in the chemical shift (Fig. 1A).
To validate the observed 1H NMR shifts in the eutectic mixtures and to draw more comprehensive conclusions around the underlying interactions between Drp and CA, we set out to compare the interaction between the eutectic mixture of both CA and Drp with the binary mixture of methyl ester caprate (MC) and Drp. MC lacks the CA hydroxy group and hence is incapable of being a hydrogen bond donor with Drp. To this end, MC and Drp binary mixtures were prepared at equimolar ratio, and their 1H NMR spectra were compared. The positive Δδ (chemical shift deviation) of the protons at 3.35 ppm (1) and 2.90 ppm (2) of Drp with CA (0.08 ppm and 0.09 ppm, respectively) as observed earlier were not observed in the case of a binary mixture of Drp and MC (Fig. 1B). The attained results confirm the essential role of the carboxylic acid group of CA in forming the eutectic mixture with Drp.
Next, we set out to study the interactions using 1H Diffusion NMR spectroscopy experiments of the binary mixtures and individual components in the CDCl3. Diffusion NMR experiments, sometimes called DOSY (Diffusion Ordered Spectroscopy), provide a means for studying the interactions between the compounds by providing a 2D spectrum in which one axis is the chemical shift while the other is the diffusion coefficient.46 The advanced data analysis module of Mestrenova was used to create two-dimensional spectra with NMR chemical shifts (X-axis) and the diffusion coefficients (Y-axis) and to calculate the diffusion coefficients (F). As shown in Fig. 2C and D, the binary mixture exhibited a related diffusion rate compared to the separate compounds (Fig. 2).
Fig. 2 DOSY 1H spectra of (A) CA, (B) Drp, (C) binary 1:1 mixture of Drp and CA and (D) binary 1:1 mixture of Drp and CA with expansion. The solvent was CDCl3. |
Previously, diffusion experiments were performed to study the strength of hydrogen bonding, where the reduction of the diffusion coefficient was directly correlated with the H-bond strength.47 Herein, following the same methodology, the change in diffusion coefficient was employed to study the strength of interactions present in our system. The change in diffusion coefficients of the binary mixture of Drp and CA was evaluated relative to the individual components in the same solvent (Table 2). A decrease in the diffusion coefficient of CA revealed the strong interactions within the eutectic mixture of the Drp and CA. Generally, diffusion NMR experiments were previously used to study the strength of intermolecular interactions.48,49 However, to the best of our knowledge, this is the first study that demonstrates the effectiveness of the diffusion NMR technique for evaluating the strength of intermolecular interactions in THEDESs.
Compound | FCA | FDrp | FMC | FDrp |
---|---|---|---|---|
(×10−6) | (×10−6) | (×10−5) | (×10−6) | |
Individual component | 9.03 | 5.72 | 1.15 | 5.72 |
Binary mixture (1:1) | 7.82 | 6.79 | 1.06 | 5.33 |
To further validate the attained results, a diffusion experiment for the equimolar mixture of Drp and MC was also conducted (Fig. 3). The diffusion coefficients were calculated for the mixture of MC and Drp (Table 2). No significant difference was observed between the diffusion coefficients of the mixture and the individual components, proposing the absence of any interactions. Altogether, the diffusion experiment results of the MC and Drp mixture signify the presence of strong interactions in the eutectic mixture of Drp and CA, as demonstrated previously.
Binary mixtures of Arp and CA were prepared with an increasing molar ratio of Arp to CA. Amongst the prepared mixtures, A1 (0.9:0.1) and A2 (0.8:0.2) were transformed into liquid form, whereas the rest of the binary mixtures were transformed into a pasty form (Table 3).
Eutectic binary mixtures | Ratios (mol) | Eutectic binary mixture form |
---|---|---|
CA:Drp | 1:0 (CA) | — |
0.9:0.1 (D1) | Liquid | |
0.8:0.2 (D2) | Liquid | |
0.7:0.3 (D3) | Liquid | |
0.6:0.4 (D4) | Paste | |
0.5:0.5 (D5) | Paste | |
0.4:0.6 (D6) | Paste | |
0.3:0.7 (D7) | Paste | |
0.2:0.8 (D8) | Paste | |
0:1 (Drp) | — | |
CA:Arp | 1:0 (CA) | — |
0.9:0.1 (A1) | Liquid | |
0.8:0.2 (A2) | Liquid | |
0.7:0.3 (A3) | Paste | |
0.6:0.4 (A4) | Paste | |
0.5:0.5 (A5) | Paste | |
0.3:0.7 (A7) | Paste | |
0.2:0.8 (A8) | Paste | |
0:1 (Arp) | — |
Interestingly, a 1H NMR study of the binary mixtures in CDCl3 revealed the interaction of the amide group with CA, mainly through hydrogen bonding. This was confirmed by the observed significant downfield chemical shift of the amide proton's peak at 8.18 ppm (1).
Similar to the Drp, the role of the piperazine ring in forming eutectic solvents with CA was evaluated. Interactions of the amine groups in the piperazine ring were predicted via the chemical shifts of the adjacent CH2 protons at 3.20 ppm (2,3). A significant downfield shift of these two –CH2 groups (2,3) was evident with increased concentration of CA in the mixture. As the piperazine ring possesses 2 amines, the lack of downfield shift of the (–CH2) protons at 2.91 ppm (4) confirms that the amine next to the aromatic ring of the piperazine is mainly involved in the CA interactions (Fig. 4A). The positive chemical shift deviation (Δδ), as calculated from the protons at 8.18 ppm (1) and 3.20 ppm (2,3) of the Arp with CA, were found to be 0.62 ppm and 0.1 ppm, respectively. The observed downfield shifts imply that these protons are deshielded due to the interactions with the carboxylic acid group of the CA.
Next, validation of the intermolecular interactions using MC, which is not capable of forming acid–base interactions with Arp, was explored. However, MC might be able to form hydrogen bonding with an amide peak at 8.18 ppm. The downfield chemical shifts observed in the Arp and CA eutectic mixture for protons (1–3) were not observed in the case of binary mixtures of MC and Arp (Fig. 4B). This observation demonstrates the essential role of the carboxylic acid protons in mediating the interactions with Arp.
A diffusion experiment confirmed the intermolecular interactions of Arp with CA, as the diffusion spectra of the equimolar eutectic mixture exhibited related diffusion behaviour (Fig. S2 A and B†). On the other hand, the binary mixture of Arp with MC (1:1) was similar to the individual components, confirming unrelated diffusion (Fig. S2C†).
Fig. 5 DSC of raw materials. (A) 1st heating cycle. (B) Cooling cycle. (C) 2nd heating cycle. (D) Observed glass transitions of Drp and Arp and Gef. |
Fig. 5B shows the DSC thermograms of the cooling cycle of the raw materials. CA was observed to recrystallise completely at 24 °C (162.12 J g−1). Arp and Drp were all observed to vitrify in the cooling cycle, with glass transitions (Tg) being visible at 25.2 °C (Arp) and 23.7 °C (Drp) (Fig. 5D).
In the second heating cycle (Fig. 5C), the melting endotherm of CA was seen at 30.61 °C. Drp and Arp were observed to recrystallise. Drp exhibited a recrystallisation exotherm at 117 °C, and the Tm of Drp was observed again at 148 °C (90.78 J g−1). Arp was observed to recrystallise at 85.3 °C, with the Tm of Arp visible at 138 °C (84.77 J g−1).
Notably, neither of the drugs recrystallised but rather exhibited vitrification during the cooling cycle. The heating/cooling rate employed during the DSC analysis was 10 °C min−1; this indicates that the drugs belong to Glass Forming Ability (GFA) Class II drugs. Drugs belonging to GFA Class I and II are notable for having weaker crystals with a higher tendency to vitrify and a general tendency to yield comparatively more stable amorphous phases than GFA Class III drugs.52
Fig. 6 DSC of Drp–CA eutectic mixtures. (A) 1st heating cycle. (B) Cooling cycle. (C) 2nd heating cycle. |
In the cooling cycle of samples D4, D5, D6, and D8 (all of which were observed to be heterogenous at room temperature), the recrystallisation of CA was not seen. However, baseline changes that likely correspond to the Tg of Drp were seen in all samples. In contrast to samples D1 and D2, the second heating cycle thermograms for samples D4, D5, D6, and D8 showed the recrystallisation and melting peaks in addition to the Tg of Drp. The recrystallisation peak was seen at 96.7 °C (D4), 102.9 °C (D5), 98.8 °C (D6), and 96.7 °C (D8). The melting endotherms were seen at 117.4 °C (23.40 J g−1) for sample D5, 127.2 °C (34.91 J g−1) for sample D5, 124.7 °C (25.37 J g−1) for sample D6, and 129.0 °C (58.71 J g−1) for sample D8, which most likely corresponds to the depressed melting of Drp. Events corresponding to CA were not observed in the samples, most likely due to the loadings being below the limits of detection of DSC.
Fig. 7 shows the DSC thermograms of the Arp–CA eutectic mixtures A3, A5, and A7. In the first heating cycle (Fig. 7A), the desolvation endotherm of form H1 of Arp was visible in all three thermograms at 58.6 °C, 69.8 °C, and 74.9 °C for samples A3, A5, and A7, respectively. Furthermore, an endothermic event was visible in the thermogram of sample A3 at 18.7 °C (21.96 J g−1), which likely represents the Tm of CA. From the enthalpy values obtained from this endotherm and the Tm of the pure CA sample, we estimate the excess precipitated CA to be approximately 13.5% of the CA content of sample A3. In the cooling cycle (Fig. 7B), no events were observed in the thermogram of sample A3. Meanwhile, the Tg of Arp was clearly visible in the thermograms of samples A5 and A7. In the second heating cycle, no events were observable in the thermogram of sample A3. In the thermograms of samples A5 and A7, the Tg of Arp was again visible. Furthermore, the peritectic melting of form H1 of Arp is again visible in the thermograms of A5 and A7, followed by the phase transformation into the metastable form III polymorph, whose melting is seen at 139 °C.51,53,54
Fig. 7 DSC of Arp–CA eutectic mixtures. (A) 1st heating cycle. (B) Cooling cycle. (C) 2nd heating cycle. |
In order to validate the hypothesis formulated by the results obtained from both DSC and NMR, a control sample using MC instead of CA was also studied. Drp was used as the model drug for this sample as Drp exhibited the widest concentration range along which a eutectic liquid could be obtained. Fig. S3† shows the DSC results of methyl caprate/Drp sample, along with the thermogram of pure. The Tm of pure methyl caprate was seen at −13 °C (90.25 J g−1). In the thermogram of the mixture (1st heating cycle), both the melting endotherm of methyl caprate and that of Drp were clearly visible at −13 °C (65.84 J g−1) and 140 °C (165.50 J g−1), respectively. In the 2nd heating cycle, the Tm of MC was not visible. However, both the recrystallisation exotherm and Tm of Drp were visible at 46.42 °C and 138.9 °C. The enthalpy of melting of Drp in the second cycle was 100.65 J g−1, corresponding to 61% of the initial amount observed in the first heating cycle, which agrees with the behaviour of pure Drp observed in the thermogram of the pure sample (Fig. 5). The remainder of Drp that was not seen in the thermogram of the second heating cycle is likely converted to the amorphous form of the drug as was seen in the pure sample. The results suggest that there is no interaction between MC and Drp, as the drug melted and recrystallised as independently as a pure sample of the drug was observed to behave, despite the presence of MC in the sample. These findings complement the earlier observations during NMR analyses of the MC–Drp mixture, validating the role of carboxylic acid groups in bringing about interactions necessary to form the DES.
As mentioned earlier, a notable observation was that both Drp and Arp are both GFA Class II drugs,52 as evidence by their failure to recrystallise on cooling at the cooling rate of 10 °C min−1 used throughout the DSC experiments. Similar observations have been previously reported,18,25,55 indicating some correlation between the ability of a drug to form a THEDES and its GFA class. More specifically, since the GFA class is a proxy for crystal strength,52 there seems to be a correlation between the crystal strength of a drug and its ability to form THEDES. The control experiment conducted using MC and Drp suggests that the THEDES formation may be brought about by the hydrogen bonding, as the substitution of the hydrogen bonding forming carboxyl group in CA with the methyl group in MC prevented the formation of a DES. The aforementioned data suggests that the underlying mechanism of DES formation is brought about by the co-former (CA) forming hydrogen bonds with the drug, and due to the weak crystals of the drug (as evidenced by its GFA), the co-former is able to significantly depress its melting,7 resulting in the formation n of a DES.
ERIS is a technique widely employed in pharmaceutical research to study the systemic availability, absorption, permeation, and interactions of drugs. By using this tool, the permeation and absorption of drugs across the intestinal epithelium of rats are studied by measuring their flux.56 DESs generally have better absorption and permeation due to their ability to improve the poor water solubility and low dissolution rate of APIs.57,58
Forming hydrophobic THEDESs with CA would increase the hydrophobic features of the drugs. Nevertheless, the mixture would exhibit improved permeability and absorption. This was clearly demonstrated in the present study (Fig. 9). All the tested samples exhibited time-dependent intestinal flux that gradually decreased over time. Notably, D1 achieved the highest flux of 1.182 mg cm−2 s−1 at 15 minutes, remarkably exceeding the other formulations. While Drp and D2 showed similar profiles, D1 consistently maintained the highest flux throughout the 3 h. This superior permeation profile observed with D1 was statistically confirmed (p-value >0.05). The results demonstrated superior drug absorption when using and maintaining green processing conditions without incorporating organic solvents.
13C NMR (126 MHz, DMSO-d6) δ 198.98, 166.31, 164.31, 153.37, 134.25, 134.23, 131.36, 131.28, 130.75, 130.33, 128.91, 123.56, 121.68, 121.06, 116.16, 115.99, 109.37, 108.76, 56.89, 51.74, 49.67, 36.16, 27.35, 22.11.
13C NMR (126 MHz, CDCl3) δ 171.62, 158.50, 138.17, 134.10, 128.71, 127.60, 127.54, 125.09, 118.82, 115.88, 108.64, 102.20, 67.64, 57.92, 53.00, 31.10, 27.03, 24.60, 22.63.
13C NMR (126 MHz, DMSO) δ 174.96, 34.12, 31.74, 29.34, 29.22, 29.12, 29.01, 24.96, 22.56, 14.40.
Footnote |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4ra01469c |
This journal is © The Royal Society of Chemistry 2024 |