Introduction

Pharmaceutical researchers have recently being devoting more and more attention to co-crystals. This interest in co-crystals is connected with the possibility of their use as a new solid form of active pharmaceutical ingredients (APIs) [1,2,3,4]. Co-crystals are defined as a connection of at least two molecules in a solid state within the same crystal lattice. In relation to pharmaceutical co-crystals, one or even two components are drugs, as opposed to the majority of cases where the second component consists of a neutral molecule or “co-former”, which should have a harmless effect on the human body [5, 6]. Non-ionisable or ionisable (only weakly acidic or basic) compounds are permitted as part of the co-crystallisation process [5]. The co-crystal components are bonded together via non-covalent bonds, predominantly via hydrogen bonding, which explains why no proton transfer occurs between drug and co-former [7]. As a result, the US Food and Drug Administration (FDA) published guidance in 2013 which stated that co-crystals are not a new entity as such but should instead be treated as a “drug product intermediate” [8]. Attempts at co-crystal formation are performed with drugs within the second or fourth group in the Biopharmaceutical Classification System (BSC), with drugs characterised by poor solubility in water, a quality which, together with dissolution rate, is the key parameter as regards the bioavailability of drugs [6]. Many studies have shown that co-crystals improve the solubility in water and dissolution rates of parent compounds [9,10,11], and moreover, co-crystals can affect other physico-chemical properties of the same drug. In the literature can be found a description of an increase in physical and chemical stability or improvement in mechanical properties [12,13,14].

Screening and preparation of co-crystals can be achieved either by solid-state or solution-based methods [4, 12, 14,15,16]. The former consists of grinding methods without (neat grinding approach) or with the addition of a small amount of solvent (solvent-assisted grinding approach) [15]. In these methods, two substances are ground together with mortar and pestle or by mechanical grinders. Solvent evaporation, slurry conversion and cooling crystallisation are examples of traditional solution-based methods [12, 14]. The recent literature provides various modifications of solution-based methods such as approaches using supercritical fluid technology [17] or an ultrasound-assisted solution process [13]. As solid-state and solution-based methods have some drawbacks [15], attempts at co-crystal formation are usually conducted using both methods together in order to obtain co-crystals with the most suitable properties, at least in the opinion of the researchers.

Regardless of the above methods, recent years have witnessed the emergence of reports about the application of differential scanning calorimetry (DSC) as a simple and rapid approach towards co-crystal screening [18,19,20]. This approach is based on the heating, by a DSC device, of binary physical mixtures of drug and co-former. Detailed interpretation of the DSC scans obtained may serve as a rapid screening approach that allows co-crystal detection. This is based on the assumption that the melting point of co-crystal differs from that of output components, since in more than 50% of cases co-crystals showed a melting point lower than that of both drug and co-former [21]. Generally, when the physical mixture capable of co-crystal formation is heated, the following thermal effects should occur on the DSC scan: an exothermic peak associated with co-crystal formation directly after an endothermic peak, together with another endothermic effect recorded at a higher temperature than the first peaks and connected with co-crystal melting, or several endothermic peaks in temperatures different to the melting points of the pure components and associated with metastable eutectic melting, eutectic melting and co-crystal melting [18, 19, 22]. Furthermore, the possibility of co-crystal screening using the DSC method has also been described, where co-crystals could be prepared at different molar ratios [23]. In this case, a DSC scan of the physical mixture shows two exothermic and three endothermic peaks. In contrast, when a physical mixture incapable of co-crystal formation is heated, only a single endothermic peak associated with eutectic melting appears on the DSC scan.

Taking all the above into consideration, the aim of this research was to verify to what extent the DSC method can be used as a rapid screening tool for detecting the co-crystallisation process in binary physical mixtures of benzodiazepines and co-formers. Benzodiazepines are psychoactive drugs with a wide range of medical applications, e.g. in treating anxiety, epileptic seizures or problems with insomnia, but are also applied as sedative drugs in premedication prior to medical or dental procedures [24]. Benzodiazepines, however, are classified as drugs very slightly or practically insoluble in water and with a low dissolution rate, a fact which has immediate effects on their poor bioavailability. As this may adversely affect efficiency or increase the occurrence of side effects in pharmacotherapy [25], many methods including physical or chemical modification of a drug have been applied to enhance the water solubility of benzodiazepines, such as particle size reduction, solubilisation by surfactants or the formation of salts and complexes [26, 27]. Because of the recent surge in interest in co-crystal formation as a promising additional way to enhance the water solubility of drugs, the preparation of co-crystals of benzodiazepines with selected co-formers, such as low molecular organic acids and other chemically neutral compounds, has become even more crucial. The co-formers chosen for this study are readily water soluble, non-toxic and with no negative effects on the human body. Nevertheless, to the best of the author’s knowledge, there have been no reports on co-crystals of benzodiazepines.

Experimental

Chemicals

Benzodiazepines—chlordiazepoxide, clonazepam, diazepam, estazolam, lorazepam, lormetazepam, oxazepam and temazepam—were supplied courtesy of Polfa Tarchomin (Warsaw, Poland). The co-formers used in the study were acquired from Sigma-Aldrich (St. Louis, USA)—p-aminobenzoic acid, caffeine, fumaric acid, glutaric acid, nicotinamide and succinic acid; Acros Organics (New Jersey, USA)—citric acid and saccharin; and POCH (Gliwice, Poland)—urea. Both benzodiazepines and co-formers were characterised by a purity level above 98% and were used without further purification.

Sample preparation

Binary physical mixtures of benzodiazepines and co-formers at 1:1 molar ratios were prepared in dust extraction hood by the gentle mixing of both ingredients in an agate mortar using plastic spatula. To achieve complete homogenisation, benzodiazepines and co-formers were thoroughly mixed with spatula over 8 min. The samples obtained were then analysed by the DSC method.

To prepare the physical mixture, fixed amount of benzodiazepine equal to 0.05 mol was mixed with 0.05 mol of co-former that amount in milligrams varied depending on molar mass of co-former used. Both ingredients were weighed using analytical balance model XA 105 Dual Range (Mettler Toledo, Schwerzenbach, Switzerland).

Differential scanning calorimetry

A heat-flux DSC instrument model 822e (Mettler Toledo, Schwerzenbach, Switzerland) was used to obtain scans of benzodiazepines, co-formers and their binary mixtures. The device was connected with STARe software. Samples of 4.00–4.05 mg were placed into flat-bottomed aluminium pans, which were then sealed with perforated lids. The samples were measured over a range from 25 to 300 °C, except for mixtures with fumaric acid, for which measurements were taken from 25 to 350 °C. All measurements were taken under nitrogen flow (flux rate of 70 mL min−1) at a heating rate 10 °C min−1. DCS instrument was calibrated using indium and zinc as standards with regard to temperature and enthalpy.

Results and discussion

In this research, a total set of 72 binary physical mixtures composed of eight benzodiazepines (diazepam, temazepam, oxazepam, lormetazepam, lorazepam, clonazepam, estazolam and chlordiazepoxide) and nine co-formers (succinic, glutaric, fumaric, citric and p-aminobenzoic acids, nicotinamide, saccharin, urea and caffeine) mixed at 1:1 molar ratios were investigated using the DSC method. This method allows the more subtle differences in heat flow between the sample under study and the reference to be measured; thus, it may also be used to illustrate the physical changes (phase transitions) that occur in the mixtures of benzodiazepines with co-formers when heated under non-isothermal conditions. Melting (endothermic process) and crystallisation (exothermic process) are the most crucial phase transitions, as knowledge of these transitions is fundamental for the detection of potential co-crystals. For this reason, each change in the DSC scan of physical mixture in comparison with those of benzodiazepine and co-former alone can be interpreted as proof of co-crystal formation. Consequently, the DSC method can serve as a reliable screening method for co-crystallisation.

Benzodiazepines and co-formers

The core structure of benzodiazepines, the 1,4-benzodiazepine ring, is shown in Fig. 1, while the substituents attached to this structure are given in Table 1. The similar chemical structure of benzodiazepines is an effect of the connection of benzene and diazepine rings [24, 25]. As benzodiazepines have different side groups attached to the core structure, the drugs used in this study were classified as 5-phenyl-1,4-benzodiazepin-2-one derivatives including 1-methyl-substituted (diazepam), 3-hydroxy-substituted (temazepam, oxazepam, lormetazepam, lorazepam) and 7-nitro-substituted (clonazepam); derivative includes the addition of fused triazolo ring at the 1,2-position (estazolam); and N4-oxide 2-methylamine derivative (chlordiazepoxide). It is worth noting that benzodiazepines also have functional groups, i.e. hydroxyl (–OH) and carbonyl (C=O), which are capable of forming intermolecular hydrogen bonds with co-formers.

Fig. 1
figure 1

1,4-Benzodiazepine ring—core structure of benzodiazepines. R1, R2, R3, R4, R5 and R6 denote different side groups attached to the core structure

Table 1 Substituents R1–R6 at the core structure of benzodiazepines—1,4-benzodiazepine ring

Inspection of the physico-chemical data given in Table 2 reveals that molecular masses of benzodiazepines are located in a very narrow range between 284.74 (diazepam) and 335.18 (lormetazepam). The drugs in question are white or pale yellow crystalline powders, very slightly soluble in water (from 1000 to 10,000 parts of solvent required for one part of solute) or practically insoluble (> 10,000 parts of solvent required for one part of solute) [28]. This merely emphasises the importance of seeking new physical or chemical ways of modifying the water solubility of benzodiazepines, for example by co-crystallisation, to improve their bioavailability. Furthermore, analysis of melting points reveals that excluding diazepam, temazepam and lorazepam, the remaining benzodiazepines melt at high temperatures in excess of 200 °C. This is in turn reflected in the high, narrow and sharp-ended DSC endothermic peaks, followed by exothermic peaks due to the decomposition of melted drugs. No exothermic peaks were observed in the case of lorazepam or estazolam. The DSC data obtained are consistent with those found in the literature [29,30,31,32], with, however, insignificant differences due to experimental conditions, such as heating rate or gas atmosphere.

Table 2 General characteristics of benzodiazepines including temperatures of peaks (T peak) and enthalpies of transitions (∆H) acquired from DSC scans

Nine co-formers were used as the second component of potential benzodiazepine co-crystals—succinic, glutaric, fumaric, citric and p-aminobenzoic acids, nicotinamide, saccharin, urea and caffeine. Their structural formulas are illustrated in Fig. 2 and their primary physico-chemical properties in Table 3. The first fundamental condition in choosing organic compounds as co-formers is that they fulfil the status of generally recognized as safe (GRAS) [33], confirmation that the co-formers are not toxic and have no adverse effects on the human body. The second condition in the choice of co-former is the existence of functional groups that in their chemical structure are able to create so-called supramolecular synthons with molecules of psychoactive drugs which leads to the formation of intermolecular hydrogen bonds [33]. The carboxylic (–COOH), amino (–NH2), carbonyl (C=O) and hydroxyl (–OH) groups are the most desirable in this respect (Fig. 2).

Fig. 2
figure 2

Structural formulas of co-formers: a succinic acid, b glutaric acid, c fumaric acid, d citric acid, e p-aminobenzoic acid, f nicotinamide, g saccharin, h urea, i caffeine

Table 3 General characteristics of co-formers including temperatures of peaks (T peak) and enthalpies of transitions (∆H) acquired from DSC scans

Table 3 shows that the co-formers used in this study are white or colourless crystalline powders, very soluble in water, over the range of 3.3 g L−1 (saccharin) to 1080 g L−1 (urea). Furthermore, they are characterised by a wide range of melting points, from 76.66 °C (glutaric acid) to 295.81 °C (fumaric acid). In the latter case, the DSC measurements for mixtures including this acid were taken up to 350 °C, whereas remaining ingredients and mixtures were heated over a range 25–300 °C. The melting of glutaric and fumaric acids, and caffeine is preceded by polymorphic transitions, which is confirmed by the small endothermic DSC peak. The melting of citric acid and urea is, in turn, followed by decomposition above 210 °C, whereas melted p-aminobenzoic acid evaporates above 240 °C. These data are consistent with those found in the literature [19, 34,35,36,37,38,39,40,41].

Rapid co-crystal screening

The results of the DSC investigation of 72 binary physical mixtures of benzodiazepines with co-formers are given in Tables 47 and presented graphically in Figs. 38. Detailed inspection of these data shows that the mixtures under study could be divided into four groups consisting of 5, 10, 48 and 9 mixtures, respectively, depending on the number of endothermic and exothermic effects detected on the DSC scans. Special attention has been paid to the appearance of new peaks due to the melting or crystallisation of co-crystals, or the disappearance of peaks due to melting of benzodiazepines and co-formers. This provides one hint to suggest the co-crystallisation of drug with co-former. Moreover, the shift in DSC peaks to lower or higher temperatures and the change in their intensities were also taken into consideration.

Table 4 Results of DSC analysis of the first group of benzodiazepines physical mixtures with co-formers
Fig. 3
figure 3

DSC scans of: (a) lorazepam, (b) nicotinamide, (c) their physical mixture

The first group of mixtures

As shown in Table 4, the first group of samples under study includes oxazepam mixtures with fumaric and citric acids, and lorazepam mixtures with fumaric acid, nicotinamide and saccharin. This group is characterised by the appearance of a new exothermic effect on the DSC scan between two endothermic peaks. In the case of three mixtures—oxazepam with fumaric acid, and lorazepam with nicotinamide and saccharin—the exothermic peak occurs immediately after the first endothermic DSC effect, as illustrated in Fig. 3. The endothermic peaks appear at 129.94 and 180.00 °C, and the new exothermic peak at 133.67 °C directly following the first endothermic effect (Table 4). The occurrence of a small exothermic peak immediately after the endothermic could suggest the possibility of co-crystal formation under mixture heating. Nevertheless, endothermic peak temperatures are close to the melting points of the individual ingredients (Tables 2, 3), which creates uncertainty over interpretation. However, it does not exclude the possibility of co-crystal formation characterised by peak temperature only slightly different to that of the melting point of lorazepam.

In the case of oxazepam mixture with fumaric acid and lorazepam with saccharin, peak temperatures are significantly different to those of the melting points of ingredients. This could suggest interaction between components and the formation of a new structure under heating, a situation presented in Fig. 4. The first endothermic effect appears at 165.35 °C and is shifted to a temperature (about ~ 45 °C, Table 4) lower than that of the melting point of oxazepam. The second endothermic peak appears at 255.50 °C, between the melting points of oxazepam and fumaric acid.

Fig. 4
figure 4

DSC scans of: (a) oxazepam, (b) fumaric acid, (c) their physical mixture

In the case of oxazepam mixture with citric acid and lorazepam with fumaric acid, a new exothermic peak appears on the DSC scan at a temperature about 20 °C higher in comparison with the first endothermic effect, which, in turn, is shifted to a lower temperature. This is untypical for mixtures capable of co-crystal formation. However, endothermic peaks (Table 4) do not overlap with the melting points of either ingredient (Tables 2, 3), but are found at significantly lower temperatures than those of the melting points of ingredients. This could suggest interaction between components, possibly co-crystallisation.

The second group of mixtures

Of 72 binary mixtures, ten (diazepam with succinic and citric acids, temazepam with glutaric acid and urea, oxazepam with glutaric and p-aminobenzoic acids and urea, lorazepam with glutaric and citric acids, and clonazepam with urea) form the second group (Table 5), distinguished by the occurrence at least three endothermic peaks on the DSC scan. These mixtures could be further divided into four subgroups, which are mixtures including glutaric acid, urea, diazepam and a final subgroup consisting of oxazepam mixture with p-aminobenzoic acid and lorazepam mixture with citric acid.

Table 5 Results of DSC analysis of the second group of benzodiazepines physical mixtures with co-formers

On the DSC scans of glutaric acid mixtures, the first peak (endothermic) appears at approximately 74 °C and could be associated with a transition from β-polymorph (the stable form at room temperature) to α-polymorph. The second endothermic effect occurs over the range of 90–99 °C and could be assigned to the melting of the α-polymorph [35]. This is illustrated in Fig. 5, which presents the DSC scan of lorazepam mixture with glutaric acid as an example. Insignificant differences in temperatures of the second peak in subsequent mixtures could be due to the influence of benzodiazepines on the melting of glutaric acid. The heats of fusion reflected by the second peak are similar for all glutaric acid mixtures (Table 5); however, the values are more than three times lower than the heat of fusion of glutaric acid (Table 3). The third peak occurs at 108.11, 123.03 and 127.24 °C for lorazepam, oxazepam and temazepam mixtures, respectively. This endothermic peak does not overlap with the peak solely due to the melting of benzodiazepine alone (Table 2). Rather it may suggest interaction between components, although additional and detailed investigation of this interaction is necessary. In all binary mixtures of glutaric acid, exothermic effects appear above 155 °C, probably due to sample decomposition.

Fig. 5
figure 5

DSC scans of: (a) lorazepam, (b) glutaric acid, (c) their physical mixture

In the case of mixtures with urea, the first endothermic effect is found at approximately 135 °C and could be associated with the melting of the co-former (Table 3). Figure 6 illustrates this through the example of a DSC scan of temazepam mixture with urea. There is an overlap with the first endothermic peak at 133.88 °C due to the melting of urea; however, the heat of fusion reflected by this peak is more than eight times lower than that for co-former alone. At 144.28 °C, the second endothermic peak appears, but the temperature differs from that for the melting of temazepam (Table 2). The last broad endothermic peak found on the DSC scan above 180 °C could not be associated with the melting point of either component, suggesting that interaction occurs under heating.

Fig. 6
figure 6

DSC scans of: (a) temazepam, (b) urea, (c) their physical mixture

The third subgroup contains diazepam mixtures with succinic and citric acids. The DSC scans of these mixtures show three endothermic effects, reflected in Fig. 7 which presents the DSC scan of diazepam mixture with succinic acid. The first and second peaks are found at temperatures approximately 10 °C lower than that of the melting points of drug and both co-formers. The third broad endothermic peak occurs above 180 °C. The situation described could imply interaction between components under heating or could be attributed to the influence of ingredients on their melting points. Again, additional study is required to clarify this issue.

Fig. 7
figure 7

DSC scans of: (a) diazepam, (b) succinic acid, (c) their physical mixture

The last subgroup consists of mixtures of oxazepam with p-aminobenzoic acid and lorazepam with citric acid. The first endothermic peak occurs at a temperature significantly lower than that for the melting of the constituent parts. The second endothermic peak, as with the first, is found at a lower temperature (about 6 °C) than that for the melting points of benzodiazepine and co-former, and the last is detected at a temperature higher than that of melting points of oxazepam and lorazepam. In the case of oxazepam mixture with p-aminobenzoic acid, the DSC scan additionally revealed an exothermic peak immediately after the third endothermic peak due to sample decomposition.

The third group of mixtures

The third group of samples under study includes 48 mixtures given in Table 6. Their DSC scans illustrate those events in which no formation of co-crystals was confirmed. In those cases, one or two endothermic peaks are found. Furthermore, in the case of sparse mixtures, an additional exothermic peak was observed at higher temperatures due to sample decomposition. The situation in which single endothermic peak occurs on the DSC scan is shown in Fig. 8. The peak temperature at 152.48 °C for lormetazepam mixture with p-aminobenzoic acid is lower than that for melting of drug and co-former alone. A single endothermic effect is evidence that both components form an eutectic mixture, which precludes the possibility of obtaining co-crystals [18]. Above 280 °C, a small exothermic effect appeared due to decomposition. On the other hand, two endothermic peaks are found in the case of estazolam mixture with p-aminobenzoic. A sharp endothermic peak at 132.13 °C on the DSC scan converting into a broad endothermic event. The temperatures of both peaks (Table 6) are lower than those for melting of estazolam and p-aminobenzoic acid (Tables 2, 3). This also suggests the formation of an eutectic mixture, which again implies the impossibility of co-crystal formation between drug and co-former.

Table 6 Results of DSC analysis of the third group of benzodiazepines physical mixtures with co-formers
Fig. 8
figure 8

DSC scans of: (a) lormetazepam, (b) p-aminobenzoic acid, (c) their physical mixture

The fourth group of mixtures

As shown in Table 7, the nine mixtures consisting of all chlordiazepoxide mixtures are gathered into the fourth group. This drug melts at 243.51 °C and then decomposes, which is reflected by two subsequent DSC peaks: endothermic and exothermic (Table 2). A detailed analysis of the DSC scans of chlordiazepoxide mixtures with co-formers reveals that the predominant exothermic peaks due to decomposition hinder the retrieval of reliable data. For this reason, the DSC scans of mixtures with chlordiazepoxide cannot be correctly interpreted. This also demonstrates the dubious effectiveness of DSC as a screening method for co-crystal formation under the heating of binary mixtures, including ingredients melted with decomposition. This result is consistent with the previous study [22].

Table 7 Results of DSC analysis of the fourth group of benzodiazepines physical mixtures with co-formers

Conclusions

The study revealed that the DSC method can be successfully applied in order to detect a co-crystallisation process occurring in a pharmaceutical mixture when heated under non-isothermal conditions. The detailed interpretation of the DSC scans of 72 binary physical mixtures including benzodiazepines (very slightly or practically insoluble in water) and selected co-formers shows that the appearance of new endothermic or exothermic DSC peaks due to melting or crystallisation of newly obtained co-crystals, or the disappearance of DSC peaks due to the melting of benzodiazepines and co-formers, provides an argument that potential co-crystals can be found in the mixture under heating.

Based on the results obtained in this study, 15 mixtures including diazepam, temazepam, oxazepam, lorazepam and clonazepam were selected as capable of forming potential co-crystals with co-formers when heated. There are mixtures of oxazepam with glutaric, fumaric, citric and p-aminobenzoic acids and urea; lorazepam with glutaric, fumaric and citric acids, nicotinamide and saccharin; diazepam with succinic and citric acids; temazepam with glutaric acid and urea; and clonazepam with urea. Furthermore, oxazepam and lorazepam were identified as the benzodiazepines which most readily undergo co-crystallisation, each drug forming at least five co-crystals. With regard to co-formers, glutaric acid, citric acid and urea participated most frequently in the co-crystallisation process, each co-former forming at least three co-crystals. On the other hand, no mixture containing lormetazepam and estazolam formed co-crystals. A similar situation has also been observed in the case of caffeine. The DSC scans confirmed that no co-crystals were obtained in any of the caffeine mixtures.

To emphasise the importance of DSC as a rapid screening tool for the detection of co-crystallisation process in binary physical mixtures under heating, the next stage of research should include other instrumental methods to confirm the findings obtained in this study. The preliminary results outlined only indicate the direction of study to enhance benzodiazepine bioavailability.