Hepatoprotective Cocrystals and Salts of Riluzole: Prediction

aSolid State Pharmaceutical Research Group (SSPRG), Department of Pharmaceutics, .... Riluzole was kindly gifted by Suven Life Sciences, Hyderabad, In...
0 downloads 0 Views 2MB Size
Subscriber access provided by READING UNIV

Article

Hepatoprotective Cocrystals and Salts of Riluzole: Prediction, Synthesis, Solid State Characterization and Evaluation Balvant Yadav, Sridhar Balasubramanian, Rahul B Chavan, Rajesh Thipparaboina, V GM Naidu, and Nalini R. Shastri Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.7b01514 • Publication Date (Web): 11 Jan 2018 Downloaded from http://pubs.acs.org on January 11, 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 free 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 accessible to all readers and 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.

Crystal Growth & Design 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 33 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

2

Hepatoprotective Cocrystals and Salts of Riluzole: Prediction, Synthesis, Solid State Characterization and Evaluation

3

Balvant Yadava, Sridhar Balasubramanianb, Rahul B Chavana, Rajesh Thipparaboinaa,

4

V.G.M. Naiduc, Nalini R Shastria*

1

5 6 7 8 9 10 11

a

Solid State Pharmaceutical Research Group (SSPRG), Department of Pharmaceutics, National Institute of Pharmaceutical Education and Research (NIPER), Hyderabad, 500037, India. b

X-ray Crystallography Division, CSIR-Indian Institute of Chemical Technology, Hyderabad, Telangana, 500007, India c

Department of Pharmacology, National Institute of Pharmaceutical Education and Research (NIPER), Hyderabad, 500037, India

12

* Corresponding author. Nalini R Shastri

13

Tel. +91-040-23423749

14

Fax. +91-040-23073751

15

E-mail: [email protected], [email protected]

16 17

Address: Department of Pharmaceutics, National Institute of Pharmaceutical Education & Research (NIPER), Balanagar, Hyderabad, India, Pin code – 500037

18 19 20 21 22 23 24 25 26 27 28 29

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

30

Abstract

31

Riluzole is a drug, used to slow the course of amyotrophic lateral sclerosis. Due to its unique

32

structure and functionalities, it is able to form both salts and cocrystals. This is a BCS class II

33

drug with poor solubility and causes hepatotoxicity which limits its application. The present

34

study aims towards development of novel solid forms of riluzole to address the said

35

limitations. Apart from this, an attempt has been made to develop a prediction model using

36

software tools to identify the appropriate synthons for formation of cocrystals. It was

37

observed that out of 33 coformers selected, prediction results were in agreement with the

38

experimental outcome for 25 coformers, which demonstrated the potential of the model

39

developed. Seven new solid forms of riluzole, five cocrystals with ferulic acid, syringic acid,

40

vanillic acid, cinnamic acid and proline; two salts with 2,4 dihydroxybenzoic acid and

41

fumaric acid were successfully developed. All the solid forms were characterized by DSC,

42

powder XRD, FTIR and single crystal XRD. Single crystal X-ray analysis of the all solid

43

form shows R22(8) motif between riluzole and coformers through N-H···O and O-H···N bond

44

except riluzole-proline zwitterionic cocrystal. In riluzole-fumaric acid, partial proton transfer

45

of O to N due to acidic H atom disorder has been observed. Dissolution profiles of all the

46

solid forms were comparable to that of plain riluzole and complete drug release was observed

47

within 60 min for all systems. In vivo hepatotoxicity study with riluzole-ferulic acid and

48

riluzole-syringic acid in mice model revealed its potential hepatoprotective effect to

49

counterattack the hepatotoxic adverse effects of riluzole.

50

Keywords: Crystal engineering, prediction model, salts, hepatotoxicity, nutraceuticals

51

ACS Paragon Plus Environment

Page 2 of 33

Page 3 of 33 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

52

1

Introduction

53

In recent times, solid state manipulation of active pharmaceutical ingredients (APIs) have

54

gained popularity to fortify the pharmaceutical and biopharmaceutical attributes such as

55

solubility, dissolution, bioavailability and stability. These solid state manipulation strategies

56

usually involve generation of new polymorphs,1 salts,2 amorphous systems,3 co amorphous

57

systems,4 and cocrystals.5 Among these, salts can effectively improve the pharmaceutical

58

attributes of ionisable API where as cocrystals are reported for alteration of physicochemical

59

properties such as solubility,6 compressibility,7 bioavailability,8 melting point,9 and

60

hygroscopicity10 of poorly ionisable or non-ionisable APIs. Designing a new cocrystal

61

requires in-depth understanding of chemistry of molecules, which can be solved by crystal

62

engineering approach. 'Crystal engineering’ as defined by Desiraju mainly deals with the

63

understanding of intermolecular interactions and utilization of such understanding in

64

designing new solids with desired physical and chemical properties.11 Synthon approach is

65

commonly used for generation of cocrystals which decides the direction and strength of bond.

66

Apart from synthon approach, various prediction models are reported in the literature to aid

67

in the selection of appropriate coformers for cocrystallization. Most of the prediction model

68

are based on the calculation of difference in lattice energies of cocrystal and its component,12

69

and calculation of gas phase molecular electrostatic potential surfaces of the individual

70

components.13 These prediction models possess some limitations. Hence, there exists a scope

71

to develop a robust and more accurate prediction tool which will speed up the process of

72

cocrystal generation by reducing number of the screening experiments.

73

Riluzole (2-amino-6-trifluoromethoxybenzothiazol) is used to slow the course of

74

amyotrophic lateral sclerosis (ALS).14 Riluzole exhibits strong neuroprotection,15

75

anticonvulsant,16 antidepressant effects17 and sedative properties.18 The molecular weight of

76

the compound is 234.20 g/mol and the empirical formula is C8H5F3N2OS. It is almost

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

77

insoluble in water at pH 7.4. Mondal P et al., had reported various cocrystals of this drug

78

with succinic acid, glutaric acid, adipic acid, sorbic acid, pimelic acid, suberic acid, azelic

79

acid, sebacic acid, 4- dimethylaminopyridine and nicotinic acid and a malonate salt for

80

increasing the solubility and dissolution.19

81

Severe hepatotoxicity was reported in clinical data for riluzole, which limits its use in

82

patients.20, 21 This necessitates careful consideration in the selection process of coformers. In

83

addition to cocrystal formation ability, the selection of coformers should entail evaluation of

84

the toxicity profile of the coformers, the potential therapeutic benefits provided and their

85

ability to tackle the hepatotoxicity problem of riluzole. Apart from the routine list of

86

pharmaceutically acceptable salt/ cocrystal formers, there is need to consider an additional

87

group of molecules as possible coformers i.e. naturally occurring compounds, nutraceuticals.

88

These naturally occurring coformers besides aiding cocrystal formation may additionally

89

provide some physiological benefits to the patients such as antioxidant, anti-inflammatory,

90

anti-aging and hepatoprotective activity. Hence, the present study was designed with the

91

objective of developing the novel solid forms i.e. cocrystals or salts using hepatoprotective

92

coformers, mainly nutraceuticals like ferulic acid, vanillic acid and syringic acid. These

93

nutraceuticals selected had already proven their hepatoprotective22-25 potential in addition to

94

their antioxidant,26 and antiepileptic27 properties. Additionally, an attempt was made to

95

develop a simple prediction model for screening of the coformers. The model was designed

96

to predict the possible synthon groups in drug and coformers which would participate in

97

hydrogen bond formation during cocrystal generation. Total 33 coformers were employed in

98

this study; structures of these coformers are represented in figure S1, supporting information.

ACS Paragon Plus Environment

Page 4 of 33

Page 5 of 33 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

99

2

Experimental section

100

2.1

Material

101

Riluzole was kindly gifted by Suven Life Sciences, Hyderabad, India. Ferulic acid, syringic

102

acid, vanillic acid, cinnamic acid, 2,4 dihydroxybenzoic acid and fumaric acid were

103

purchased from Alfa Aesar, India. Proline was purchased from Avra Synthesis Private Ltd.,

104

India. All other solvents and chemicals were of analytical grade. HPLC grade acetonitrile was

105

purchased from Merck, India. Double distilled water was generated in house.

106

2.2

107

Hydrogen bond formation between riluzole and coformers were predicted using Materials

108

Studio 7.0 (Accelrys Inc., USA) software tools. Amorphous cells were constructed to identify

109

the interaction sites and possible synthons. Briefly, MOL files were downloaded from

110

Chemical Entities of Biological Interest (CHEBI). Amorphous cell was constructed, after

111

structure clean up and adjustment of hydrogen atoms with the following parameters for

112

temperature, cell type and configurations at 298 K, periodic and 1-100 respectively. Post

113

processing involved calculation of hydrogen bonds and assignment of bond length to identify

114

the possible synthons.

115

2.3

116

2.3.1 Liquid assisted grinding

117

Preliminary screening for novel forms was carried out by grinding drug and coformers in 1:1,

118

1:2 and 2:1 molar ratio in an agate mortar in the presence of 120 µL of methanol for 30 min.

119

Riluzole-ferulic acid, riluzole-syringic acid, riluzole-vanillic acid, riluzole-cinnamic acid,

120

riluzole-proline, riluzole-2,4-dihydroxybenzoic acid and riluzole-fumaric acid were

121

successfully synthesized by liquid assisted grinding.

Prediction and Modelling

Solid form screening

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

122

2.3.2 Recrystallization experiments

123

Recrystallization experiments were carried out by dissolving the drug and coformer (1:1

124

molar ratio) in methanol or ethanol in a beaker at 55-65 °C. The beaker was covered with

125

pin-holed aluminium foil and left at room temperature for slow evaporation. Single crystals

126

obtained after 2-3 days were filtered, air-dried, and stored in a glass vial until further analysis.

127

Information regarding the crystallization experiments is provided in supporting information

128

table S1. Single crystals were obtained from all recrystallization systems except riluzole-

129

syringic acid, wherein single crystals suitable for structural analysis were not obtained.

130

2.4

131

2.4.1 DSC

132

Thermal analysis of plain riluzole and novel solid forms were carried out using Mettler

133

Toledo DSC system operating with Stare software to determine the melting point and

134

enthalpy. DSC instrument was calibrated for temperature and heat flow using high purity

135

indium standard. Accurately weighed samples (5−10 mg) were scanned at a heating rate of 10

136

°C /min from 20 to 250 °C in aluminium crimped pans with pinhole. The measurements were

137

conducted in nitrogen gas environment with a purging rate of 60 mL/min.

138

2.4.2 Powder X-ray Diffraction (P-XRD)

139

P-XRD patterns of samples were recorded at room temperature using PANalytical X’Pert

140

PRO X-ray Powder Diffractometer (Eindhoven, Netherlands), using Ni-filtered Cu Kα

141

radiation (λ = 1.5406 Å). The data were recorded over a scanning 2θ range of 2° to 50° at

142

step time of 0.045 steps/0.5 s.

143

2.4.3 FTIR

144

IR discs were prepared by using 100 mg potassium bromide and 2 mg of samples with

145

blending. The blend was subjected to compression under vacuum at a pressure of 12 ρsi for 3

Characterization

ACS Paragon Plus Environment

Page 6 of 33

Page 7 of 33 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

146

min. FTIR was performed using a PerkinElmer IR spectrophotometer with a suitable holder

147

for the disc. FTIR spectra were obtained in the range of 4000 to 400 cm−1.

148

2.4.4 Thermogravimetric Analysis (TGA)

149

TGA of riluzole-2,4 dihydroxy benzoic acid methanol solvate was carried out using ExStar

150

TGA/DTA 7200 instrument operating with Muse software. Sample was loaded in alumina

151

crucibles and heated at a rate of 20 °C/min over a temperature range of 30 to 300 °C under a

152

nitrogen purge of 60 mL/min.

153

2.4.5 Single Crystal Analysis

154

The intensity data for riluzole-ferulic acid was collected at room temperature using a Bruker

155

Smart Apex CCD diffractometer with graphite monochromated MoKα radiation

156

(λ=0.71073Å) by the ω-scan method. Preliminary lattice parameters and orientation matrices

157

were obtained from four sets of frames. Integration and scaling of intensity data were

158

accomplished using the program SAINT.28 The structures were solved by direct methods

159

using SHELXS and refinement was carried out by full-matrix least-squares technique using

160

SHELXL.29 Anisotropic displacement parameters were calculated for all non-hydrogen

161

atoms. X-ray data of riluzole-vanillic acid, riluzole-cinnamic acid, riluzole-proline, riluzole-

162

2,4 dihydroxy benzoic acid methanolate and riluzole-fumaric acid were collected at room

163

temperature on a Bruker D8 QUEST instrument with an IµS Mo micro source (λ = 0.7107 A)

164

and a PHOTON-100 detector. The raw data frames were reduced and corrected for absorption

165

effects using the Bruker Apex 3 software suite programs. The structure was solved using

166

intrinsic phasing method and further refined with the SHELXL program and expanded using

167

Fourier techniques. The H atoms bound to the N and O atoms were located in difference

168

Fourier maps, and their positions and isotropic displacement parameters were refined. C-

169

bound H atoms were located in a difference density map but were positioned geometrically

170

and included as riding atoms, with C—H distance = 0.93 -0.96 Å and with Uiso(H) values of

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

171

1.2Ueq(C). Atoms O1 and F1 of riluzole-2,4 dihydroxy benzoic acid methanolate were

172

disordered and the site-occupancy factors of the disordered atoms O1/O1D/F1/F1D were

173

refined to 0.688(7) and 0.312(8). The anisotropic displacement parameters for the disordered

174

atoms were restrained to be similar. The C1-O1 distance was constrained to 1.32(1) Å. The H

175

atom bound to the O atom of the fumaric acid of riluzole-fumaric acid during the initial

176

refinement shows high Uiso(H) value. Contoured difference Fourier maps show significant

177

electron density at the location of the potential H-atom site and the electron density is

178

smeared along the O···N axis in riluzole-fumaric acid. However, the electron density has a

179

maximum at atom O. Hence the H atoms were treated as disorder over two sites (H1O and

180

H1N) and refined with their site-occupancy factors to 0.67 (3) and 0.33 (3). The disordered H

181

atoms were located in difference Fourier and refined as riding, with Uiso(H) =

182

1.2Ueq(N)/1.5Ueq(O), using constrained distance of 0.82 and 0.86 Å for O—H and N—H,

183

respectively.

184

2.4.6 In vitro Dissolution Study

185

Dissolution study was performed for the plain riluzole and novel solid forms of riluzole. USP

186

apparatus II (paddle type) was selected for the studies. The dissolution studies were

187

performed in 1000 ml of media maintained at 37± 0.5 °C, at 50 rpm. Phosphate buffer (pH

188

6.8) was used as a dissolution media.20 Plain drug and all novel solid forms equivalent to 20

189

mg of drug were added to the dissolution media. All samples were passed through BSS

190

#sieve no. 44 and retained on BSS #sieve no. 60 before addition of dissolution media.

191

Aliquots of 5 ml were withdrawn at predetermined time points (5, 10, 15, 30, 45, 60, 120,

192

180, and 240 min) substituting the same with equal quantity of fresh dissolution media.

193

Samples were filtered using 0.22 µ filters and analyzed by a validated HPLC method after

194

suitable dilution with the mobile phase. Dissolution profiles of plain drug and novel solid

ACS Paragon Plus Environment

Page 8 of 33

Page 9 of 33 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

195

forms were evaluated for the amount of drug released at 15 min (Q15), % dissolution

196

efficiency (DE) at 15 min and similarity factor using DDsolver.30

197

2.4.7 HPLC Method

198

Drug was quantified using a HPLC system (e2695 Waters) consisting of a HPLC pump, an

199

automated injector equipped with a UV detector (2998 PDA) and an auto sampler. Formic

200

acid (0.005 % v/v) in water and acetonitrile were used as mobile phase at a flow rate of 1

201

ml/min in an isocratic mode on an Inertsustain C18 (5µ, 4.6 x 100 mm) column. After

202

suitable dilution, 20 µL of sample was injected, and the absorbance of elute was recorded at

203

222 nm. Details of the HPLC method used are given in supporting information table S2.

204

2.4.8 In vivo hepatotoxicity study

205

In vivo hepatotoxicity study for plain drug and novel solid forms containing hepatoprotective

206

coformers were carried out in 6-7 weeks old, male BALB/c mice weighing 20-25 gm.

207

Animals were obtained from the central animal facility, National Institute of Pharmaceutical

208

Education and Research (NIPER), Hyderabad, India. Animals were housed at 22 ± 2 °C, 50-

209

60 % relative humidity (RH), under 12 hr natural light/dark cycle. Standard pellet diet and

210

water was given ad libitum. The study protocol was duly approved by the Institutional

211

Animal Ethics Committee (IAEC), NIPER, Hyderabad, India (NIP/12/2016/PE/220). The

212

animals were divided into six groups (control group, riluzole treated group, ferulic acid

213

treated group, syringic acid treated group, riluzole-ferulic acid cocrystal treated group and

214

riluzole-syringic acid cocrystal treated group). All treatment doses were prepared in 0.5 %

215

sodium carboxy methyl cellulose to form a suspension. Plain drug riluzole and cocrystals

216

(riluzole-ferulic acid and riluzole-syringic acid were administered orally to animals in a dose

217

equivalent to 25 mg/kg of body weight. While, ferulic acid and syringic acid were given at an

218

equivalent dose of ferulic acid and syringic acid present in riluzole-ferulic acid and riluzole-

219

syringic acid cocrystal respectively. The blood samples were collected at 24 hr after dosing

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 10 of 33

220

via retro orbital route in vials prefilled with heparin solution. The blood samples were

221

centrifuged at 3000 rpm for 10 min, plasma was collected and analysed for enzyme levels.

222

The enzyme levels of aspartate transaminase (AST), alanine transaminase (ALT) and alkaline

223

phosphatase (ALP) in plasma were estimated using commercial kits (Accurex, India) as per

224

manufacturer’s protocol.

225

2.4.9 Statistical analysis

226

The data obtained in animal study were expressed as the mean of replicate determinations

227

(n=5) ± standard error of mean (S.E.M). Statistical comparisons were made using one-way

228

analysis of variance (ANOVA). The intergroup variations were analyzed by “Tukey’s

229

multiple comparison test” using the Graph Pad Prism software, Version 5.0.

230

3

231

With its unique fluorine containing chemical structure and easily accessible functional

232

groups, riluzole is found to have propensities to form salt or cocrystals through hydrogen or

233

ionic interactions. Hence, to screen the appropriate coformers for generation of novel

234

cocrystal or salt of riluzole, prediction studies were carried out using Material Studio

235

software. Outcome of prediction study are depicted in supporting information table S3, along

236

with the predicted bond length and rank orders. From the

237

concluded that some of the coformers such as ferulic acid, fumaric acid, malonic acid,

238

succinic acid, ascorbic acid, azelic acid, cinnamic acid and 2,4 dihydroxybenzoic acid form

239

intermolecular (amine- carboxylic acid) hydrogen bond with riluzole (figure 1 and table S3

240

supporting information). Other coformers namely glutaric acid and sebacic acid showed the

241

presence of intermolecular interactions between carboxylic acid-halogen. In contrast,

242

phenolic –OH of adipic acid showed molecular interactions with the sulphur of riluzole.

243

Experimental outcome of the screening experiments of this current study with 24 coformers

244

and 9 coformers (adipic acid, azelaic acid, 4-dimethylaminopyridine, glutaric acid, maleic

Result and discussion

ACS Paragon Plus Environment

prediction study, it can be

Page 11 of 33 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

245

acid, nicotinic acid, sorbic acid, suberic acid and succinic acid) already reported for this drug

246

by Mondal P et al,19 was compared with the prediction outcome (table 1). It was observed

247

that out of 33 coformers reported in table 1, prediction outcome showed similarity with the

248

outcome of crystallization experiments of riluzole with the 25 coformers. Coformers in the

249

same analogue series have not shown similar pattern for cocrystal formation in the prediction

250

study. Ferulic and cinnamic acids have given cocrystals but not related caffeic and p-

251

coumaric acids. This could be attributed to the fact that this prediction model was based on

252

molecular dynamics, wherein the orientation of the drug and coformers during amorphous

253

cell construction may vary from co-former to co-former despite their analogous similarity.

254

However, it is interesting to note that the prediction behaviour matched with the experimental

255

behaviour in these molecules demonstrating the potential utility of this model. Additionally,

256

single crystal analysis of few novel solid forms of riluzole with (ferulic acid, vanillic acid,

257

fumaric acid, cinnamic acid and 2,4-dihyroxybenzoic acid) depicted similar type of

258

intermolecular interactions between drug and coformers to that of predicted results for

259

molecular interactions.

260

Furthermore, the difference in pKa of acid and base (∆pKa) was used to predict the formation

261

of cocrystal or salt. ∆pKa < 0 often gives a cocrystal, ∆pKa > 3 may form a salt and ∆pKa in

262

the range of 0 to 3 may result in salt-cocrystal hybrids31-34. ∆pKa value for the screened drug-

263

coformer systems were calculated and depicted in table 1, and S4 supporting information,

264

wherein novel solid systems like riluzole-ferulic acid, riluzole-syringic acid, riluzole-vanillic

265

acid, riluzole-cinnamic acid and riluzole-proline showed negative value of ∆pKa which

266

correspond to the cocrystal formation, while ∆pKa values for riluzole-2,4-dihydroxybenzoic

267

acid methanolate and riluzole-fumaric acid were found to be 0.69 and 0.77, respectively,

268

which led to the assumption that riluzole-2,4-dihydroxybenzoic acid methanolate and

269

riluzole-fumaric acid may form cocrystal salt continuum zone.

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

270 271

Page 12 of 33

Table 1 Comparative analysis of prediction and experimental outcome (# cocrystals and salts reported by Mondal P et al.)19 Sr. No.

System

1 2 3 4

Riluzole-adipic acid# Riluzole- ascorbic acid Riluzole-azelaic acid# Riluzole-benzoic acid

0.63,-1.6 -0.37 -0.75 -0.40

5 6 7 8 9

-0.82 -3.0 -0.66 1.01 0.69

11 12

Riluzole-caffeic acid Riluzole-chrysin Riluzole-cinnamic acid Riluzole-citric acid Riluzole2,4dihydroxybenzoic acid Riluzole- 4 dimethylaminopyridine# Riluzole-ferulic acid Riluzole-fumaric acid

13 14 15 16 17 18

Riluzole-gallic acid Riluzole-glutaric acid# Riluzole-hesperetin Riluzole-maleic acid Riluzole-malic acid Riluzole-malonic acid#

19

Riluzole-nicotinic acid#

20 21

Riluzole-p-coumaric acid Riluzole-proline

22

Riluzole-quercetin

23 24 25 26 27 28

Riluzole-salicylamide Riluzole-salicylic acid Riluzole-sebacic acid Riluzole-silibinine Riluzole-sorbic acid# Riluzole-suberic acid#

29

Riluzole-succinic acid#

30 31

Riluzole-syringic acid Riluzole-tartaric acid

32 33

Riluzole-vanillic acid Riluzole-vanillin

10

∆pKa

-5.9 -0.78 0.77, 0.64 -0.64 -0.54 -4.12 1.97,-2.4 0.4,-1.3 0.97,1.89 -0.95 -0.2 1.86,0.74 -2.07,4.68 -4.57 0.83 -0.92 -3.95 -0.96 -0.72 -0.46,1.92 -0.13 0.82,0.54 -0.71 -3.58

Prediction model intermolecular interactions Intermolecular Intermolecular Intermolecular Absence of interaction Intramolecular Intramolecular Intermolecular Intramolecular Intermolecular

Cocrystal or salt formation

Melting point (°C)

Predicted outcome

Experimental outcome

Yes Yes Yes No

Yes No Yes No

115 97 110 101

No No Yes No Yes

No No Yes No Yes

102 110 145 72 157

Absence of interaction Intermolecular Intermolecular

No

Yes

96

Yes Yes

Yes Yes

159 163

Intramolecular Intermolecular Intramolecular Intramolecular Intramolecular Intermolecular

No Yes No No No Yes

No Yes No Yes No Yes

86 120 105 153 78 139

Absence of interaction Intramolecular

No

Yes

118

No

No

92

Absence of interaction Intramolecular

No

Yes

206

No

No

101

Intramolecular Intramolecular Intermolecular Intramolecular Intramolecular Absence of interaction Intermolecular

No No Yes No No No

No No Yes No Yes Yes

130 110 125 109 91 136

Yes

Yes

133

Intermolecular Intermolecular

Yes Yes

Yes Yes

142 158

Intramolecular Intramolecular

No No

Yes No

139 78

272

ACS Paragon Plus Environment

Page 13 of 33 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

273 274 275 276 277 278

Crystal Growth & Design

Figure 1 Prediction study outcome: (a) riluzole-adipic acid, (b) riluzole-ascorbic acid, (c) riluzoleazelic acid, (d) riluzole-benzoic acid, (e) riluzole-caffeic acid, (f) riluzole-chrysin, (g) riluzolecinnamic acid, (h) riluzole-citric acid, (i) riluzole-p-coumaric acid, (j) riluzole-2,4-dihydroxybenzoic acid, (k) riluzole-2,4 -dimethylaminopyridine, (l) riluzole-ferulic acid, (m) riluzole-fumaric acid, (n) riluzole-gallic acid, (o) riluzole-glutaric acid, (p) riluzole-hesperetin, (q) riluzole-maleic acid, (r)

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

279 280 281 282

riluzole-malic acid, (s) riluzole-malonic acid, (t) riluzole-nicotinic acid, (u) riluzole-proline, (v) riluzole-quercetin, (w) riluzole-salicylic acid, (x) riluzole-salicylamide, (y) riluzole-sebacic acid, (z) riluzole-silibinin, (aa) riluzole-sorbic acid, (ab) riluzole-suberic acid, (ac) riluzole-succinic acid, (ad) riluzole-syringic acid, (ae) riluzole- tartaric acid, (af) riluzole-vanillic acid and (ag) riluzole-vanillin

283

3.1

284

Single endotherm at 119.2 °C corresponding to its melting appeared in the DSC curve of

285

riluzole (figure 2). The melting endotherm for coformers such as ferulic acid, syringic acid,

286

fumaric acid, cinnamic acid, 2,4 dihydroxybenzoic acid, proline, and vanillic acid were

287

observed at 172.8, 208.4, 207.1, 133.8, 194.9, 219.4 and 211.5 °C respectively. The

288

generated novel solid forms such as riluzole-ferulic acid, riluzole-syringic acid, riluzole-

289

fumaric acid, riluzole-cinnamic acid, riluzole-2,4-dihydroxybenzoic acid methanolate,

290

riluzole-proline and riluzole-vanillic acid gave single melting endotherms at 159.4, 142.4,

291

163.1, 145.5, 157.3, 206.9 and 139.7 °C respectively (figure 2). Riluzole-2,4-

292

dihydroxybenzoic acid methanolate gave a small endotherm at 86 °C indicating a presence

293

of solvate form, which complemented the single crystal structural analysis. The onset of

294

melting and ∆H values for all novel forms along with coformers, and respective thermograms

295

along with controls are given in supporting information (table S5). Changes in melting point

296

and enthalpy indicated the formation of new solid forms, which were additionally confirmed

297

by P-XRD studies.

DSC

ACS Paragon Plus Environment

Page 14 of 33

Page 15 of 33 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

298 299 300 301 302

Figure 2 Overlay of DSC thermogram of R- riluzole, RFRA-riluzole-ferulic acid, RSYRA-riluzolesyringic acid, RVNLA– riluzole-vanillic acid, RCA-riluzole-cinnamic acid, RPRO-riluzole-proline, RDHBA methanolate riluzole-2, 4-dihydroxybenzoic acid methanolate and RFUM-riluzole-fumaric acid

303

3.2

304

FTIR analysis of plain riluzole and novel solid forms were performed to detect the presence

305

of hydrogen bond between drug and coformers. Plain riluzole showed characteristic N-H

306

stretching band at 3371.5, 3279.4 cm−1 which confirms the presence of primary amine group

307

in riluzole. Apart from this, presence of peak at 1606.6 and 1327.1 cm−1 in the spectra

308

correspond to N-H bending and C-N stretching, respectively. O-H and -C=O stretching

309

corresponding to the carboxylic acids appear in the region of 3400-2400 cm−1 and 1720-1680

310

cm−1, respectively (figure S4, supporting information). It is expected that due to hydrogen

311

bond formation in novel solid forms, there will be shift in these stretching peaks. Specifically,

312

-N-H stretching and bending peaks at 3279.0 cm−1 and 1606.5 cm−1 showed significant shift

313

in IR spectra of novel solid forms. Significant shift in the characteristic peaks of riluzole

314

confirms the formation of novel solid forms. To differentiate formation of cocrystal and salt

FTIR

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

315

by FTIR, presence or absence of carboxylic acid O-H stretching peak (2500-3200 cm−1) in

316

FTIR spectra of solid forms is considered as an indicator.19 As highlighted in the figure S4

317

supporting information, generated cocrystals showed carboxylic acid -O-H stretching peak

318

that confirmed the presence of carboxylic acid group, while this peak was absent in fumaric

319

acid and DHBA based solid forms confirming formation of salt for these two coformers.

320

3.3

321

As shown in figure 3, characteristic peaks of riluzole-ferulic acid, riluzole-syringic acid,

322

riluzole-vanillic acid, riluzole-cinnamic acid, riluzole-proline, riluzole-2, 4-dihydroxybenzoic

323

acid methanolate and riluzole-fumaric acid in PXRD patterns differ from riluzole and its

324

corresponding coformer. Riluzole has shown characteristic peaks at 2θ values 9.0, 13.5, 18.1

325

and 26.4 whereas new characteristic peaks appeared in riluzole-ferulic acid (13.8, 14.9, 27.3

326

and 27.9); riluzole-syringic acid (13.2, 14.9, 18.4 and 26.2), riluzole-vanillic acid (14.3, 15.2,

327

20.3, 26.5 and 26.9), riluzole-cinnamic acid (11.7, 18.8, 23.9 and 27.3), riluzole-proline

328

(11.1, 16.7, 20.9 and 22.5), riluzole-2,4 dihydroxybenzoic acid (8.1, 24.8, 26.9, and 27.5) and

329

riluzole-fumaric acid (12.2, 22.6, 25.0 and 26.9). The appearance of new peaks in PXRD

330

pattern when compared with R and respective coformers, collaborate the generation of novel

331

solid forms. Bulk purity of all generated system was confirmed by superimposable

332

experimental and simulated diffraction patterns overlay (supporting information figure S5).

333

Overall, outcome of DSC, PXRD, FTIR and single crystal XRD confirmed the occurrence of

334

the new solid forms.

335

Plain salicylamide and riluzole-salicylamide combination from liquid assisted grinding

336

showed a melting point at 141.0 and 138.7 ⁰C, respectively in DSC (data not shown).

337

Although the riluzole-salicylamide combination gave a melting point which was intermediate

338

to the starting components, the PXRD result of the combination showed overlapping peaks as

PXRD

ACS Paragon Plus Environment

Page 16 of 33

Page 17 of 33 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

339

that of plain drug and salicylamide indicating no cocrystal formation (figure S6 supporting

340

information).

341 342 343 344 345 346 347

Figure 3 Overlay of diffractograms of R-riluzole, FRA- ferulic acid, RFRA-riluzole-ferulic acid, VNLA-vanillic acid, RVNLA-riluzole-vanillic acid, SYRA-syringic acid and RSYRA-riluzolesyringic acid, CA- cinnamic acid, RCA- riluzole-cinnamic acid, PRO-proline, RPRO-riluzole-proline, DHBA-2,4-dihydroxybenzoic acid, RDHBA methanolate- riluzole-2,4-dihydroxybenzoic acid methanolate, FUM-fumaric acid and RFUM-riluzole-fumaric acid. (* indicates new peaks in solid forms).

348

3.4 TGA

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

349

Riluzole-2, 4-dihydroxybenzoic acid methanolate showed a weight loss of 3.248 % (0.0819

350

mg) in the region of 86-100 °C in TGA thermogram (figure S2, supporting information). This

351

finding is in line with the DSC outcome and the single crystal structural analysis. The

352

outcome hence confirms the presence of solvate in novel solid form, riluzole-2,4-

353

dihydroxybenzoic acid methanolate. However, the discrepancy in weight loss during TGA

354

analysis from the theoretical mass loss of one molecule of methanol might have occurred due

355

to slipping of methanol from the crystal lattice at lower temperature.35 This was confirmed by

356

integration of results from 25 ⁰C up to 100 and 154 ⁰C, which showed enhancement in

357

percentage weight loss (figure S2 and S3, supporting information).

358

3.5

359

All six novel solid forms (riluzole-ferulic acid, riluzole-vanillic acid, riluzole-cinnamic acid,

360

riluzole-proline, riluzole-2,4-dihydroxybenzoic acid methanolate and riluzole-fumaric acid)

361

crystallized in monoclinic system. The state of the carboxylic acid group (whether neutral or

362

anionic) can be determined from the C-O and C = O bond lengths. For a neutral carboxylic

363

acid group, the C-O and C=O distances are around 1.2 and 1.3A˚, respectively. However,

364

when deprotonation occurs, both C-O bond lengths will be similar, around 1.25 A˚.36 Based

365

on the bond distances listed in table 3, riluzole-ferulic acid, riluzole-vanillic acid and riluzole-

366

cinnamic acid exist as cocrystals. In case of riluzole-proline, internal proton transfer takes

367

place from carboxylic acid to the N atom of the pyrrolidine ring and exists as zwitter ion. As

368

a result, the C-O distances of the carboxylic group of proline are similar. However, riluzole-

369

proline as a whole exists as zwitterionic cocrystal.37 The near equality of C-O distances of

370

carboxylic group of riluzole-2,4-dihydroxybenzoic acid methanolate, show deprotonation

371

from O to N and confirm the presence of salt. However, in riluzole-fumaric acid, the C-O

372

bond distances are intermediate between the neutral and ionic because of partial proton

373

transfer between O to N. The geometry of bond angles of carboxylic acid group also show

Single crystal structure analysis

ACS Paragon Plus Environment

Page 18 of 33

Page 19 of 33 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

374

similar trend for neutral, ionic and intermediate state. The crystal structures of riluzole-ferulic

375

acid,

376

dihydroxybenzoic acid methanolate and riluzole-fumaric acid are mainly determined by

377

N−H···O and O−H···O hydrogen bonds. Hydrogen bond geometry of crystal had been

378

explained in table 4 for all systems.

riluzole-vanillic

acid,

riluzole-cinnamic

acid,

riluzole-proline,

riluzole-2,4-

Table 2 Crystallographic parameter

379 Identification code

Riluzoleferulic acid

Riluzolevanillic acid

Riluzole cinnamic acid

Riluzole-proline

Riluzolefumaric acid

1557240

Riluzole-2,4dihydroxybenzo ic acid methanolate 1557241

CCDC deposition no. Empirical formula Formula weight

1538080

1557243

1557244

C8H5F3N2OS, C10H10O4 428.38

C8H5F3N2OS,C8 H8O4 402.34

C8H5F3N2OS+, C9H8O2382.35

C8H5F3N2OS C5H9NO2 349.33

C8H6F3N2OS+ C7H5O4-,CH4O 420.36

C8H5.33F3N2OS,0 .5(C4H3.34O4) 292.24

Temperature K Wavelength Å

294(2) 0.71073

294(2) 0.71073

294(2) 0.71073

294(2) 0.71073

294(2) 0.71073

294(2) 0.71073

Crystal system Space group a (Å) b (Å) c (Å) α (deg) β (deg) γ (deg) Volume Å3

Monoclinic P21/c 6.366 (8) 7.5981(9) 38.608 (5) 90 93.901(2) 90 1863.1(4)

Monoclinic P21/n 6.655 (6) 35.02 (3) 7.437 (8) 90 95.61(2) 90 1725 (3)

Monoclinic P21/c 7.602 (4) 6.397 (3) 35.158 (18) 90 90.946 (14) 90 1709.5 (14)

Monoclinic P21 8.457 (14) 5.425 (9) 16.02 (3) 90 97.54 (2) 90 729 (2)

Monoclinic P21/n 13.874 (13) 8.694 (8) 15.219 (12) 90 97.102 (13) 90 1822 (3)

Monoclinic P21/C 8.140 (2) 15.658 (3) 9.508 (2) 90 94.46 (3) 90 1208.1 (5)

1557242

Z

4

4

4

2

4

4

Density (calculated) Mg/m3

1.527

1.549

1.486

1.592

1.533

1.607

Absorption

0.237

0.250

0.240

0.275

0.245

0.312

coefficient mm-1 F(000)

880

824

784

360

864

592

Crystal size mm

0.480x0.320x

0.420x0.380x0.2 80

0.380x0.360x0.2 80

0.400x0.360x0.2 20

0.420x0.180x0.1 60

0.360x0.220x0.1 80

0.270 θ range

1.05-28.364

2.36-28.345

2.317-28.356

2.429-28.38

2.697-28.470

2.510-28.382

Range h Range k Range l Reflections collected Observed reflections Data / restraints / parameters Goodness-of-fit on F2

-8 to 8 -10 to 10 -51 to 49 20608

-8 to 8 -46 to 46 -9 to 9 43639

-9 to 10 -8 to 8 -46 to 46 17852

-10 to 11 -7 to 7 -21 to 21 17399

-18 to 18 -11 to 11 -20 to 19 56736

-10 to 10 -20 to 20 -12 to 11 28755

4514

4280

4256

3644

4596

3022

4514/0/279

4280/0/262

4256/0/248

3644/1/224

4596/81/297

3022/0/183

1.087

1.150

1.075

1.076

1.051

1.029

Final R indices [I>2σ(I)] R1, wR2 R indices (all data) R1, wR2

0.0547, 0.1395

0.0623, 0.1661

0.0511, 0.1389

0.0394, 0.0900

0.0705, 0.2034

0.0381, 0.1040

0.0600, 0.1436

0.0722, 0.1723

0.0562, 0.1427

0.0497, 0.0942

0.0845, 0.2197

0.0420, 0.1088

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 33

380 Table 3 C-O bond distances and angles value of novel solid forms

381

Compounds Riluzole-ferulic acid Riluzole-vanillic acid Riluzole-cinnamic acid Riluzole-proline Riluzole-2,4-dihydroxybenzoic acid methanolate Riluzole-fumaric acid

Bond distances (Ǻ) 1.218 (2) 1.308 (2) 1.219 (2) 1.311 (2) 1.219 (2) 1.305 (2) 1.238 (4) 1.243 (4) 1.255 (3) 1.259 (3)

Bond angles (°) 122.90 (18) 114.03 (17) 122.50 (2) 114.40 (2) 122.91 (17) 113.89 (17) 118.3 (3) 116.2 (2) 119.3 (2) 117.9 (2)

1.2228 (16)

120.43 (12)

1.2911 (16)

114.88 (11)

Table 4 Hydrogen bond geometry in crystal structures

382 Compounds Riluzoleferulic acid

Riluzolevanillic acid

Riluzolecinnamic acid Riluzoleproline

Riluzole2,4dihydroxy benzoic acid methanolate

Interaction N(2)-H(1N)···O(1)#1 N(2)-H(2N)···O(1) O(3)-H(3O)···O(4) O(2)-H(2O)···N(1) N(2)-H(1N)···O(3) N(2)-H(2N)···O(3)#1 O(2)-H(1O)···N(1) O(4)-H(2O)···F(2)#2 O(4)-H(2O)···O(5) N(2)-H(1N)···O(3) N(2)-H(2N)···O(3)#1 O(2)-H(1O)···N(1) N(3)-H(1N)···N(2)#1

d(D-H) 0.88(3) 0.90(3) 0.82(3) 0.96(4) 0.85(4) 0.88(4) 0.88(4) 0.74(4) 0.74(4) 0.93(3) 0.85(3) 0.92(3) 0.98(5)

d(H···A) 2.11(3) 1.95(3) 2.09(3) 1.74(3) 2.07(4) 2.18(4) 1.80(4) 2.50(4) 2.13(4) 1.92(3) 2.16(3) 1.75(3) 2.59(4)

d(D···A) 2.925(3) 2.840(3) 2.611(3) 2.690(2) 2.908(4) 2.988(3) 2.671(3) 3.126(4) 2.602(3) 2.846(3) 2.921(3) 2.661(2) 3.163(6)