Cocrystals of a coumarin derivative: an efficient approach towards anti-leishmanial cocrystals against MIL-resistant Leishmania tropica

This study demonstrates the synthesis of non-cytotoxic active candidates (co-crystals) of coumarin-3-carboxylic acid with various coformers to target the MIL-resistant Leishmania tropica. These promising anti-leishmanial results indicate the importance of crystal engineering by highlighting that manipulation of supramolecular architecture in the solid state can impact the biological response.


Introduction
Leishmaniasis is one of the several neglected tropical diseases (NTDs) and the ninth most burdened among infectious diseases (Berry & Berrang-Ford, 2016).It is caused by the parasitic protozoan leishmanial and transferred to mammals by the bite of a female phlebotomine sandfly vector (Desjeux, 1992).Every year, approximately 700 000 to 1 000 000 new cases of leishmaniasis are reported (WHO, 2022).In the last 70 years, no major modifications have been made in the treatment (oral or topical application) or prevention (vaccination) of leishmaniasis (Rezvan & Moafi, 2015).The available drugs have associated limitations involving toxicity, hydrophobicity and high cost (Castro-Gomes et al., 2009;Haldar et al., 2011).
The situation is motivation for many chemists to design structurally diverse libraries to identify promising hits against leishmaniasis (Hussain et al., 2014).In addition, the increase in the number of outbreaks and scarcity of safe and effective therapies against cutaneous leishmaniasis (CL) demand urgent action to develop new anti-leishmanial agents.Several natural products (Boluda et al., 2007), sulfonamide analogs (Pinheiro et al., 2019) and nitrogen-containing heterocycles (Hussain et al., 2014;Chen et al., 2008) are being explored.The toxicity and high cost of existing anti-leishmanial drugs (miltefosine, pentamidine, antimonials and amphotericin B) not only contribute towards economic burden but also prompted an energetic search for more effective treatment from a range of resources available, such as naturally occurring compounds, repurposing of the current drugs and structural modification of drug candidates.
The fundamental physicochemical and biological properties of compounds are associated with the structural features and molecular configuration in the solid state.A change in the functional group, molecular arrangement or interactions exerts a direct effect on the properties of a solid (Seddon & Zaworotko, 1999).Cocrystals are composed of two or more neutral molecules in a crystal structure with a distinct stoichiometry.They are formed via non-covalent interactions such as hydrogen bonding, �-stacking and van der Waals forces, as well as halogen bonding (Mandal et al., 2019;Bauza ´et al., 2016).Crystal engineering is a well established approach for designing organic solids with a wide range of applications in the field of pharmaceutical sciences (Moorthy et al., 2010).Pharmaceutical cocrystals are often designed based on crystal engineering approaches that are effective in improving physicochemical properties of clinical relevance (Desiraju & Parshall, 1989;Bolla & Nangia, 2016;Malamatari et al., 2017).Moreover, the US FDA's consideration of cocrystals as new and legitimate forms of active pharmaceutical ingredients (APIs) further contributed to a rise in the interest of pharmaceutical manufacturers towards the development of certain cocrystals as new drug leads (Kale et al., 2017).In the literature there are many studies that report improvement in the biological activity of pharmaceutical ingredients via cocrystallization (Aakero ¨y et al., 2011;Nascimento et al., 2021).Pharmaceutical cocrystals have been widely used in industries and academia in the last two decades to improve the ADME (absorption, distribution, metabolism and excretion) properties of APIs, as well as bioavailability, solubility, chemical stability, hygrostability, dissolution rate, tabletability etc. (Thakuria et al., 2013;Kumari & Ghosh, 2020;Box et al., 2016).As a result, numerous studies covering the fundamental aspects of cocrystallization have been published.The literature provides several examples of cocrystals, such as ertuglifozin L-pyroglutamic acid, sacubitril-valsartan, escitalopram oxalate-oxalic acid and termidol-celecoxib that are currently on the market or in various clinical trial phases (Kaduk et al., 2021;Videla et al., 2017).These indicate that cocrystallization is an effective approach for enhancing the physicochemical properties of APIs.

Single-crystal X-ray diffraction
SCXRD analyses of all single-crystals were carried out on a Bruker D8 venture (Germany), fitted with a photon detector with CMOS 100 technology.The crystals were irradiated by graphite-monochromated Cu K� radiation (� = 1.54178A ˚) at 100 (2) to 300 (2) K. Integration and reduction of data were completed using the Bruker SAINT software (Bruker, 2016).The structures were solved by direct methods and Fourier transformation techniques using the SHELXL program (Sheldrick, 2015).Structures were refined by full-matrix leastsquares calculations on F 2 .All non-hydrogen atoms were refined with anisotropic displacement parameters and placed at geometrically idealized positions, and all hydrogen atoms were located by difference maps and refined isotropically.The inter-molecular interactions between the molecules were calculated using PLATON (Spek, 2003).The crystal-packing diagrams and 3D structures were drawn using Mercury (Macrae et al., 2008) and ORTEP (Farrugia, 1997).Crystallographic and refinement data are summarized in Table 1.

Powder X-ray diffraction
The bulk samples of all synthesized cocrystals were characterized via PXRD analysis on a Bruker D8 Advance diffractometer equipped with a LynxEye detector and monochromatic Cu K� radiation (� = 1.54060A ˚) sources at 25 � .The powdered samples were placed in an acrylic sample holder.The data were collected initially within the range 5 to 65 � (2�) with a step size of 0.036 � .In order to determine the full structure, a continuous scan mode was used.

Hirshfeld surface analysis
Hirshfeld surfaces and 2D fingerprint plots were generated with Crystal Explorer (version 17.5; Spackman & Jayatilaka, 2009;Mackenzie et al., 2017) using the automatic procedures implemented in the software.These surfaces were mapped with a normalized contact distance (d norm ), shape-index, curvedness and ab initio electrostatics surface parameters, with automatic values.

Fourier transform infrared studies
FTIR spectra were recorded for CU, the coformers and its cocrystals on a Bruker Vector 22 FTIR spectrometer (Germany).All samples were analyzed via the KBR disk technique, and a spectrum was collected under identical conditions, the spectrum scan range was 400 to 4000 cm À 1 with a resolution of 2 cm À 1 and an accumulation of 10 scans.

Thermal analysis
DSC and TGA were performed on a LINSEIS STA PT1600 with heating rate of 10 � C min À 1 .About 20-25 mg of samples were crimped in a ceramic pan and scanned from 30 to 600 � C under dry N 2 gas purging.The Linseis TA software (version 1.0; Linseis Messgeraete GmbH) was employed for collecting data.

Biological screening
In vitro biological activities of CU, coformers and synthesized cocrystals 2-6 were evaluated for their MIL-resistant L. tropica promastigotes and cytotoxicity against T3 normal mouse fibroblast cell line.Detailed methodologies of the biological assays are provided in the supporting information.

Statistical analysis
Three replicates were used in each experiment, unless otherwise stated.All results were presented as mean standard deviations.A one-way ANOVA was used to analyze statistical differences at P < 0.05 (95% confidence interval) in conjugation with Tukey's Multiple Comparison Test using the Graph Pad Prism software (version 5; California, USA; https://www.graphpad.com).

Selection of coformers
Based on the literature review, functional groups capable of forming supramolecular synthons via hydrogen bonds such as acid� � �acid (COOH� � �COOH) and acid� � �amino (COOH� � �NH 2 ) are the essential structural features known to facilitate the formation of cocrystals (Nugrahani & Jessica, 2021;Desiraju, 1995).The coformers in the present study evidenced the above statement as 1b, 1c, 1d and 1f possess an amino-pyridine functionality and are well known for forming dimeric hetrosynthons in crystal structures.On the other hand, coformer 1e showed the carboxylic acid homosynthons motif.Moreover, the present study revealed that all the coformers  (1b-1f) demonstrate non-cytotoxicity against the normal 3T3 fibroblast cell line.Therefore, it was considered worthwhile to explore the potential of conformers not only as supramolecular synthons, but also as coformers of bioactive cocrystals.

Hirshfeld surface analysis
The Hirshfeld surface analysis was used to quantify the nature, regions and types of inter-molecular interactions in the crystal structure via mapping their properties in various modes, such as d norm , shape-index, curvedness, electrostatic potential surface and 2D plots.The dark-red and blue regions indicate the shorter (close contacts) and longer (distant contacts) distances in comparison with the van der Waals radii, respectively, and the white regions reflect a distance equal to the sum of the van der Waals radii (Venkatesan et al., 2016).The darkest red spots on the Hirshfeld surface exhibit the O-H� � �O, O-H� � �N and N-H� � �O contacts.These strong inter-molecular interactions facilitate the formation of cocrystals (Fig. S1).The Hirshfeld surfaces, mapped over the shape-index and curvedness surface, are depicted in Figs.S2 and S3.These surfaces were used to present weak inter-molecular interactions in the cocrystal and the overall packing in the crystal structure.The presence of blue and red triangles on shape-index surfaces and flat green regions on the curvedness indicate the C-H� � �� or �-stacking in cocrystals.Another Hirshfeld surface was mapped over the calculated ab initio electrostatic potential on the Hartree-Fock (HF) level of theory using the 6-311G(d.p)basis set.Fig. S4 shows that the positive electrostatic potential over the surface are hydrogen-bond donors (blue regions) and the negative electrostatic potential are the hydrogen-bond acceptors (red regions). The

FTIR analysis
FTIR analyses of CU (1a), 1b, 1c, 1d, 1e, 1f and their cocrystals 2-6 were performed and are presented in the supporting information (Figs.S10-S14).In cocrystal 2 (CU:1b), the -NH 2 stretching vibrations appear to be red-shifted at 3363 cm À 1 compared with 3459 cm À 1 observed for 1b.CU (1a) revealed a stretching C O bond of the lactone carbonyl at 1745 cm À 1 , whereas the red-shifted absorption band appeared at 1723 cm À 1 in cocrystal 2. The slight red-shift in the C O bond of the acid carbonyl from 1683 to 1680 cm À 1 in the cocrystal and blueshift in C-O from 1225 to 1250 cm À 1 clearly support the involvement of hydrogen bonding in the formation of cocrystal 2 (Fig. S10).Similarly, in cocrystal 3 (CU:1c), the -NH 2 -stretching vibration appeared as a strong band at 3386 cm À 1 which was found to be red-shifted in comparison with the -NH 2 stretching vibration (3504 cm À 1 ) in 1c.Similarly, in cocrystal 3, a strong absorption band at 1759 cm À 1 appeared due to the C O of the lactone moiety and showed blue-shifting from 1745 cm À 1 , the stretching frequency of the lactone carbonyl in CU.In addition, the blue-shift in C-F from 1329 to 1334 cm À 1 indicates the involvement of the -CF 3 functionality of the coformer in hydrogen bonding (Fig. S11).In cocrystal 4 (CU:1d), the C O (lactone carbonyl) absorption band appeared at 1731 cm À 1 and revealed a red-shift compared with the stretching frequency observed for CU (1745 cm À 1 ).The blue-shifted olefinic C C bond-stretching frequency from 1610 to 1625 cm À 1 and red-shift in stretching frequency of C O (acid carbonyl) from 1683 to 1663 cm À 1 further support the involvement of the carboxylic acid moiety of CU in hydrogen bonding with conformer 1d.Furthermore, the broadening of -NH 2 and carboxylic -OH absorption bands (3500 to 2756 cm À 1 ) in cocrystal 4 indicates strong hydrogen bonding between the acid and amine groups (Fig. S12).In the IR spectrum of cocrystal 5 (CU:1e), the stretching frequencies signify that both the -COOH group of CU and 1e do not get deprotonated.The IR spectra showed three intense bands at 1748, 1675 and 1632 cm À 1 .The strong band at 1675 cm À 1 of the cocrystal may be due to overlapping of the acid carbonyl group of CU (1683 cm À 1 ) and that of coformer 1e (1667 cm À 1 ).In 1e, stretching of the N-H hydrogen bond at 3462 cm À 1 was observed, whereas in the cocrystal the N-H stretching vibration was observed at 3460 cm À 1 with a redshift (Fig. S13).In cocrystal 6 (CU:1f) the C O carbonyl absorption band of lactone appeared at 1738 cm À 1 , whereas in CU it appears at 1745 cm À 1 , this red-shift indicates the involvement of C O (lactone carbonyl) in hydrogen bonding with the 1f coformer.The C O of the acid carbonyl from 1683 to 1703 cm À 1 was attributed to red-and blue-shifts in the olefinic bond-stretching frequency from 1610 to 1613 cm À 1 .The -NH 2 stretching was observed as a sharp band at 3431 cm À 1 in 1f, whereas in cocrystal 6 the -NH 2 stretching was observed with a blue shift at 3330 cm À 1 (Fig. S14).In conclusion, the red and blue shifts of the characteristic functional group stretching frequencies in the IR spectra of synthesized cocrystals 2-6 clearly demonstrate the role of hydrogen bonding in cocrystallization.

Thermal analysis
To analyze the thermal behavior of newly synthesized cocrystals, DSC and TGA measurements were performed.The thermal properties of the synthesized cocrystals were significantly different from those of the pure APIs and coformers.DSC and TGA thermograms of pure CU (1a), coformers 1b-1f and their synthesized cocrystals 2-6 are presented in the supporting information (Fig. S15-S19).The DSC curve of CU (1a) revealed a eutectic endotherm at 191.78 � C whereas 1b revealed two endotherm peaks at 67.29 and 199.8 � C.However, the DSC spectra of cocrystal 2 exhibited a small endotherm at 113.6 � C, followed by a second larger endotherm at 189.1 � C, indicating the development of a new solid phase.The TGA analysis of cocrystal 2 showed it is thermally stable up to 113.8 � C, followed by a percentage mass loss of 7.2% with a temperature increase up to 189.1 � C; 99.9% mass loss occured at 491.8 � C, compared with CU and 1b, which showed thermal stability up to 191.78 and 67.29 � C, respectively (Fig. S15).The DSC spectrum of cocrystal 3 was found to have an exothermic peak at 160.  increased temperature up to 310.30 � C (Fig. S17).The DSC thermogram of cocrystal 5 showed an endothermic peak at 170.5 � C, different from CU (191.78 � C) and 1e (186.0 � C).This clearly indicates the development of a new solid state (i.e.cocrystal 5), which was found to be stable up to 170.5 � C with a mass loss of 7.20%.The TGA curves of synthesized cocrystal 5 revealed noticeable changes in the thermal decomposition pattern (Fig. S18).The DSC spectra of cocrystal 6 exhibited an endotherm peak at 149.0 � C which indicates the development of a new crystalline phase.Moreover, TGA of cocrystal 6 revealed that it is thermally stable up to 123.58 � C, with a loss of 5.9%.In cocrystal 6, we noted that, on the basis of TGA curve analysis, the decomposition of the cocrystal appears to pass through three stages, although further work is required to better understand this mechanism (Fig. S19).

MIL-resistant L. tropica
The alkylphosphocholine drug MIL is a broad-spectrum drug that is active against various parasitic species, cancer cells, as well as against a number of pathogenic fungi and bacteria (Dorlo et al., 2012).Knowledge about MIL resistance in L. tropica is limited to defects in drug internalization (defected inner translocation of MIL) and increased drug efflux (Pe ´rez-Victoria et al., 2006).According to Hendrickx et al. (2014Hendrickx et al. ( , 2012)), when emergence of any degree of resistance occurs in the MIL-resistant culture, the resistance does not revert back to wild-type (WT) phenotype, despite the removal of MIL-selective pressure.
L. tropica MIL-unresponsive/resistant parasites were developed and maintained using a step-wise selection of the drug MIL up to a concentration of 196 mM.No significant differences in growth patterns were observed between WTand MIL-resistant strains.Fig. 9 shows that no inhibitory effects of MIL on cell proliferation of L. tropica were observed, demonstrating successful emergence of resistance via a dose-dependent increase of MIL.
A potential disadvantage in the use of MIL in leishmanial assay is the emergence of in vitro drug resistance (Varela-M et al., 2012).Furthermore, this drug is found to be potentially teratogenic, and is not recommended for pregnant women (Committee, 2010; Murray et al., 2005).Mechanisms that are responsible for the resistance acquisition in the L. tropica parasite against MIL include reduction in drug uptake, increased efflux and alteration in permeability of the plasma membrane (Pe ´rez- Victoria et al., 2003;Seifert et al., 2003;Kulshrestha et al., 2014;Sa ´nchez-Can ˜ete et al., 2009;Mondelaers et al., 2016).
Hence it is a necessity to find alternative therapeutic options for leishmaniasis.During the current study, the resistant strain was generated and a series of cocrystals were evaluated against the parasitic line.

Inhibitory potential of the synthesized cocrystal against MIL-resistant L. tropica in vitro
Understanding the role of mixtures of molecules in any biosystem is very complicated and difficult to understand.The synergistic effect on biological activities due to a multi molecular system is well reported in the literature, perhaps the best explanation for biological activities of plant extracts (the complex mixture of natural products).The co-crystals are systematically designed multi-component molecules in a stoichiometric ratio.Therefore, the orientation of cocrystal components due to hydrogen bonding is responsible for the change in physicochemical and biological properties compared with the individual components.The present study demonstrates the susceptibility of a synthesized series of cocrystals (2-6) of CU (1a) with coformers (1b and 1c) against the MIL-resistant L. tropica.Cocrystal CU:1d (4) appeared as a potent (<0.05) anti-leishmanial agent against resistant promastigotes with a IC 50 value of 48.71 � 0.75 mM against the tested standard drug (IC 50 = 169.55� 0.078 mM).Cocrystal CU:1b (2) appeared to be the second most potent (<0.05) antileishmanial agent with a IC 50 value of 61.83 � 0.59 mM, followed by cocrystal CU:1c (3, IC 50 = 125.7 � 1.15 mM).Among all coformers, only 1f showed potent anti-leishmanial effects (IC 50 = 78.0� 0.096 mM); however, the synthesized cocrystal of CU with 1f (6) appeared to be inactive and therefore demonstrated the role of supramolecular features in the modification of the orientation of the API and coformer and finally the molecular properties in a way to make the molecule inactive.All 1:1 physical mixtures of APIs and coformers also appeared to be inactive.In the case of mixtures, the reason for the complete loss of anti-leishmanial activity cannot be explained clearly; however, both the role of concentration and the free dispersion of CU and coformers (1b-1e) in the system could be possible reasons (Table 3).

Conclusions
Five new non-cytotoxic cocrystals of coumarin-3-carboxylic acid with pharmaceutically acceptable coformers were successfully synthesized via a neat grinding approach in a 1:1 stoichiometric ratio.Hirshfeld surface analysis demonstrated the impact of various non-covalent interactions towards the stability of the cocrystal in the solid state.Importantly, the anti-leishmanial activity evaluation against the MIL-resistant L. tropica revealed that synthesized cocrystals are more effective and non-toxic anti-leishmanial candidates compared with tested standard miltefosine against the resistant lines of clinical isolates of cutaneous leishmaniasis.Evidence that modification of supramolecular features via co-crystallization contributed towards anti-leishmanial activity is further supported by the fact that the physical mixtures (1:1) of API and amitrole were found to be inactive.Although further studies are required, the current work emphasizes the importance of cocrystallization of commercially available candidates with suitable coformers to enhance their therapeutic potential.

Figure 1
Figure 1Schematic representation of the synthesis of cocrystals 2-6 by neat grinding in a mixer mill.
overall 2D fingerprint plots resolved into all types of contacts (H� � �H, O� � �H, C� � �H, C� � �N, N� � �H, C� � �O, C� � �C, H� � �F, F� � �F, H� � �N, O� � �O, Br� � �H, Br� � �C and N� � �O) and their relative percentage populations are shown in the bar graph (presented in Fig. 7).The main contacts (H� � �H, O� � �H, N� � �H and H� � �F) are the major contributors towards the formation of the Hirshfeld surface.The O� � �H interactions are depicted in Figs.S5-S9 by inside sharp spikes, the N� � �H contacts are revealed by sharp edge spikes and H� � �H contacts are indicated by the main body of the fingerprint plots.The O� � �H interactions make up the largest proportion, indicating that they are the main contributors to the stabilization of cocrystals 2-6.

Figure 7
Figure 7Bar plot representing the 2D fingerprint plots, showing the percentage contributions of the contents of cocrystals 2-6.
4 � C, whereas a sharp endotherm appeared at 180.04 � C, distinctly different from CU (eutectic melting endotherm at 191.78 � C).The coformer 1c exhibited two sharp endotherms at 48.38 and 175.56 � C, demonstrating the new crystalline phase, i.e. cocrystal 3. The TGA profile of cocrystal 3 exhibited a thermal stability up to 160.4 � C with a percentage mass loss of 17.3% and complete mass loss of 99.68% with an increased temperature up to 250.10 � C. The results indicate that the synthesized cocrystal was stable up to 160.4 � C, compared with CU (191.78 � C) and 1c (48.38 � C) (Fig. S16).The DSC spectrum of cocrystal 4 exhibited a eutectic endotherm at 178.73 � C, which differed from pure CU (191.78 � C), demonstrating the new cocrystal phase.TGA of cocrystal 4 revealed thermal stability up to 178.73 � C with a mass loss of 18.15% and with complete mass loss of 99.89% with an

Figure 9 MIL
Figure 9 MIL-resistant L. tropica line generated by step-wise selection.1 � 10 6 log-phase promastigotes were incubated in the presence of a range of drug concentrations.The surviving cells were quantified with trypan blue dye.Populations of parasites were grown in increasing concentrations of MIL, showing increased resistance to MIL.Bars of both resistant and WT parasites represent the more or less similar growth patterns.