First Comparative Study of the Three Polymorphs of Bis

Publication Date (Web): June 5, 2018. Copyright © 2018 American Chemical Society. *E-mail: ... Crystal Growth & Design. Freitas, de Melo, Viana, Ferr...
0 downloads 0 Views 2MB Size
Subscriber access provided by Kaohsiung Medical University

Article

First Comparative Study of the Three Polymorphs of Bis(isonicotinamide) Citric Acid Cocrystals and the Concomitant Salt 4-carbamoylpyridinium Citrate Isonicotinamide Paul Stainton, Tudor Grecu, James McCabe, Tasnim Munshi, Elisa Nauha, Ian J. Scowen, and Nicholas Blagden Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.8b00597 • Publication Date (Web): 05 Jun 2018 Downloaded from http://pubs.acs.org on June 8, 2018

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

First Comparative Study of the Three Polymorphs of Bis(isonicotinamide) Citric Acid Cocrystals and the Concomitant Salt 4-carbamoylpyridinium Citrate Isonicotinamide. Paul Stainton§, Tudor Grecu†, Jim McCabe†, Tasnim Munshi‡, Elisa Nauha‡, Ian J. Scowen‡, Nicholas Blagden§*.

Abstract This contribution presents a previously unreported example and comparative study of a binary crystalline adduct which exhibits concomitant salt and cocrystal forms. Specifically, three new polymorphs (α, β, γ) of a novel 2:1 cocrystal of bis(isonicotinamide) citric acid are reported and the route to isolating each polymorph from aqueous solution is presented and the crystal chemistry of the system characterised. In addition, the metastable 2:1 salt (ionic adduct), was isolated as a transient (overlapping) phase under identical crystallisation conditions to the polymorphs; is also presented. As far as the authors can ascertain, this may be the first reported example of a salt-cocrystal-polymorphic concomitant system. A comparative study of the solid form landscape is presented, wherein mapping reveals the unique structural complexity of this multiple form salt – cocrystal concomitant system. The α and β forms are structurally very similar, where comparisons of packing behaviour are shown to only be different outside the glide plane. The γ form is very different in structure, where supramolecular chirality is present. In addition, further characterization of the isolated solid-state forms includes thermal, spectroscopic and computational analysis, all of which are employed to further verify the solid form landscape of the isonicotinamide - citric acid adduct and were found to have an order of stability of β, α, γ, salt; the salt being the metastable form.

§

School of Pharmacy, University of Lincoln, LN6 7TS, UK. Pharmaceutical Sciences, IMED, AstraZeneca, Macclesfield, SK10 2NA, UK. ‡ School of Chemistry, University of Lincoln, LN6 7TS, UK. * email: [email protected]

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Introduction Polymorphs are different crystal structures of the same compound composition (molecule or ion), arising from different crystallisation conditions such as temperature, solvent, and pressure. The concomitant aspects of multiple crystal forms have been extensively reviewed1,2. Screening efforts are undertaken to optimise and probe for all possible phases. A prerequisite for regulation and intellectual property of drug products includes awareness and knowledge of multiple crystal forms3,4, as the tendency to convert to different forms may be detrimental to these due diligence steps5. This is particularly important where polymorphs exhibit differing bioavaibilities6. Cocrystals, i.e. multi-component molecular crystals7 which contain neutral non-solvent molecules, are becoming an increasingly important dosage option for pharmaceuticals comparable to the often used salt forms8. Within the context of dosage form design, examples have been presented in the literature as potential crystal engineering strategies to improve the solubility9, bioavailability10 and other properties of neutral active pharmaceutical ingredients (APIs). The Cambridge Structural Database (CSD) cites several examples of polymorphism exhibited by cocrystals11.

Gadade et al.12 recently compiled a table of representative

polymorphic cocrystals. Overall, citations and the database reveal a wide interest in cocrystal phases and an emerging recognition of the polymorphic behaviour of cocrystals13,14. Consequently, a screen for polymorphism in cocrystals is a necessary step in dosage form design. The aim of this contribution was to explore the crystal phases of citric acid (CA) to investigate reducing the solubility of adduct forms through recrystallization with co-formers. Reports of cocrystals exist for a range of carboxylic acids, for example; malonic acid or oxalic acid15, with nicotinamides.

However, it was noted that citric acid and isonicotinamide (IN) was one

combination for which a structure is not found in the CSD to date.

ACS Paragon Plus Environment

Page 2 of 32

Page 3 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

CA is reported to form many salts and cocrystals. Specifically, CA does form cocrystals with the isomer of IN, nicotinamide (NIC) which has been reported by Lemmerer et al.16 to exist as a 2:1 ratio NIC:CA (CSD code CUYXUQ). Another relevant system is the cocrystal formation between CA and isonicotinic acid (INA), a very similar hydrogen bond motif structure to IN where an acid group replaces the amide group. In this CA:INA example the isonicotinic acid presented as a zwitterion in the lattice, which may have contributed to the stabilization of CA preventing proton donation and consequent salt formation (CSD code RUWGAS)17. It is important to note IN is itself polymorphic. Li et al. succeeded in isolating three different forms18. Many salts for IN have been reported with various proton donors as it readily accepts a proton at the pyridine nitrogen. For example with picric acid,19 in a 1:1 ratio where picric acid donates its single OH proton to one IN (CSD code GAHJAB). Another salt forms with the diacid squarate20 which donates two protons, one to each molecule of IN in the asymmetric unit (CSD code DAYFEP). A similar system was found with oxalic acid (CSD code ULAWAF)21. IN does not always form salts as in the case with another of Lemmerer et al.’s works22 where cocrystals were produced with several cycloalkanecarboxylic acids, some of which were polymorphic. The investigations reported herein yielded four phases containing citric acid and isonicotinamide. These were found to be 3 polymorphs of a 2:1 (IN:CA) cocrystal and a 2:1 salt. This study adds an example of a trimorphic cocrystal, and contributes, as far as the authors are aware, the first case of a salt form identified as the metastable phase transient to the polymorphs of the cocrystal forms.

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 32

Experimental Part of the analysis of the paper is to connect the crystallographic view with that of the observed crystallization outcome. This requires motif analysis; crystal growth identification and structure refinement checks to mitigate issues pertaining to complex/co crystal and salt form employed in the paper. See supplementary information.

Materials Citric acid anhydrous (99% ACROS), Isonicotinamide (99% Aldrich), Absolute Ethanol (Fisher), were purchased from Fisher Scientific or Sigma Aldrich and used without further purification. Water was used from a Millipore 18.2MΩ water purification system. Cocrystal Synthesis Cocrystals were synthesised by recrystallization in water. CA and IN were weighed in beakers at molar ratio 1:2, weighing approximately 0.75g IN to 1g CA adding 5ml water.

As the

dissolution of the coformer mixture was taking place, the solution rapidly became highly supersaturated with regards to the salt form, which began to nucleate before the coformers were fully dissolved at room temperature. To produce the three polymorphic cocrystals, dissolution of the salt form was required. The salt was dissolved by heating to 80°C using a hotplate. The samples were then removed from the hotplate and allowed to cool on the bench, where crystallization occurred for α and β within 2 hours, and γ within 24 hours. Unique to the β form, the solution appeared to turn bright yellow although the crystals remained white. PXRD Analysis Samples were analysed using a Bruker D8 Discover Powder Diffractometer. Samples were ground before being placed on a transmission well plate where they were subjected to X-ray

ACS Paragon Plus Environment

Page 5 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

radiation emitted from a Cu source (Kα, λ = 1.54056 Å) and measured at a 2θ range of 3-40°. Powder patterns were compared to simulated powder patterns generated using the cif files from single crystal analysis on the Mercury Software. The indices were assigned and presented in Figures S1-4 in the supplementary information.

SCXRD Analysis Single crystals suitable for X-ray diffraction measurements were mounted on MiTeGen DualThickness MicroMounts and analysed using a Bruker D8 Venture diffractometer with a Photon detection system. Unit cell measurements and data collections were performed at 173 K using Cu Kα radiation (λ = 1.54056 Å). Crystal data and refinement parameters are presented in Table 1 (CCDC 1837701-1837704 and Tables S5-8 in the supplementary information). Structure solutions were carried out by direct methods and refinement with SHELXL23 was finished using the ShelXle24 software. All non-hydrogen atoms were refined anisotropically, the C-H hydrogen atoms placed in idealised positions and the N-H and O-H hydrogen atoms found from the electron density map and refined with fixed bond distances and thermal parameters riding on the parent atom. In form α the CA is disordered over two positions regarding the arrangement of the acid and hydroxyl groups in the central carbon atom. Two positions were identified with occupancies of approximately 80% and 20%. OH hydrogen atoms for the 20% occupancy part were fixed with AFIX 147, while the 80% occupancy hydrogen atoms were identified from the electron density map. Because this disorder affects the hydrogen bonding, only the higher occupancy part is used in the discussion.

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The γ form was originally solved in a centrosymmetric space group (C2/c) with CA disordered over a two-fold rotation but it would not refine well (R=14.6). Changing to Cc made the refinement better (R=4.5) and the disorder was no longer present. The structure was refined as an inversion twin (Flack parameter = 0.4(3)). This structure seems to be a case of supramolecular chirality of the citric acid molecule25. The salt form has a disorder on one of the IN-amide groups where a slight twisting is seen. Two positions were identified with occupancies of approximately 66% and 34%. This disorder does not affect the hydrogen bonding.

Thermal Analysis Thermomicroscopic experiments were performed on a Linkam THMS600 stage mounted on a Zeiss microscope. Samples were heated at a rate of 10°C/min until the melting point was attained. Differential Scanning Calorimetry (DSC) and Thermo Gravimetric Analysis (TGA) was measured on a Netzsch STA449F3 from 30 to 180°C in three rates of 5°C/min, 10°C/min and 20°C/min under N2 flow at 50ml/min. Computational Analysis Morphology and attachment energy (Eatt) were calculated using Accelrys Materials studio by calculating the growth morphology using COMPASS forcefield26 using coordinates from the single crystal structures. The CSD-Materials 2017 release implementation of the UNI potentials Gavezzotti et al. was used to calculate intermolecular potentials. This approach deconvolutes the intermolecular potential contributions and illustrates the difference between the total energy of a pair of molecules in a given relative orientation and the molecules in isolation27. The calculations for each phase were done using UNI calculations28, 29 in Mercury software from the single crystal structures.

ACS Paragon Plus Environment

Page 6 of 32

Page 7 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Results and Discussion Cocrystal Synthesis and Crystal Morphology Three polymorphs: α, β, and γ along with a salt were discovered, forms α, β, and the salt were found to have a similar block morphology whereas form γ appeared as needles/plates. Images of crystal habit found in Figure 1.

Figure 1 Images of the crystal habits of forms α (a), β (b), γ (c) and the salt (d).

ACS Paragon Plus Environment

Page 8 of 32

Page 9 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

A schematic for the formation of each form is found in Figure 2. The salt is the kinetic form and is easily isolated due to rapidly forming before the coformers had fully dissolved. All forms α, β, and γ crystallised on cooling, growing at the following rates: Form α crystallised first within a few hours, due to it being the second kinetic product after the salt. β formed at very similar rates as α, though at much lower yields. γ only formed after at least 24 hours but rapidly nucleated at higher supersaturation.

Figure 2 Schematic for the formation of each phase, morphologies calculated using Accelrys Materials Studio from the CIF files.

An explanation for the rate of formation of each can be offered by examination of the induction times. Which is to say, the elapsed time after supersaturation was achieved and the onset of primary nucleation. Primary nucleation is known to occur randomly30, especially in highly viscous solutions where mass transfer is low, as in this system. Once the solution becomes supersaturated after the salt is dissolved, the α form is preferred at lower supersaturations (higher temperatures). The low supersaturation favours both α and β forms but

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

the α form was found in approximately 9 out of 10 samples formed within 2 hours, the other being the β form, suggesting α being more kinetically favoured over β. This was observed in hot stage microscopy (HSM) and differential scanning calorimetry (DSC) where melts of γ and the salt may prefer to recrystallise into the α form rather than the β form at the onset of melting. If primary nucleation had not occurred within 24 hours, the supersaturation had continued to rise as the temperature decreases and the solvent evaporates, moving further into the metastable zone where the γ form is preferred. Melting points were determined using HSM and are discussed further in the thermal analysis section but stated here from lowest to highest: Salt: 146.7°C, γ: 150.0°C, α: 154.1°C, β: 156.0°C. PXRD Analysis Results The powder patterns for each phase are shown in Figure 3. When the patterns of the α and β forms are compared there are many similarities. However, β exhibits noticeable unique peaks at 2θ 11.1° and especially the large peak at 25.0° annotated with coloured marks on Figure 3. Initially, these peaks were disregarded as an impurity or mixture of cocrystal and coformer, however, SCXRD analysis removed these possibilities. The γ form has small unique peaks at 2θ 5.5° and 11.5° with especially large peaks at 18.7° and 28.2°. The salt exhibits high crystallinity with many unique peaks notably: 13.5°, 21.2°, 24.0° and 27.3°.

ACS Paragon Plus Environment

Page 10 of 32

Page 11 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

Figure 3 Powder XRD patterns of four forms; α (blue), β (red), γ (green), salt (orange). SCXRD Structures The overall crystal chemistry, the key features of motif, conformation and packing networks are presented. Along with the appropriate packing analysis between the isolated forms. The crystal data and refinement parameters of the structures are summarised in Table 1.

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 32

Table 1 Crystal Data and Refinement Parameters of Four Forms: α, β, γ, and Salt. β C6H8O7·2(C6H6N2O)

γ C6H8O7·2(C6H6N2O)

436.38 Triclinic, P-1

436.38 Monoclinic, P21/c

436.38 Monoclinic, Cc

Salt C6H7O7·C6H6N2O· C6H7N2O 436.38 Triclinic P-1

10.4670 (4) 10.4844 (4) 10.5487 (4) 98.7808 (11) 112.535 (1) 109.0792 (12) 957.90 (6) 2 1.06 0.36 x 0.24 x 0.12 0.649, 0.754 14673 3746 3598 0.023 0.040 0.101 1.15 362 16 0.28, -0.19

10.3854 (4) 21.0202 (8) 17.6219 (7) 90 94.0777 (16) 90 3837.2 (3) 8 1.06 0.36 x 0.15 x 0.11 0.651, 0.754 40494 7534 6685 0.047 0.050 0.114 1.13 607 16 0.43, -0.25

31.4521 (12) 5.3196 (2) 11.7729 (5) 90 101.4741 (19) 90 1930.39 (13) 4 1.05 0.46 x 0.11 x 0.05 0.627, 0.754 9836 3610 3386 0.046 0.045 0.120 1.05 305 12 0.44, -0.27

9.5658 (3) 10.6257 (4) 11.3057 (4) 111.3323 (10) 102.8329 (11) 108.2681 (10) 940.21 (6) 2 1.08 0.46 × 0.20 × 0.17 0.603, 0.754 18147 3691 3473 0.035 0.035 0.093 1.05 329 10 0.32, -0.26

α Chemical formula C6H8O7·2(C6H6N2O) Mr Crystal system, space group a (Å) b (Å) c (Å) α (°) β (°) γ (°) 3 V (Å ) Z -1 µ (mm ) Crystal size (mm) Tmin, Tmax Measured refl. Independent refl. Refl. [I > 2σ σ(I)] Rint 2 2 R[F > 2σ σ(F )] 2 wR(F ) S Parameters Restraints -3 ∆ρmax, ∆ρmin (e Å )

The differences in the relative stability of each phase in the first instance can partly be attributed to the number of degrees of freedom of CA related to conformation. As CA is a triacid where two acid groups are present at each end of the molecule and one at the centre, each acid group can rotate freely and remain fixed in cisoid or transoid positions, depending on the relaxation period. CA would ideally prefer to relax into its lowest energy state, where torsion angles are largest with little steric strain as per its transoidal conformation in the α and β forms (Figure 4a). A slightly modified cisoid conformation is present in the γ form (Figure 4b). This conformation exhibits larger strain where the end acids have rotated with the OH groups and the central alcohol and acid groups are both facing upwards. This leaves the lower part of the

ACS Paragon Plus Environment

Page 13 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

molecule less able to form H-bonds, limiting integration sites. Other conformations of CA are possible such as with the salt where only one end acid group is rotated giving a twisted molecule (Figure 4c). The salt forms first with lowest thermal stability suggesting that this conformation is the least relaxed.

Figure 4 Three possible conformations of CA. a) Transoidal as in α and β. b) Cisoidal as in γ. c) A mixture of trans and cis arrangements as in the ionic form.

α and β Forms – Crystal Packing The α and β forms are considered together due to the similarity of the structures. The α form crystallised in spacegroup P-1 with two molecules of IN and one disordered molecule of CA in the asymmetric unit (Figure 5a). The β form crystallised in spacegroup P21/c with 2 molecules

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 32

of CA and four molecules of IN in the asymmetric unit (Figure 5b). In β there is no disorder of the CA – the conformations of the two CA molecules match with the main CA in α.

a

b CA

IN1 CA2

IN3

IN2 IN2 IN4

CA

IN1 d c

Figure 5 a) Asymmetric unit of form α; b) Asymmetric unit of form β; c) H-bonding dimer motif of IN; d) H-bonding dimer motif of CA. A hydrogen bonding dimer motif exists in both forms α and β with the pairing of IN and of CA respectively. IN is paired as a dimer in all four forms with IN bonded at the amide end (Figure 5c). In α and β, CA is paired as a dimer with a CA H-bonded to another CA at two points (Figure 5d). For both forms, the difference in symmetry of the two IN molecules within the asymmetric unit is evident and will be labelled IN1, IN2 etc. One-dimensional ribbons form in the order of CA, IN1, IN1, CA, IN2, IN2 where IN1 is flat and IN2 angles upwards creating a

ACS Paragon Plus Environment

Page 15 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

slight stepwise structure in 2D (Figure 6a). The same stepwise structure is evident in the β form (Figure 6b).

Figure 6 a) Hydrogen bonding in the α form showing the paring of symmetrically equivalent IN molecules and CAs. b) Hydrogen bonding in the β form showing the pairing of symmetry in equivalent IN molecules and CA. c) Crystal packing similarity of α (grey) and β showing the mismatch when extended outside the pink glide plane. CA disorder and hanging contacts hidden for clarity.

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Within the CA dimer, each molecule is hydrogen bonded to the carboxylic acid OH of CA1 to the carbonyl at the end of CA2, and vice versa. One CA is flipped so equivalent bonds can be made from one CA to the other. The remaining two H-bonds on the CA are bonded with another CA in the same way as the α form where the molecules face each other and are bonded from the central methyl acid OH to the end carbonyl. As well as being paired via their amide groups, the IN-amide nitrogens are bonded to the remaining oxygens on the CA, IN2 bonding to the carbonyl on the central methyl acid and IN1 bonding to the central alcohol group. The remaining IN is bonded to each CA from the nitrogen on the amide group to the alcohol group on the central methyl acid. Ergo, for CA1, this occurs with IN3 and for CA2 with IN1. When the α and β forms are compared, the packing of β is almost identical to that in α. Structure similarity search in Mercury31 yields a cluster of 13/15 molecules matching but when larger numbers are compared this ratio is reduced to 21/40 molecules matching. It was found that in one layer (two molecules thick between the glide planes in β) the packing is identical (Figure 6c). But when another layer above or below is added there are clear differences. The β form exhibits reflection then translational movement outside the glide plane whereas the α form only exhibits translational movement. In this way, the layering is the same in α, but different in the β form. Overall packing trends are an easy way to distinguish between the cocrystal polymorphs. When packing is observed down the a and b-axis (Figures S5 and S6 in the SI) layers are not observed. Instead, pairs of CA are observed at alternate positions between groups of mixed IN. Down the b-axis, no layering is observed but CA appears to be distributed between IN. At this angle, the difference in angle of each IN can be observed, where IN1 lays flat and IN2 is angled

ACS Paragon Plus Environment

Page 16 of 32

Page 17 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

upwards. Down the c-axis, a similar pattern is observed where CA is paired between groups of IN. Differences can be observed in the layering patterns. The α form exhibits layers that consist of rounder bumps whereas the β form layers are wavy. γ Form - Crystal Packing The γ form crystallises in space group Cc with one molecule of CA and two IN in the asymmetric unit (Figure 7a). The structure is clearly different than the α and β forms. Similar dimers exist with symmetrically inequivalent IN as in Figure 5c, although CA does not form dimers, but rather a chain of CA.

Figure 7 a) Asymmetric unit of form γ showing cisoid conformation from the side of CA. b) γ form showing H-bonding, rotated viewing CA from above.

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The central methyl acid OH is bonded to the next CA at an outer acid OH group and the inner OH (Figure 7b). Clear layers of IN1, IN2, CA are formed. The pyridine nitrogen is bonded to an outer acid OH as per both α and β forms. As well as bonding in its dimer motif, the NH2 on the amide is bonded to the carbonyls at the ends of CA. In this way, one IN is bonded to two CA. Overall CA is bonded to four IN, two symmetrically equivalent IN molecules on one side and two on the other. The central alcohol OH is bonded to the central acid OH. When packing is observed down the a-axis layers are easily distinguished (SI Figure S7). Clear layers of IN1, IN2, CA are shown. The layering creates a perforated sheet. Down the baxis layering is easily observed. The layers are equally spaced and are parallel. Down the c-axis layering is again observed. Salt Form - Crystal Packing The salt form crystallises in space group P-1 with one molecule of CA and two of IN in the asymmetric unit (Figure 8a). Although technically this is a 1:1:1 ternary system where a 4carbamoylpyridinium (4-CP) citrate salt forms with the addition of another IN molecule. An ionic bond is produced when proton transfer occurs between the central methyl acid OH and the nitrogen on the pyridine of IN1. In all the previous structures the IN-pyridine nitrogen has bonded to the carboxyl groups at the end of the citrate molecule.

ACS Paragon Plus Environment

Page 18 of 32

Page 19 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

Figure 8 a) Asymmetric unit of the salt form. b) Ring shape formed in the salt. c) Hydrogen bonding in the salt. d) Citrate 2 H-bonding arrangement. IN and 4-CP form dimers as in Figure 5c. IN bonds from the nitrogen on the pyridine ring to the carboxyl OH at the end of the citrate. This bond also shows some limited evidence of a salt formation with the slightly longer bond length between the proton and the oxygen on the end acid OH group. There is also some electron density where the proton would be if it were a true salt as with 4-CP. The arrangement of the pairing of IN and 4-CP and the ionic bond between 4-CP and the central acid OH of the citrate gives the extended ring shape shown in Figure 8b. IN is additionally bonded from the amide N-H to the central citrate hydroxyl group. 4-CP is bonded from the nitrogen on the amide to the carbonyl on the end acid group on the citrate. The citrate is consequently bonded to two IN and two 4-CP molecules. π – π stacking is observed in Figure 8c. As the structure is extended a chain of citrate molecules is observed built by two hydrogen bonding arrangements. One involves the central OH and an ending carbonyl and the other the central carboxylate and an ending acid OH (Figure 8d).

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 32

When packing is observed down the a-axis layers are observed (SI Figure S8). A layer of citrate molecules is flanked by alternate pairs of each IN which are arranged in a zigzag shape. Down the b-axis, no layering is observed. Down the c-axis layering is again observed as described in the bonding section with layers of IN intercalated.

Thermal Analysis - Differential Scanning Calorimetry (DSC) and Thermogravimetric Analysis (TGA) The DSC and TGA were measured for each sample at three rates (5°C/min, 10°C/min and 20°C/min) to determine any kinetic effects for the thermal behaviour of each phase.

The

10°C/min rate is presented in Figure 9 and the other two rates contained in Figures S9 and S10 in the supplementary information.

Figure 9 DSC and TGA of all forms Ramped at 10°C/min.

ACS Paragon Plus Environment

Page 21 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

The thermogravimetric analysis (TGA) of each sample suggests the start of decomposition after each sample melts.

The melting energy cannot be accurately calculated due to the

decomposition process interfering with the melting energies.

Table 2 Summary of DSC Data.

Sample

β

α

γ

Salt

Ramp Rate Onset 1 (°C) (°C/min)

Onset 2 (°C)

5

156.0

-

10

156.8

-

20

157.5

-

5

153.5

-

10

154.1

-

20

153.8

-

5

145.8

153.3

10

146.6

154.7

20

147.3

156.7

5

142.7

150.7

10

138.6

154.9

20

139.6

154.8

Form γ and the salt both undergo two phase changes which may be a melt and recrystallization into a form consistent with the MP of the α form and sometimes the β form. This, however, was noted as requiring further analysis to confirm the transition from salt to cocrystal by melting such as variable temperature PXRD which should be a concern for future work.

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Hot Stage Microscopy (HSM) and thermal stability The melting points were determined for all forms and presented in Table 3. Images are available in Figures S11-14 in the supplementary information. The order of melting points at high temperatures, close to the melting point but not necessarily at normal ambient temperatures suggests the order of stability of forms to be from lowest to highest; salt, γ, α and β.

Table 3 Hot Stage Microscopy Data Sample

Melting Point (°C)

β

156.0

α

154.1

γ

150.0

Salt

146.7

Computational Analysis – Morphology, Attachment Energy, and Intermolecular Potentials An analysis of the respective forms crystal habit was conducted using Accelrys Materials Studio 6 software32 implementation of HABIT. These calculations utilised all experimental structural data to further evaluate the order of stability of the different forms from a growth morphology perspective. Here, the growth rate is related to Eatt, examined for each phase, and a comparison undertaken pertaining to the morphology of each phase, which further contributes gaining insight into the overall stability of the phases. Morphology and Attachment Energy The predicted morphology and attachment energy (Eatt) was calculated for all the forms isolated, and the projection with faces indices are presented in Figure 10. A full summary table of all Eatt values is given in Table S10 in the supplementary information. The subsequent predicted morphology utilising the Eatt ,

illustrates which faces accept the integration of

ACS Paragon Plus Environment

Page 22 of 32

Page 23 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

molecules into the bulk lattice more favourably than others; thus, depicts a prediction of the slow and fast-growing faces allied to each form, which contributes to the overall morphology observed. This method of analysis develops an understanding of the relative growth of faces and was developed to improve the molecular basis of crystal morphology.

Figure 10 Calculated morphology of α (a), β (b), γ (c), Salt (d).

ACS Paragon Plus Environment

Crystal Growth & Design 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 32

α, β and the salt all have predicted block morphologies whereas γ exhibits a more needle/platelike habit, which does reflect experimental observations, (see figure 1) and future work would examine the role of solvent on morphology in greater detail.

For this contribution, these

morphology calculations were employed to explore in silico the relationship between fastest growing faces of the isolated forms. By definition, corner faces display a high attachment energy, and the slowest growing face corresponds to the face with the lowest attachment energy and are the largest expressed faces. Firstly, the attachment energies of all faces of α, β, and the salt are not very disperse suggesting similar growth rates in all directions and suggest the block type habit for these forms. In order to explore in silico the relative growth rates of each form, the key relationship is Eatt ∝ growth rate of the face, for each possible face; consequently the calculation are indicative of the slowest and fastest growing faces of the crystal forms. Overall, γ exhibit the largest energy difference of all forms which gives rise to needle-like morphology. The β form has the highest Eatt of all forms; this suggests in growth terms is the quickest growing form, implying the slowest faces of β are relatively faster growing than the fastest faces expressed on α form and the salt. The salt has the lowest overall Eatt indicating the slowest growing form. Therefore, the order of growth rate from slowest to fastest is salt, α, γ, and β. A comparison between the evaluated Elatt, and Eatt for the predicted fastest and slowest faces for each phase is summarised in Table 4 along with the relevant experimental observations. Overall, the α and salt forms retain the order of stability when both the order of appearance and energies are considered together, whereas according to the energy calculations, the β form grows the fastest but appears third, the opposite being true for the γ form.

ACS Paragon Plus Environment

Page 25 of 32 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Crystal Growth & Design

Table 4 Comparison Between Phases: Calculated Energy Parameters, Phase Appearance and Molecular Volume within the Lattice. Property The appearance of phase (grey) and calculated crystal properties

Initial or Lowest Ranking property

Final or Highest Ranking property

Order of appearance of phase sampled by PXRD

Salt (