Next Article in Journal
Dependence of the Nanoscale Composite Morphology of Fe3O4 Nanoparticle-Infused Lysozyme Amyloid Fibrils on Timing of Infusion: A Combined SAXS and AFM Study
Next Article in Special Issue
First and Second Dissociation Enthalpies in Bi-Component Crystals Consisting of Maleic Acid and L-Phenylalanine
Previous Article in Journal
Beneficial Effect of Quercetin on Erythrocyte Properties in Type 2 Diabetic Rats
Previous Article in Special Issue
Diagnosis of Agglomeration and Crystallinity of Active Pharmaceutical Ingredients in Over the Counter Headache Medication by Electrospray Laser Desorption Ionization Mass Spectrometry Imaging
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Solvatomorphism of Moxidectin

by
Toni Grell
1,
Mauro Barbero
2,
Franco Pattarino
2,
Giovanni Battista Giovenzana
2 and
Valentina Colombo
1,*
1
Dipartimento di Chimica, Università degli Studi di Milano, Via Golgi 19, I-20133 Milano, Italy
2
Dipartimento di Scienze del Farmaco, Università del Piemonte Orientale A. Avogadro, Largo Donegani 2/3, I-28100 Novara, Italy
*
Author to whom correspondence should be addressed.
Molecules 2021, 26(16), 4869; https://doi.org/10.3390/molecules26164869
Submission received: 8 July 2021 / Revised: 2 August 2021 / Accepted: 7 August 2021 / Published: 11 August 2021
(This article belongs to the Special Issue Solid-State of Organic Pharmaceutical Compounds)

Abstract

:
The solvatomorphism of the anthelmintic drug moxidectin is investigated, and a new solvatomorph with nitromethane is reported. Moreover, the hitherto unknown crystal structures of the solvatomorphs with ethanol and 2-propanol are reported and discussed. The thermal characterization of these solvatomorphs through variable-temperature powder X-ray diffraction analysis (VT-PXRD) is also described, providing new insights into the crystallochemistry of this active pharmaceutical ingredient.

1. Introduction

Polymorphism and solvatomorphism are highly important phenomena in the pharmaceutical industry [1,2]. Different arrangements of the molecules of an active pharmaceutical ingredient (API) in the solid-state give crystal structures with properties which have impact on the manufacturing process (API melting temperature, hygroscopicity and compactability), pharmaceutical behavior (solubility and dissolution rate of the API) and shelf life (stability) of medicinal products. As a consequence, the crystal form of the API may significantly affect the efficacy and safety of a medicinal product [3]. Therefore, the identification of the solid-state forms of an API and the investigation of their properties is of paramount importance for the marketing authorization of the medicine. The well-known case of the antiviral drug ritonavir [4] illustrates this relevance and the reasons why nowadays the presence of polymorphs and solvatomorphs are thoroughly studied. In contrast, the detailed investigation of their solid-state structures is often paid little attention even though this can reveal further insights that might allow for a deeper understanding and potentially even predictions [5].
This applies to moxidectin (hereafter Moxi, Scheme 1), a bioactive substance used in long-acting formulations such as tablet, oral drench, oral gel/paste, injectable, pour-on and topical solution, in the prevention of heartworm in dogs (e.g., ProHearth®) and for the treatment of internal and external parasites in cattle, sheep, horses and deer (e.g., Cydectin®, Quest®, Equest®) [6,7,8,9]. Moxidectin has recently been approved by the FDA for the treatment of river blindness [10] and is nowadays under investigation for pest control on crops. The molecular structure of moxidectin is based on a 16-membered macrocyclic lactone ring, a feature shared with the other members of the milbemycin group, as fermentation products of Streptomyces [11].
While the overwhelming part of the research on moxidectin is focused on its biological activity, preparation and formulation, investigations on the poly- or solvatomorphism are scarce. To date, one polymorph of the pure form [12] (usually referred to as Form I), an amorphous form and three solvatomorphs with ethanol (EtOH), 2-propanol (iPrOH) and 1-butanol (nBuOH) have been reported [13]. However, structural data only exist for the pure Form I (summarized in Table 1). This is not unexpected for pharmaceutical compounds, for which the complete structural information is often lacking in the scientific literature and in the patent landscape, mostly because of the difficulty in growing single crystals of metastable phases from solution (thus preventing structure determination by conventional X-ray diffraction methods).
Our research group is interested in the characterization of polymorphic drugs by single-crystal XRD [14] and, when single-crystals are not available, with ab initio methods applied to laboratory-scale powder X-ray diffraction data (PXRD), collected from microcrystalline powders [15,16]. Regarding solid moxidectin, the scientific literature and the patents deposited in the last years, witness the occurrence of different crystal forms, for which very limited structural knowledge is available. Indeed, given the lack of structural information on this drug, we have undertaken a systematic screening for crystallization solvents and conditions, with the aim of reaching a deeper knowledge of its crystal chemistry. Table 1 reports a synoptic collection of the crystal forms described in the literature and patents, while this paper reports on the understanding of the solid-state forms of moxidectin, with the comprehensive structural characterization of its solvatomorphs and the discovery of a novel one.

2. Results and Discussion

2.1. Preliminary Considerations on Moxidectin Reports

As already anticipated in the introduction, despite the recent FDA approval, very few reports have been published on the solid-state chemistry of moxidectin [12,13,17,18]. Crystalline form I was identified in 1988 and its structure revealed by single-crystal X-ray diffraction analysis (deposited in the CSD [19] with reference code GETBOW [12]). The calculated powder diffraction pattern is reported in Figure 1. Later, in 2013, a patent [13] was published which reports on the crystallization of the amorphous form of moxidectin into Form I, exploiting heating/cooling procedures in an oil bath. Moreover, the existence of three moxidectin solvates was stated, namely with ethanol, 2-propanol and 1-butanol. In this patent, these solvates are compared by DSC (Differential Scanning Calorimetry) analysis and visual comparison of the PXRD patterns. In addition to that, a graphical representation of the crystal structure, together with crystallographic space group, and unit cell metrics is reported for the ethanol solvate. However, this crystal structure has never been deposited and thus the quality and origin of the crystallographic data remains doubtful. More recently, two additional works appeared on the crystal growth of the moxidectin-ethanol solvate. In 2015, Shi et al. described the solid-state transformation of the ethanol solvate to the amorphous form, indicating that the ethanol solvate may be not stable during drying, storage and transport of the API or of solid pharmaceutical forms [18]. Moreover, in the same paper, the reversible transformation between the ethanol solvate and form I in water/ethanol mixtures was reported. Two years later, the same group studied the mechanism of formation of a gel and its subsequent transformation to form I [17]. Moxidectin shows indeed a promising rich crystal chemistry, which needs to be clarified and accurately described. Here we intend to give a comprehensive characterization of the known crystal forms of moxidectin as well as a new one, reporting detailed structural data, as a benchmark on its solid-state chemistry for the stakeholders dealing with polymorphism and solvatomorphism of this drug.

2.2. Ethanol, 2-Propanol, Methanol and 1-Butanol Solvates

Systematic crystallization trials have been performed to verify the existence of the already claimed solvatomorphs and furthermore explore yet unknown ones. From the respective solvents we could grow single crystals suitable for SC-XRD experiments for the solvatomorphs containing ethanol and 2-propanol. All our attempts to crystallize the 1-butanol solvatomorph instead failed and yielded only powders of the amorphous form, which showed the typical amorphous trace when analyzed with PXRD. Varying crystallization conditions and solvents, we furthermore discovered a hitherto unknown solvatomorph containing nitromethane (Section 2.3).
Results from SC-XRD analyses showed that the crystal structures of the ethanol and 2-propanol solvatomorphs are isostructural. The compounds crystallize in the space group P21 with one moxidectin and two solvent molecules in the asymmetric unit (here reported as Moxi·2EtOH, Moxi·2iPrOH), all lying in general position (Figure 2). Each moxidectin molecule is connected to two neighboring ones through intermolecular hydrogen bonds (HBs) between the hydroxyl groups of the hexahydrobenzofuran ring (O2 and O4, see Scheme 1 for numbering scheme). These HBs contacts (O4–H···O2, 2.760(2) Å) generate an alignment of the moxidectin molecules creating chains along the crystallographic [010] direction (Figure 3). Interestingly, the two independent solvent molecules are paired by hydrogen bonds and furthermore connected to the moxidectin molecules, thus forming a four-membered cyclic HB pattern containing two solvent (ROH) and two moxidectin molecules (Figure 4). This cyclic arrangement is based on intermolecular hydrogen bonds of the type O–H···O, which range between 2.700(3) and 2.782(4) Å. The fact that the hydrogen atom of the hydroxyl group of one moxidectin molecules (O4) is engaged in this HB pattern, and that the OH group of the involved ethanol molecule (O9) forms a bifurcated HB interacting with a second oxygen (O1) of the same moxidectin molecule (Figure 4), leads to the absence of the intramolecular hydrogen bond, which was observed in form I, Moxi·1.5MeNO2 (see below) and in other milbemycins [20]. It should be noted here that the presence of many electronegative oxygen atoms that can behave as HB acceptors, may generate weak HBs of the type C–H···O that help stabilizing the structure. As for form I, one such interaction is also observed in these structures with distances of 3.377(3) Å (Moxi·2EtOH) and 3.429(3) Å (Moxi·2iPrOH) (see discussion below and Section S3 in the Supplementary Material).
In addition to the crystal structures of Moxi·2EtOH and Moxi·2iPrOH, observations on the crystallization behavior of moxidectin in other alcohols, namely methanol (MeOH) and 1-butanol (nBuOH), have been reported in the previously cited patent [13]. The authors describe that crystallization from methanol solutions exclusively gives rise to form I. Most likely the methanol molecule is too small to stabilize the formation of the isostructural solvatomorph, as in the case of bulkier EtOH and iPrOH alcohols. Beyond that, the same patent claims the existence of a 1-butanol solvatomorph, containing 1.5 solvent molecules per moxidectin unit. Unfortunately, we failed to reproduce a sample of this 1-butanol solvatomorph, instead obtaining merely amorphous powders. However, given our findings on the isostructural nature of ethanol and 2-propanol solvates, we considered the possibility that this also applies to the chemically similar, yet bulkier, nBuOH. We assumed that nBuOH would stabilize the isomorphous crystal structure containing 2 instead of 1.5 solvent molecules, as previously claimed. In fact, the authors based this conclusion on DSC and TGA measurements, reporting that the nBuOH solvated sample was air-dried prior to analysis. This treatment could thus have induced a partial desolvation of the clathrated solvent molecules [13]. This effect was already observed for the ethanol solvatomorph, where a slow loss of ethanol molecules has been reported at relatively low temperatures (303.15 K) [18]. We confirm the same observation, as Moxi·2EtOH single crystals progressively lose their crystallinity during a room temperature overnight XRD data collection, owed to slow solvent loss. Indeed, the loss of some solvent molecules from the bulk of the nBuOH solvatomorph, at room temperature, is plausible. To support this assumption, quantum mechanical calculations (details in Section 3.4 and Section S3 in the Supplementary Material) were carried out. We performed a geometry optimization on the hypothetical isostructural solvatomorph Moxi·2nBuOH using DFT calculations. For that purpose, we modified the crystal structure of the ethanol solvatomorphs by replacing hydrogen atoms in the solvent molecules by ethyl groups. This input was subsequently optimized, and a structure was obtained which was practically identical to the given input and isomorphous to the structures of ethanol and 2-propanol. This shows that the claimed Moxi·2nBuOH structure is theoretically possible and that also nBuOH could stabilize its formation. Nevertheless, further investigations and crystal growth experiments are indispensable for the understanding of the crystallization kinetics as well as dynamics in alcohols. Furthermore, powder or SC-XRD analysis will be essential to clarify and confirm the identity of the nBuOH form.

2.3. The Nitromethane Solvatomorph, Moxi·1.5MeNO2

By keeping a nitromethane solution of moxidectin at 4 °C, crystals of a new solvatomorph have been obtained. The SC-XRD analysis revealed that in this case moxidectin has crystallized in the space group P21212 with a unit cell containing four moxidectin molecules as well as 6 nitromethane molecules, corresponding to the chemical formula Moxi·1.5MeNO2.
In the crystal structure of Moxi·1.5MeNO2, groups of three nitromethane molecules are located in a cavity and disordered around the twofold axis (Figure 5, in the crystal structure an overlap of the two positions is observed). In contrast to the already described alcoholate solvates, there are no strong HBs between the solvent molecules and the moxidectin. Instead, two intermolecular HBs are observed between adjacent moxidectin molecules (O2–H···O3, 2.869(4) Å) forming a dimer located on the twofold axis (Figure 6). In addition, the same intramolecular HB observed in crystal form I is present here (O4–H53···O1, 2.763(5)). Interestingly, six weak HBs of the type C–H···O can be observed, with oxygen being acceptor atom in both moxidectin and the solvent molecules. The distances observed for these contacts fall in the range of 3.254(6)-3.58(1) Å. This is much more compared to the solvatomorphs with ethanol and 2-propanol, where only one of these interactions is observed.

2.4. Comparative Analysis of Moxidectin Crystal Forms

A deeper understanding of the inter- and intramolecular interactions in the crystal structures of the solvatomorphs can be gained by inspecting their molecular structures. Figure 7 shows a comparison between all conformations of the moxidectin molecules in the different solvatomorphs. The conformation of the macrocycle in the different forms is generally conserved, as expected from the many structural features contributing to the reduction of its conformational freedom, represented by: (i) the ring fusion with the bicyclic hexahydrobenzofuran system, (ii) the ring fusion with the spirobipyran system, (iii) the endocyclic olefinic C=C bonds (two of them conjugated). Nevertheless, even if small, significant differences can be found. Within the four crystal structures, two conformations can be distinguished: the first is found in the alcoholate solvatomorphs (Moxi·2ROH, R = Et, iPr), whose conformations (depicted in yellow and red in Figure 7) are practically identical, further underlining their isomorphism. A different conformation is found in form I and Moxi·1.5MeNO2 (depicted in blue and green in Figure 7). Here, besides a negligible variation of the torsion angles of the terminal groups (C26–C30), a different orientation of the lactone group C1(=O1)–O6 (highlighted with dotted circle) is clearly observed. Moreover, in form I and Moxi·1.5MeNO2, the exocyclic lactone carbonyl group C1=O1 points inwards the macrocycle, thereby reducing its distance to the hexahydrobenzofuran ring. In this way, the previously described intramolecular HB with the respective hydroxyl group (O4–H53···O1, solid circles in Figure 7), can be established. In contrast to that, in the alcohol-based solvatomorphs (Moxi·2ROH, R = Et, iPr) more attractive intermolecular interactions are possible as suitable HB-donor solvent molecules are present: the hydroxyl group O4–H53 here is engaged in the HB pattern described above (Figure 4), not allowing the formation of the intramolecular HB with O1. To achieve this, the carbonyl group C1=O1 alters its torsion to point outwards from the macrocycle, consequently inducing a change in the orientation of the hexahydrobenzofuran ring. These effects sum up to have noticeable influence on the mutual orientation of the 16-membered ring (labeled A) and the spirocycle (labeled B) as visualized in Figure 7, leading to a significant increase in the puckering of the moxidectin molecule, with the latter showing a more pronounced concavity.
These considerations demonstrate the substantial relationship between the conformation and the crystal structure of the moxidectin solvatomorphs: in the presence of suitable protic solvent molecules, a HB network is established between the two co-formers, stabilizing the crystal structure (Moxi·2ROH, R = Et, iPr). Instead, in the absence of protic solvent molecules fewer strong HBs are established, and their absence is compensated by a larger number of weak intermolecular interactions of the type C–H···O (form I, Moxi·1.5MeNO2), contributing to the stabilization of the crystal structure. A list of all HB interactions is given in Table S1 (Supplementary Material).

2.5. Variable-Temperature Powder X-ray Diffraction (VT-XRPD) Analysis

In order to answer the question whether the observed solvatomorphs discussed in this work show possible additional phases or transitions, we investigated them using variable-temperature powder X-ray diffraction (VT-PXRD) analysis. This analytical tool is very useful for pharmaceutical compounds [15,16,21] or coordination polymers [22,23]. While classical thermal analyses can detect solid-state phenomena without giving detailed information about the crystallinity of the species under investigations, VT-PXRD has the power to detect all irreversible and reversible phenomena which occur at specific temperatures, giving vital information on the crystal chemistry of the investigated samples. The observable phenomena are numerous, and can be related, for example, to a solid-solid phase transition; the loss of solvent molecules with subsequent formation of new (amorphous or crystalline) phase; melting and subsequent amorphization or recrystallization toward a new crystal phase, etc. For moxidectin, we performed experiments on Moxi·2iPrOH and Moxi·1.5MeNO2. The crystalline samples were heated from room temperature to 200 °C and diffractograms taken at 5 °C intervals. For both samples, amorphization was observed. After amorphization, no further changes were observed until melting of the sample. Figure 8 and Figure 9 report the two-dimensional plots obtained for Moxi·2iPrOH and Moxi·1.5MeNO2, respectively. In the literature, the slow desolvation with concomitant recrystallization to form I has been reported for the EtOH solvate [18] in water/ethanol mixture, concluding that water plays a crucial role in this transformation. We did not observe the appearance of form I while performing the VT-PXRD experiment with Moxi·2iPrOH and Moxi·1.5MeNO2, drawing the same conclusion.
Moreover, slow recrystallization of amorphous moxidectin at high temperature has also been reported in the literature [13]. To verify this, we performed an additional experiment heating the amorphous moxidectin for a prolonged time (5 h) between 170 and 190 °C (Figure S1). However, no formation of a new crystalline phase could be detected in these conditions.

3. Experimental Section

3.1. Materials and Methods—Crystallization

Moxidectin was purchased from Merck-Sigma Aldrich. The solvents were purchased from Carlo Erba Reagent and used without further purification.
The search for experimental conditions for the crystallization of moxidectin was undertaken by systematic variation of solvents and temperature. Protic and aprotic polar solvents such as water, methanol, 1-propanol, 2-propanol, 1-butanol, ethanol, nitromethane (MeNO2), acetone, acetonitrile, THF, as well as nonpolar solvents such as toluene, cyclohexane, diisopropylether, MTBE and dichloromethane were tested.
Moxidectin was suspended in the respective solvent at different concentrations and the mixture was gently heated until complete dissolution. The solutions were placed in a closed vial to avoid evaporation of the solvent and stored at 4 °C or 25 °C. This operation was carried out in the dark to avoid photodegradation of moxidectin [24]. Crystals were obtained from the ethanol, 2-propanol and nitromethane solutions stored at 4 °C. No crystals were obtained from solutions stored at 25 °C except for the ones in 2-propanol.

3.2. Single Crystal X-ray Diffraction

Single crystal X-ray diffraction (SC-XRD) data were collected with an AXS APEXII CCD area-detector diffractometer (BRUKER). The radiation source was a molybdenum anode (Mo-Kα, λ = 0.71073 Å) with the generator working at 50 kV and 30 mA. The data reduction was carried out with CrysAlis Pro [25] version 1.171.40.53 using an empirical absorption correction with spherical harmonics (SCALE3 ABSPACK). The structure was solved by dual space methods with SHELXT-2017 [26] and refined with SHELXL-2018 [27] using the WinGX program suite [28]. Structure refinement was done using full-matrix least-square routines against F2. All hydrogen atoms with the exception of Hydroxyl-H atoms were calculated on idealized positions. The pictures were generated with the programs ORTEP-3 [28] (representation with anisotropic displacement parameters), MERCURY [29] (comparison of molecular conformations) and DIAMOND [30] (all other images). Hydrogen atoms besides for OH groups were omitted for clarity in these representations. CCDC 2089844 (Moxi·2EtOH), CCDC 2089842 (Moxi·2iPrOH), CCDC 2089843 (Moxi·2iPrOH, room temperature), CCDC 2089841 (Moxi·1.5MeNO2) contain the supplementary crystallographic data in the Supplementary Materials for this paper. These data and additional information can be obtained free of charge via https://summary.ccdc.cam.ac.uk/structure-summary-form (or from the Cambridge Crystallographic Data Centre, 12 Union Road, Cambridge CB2 1EZ, UK; Fax: (+44)1223-336-033; or [email protected]). Table 2 contains the summary of the crystallographic data for each structure. The calculation of the theoretical powder pattern from the single crystal structures was performed using the powder diffraction simulation tool implemented in Mercury [29].

3.3. Variable-Temperature X-ray Powder Diffraction (VT-XRPD) Analysis

Diffraction experiments were performed using Cu-Kα radiation (λ = 1.5418 Å) on a vertical-scan Bruker AXS D8 Advance diffractometer in θ:θ mode, equipped with a Goebel Mirror and a Bruker Lynxeye linear Position Sensitive Detector (PSD), with the following optics: primary and secondary Soller slits, 2.3° and 2.5°, respectively; divergence slit, 0.1°; receiving slit, 2.82°. Generator setting: 40 kV, 40 mA. The nominal resolution for the present set-up is 0.08° 2θ (FWHM of the α1 component) for the LaB6 peak at about 21.3° (2θ). VT-XRPD experiments were carried out by coupling a custom-made sample heater, assembled by Officina Elettrotecnica di Tenno, Ponte Arche, Italy, to the instrumental set-up described above. For that, a powdered microcrystalline sample was ground in an agate mortar and was deposited in the hollow of a quartz zero-background plate framed by an aluminum skeleton. The data were acquired within a sensible, low-angle 2θ range. The samples were heated in situ in the temperature range 20–130 °C (steps of 10 °C, Moxi·1.5MeNO2), 20–110 °C (steps of 5 °C, Moxi·2iPrOH) and 170–200 °C (over the course of 5 h, amorphous Moxi). Before each measurement we waited 5 minutes for temperature equilibration, while each measurement took about 8 minutes to be completed (0.02°/step and 0.5 s/step, as experimental conditions for the data collection). The VT diffractograms are depicted in Figure 8 and Figure 9 and Figure S1. The diffraction patterns were elaborated with the software EVA [31]. Performing the experiments both in air and in nitrogen atmosphere gave identical results.

3.4. Quantum Chemical Calculations

The calculations were carried with DFT using the software suite CRYSTAL-17 [32], applying the B3LYP functional [33,34] in combination with the basis set 6–31 G [35,36,37]. Furthermore, the atom-pairwise dispersion correction D3 [38] was used and the convergence criterion set to 10−8 Hartree.
The input for the optimization of the hypothetical solvatomorph of 1-butanol was created by modifying the crystal structure of the ethanol solvatomorph through replacing hydrogen atoms in the solvents by ethyl groups. The optimized crystal structure is isostructural to the solvatomorphs of ethanol and 2-propanol. The obtained unit cell parameters were a = 11.4065, b = 8.9996, c = 21.5342, β = 93.7797°; the theoretical temperature corresponds to 0 K. Further details, including the structure model, are listed in the Supplementary Material.

4. Conclusions

In this paper we clarified the solid-state chemistry of moxidectin API. Through SC-XRD analysis we were able to report a crystallographic characterization of moxidectin solvatomorphs containing two equivalents of ethanol and 2-propanol. Furthermore, a new solvatomorph with nitromethane was obtained and structurally characterized. Our structural analysis confirmed the isostructural nature of the alcoholate solvatomorphs (Moxi·2ROH, R = Et, iPr). Moreover, a deep analysis on their crystal structures and literature reports, together with the aid of quantum chemical calculations, lead to the conclusion that the claimed solvatomorph with 1-butanol likely possesses 2 instead of 1.5 solvent molecules in its crystal structure, being isostructural to the other alcohol-based solvatomorphs. Careful analysis of the intermolecular interactions in the reported crystal structures revealed a substantial relationship between the molecular conformation of moxidectin and the hydrogen bonds present in the solid-state of its solvatomorphs. This work offers a comprehensive view of the solvatomorphism of moxidectin and provides important information for the manufacture of medicinal products of this active pharmaceutical ingredient, which are useful to meet the strict requisites of efficacy, safety, quality mandatory for the marketing authorization.

Supplementary Materials

The following are available online. Variable-temperature X-ray powder diffraction (VT-XRPD) of amorphous moxidectin (Figure S1).; distances [Å] and type for strong and weak hydrogen bonds observed in form I and in the solvatomorphs discussed in this work (Table S1) and additional information on quantum chemical calculations.

Author Contributions

Conceptualization, G.B.G., V.C. and F.P.; preparation of crystalline samples, T.G., M.B. and G.B.G.; diffractometric analyses, T.G. and V.C.; data analyses and software, T.G. and V.C.; resources, V.C., F.P. and G.B.G.; writing—original draft preparation, T.G., M.B. and V.C.; writing—review and editing, V.C, F.P. and G.B.G.; funding acquisition, V.C. All authors have read and agreed to the published version of the manuscript.

Funding

Università degli Studi di Milano is gratefully acknowledged for support through the PSR2018 program (grant No. PSR2018_DIP_005_COLOMBO_VALENTINA to VC) and the APC initiative.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are available from the authors upon reasonable request.

Acknowledgments

Andrej Ckatroski is acknowledged for experimental assistance.

Conflicts of Interest

The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Sample Availability

Samples are not available from the authors.

References

  1. Brittain, H.G. Polymorphism in Pharmaceutical Solids; CRC Press: Boca Raton, FL, USA, 2018. [Google Scholar]
  2. Bernstein, J. Polymorphism in Molecular Crystals; Clarendon Press: Oxford, UK, 2007. [Google Scholar]
  3. Hilfiker, R.; von Raumer, M. Polymorphism in the Pharmaceutical Industry. In Solid Form and Drug Development; John Wiley & Sons Incorporated: Newark, NJ, USA, 2018. [Google Scholar]
  4. Bauer, J.; Spanton, S.; Henry, R.; Quick, J.; Dziki, W.; Porter, W.; Morris, J. Ritonavir: An extraordinary example of conformational polymorphism. Pharm. Res. 2001, 18, 859–866. [Google Scholar] [CrossRef] [PubMed]
  5. Price, S.L. The computational prediction of pharmaceutical crystal structures and polymorphism. Adv. Drug Deliv. Rev. 2004, 56, 301–319. [Google Scholar] [CrossRef]
  6. Cobb, R.; Boeckh, A. Moxidectin: A review of chemistry, pharmacokinetics and use in horses. Parasites Vectors 2009, 2, S5. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  7. Prichard, R.; Ménez, C.; Lespine, A. Moxidectin and the avermectins: Consanguinity but not identity. Int. J. Parasitol. Drugs Drug Resist. 2012, 2, 134–153. [Google Scholar] [CrossRef] [PubMed]
  8. Paul, A.J.; Tranquilli, W.J.; Hutchens, D.E. Safety of moxidectin in avermectin-sensitive collies. Am. J. Vet. Res. 2000, 61, 482–483. [Google Scholar] [CrossRef] [Green Version]
  9. El-Saber Batiha, G.; Alqahtani, A.; Ilesanmi, O.B.; Saati, A.A.; El-Mleeh, A.; Hetta, H.F.; Magdy Beshbishy, A. Avermectin Derivatives, Pharmacokinetics, Therapeutic and Toxic Dosages, Mechanism of Action, and Their Biological Effects. Pharmaceuticals 2020, 13, 196. [Google Scholar] [CrossRef]
  10. Boussinesq, M. A new powerful drug to combat river blindness. Lancet 2018, 392, 1170–1172. [Google Scholar] [CrossRef] [Green Version]
  11. Hood, J.D.; Banks, R.M.; Brewer, M.D.; Fish, J.P.; Manger, B.R.; Poulton, M.E. A novel series of milbemycin antibiotics from Streptomyces strain E225. I. Discovery, fermentation and anthelmintic activity. J. Antibiot. 1989, 42, 1593–1598. [Google Scholar] [CrossRef] [Green Version]
  12. Beddall, N.E.; Howes, P.D.; Ramsay, M.V.J.; Roberts, S.M.; Slawin, A.M.Z.; Sutherland, D.R.; Tiley, E.P.; Williams, D.J. Chemical transformations of S541 factors (A)-(D): Preparation and reactions of the 23-ketones. Tetrahedron Lett. 1988, 29, 2595–2598. [Google Scholar] [CrossRef]
  13. Cady, S.M.; Ma, B.; Chapman, R.C.; Yang, C.; Jain, U. Long-Acting Injectable Moxidectin Formulations and Novel Moxidectin Crystal Forms. Patent Number WO2013082373-A1, 6 June 2013. [Google Scholar]
  14. Barbero, M.; Mossotti, M.; Sironi, A.; Giovenzana, G.B.; Colombo, V. Crystal structure of pirfenidone (5-methyl-1-phenyl-1H-pyridin-2-one): An active pharmaceutical ingredient (API). Acta Crystallogr. Sect. E Crystallogr. Commun. 2019, 75, 984–986. [Google Scholar] [CrossRef] [PubMed]
  15. Colombo, V.; Masciocchi, N.; Palmisano, G. Crystal chemistry of the antibiotic doripenem. J. Pharm. Sci. 2014, 103, 3641–3647. [Google Scholar] [CrossRef]
  16. Tasso, B.; Pirisino, G.; Novelli, F.; Garzon, D.; Fruttero, R.; Sparatore, F.; Colombo, V.; Sironi, A. On the self-condensation of aminoguanidine leading to 1,1,4,10,10-pentaamino-2,3,5,6,8,9-hexaazadeca-1,3,5,7,9-pentaene (structure elucidation through X-ray powder diffraction). Tetrahedron 2014, 70, 8056–8061. [Google Scholar] [CrossRef]
  17. Shi, Z.; Hao, L.; Zhang, M.; Dang, L.; Wei, H. Gel formation and transformation of Moxidectin during the anti-solvent crystallization. J. Cryst. Growth 2017, 469, 8–12. [Google Scholar] [CrossRef]
  18. Shi, Z.; Wang, Z.; Zhang, T.; Dang, L.; Wei, H. The in situ monitoring of the transformation of moxidectin ethanol solvate to form I in an ethanol–water mixture. RSC Adv. 2015, 5, 98050–98056. [Google Scholar] [CrossRef]
  19. Groom, C.R.; Bruno, I.J.; Lightfoot, M.P.; Ward, S.C. The Cambridge Structural Database. Acta Crystallogr. Sect. B Struct. Sci. Cryst. Eng. Mater. 2016, 72, 171–179. [Google Scholar] [CrossRef] [PubMed]
  20. Bizdena, E.; Belyakov, S.; Jure, M.; Grinsteine, I.; Kumpiņš, V.; Turks, M. X-ray structure analysis of a solid solution of milbemycins A3 and A4. Nat. Prod. Res. 2013, 27, 1936–1939. [Google Scholar] [CrossRef]
  21. Terruzzi, S.; Bellomi, S.; Marras, G.; Barreca, G.; Ventimiglia, G.; Cervellino, A.; Masciocchi, N. Disclosing the Rich Crystal Chemistry of Lesinurad by ab Initio Laboratory X-ray Powder Diffraction Methods. Cryst. Growth Des. 2018, 18, 6863–6872. [Google Scholar] [CrossRef]
  22. Aprea, A.; Colombo, V.; Galli, S.; Masciocchi, N.; Maspero, A.; Palmisano, G. Thiazolo[5,4-d]thiazole-2,5-dicarboxylic acid, C6H2N2O4S2, and its coordination polymers. Solid State Sci. 2010, 12, 795–802. [Google Scholar] [CrossRef] [Green Version]
  23. Chierotti, M.R.; Gobetto, R.; Nervi, C.; Bacchi, A.; Pelagatti, P.; Colombo, V.; Sironi, A. Probing hydrogen bond networks in half-sandwich Ru(II) building blocks by a combined 1H DQ CRAMPS solid-state NMR, XRPD, and DFT approach. Inorg. Chem. 2014, 53, 139–146. [Google Scholar] [CrossRef]
  24. Awasthi, A.; Razzak, M.; Al-Kassas, R.; Harvey, J.; Garg, S. Analytical profile of moxidectin. Profiles Drug Subst. Excip. Relat. Methodol. 2013, 38, 315–366. [Google Scholar] [PubMed]
  25. Agilent, CrysAlis PRO; Agilent Technologies Ltd.: Yarnton, Oxfordshire, UK, 2019.
  26. Sheldrick, G.M. SHELXT—Integrated space-group and crystal-structure determination. Acta Crystallogr. A Found. Adv. 2015, 71, 3–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Sheldrick, G.M. Crystal structure refinement with SHELXL. Acta Crystallogr. C Struct. Chem. 2015, 71, 3–8. [Google Scholar] [CrossRef]
  28. Farrugia, L.J. WinGX and ORTEP for Windows: An update. J. Appl. Crystallogr. 2012, 45, 849–854. [Google Scholar] [CrossRef]
  29. Macrae, C.F.; Bruno, I.J.; Chisholm, J.A.; Edgington, P.R.; McCabe, P.; Pidcock, E.; Rodriguez-Monge, L.; Taylor, R.; van de Streek, J.; Wood, P.A. Mercury CSD 2.0—new features for the visualization and investigation of crystal structures. J. Appl. Crystallogr. 2008, 41, 466–470. [Google Scholar] [CrossRef]
  30. Brandenburg, K.; Putz, H. DIAMOND; Crystal Impact GbR: Bonn, Germany, 1999. [Google Scholar]
  31. Bruker AXS Inc. EVA; Bruker AXS Inc.: Madison, WI, USA, 2012. [Google Scholar]
  32. Dovesi, R.; Erba, A.; Orlando, R.; Zicovich-Wilson, C.M.; Civalleri, B.; Maschio, L.; Rérat, M.; Casassa, S.; Baima, J.; Salustro, S.; et al. Quantum-mechanical condensed matter simulations with CRYSTAL. WIREs Comput. Mol. Sci. 2018, 8, e1360. [Google Scholar] [CrossRef]
  33. Becke, A.D. Density-functional exchange-energy approximation with correct asymptotic behavior. Phys. Rev. A 1988, 38, 3098–3100. [Google Scholar] [CrossRef]
  34. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-Salvetti correlation-energy formula into a functional of the electron density. Phys. Rev. B Condens. Matter 1988, 37, 785–789. [Google Scholar] [CrossRef] [Green Version]
  35. Clark, T.; Chandrasekhar, J.; Spitznagel, G.W.; Schleyer, P.V.R. Efficient diffuse function-augmented basis sets for anion calculations. III. The 3-21+G basis set for first-row elements, Li-F. J. Comput. Chem. 1983, 4, 294–301. [Google Scholar] [CrossRef]
  36. Ditchfield, R.; Hehre, W.J.; Pople, J.A. Self-Consistent Molecular-Orbital Methods. IX. An Extended Gaussian-Type Basis for Molecular-Orbital Studies of Organic Molecules. J. Chem. Phys. 1971, 54, 724–728. [Google Scholar] [CrossRef]
  37. Hehre, W.J.; Ditchfield, R.; Pople, J.A. Self—Consistent Molecular Orbital Methods. XII. Further Extensions of Gaussian—Type Basis Sets for Use in Molecular Orbital Studies of Organic Molecules. J. Chem. Phys. 1972, 56, 2257–2261. [Google Scholar] [CrossRef]
  38. Grimme, S.; Antony, J.; Ehrlich, S.; Krieg, H. A consistent and accurate ab initio parametrization of density functional dispersion correction (DFT-D) for the 94 elements H-Pu. J. Chem. Phys. 2010, 132, 154104. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Scheme 1. Structural formula of moxidectin including the labels used in this paper for SC-XRD structure description (numbers in black for carbon and red for oxygen atoms).
Scheme 1. Structural formula of moxidectin including the labels used in this paper for SC-XRD structure description (numbers in black for carbon and red for oxygen atoms).
Molecules 26 04869 sch001
Figure 1. Calculated PXRD patterns (radiation Cu-Kα) for moxidectin known crystal forms: bottom to top: Form I; Moxi·2EtOH; Moxi·2iPrOH; Moxi·1.5MeNO2. Patterns were calculated from single crystal XRD using the Mercury (CCDC) powder diffraction simulation tool.
Figure 1. Calculated PXRD patterns (radiation Cu-Kα) for moxidectin known crystal forms: bottom to top: Form I; Moxi·2EtOH; Moxi·2iPrOH; Moxi·1.5MeNO2. Patterns were calculated from single crystal XRD using the Mercury (CCDC) powder diffraction simulation tool.
Molecules 26 04869 g001
Figure 2. Asymmetric unit in the crystal structure of the ethanol solvate (Moxi·2EtOH) containing the molecular structures of moxidectin and two independent ethanol molecules shown with anisotropic displacement parameters (ADPs). ADPs are shown with 50% probability. Color code: C, grey; N, blue; O, red; H, white.
Figure 2. Asymmetric unit in the crystal structure of the ethanol solvate (Moxi·2EtOH) containing the molecular structures of moxidectin and two independent ethanol molecules shown with anisotropic displacement parameters (ADPs). ADPs are shown with 50% probability. Color code: C, grey; N, blue; O, red; H, white.
Molecules 26 04869 g002
Figure 3. Representation of the chains along the crystallographic [010] direction formed by hydrogen bonds between the molecules in Moxi·2ROH (R = Et, iPr; exemplarily shown for the EtOH solvate). Color code: C, grey; N, blue; O, red; H, white.
Figure 3. Representation of the chains along the crystallographic [010] direction formed by hydrogen bonds between the molecules in Moxi·2ROH (R = Et, iPr; exemplarily shown for the EtOH solvate). Color code: C, grey; N, blue; O, red; H, white.
Molecules 26 04869 g003
Figure 4. Four-membered cyclic hydrogen bond motif between two moxidectin molecules (only fragments shown) and two ethanol solvent molecules. Selected distances O1–O9 3.255(4) Å, O2–O4 2.760(2) Å, O2–O10 2.700(3) Å, O4–O9 2.782(4) Å, O9–O10 2.719(5) Å. Colour code: C, grey; N, blue; O, red; H, white.
Figure 4. Four-membered cyclic hydrogen bond motif between two moxidectin molecules (only fragments shown) and two ethanol solvent molecules. Selected distances O1–O9 3.255(4) Å, O2–O4 2.760(2) Å, O2–O10 2.700(3) Å, O4–O9 2.782(4) Å, O9–O10 2.719(5) Å. Colour code: C, grey; N, blue; O, red; H, white.
Molecules 26 04869 g004
Figure 5. Asymmetric unit of the crystal structure of the nitromethane solvate Moxi·1.5MeNO2 showing a one moxidectin and three nitromethane molecules. For the disordered nitromethane molecules only one position is shown. Anisotropic displacement parameters are shown with 30% probability.
Figure 5. Asymmetric unit of the crystal structure of the nitromethane solvate Moxi·1.5MeNO2 showing a one moxidectin and three nitromethane molecules. For the disordered nitromethane molecules only one position is shown. Anisotropic displacement parameters are shown with 30% probability.
Molecules 26 04869 g005
Figure 6. Crystal structure of Moxi·1.5MeNO2 showing a dimer of two moxidectin molecules. Selected distances O2–O3′ 2.863(5) Å, O1–O4 2.759(6) Å. Color code: C, grey; N, blue; O, red; H, white.
Figure 6. Crystal structure of Moxi·1.5MeNO2 showing a dimer of two moxidectin molecules. Selected distances O2–O3′ 2.863(5) Å, O1–O4 2.759(6) Å. Color code: C, grey; N, blue; O, red; H, white.
Molecules 26 04869 g006
Figure 7. Comparison of the conformations of the moxidectin molecule in the crystal forms described in this paper. The graphical representation shows the different conformations of moxidectin solvatomorphs. Color code: yellow (Moxi·2EtOH), red (Moxi·2iPrOH), green (Moxi·1.5MeNO2) and blue (form I, CSD reference code GETBOW [12]). The orientation of the 16-membered ring (A) and the spirocycle (B) are indicated.
Figure 7. Comparison of the conformations of the moxidectin molecule in the crystal forms described in this paper. The graphical representation shows the different conformations of moxidectin solvatomorphs. Color code: yellow (Moxi·2EtOH), red (Moxi·2iPrOH), green (Moxi·1.5MeNO2) and blue (form I, CSD reference code GETBOW [12]). The orientation of the 16-membered ring (A) and the spirocycle (B) are indicated.
Molecules 26 04869 g007
Figure 8. Two-dimensional contour plot as a function of 2θ and temperature for the collection of PXRD patterns measured at elevating temperatures in the range of 25–100 °C for Moxi·2iPrOH. The collection of PXRD patterns shows that, by raising the temperature, the sample loses crystallinity at about 70 °C, and becomes amorphous. No peak shifts are observed. As a sign of amorphization, the lowering of peak intensities, starting from 65 °C, is observed. The sample can be considered fully amorphous at about 85 °C.
Figure 8. Two-dimensional contour plot as a function of 2θ and temperature for the collection of PXRD patterns measured at elevating temperatures in the range of 25–100 °C for Moxi·2iPrOH. The collection of PXRD patterns shows that, by raising the temperature, the sample loses crystallinity at about 70 °C, and becomes amorphous. No peak shifts are observed. As a sign of amorphization, the lowering of peak intensities, starting from 65 °C, is observed. The sample can be considered fully amorphous at about 85 °C.
Molecules 26 04869 g008
Figure 9. Two-dimensional contour plot as a function of 2θ and temperature for the collection of PXRD patterns measured at elevating temperatures in the range of 25–130 °C for Moxi·1.5MeNO2. The collection of PXRD patterns shows that the sample loses crystallinity at about 50 °C and becomes amorphous by raising the temperature. No peak shifts are observed. As a sign of amorphization, the lowering of peak intensities, starting from 50 °C, is observed. The sample can be considered fully amorphous at about 100 °C.
Figure 9. Two-dimensional contour plot as a function of 2θ and temperature for the collection of PXRD patterns measured at elevating temperatures in the range of 25–130 °C for Moxi·1.5MeNO2. The collection of PXRD patterns shows that the sample loses crystallinity at about 50 °C and becomes amorphous by raising the temperature. No peak shifts are observed. As a sign of amorphization, the lowering of peak intensities, starting from 50 °C, is observed. The sample can be considered fully amorphous at about 100 °C.
Molecules 26 04869 g009
Table 1. Synoptic collection of moxidectin crystal forms as well as structural data and literature reports.
Table 1. Synoptic collection of moxidectin crystal forms as well as structural data and literature reports.
Crystal FormFirst ReportStructure Report
Form I (Moxi, unsolvated)Ref. [12]Ref. [12], CSD: GETBOW
Amorphous (Moxi, unsolvated)Ref. [13]
Moxi·2EtOHRef. [13]this work
Mox·2iPrOHRef. [13]this work
Moxi·2nBuOHRef. [13] (reported as Moxi·1.5nBuOH)
Moxi·1.5MeNO2this workthis work
Table 2. Summary of crystallographic data.
Table 2. Summary of crystallographic data.
CompoundMoxi·2EtOHMoxi·2iPrOHMoxi·1.5MeNO2
FormulaC40H65NO10C43H69NO10C38.23H56.70N2.23O10.47
FW731.94759.99715.05
Crystal SystemMonoclinicMonoclinicOrthorhombic
Space groupP21P21P21212
a (Å)11.1913(4)11.4417(3)13.9453(5)
b (Å)8.8440(3)8.9507(2)19.8096(6)
c (Å)21.9594(8)21.7210(6)15.2680(6)
α (°)909090
β (°)93.747(3)92.301(3)90
γ (°)909090
V (Å)32168.80(13)2222.68(10)4217.8(3)
Z224
Temperature (K)120(2)120(2)100(2)
ρ(calcd) (g·cm–3)1.1211.1361.126
μ (mm–1)0.0790.0790.081
F(000)7968281542
Crystal size (mm3)0.150 × 0.120 × 0.1000.200 × 0.100 × 0.1000.120 × 0.120 × 0.030
θ range (°)1.824 to 31.5711.781 to 31.7191.684 to 26.372
Limiting indices−16 ≤ h ≤ 16,
−12 ≤ k ≤ 13,
−32 ≤ l ≤ 32
−16 ≤ h ≤ 16,
−13 ≤ k ≤ 13,
−31 ≤ l ≤ 30
−17 ≤ h ≤ 17,
−24 ≤ k ≤ 24,
−19 ≤ l ≤ 19
Reflns collected/unique a25503/13011 [Rint = 0.0249]26404/13370 [Rint = 0.0180]36172/8640 [Rint = 0.0244]
Data/restraints/param13011/9/48513370/9/5388638/2/508
Completeness to
θ = 25.242° (%)
100.0100.0100.00
Max. and min. transmission1.00000 and
0.96333
1.00000 and
0.90598
1.00000 and
0.61508
Final R indices
[I > 2σ(I)] b
R1 = 0.0525, wR2 = 0.1160R1 = 0.0448, wR2 = 0.1112R1 = 0.0536, wR2 = 0.1326
R indices (all data)R1 = 0.0720, wR2 = 0.1254R1 = 0.0549, wR2 = 0.1169R1 = 0.0584, wR2 = 0.1357
GooF c on F21.0131.0211.100
Largest diff peak
and hole (e·Å−3)
0.652 and −0.4260.355 and −0.3440.601 and −0.435
a Rint = Σ|Fo2 − Fo2(mean)|/ΣFo2, b R1 = Σ||Fo| − |Fc||/Σ|Fo|, wR2 = {Σ[w(Fo2 − Fc2)2]/Σ[w(Fo2)2]}1/2, c GooF = {S/(n −p)}1/2 = {Σ[w(Fo2 − Fc2)2]/(n − p)}1/2.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Grell, T.; Barbero, M.; Pattarino, F.; Giovenzana, G.B.; Colombo, V. Solvatomorphism of Moxidectin. Molecules 2021, 26, 4869. https://doi.org/10.3390/molecules26164869

AMA Style

Grell T, Barbero M, Pattarino F, Giovenzana GB, Colombo V. Solvatomorphism of Moxidectin. Molecules. 2021; 26(16):4869. https://doi.org/10.3390/molecules26164869

Chicago/Turabian Style

Grell, Toni, Mauro Barbero, Franco Pattarino, Giovanni Battista Giovenzana, and Valentina Colombo. 2021. "Solvatomorphism of Moxidectin" Molecules 26, no. 16: 4869. https://doi.org/10.3390/molecules26164869

APA Style

Grell, T., Barbero, M., Pattarino, F., Giovenzana, G. B., & Colombo, V. (2021). Solvatomorphism of Moxidectin. Molecules, 26(16), 4869. https://doi.org/10.3390/molecules26164869

Article Metrics

Back to TopTop