Jovana Panića,
Aleksandar Tot*a,
Nenad Jankovićb,
Patrik Dridc,
Slobodan Gadžurića and
Milan Vraneša
aFaculty of Science, Department of Chemistry, Biochemistry and Environmental Protection, University of Novi Sad, Trg Dositeja Obradovića 3, 21000 Novi Sad, Serbia. E-mail: aleksandar.tot@dh.uns.ac.rs
bDepartment of Chemistry, Faculty of Science, University of Kragujevac, Radoja Domanovića 12, 34000 Kragujevac, Serbia
cFaculty of Sport and Physical Education, University of Novi Sad, Lovćenska 16, 21000 Novi Sad, Serbia
First published on 6th April 2020
The purpose of this paper was to examine the density, viscosity and electrical conductivity at different temperatures, as well as the thermal stability and structural properties of previously reported ionic liquids based on active pharmaceutical ingredients. Lidocaine-based ionic liquids, with ibuprofen and salicylate as counterion, were prepared first. Their structures were confirmed by infrared, mass and 1H and 13C nuclear magnetic resonance spectroscopy. The Newtonian behaviour of lidocaine ibuprofenate was confirmed from viscosity measurement results, while it was impossible to determine for lidocaine salicylate. The interactions and structures of the studied ionic liquids were analyzed based on the measured density values, viscosity, electrical conductivity, and calculated values of thermal expansion coefficients and activation energy of viscous flow. From a theoretical aspect, DFT and MD calculations were performed. The obtained descriptors and radial distribution, as well as structural functions, were used to understand the structural organization of the synthesized ionic liquids.
In this work, attention was dedicated to two API-ILs, lidocaine ibuprofenate ([Lid][Ibp]) and lidocaine salicylate ([Lid][Sal]). It is generally known that lidocaine has been using for local anaesthesia, while ibuprofen and salicylate are nonsteroidal anti-inflammatory drugs (NSAID). The combination of local anaesthetic and NSAID drugs has shown positive effects in the treatment of analgesia in animals.5 In human practice, such formulations would find the most extensive use in the local treatment of acute pain and inflammation, with the most effortless topical application. To achieve these dual drug function, necessary is to modify or improve some of the physicochemical properties that these drug substances have while in solid form. By combining these drugs into the ionic liquid, certain problems induced by their structure properties, like low skin permeability for some NSAIDs, which limit their transdermal use, are overcome. Recently, studies about novel ILs with the potential anti-inflammatory application reported that transformation to an ionic liquid leads to the improved lipophilicity/hydrophilicity of drugs, and thus better permeability through a cell membrane compared with solid salts, lidocaine and ibuprofen.6–9 In the literature, is still arguing if lidocaine ibuprofenate is an ionic liquid or deep eutectic salt. Although lidocaine ibuprofenate was first labelled as an ionic liquid and that term further expanded,8,9 in the work of Berton et al.7 it has been found as a deep eutectic liquid cocrystal with strong hydrogen bonds or interactions between partially ionized species. Lidocaine salicylate and acetylsalicylate have been prepared as dual-functional ILs also, but lidocaine acetylsalicylate showed limited stability due to conversion to salicylate.10
Altogether, both experimental and theoretical studies about structure, physicochemical properties, as well as studies on contacts between ILs and lipid membranes, remain rather popular. Although the number of new ionic liquids based on active pharmaceutical ingredients is rapidly rising, there is a lack of fundamental information about density, viscosity, electrical conductivity and thermal properties. Thus, this paper aims is to determine the physicochemical properties of earlier reported ionic liquids based on APIs, lidocaine ibuprofenate and lidocaine salicylate. Such study is significant when combined with complementary investigations, for developing and finding new active ingredients with improved pharmacokinetic characteristics.
The IR spectrum data of synthesized lidocaine ibuprofenate and lidocaine salicylate, as well as spectra of lidocaine, salicylic acids and ibuprofen, was used to estimate the level of the conversion to ionic liquid (Fig. 1 in manuscript and Fig. S4–S6 in ESI†). This approach was successfully applied in determination of starting compounds conversion to lidocaine etodolac ionic liquid.6 The wave numbers, v (cm−1) of the carbonyl groups in the carboxylic group of salicylic acid and ibuprofen, the carbonyl groups of amides in the lidocaine, the carbonyl groups of carboxylate anions (–COO−) and the amides in the lidocaine salicylate and lidocaine ibuprofenate were compared. Thereby, the presence of obtained products in form of ionic liquids was concluded. The measured spectra were presented in Fig. S4–S6 in ESI,† while characteristic bands were presented in Table 1.
Compound | Assignment | Wave number (cm−1) |
---|---|---|
Salicylic acid | –COOH | 1658 |
Ibuprofene | –COOH | 1720 |
Lidocaine salicylate | –COO− | 1589 |
Lidocaine ibuprofenate | –COO− | 1688 |
As can be seen from Table 1 and Fig. 1a, in lidocaine salicylate, the absorption bands characteristic to CO stretching of the carboxyl group of salicylic acid were diminished (1658 cm−1) and a new band appeared at about 1589 cm−1. This provided evidence that salicylic acid is completely deprotonated in lidocaine salicylate, proving completely conversion to ionic liquid.
On the other hand, the CO stretch in lidocaine ibuprofenate is located at around 1688 cm−1, and present combination of three peaks (Fig. 1b). These peaks are located at 1664 cm−1 (CO stretch of lidocaine amide group), 1720 cm−1 (CO of ibuprofen) and a new broad peak at 1690 cm−1 which indicate a dynamic equilibrium between ibuprofene, lidocaine and H-bonded lidocaine ibuprofenate.7 These observation suggest that ibuprofene is not completely in deprotonated state, indicating that conversion of lidocaine ibuprofenate in form of ionic liquid is not complete.
The thermal expansion coefficients for the pure ILs, αp, was calculated from obtained experimental densities:
The thermal expansion coefficient was determined in the same temperature range as measured densities. The obtained results are tabulated in Tables 2 and 3 and graphically shown as a variation with temperature, in Fig. 2 for d and Fig. 3 for αp.
T/K | da/g cm−3 | ηb/mPa s | kb × 104/mS cm−1 | λm × 104/S cm2 mol−1 | αp × 104/K−1 |
---|---|---|---|---|---|
a Standard uncertainty: u(d) = 3.0 × 10−3g cm−3, u(T) = 0.015 K.b Relative standard uncertainty: ur(η) = 0.02, ur(k) = 0.01. | |||||
293.15 | 1.02546 | 20712 | 1.44 | 0.62 | 8.28 |
298.15 | 1.02145 | 5606.9 | 2.51 | 1.08 | 8.32 |
303.15 | 1.01715 | 2806.0 | 4.75 | 2.06 | 8.35 |
308.15 | 1.01278 | 1467.6 | 7.48 | 3.26 | 8.39 |
313.15 | 1.00860 | 832.08 | 12.5 | 5.47 | 8.42 |
318.15 | 1.00436 | 468.49 | 18.1 | 7.95 | 8.46 |
323.15 | 1.00006 | 298.75 | 28.1 | 12.4 | 8.49 |
T/K | da/g cm−3 | ηb × 10−4/mPa s | kb × 104/mS cm−1 | λm × 104/S cm2 mol−1 | αp × 104/K−1 |
---|---|---|---|---|---|
a Standard uncertainty: u(d) = 3.0 × 10−3g cm−3, u(T) = 0.015 K.b Relative standard uncertainty: ur(η) = 0.02, ur(k) = 0.01. | |||||
298.15 | — | 5836.9 | — | — | — |
303.15 | — | 1483.1 | — | — | — |
308.15 | — | 290.52 | — | — | — |
313.15 | — | 89.149 | 0.40 | — | — |
318.15 | 1.12741 | 29.004 | 1.75 | 0.57 | 8.10 |
323.15 | 1.12311 | 10.226 | 5.62 | 1.84 | 8.13 |
328.15 | 1.11885 | 5.0573 | 12.6 | 4.16 | 8.17 |
333.15 | 1.11468 | 2.2890 | 29.0 | 9.57 | 8.20 |
338.15 | 1.11025 | 0.9496 | 61.8 | 20.5 | 8.23 |
343.15 | 1.10558 | 0.4308 | 116 | 38.6 | 8.26 |
348.15 | 1.10065 | 0.1809 | 209 | 70.0 | 8.30 |
353.15 | 1.09597 | 0.0970 | 349 | 118 | 8.34 |
Fig. 3 Variation of thermal expansion coefficients of ILs with temperature: ■ lidocaine ibuprofenate, and ○ lidocaine salicylate. |
By virtue of the fact that the density of ILs dependent upon directly on cation–anion interaction, it can be concluded that denser lidocaine salicylate has a better geometric fit of planar aromatic systems and stronger stabilization π–π interactions existence between lidocaine and salicylate than between lidocaine and ibuprofenate.13 Additionally, a higher thermal expansion coefficient of lidocaine ibuprofenate evidenced a faster movement of the cation from the anion, due to weaker interactions between them.
Variation of viscosity with temperature was fitted using the Vogel–Fulcher–Tammann (VFT) equation14–16 (Fig. 4):
Fig. 4 Variation of viscosity with temperature for pure ionic liquids: ■ lidocaine ibuprofenate, and ○ lidocaine salicylate. |
The determination of Newtonian or non-Newtonian behaviour was applied for lidocaine ibuprofenate, while for lidocaine salicylate it was impossible. Namely, for determination of Newtonian behaviour the measurements with a wide range of RPM are necessary. Unfortunately, due to the high viscosity of lidocaine salicylate and the limitation of equipment at higher RPM measurements with LV-4 spindle were impossible. The viscosity of fluids with Newtonian behaviour at a given temperature remains constant regardless of which Viscometer model, spindle or speed was used for the measurements. A linear dependency of shear stress from shear rate presented in Fig. 5 (Table S4†), together with the constant viscosity with shear rate values indicate Newtonian behaviour of lidocaine ibuprofenate.
Fig. 5 (a) Flow curve, showing shear stress as a function of shear rate for pure lidocaine ibuprofenate; (b) changes of viscosity values at 298.15 K with shear rate along with error bars. |
Fig. 6 Variation of electrical conductivity of pure ionic liquids with temperature: ■ lidocaine ibuprofenate, and ○ lidocaine salicylate. |
To predict ionicity for lidocaine salicylate and lidocaine ibuprofenate, Walden plot was applied.
The relation between molar conductivity and viscosity can be demonstrated by the equation:
logλm = logC + αlogη−1 |
Fig. 7 Walden plot for lidocaine ibuprofenate, ■, in temperature range 293.15–323.15 K, and lidocaine salicylate, ○, in temperature range 318.15–353.15 K. |
In the comparison of ideal KCl line, it can be seen that lidocaine salicylate and lidocaine ibuprofenate lie below this line. Further, the lidocaine salicylate is closer to KCl line, indicating better ionicity compared to lidocaine ibuprofenate. Typically, data for ionic liquids fall below this ideal line because of the strong interaction between ions in the pure ionic medium, which reduce the movement of most ions that are partially correlated with their closest neighbours. The traditional Walden plot shows three distinct groups of ionic liquids that are distinguished based on vertical distance to the KCl line, denoting as ΔW.17 The first group lies quite close to the KCl line (0.1 < ΔW < 0.5) suggesting that the ionic liquid is made up of almost independently mobile ions and can be classified as good ionic liquids.17–19 The second groups lies lower on the plot (0.5 < ΔW < 1.0) and is named “poor” ionic liquids. In this group existence of H-bonds and other specific interaction between ions in the pure state is more pronounced, leading to a significant reduction of ions mobility. The third group represents the liquids that are at least an order of magnitude below the ideal line, also described as liquid ion pairs or subionic liquids (ΔW > 1.0).17–19 In these liquids, ion conductivity is substantially less than would be expected based on their viscosity. The simplest example is the ion pair, which is neutral and therefore, does not contribute at all to the measured conductivity.
From Table 3 it can be seen that lidocaine salicylate, although being more viscous (290 Pa s) than lidocaine ibuprofenate (0.5 Pa s) nonetheless has similar conductivity, presumably because the degree of ion association is lower for lidocaine salicylate. This assumption is confirmed using Walden plot, where lidocaine salicylate is closer to idel KCl line (ΔW = 0.71) and can be classified as “poor” ionic liquid. On the other side, lidocaine ibuprofenate lies in region of subionic liquids (ΔW = 1.34) proving the predominance of ion pairs in its structure. This result is in accordance with conclusions driven from FTIR measurements, proving that conversion of lidocaine ibuprofenate in classical ionic liquid is not fully accomplished.
Fig. 9 Representation of optimized structures along with visualization of non-covalent interactions: (a) lidocaine ibuprofenate, (b) lidocaine salicylate. |
In addition, from Fig. 9 a potential intramolecular interaction between the keto group of lidocaine and the methyl group from the ortho position can be seen, which occurs as a result of an increase in the acidity of hydrogen atoms from the methyl group due to the expressed positive inductive and hyperconjugation effect, indicating the potential cation–cation interactions.
The next step was an analysis of frontier molecular orbitals (HOMO and LUMO) as a convenient tool for visualization of the global stability and reactivity of investigated compounds. The atomic contribution to HOMO–LUMO orbitals were generated and presented in Fig. 10 and 11 (positive electron density represents red colour and negative is in blue colour). From mentioned figures, in both ionic liquids, the HOMO and LUMO orbitals are mostly located above anions. The HOMO and LUMO orbitals shows which part of molecule had potential to interact with other structures, while the HOMO–LUMO difference represents the descriptor of this molecule stability. The HOMO and LUMO energy gap, ΔEgap (Table 4), quantify the possible charge transfer interactions within the molecule. The smaller value of band ΔEgap is connected to increased reactivity of the molecule. As can be seen in Table 4, the ΔEgap is lower for lidocaine salicylate. Based on this observation lidocaine ibuprofenate is more stable IL implying lower chemical reactivity because it hinders electrons to high-lying LUMO or extracting electrons from a low-lying HOMO inhibiting the formation of an activated complex for any further reaction.
Ionic liquid | ΔEgap/eV |
---|---|
Lidocaine ibuprofenate | 5.76 |
Lidocaine salicylate | 4.95 |
Further, hydrophobic/hydrophilic surfaces for both ionic liquids were calculated and presented in Fig. 12. Moreover, the percentage of hydrophilic surfaces, as well as lipophilic descriptor AlogP, were calculated and shown in Table 5.
Fig. 12 Representation of hydrophilic (turquoise colour) and hydrophobic (orange colour) for investigated ionic liquids: (a) lidocaine ibuprofenate, (b) lidocaine salicylate. |
Ionic liquid | AlogP | Percentage of hydrophilic surfaces (%) |
---|---|---|
Lidocaine ibuprofenate | 3.97 | 36.96 |
Lidocaine salicylate | 1.55 | 46.83 |
As can be seen from Table 5, both descriptors indicate that lidocaine salicylate is more hydrophilic than lidocaine ibuprofenate.
The molecular dynamic simulations were applied in order to obtain detailed information about global structure and interactions of pure ionic liquids. The using of MD simulations offers the opportunity to predict physical properties directly from a molecular interaction model and also to visualize the organization of molecules and their structure at a molecular level. This method was refined to present the position of a given ion corresponding to the position of the selected counterion, describing the local structure around ions using the structural analysis tools like radial distribution functions (RDFs) and structural factor calculations (SDFs).
The information about the distance between cations and anions were obtained from MD simulations using radial distribution function (RDFs). RDFs for anion (oxygens from COO−) and hydrogens from N1 and N2 in all ionic liquids are presented in Fig. 13. From Fig. 13, it can be seen higher values of g(r) at a smaller distance for lidocaine salicylate. This behaviour indicates that salicylate, due to a smaller size, can approach closer to lidocainium cation, than more bulk ibuprofenate. This result, indicate better, denser packaging of smaller salicylate between lidocainium cations, causing higher values of density and viscosity of lidocaine salicylate comparing to lidocaine ibuprofenate. The lidocaine ibuprofenate with branched, more hydrophobic tail prevents packaging of the ions and decrease the efficiency of π–π interactions more than salicylate. These results are in accordance with lower Tg values and higher conductivity of lidocaine ibuprofenate.
Further, comparing results for hydrogen attached to N1 and N2 atoms, it can be seen that hydrogen connected to N2 is closer to both anions. On the other hand, this makes N1 atom of lidocaine more available for additional interaction. With this purpose, the RDFs between N1 and O1 atom from lidocaine cation were calculated and shown in Fig. 14, to present the possibility of self-aggregation of lidocaine cation. From Fig. 14, the significant peak at 2.8 Å is noted, suggesting the self-aggregation tendency of lidocaine due to possible formation of H-bond between N1 and O1 atoms of lidocaine.
To separate anionic and cationic contributions to RDFs, structure functions were calculated for the representative centre of cation's, and anion's center of mass and results are presented in Fig. 15. As can be seen from Fig. 15, the same contribution of cation–cation and cation–anion to interaction in pure ionic liquid were obtained. The only differences are in anion–anion interactions, where the negative contribution to structure factor is noted in lidocaine salicylate, which is in contrary to lidocaine ibuprofenate. All obtained theoretical results indicate that ibuprofenate have a tendency to interact with itself, creating more cluster structures comparing to salicylates where this type of interactions are absent. The existence of anion segregation in lidocaine ibuprofenate originates from the partially deprotonated carboxylic group, which diminish electrostatic repulsion between negative charged anions, allowing the formation of H-bond between protonated and deprotonated –COO− group of ibuprofen. As the last evidence of stronger interactions between cation and anion in lidocaine salicylate, the ion pair binding energy (ΔGbin) was calculated. This is important parameter that can be correlated with physico-chemical properties of ILs such as melting point, thermal stability and transport properties.20 The significantly more negative energy was obtained for lidocaine salicylate (Ebin = −116.21 kJ mol−1) than for lidocaine ibuprofenate (Ebin = −71.44 kJ mol−1). From work of Bernard et al. and from obtained experimental results in this work, the dependence between low value of ΔGbin and higher values of viscosity and density were obtained.21 This confirms that ion pair binding energy is correlated with the strength of the intermolecular interaction between ions. Lower value of binding energy points to stronger molecular packaging, suggesting higher value of density and viscosity, as in the case of investigated ionic liquids.
Fig. 15 Partitioning of S(q) to cation–cation (red line), cation–anion (black line) and anion–anion (blue line) correlations: (a) lidocaine ibuprofenate, (b) lidocaine salicylate. |
Footnote |
† Electronic supplementary information (ESI) available. See DOI: 10.1039/c9ra08815f |
This journal is © The Royal Society of Chemistry 2020 |