Distinguishing of cocrystals from simple eutectic mixtures: phenolic

reduce experimental effort by narrowing down the list of potential coformers. Based on different methodologies .... software.46 The values of all 1444...
0 downloads 3 Views 1MB Size
Subscriber access provided by UNIV OF NEW ENGLAND ARMIDALE

Distinguishing of cocrystals from simple eutectic mixtures: phenolic acids as potential pharmaceutical coformers Maciej Przyby#ek, and Piotr Cysewski Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.8b00335 • Publication Date (Web): 23 Apr 2018 Downloaded from http://pubs.acs.org on April 24, 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 39 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

Distinguishing of cocrystals from simple eutectic mixtures: phenolic acids as potential pharmaceutical coformers Maciej Przybyłek* and Piotr Cysewski

Chair and Department of Physical Chemistry, Faculty of Pharmacy, Collegium Medicum of Bydgoszcz, Nicolaus Copernicus University in Toruń, Kurpińskiego 5, 85-950 Bydgoszcz, Poland

ABSTRACT The multiparameter model comprising 1D and 2D QSPR/QSAR descriptors was proposed and validated for phenolic acids binary systems. This approach is based on the optimization of regression coefficients for maximization of true positives percentage in the pool of systems comprising either simple binary eutectics or cocrystals. The training set consisted of 58 eutectics and 168 cocrystals. The solid dispersions collection used for models generation comprised literature data enriched with our new experimental results. From all 1445 descriptors computable in PaDEL only 13 orthogonal descriptors with the highest predicting power were taken into account. The analysis revealed importance of the parameters characterizing atom types (naaN, SHsOH, SsssN, nHeteroRing, maxHBint6, C1SP2), autocorrelation functions (ATSC1i, AATSC1v, MATS8m, GATS1i) and also other molecule structure measures (WTPT-5, MLFER_A, MDEN-22). The proposed approach is very simple and requires only information about the structure encoded in canonical SMILES string. The inverting of the problem of cocrystals screening and focusing on the homogeneous group of 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 2 of 39

coformers for cocrystallization with variety of drugs rather than seeking coformers for particular active pharmaceutical ingredient proved to be very efficient. This led to very valuable clues for selection of pairs for cocrystallization with probability about 80%.

Keywords: Pharmaceutical cocrystals, Phenolic acids, Cocrystal screening, Molecular descriptors, Nutraceutical cocrystals

Introduction It is well known that the pharmacokinetic properties of drugs are very closely related to the crystal form. Recently, extensive studies have been conducted to improve active pharmaceutical ingredients’ (APIs) relevant properties via crystal engineering including polymorphism and crystal morphology control1,2 and cocrystallization.3–6 Cocrystals deserve particular attention, due to the wide range of potential coformers, which can be utilized for tuning the physicochemical properties of APIs. Apart from bioavailability improvement and related aspects, there are several other features important from the pharmaceutical viewpoint which can be modified through cocrystallization, such as mechanical properties and stability.7,8 The latter feature seems to be interesting in the context of preserving properties of certain phytochemicals such as phenolic acids. These compounds are typical strong antioxidants, since they exhibit high radical scavenging activity,9–11 lipid peroxidation diminishing properties5,9,12 and endogenous defense systems activation ability.9,13 Phenolic acids can be divided into two major subclasses, namely hydroxybenzoic acids and hydroxycinnamic acids. Compounds belonging to the former subclass can be found in many fruits, teas, and vegetables.12,13 The latter subclass, cinnamic acid derivatives, are stronger antioxidants than hydroxylated benzoic acid analogues.14 Similarly to other phenolic acids,

2 ACS Paragon Plus Environment

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

they are also effective anti-inflammatory agents.15 Due to their properties, these nutraceuticals exhibit promising therapeutic effects for the treatment of diabetes and obesity.15 Phenolic acids cocrystals are interesting from both structural and practical aspects. Due to the hydroxyl group attached to the benzene ring hydroxybenzoic acids can form various intermolecular interactions. This diversity is manifested by the variety of supramolecular synthons or even synthon polymorphism.16–18 From the pharmaceutical viewpoint APIs cocrystallization with phenolic acids seems to be advantageous for several reasons. First, it has been demonstrated that these coformers are solubility and dissolution rate enhancers.19–23 Second, recent studies7 showed that phenolic acids cocrystals exhibit antimicrobial and anticancer activity improvement, which is probably caused by their properties (high antioxidant potential, antifungal and antibacterial activity). This shows that phenolic acids can not only play the role of excipient, but can also act as active ingredients in drug-drug cocrystals. Third, it has been demonstrated that phenolic acids cocrystals are more stable than pure API. This is understandable, since these compounds are commonly used as preservatives in food industry, pharmacy and cosmetics.24 Some interesting examples of stability improvement via cocrystallization are molecular complexes of popular antituberculosis agent, isoniazid with ferulic and caffeic acids which exhibited significant stability under accelerated ICH conditions (relative humidity of 75% and 40°C).18 Cocrystals screening is one of many important strategies for seeking new forms of drugs.25 Since numerous documented examples of cocrystals and quite a lot of simple eutectics can be found in the literature, it is possible to create theoretical models that allow to reduce experimental effort by narrowing down the list of potential coformers. Based on different methodologies, several approaches enabling new multicomponent crystals designing have been raised including hypothetical cocrystal lattice energy calculation,26 Hansen solubility parameters analysis,27,28 Conductor like Screening Model for Real Solvents

3 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 39

(COSMO-RS) based enthalpy of mixing calculation29–34 and various theoretical models utilizing structural parameters for thermodynamic characteristics of solids including melting point prediction.35–39 The last mentioned approaches deserves special attention due to the application of molecular descriptors, which nowadays can be easily calculated using many available on the market or freeware software. The main goal of this study is to develop effective and accurate phenolic acids cocrystals screening model based on the common 1D and 2D QSPR/QSAR descriptors. Although these parameters are very fast to compute, they give simplified characteristics of molecules’ structural features compared to the quantumchemical indices. Nevertheless, constitutional and topological descriptors were utilized in many QSPR/QSAR models and their significant contribution in predicting physicochemical properties and biological activity was confirmed.40 Our recent studies showed that hydrophilicity indices can be successfully used for classification of dissolution rate enhancing/reducing cocrystal formers.31 Therefore, theoretical predictive models based on the 1D and 2D molecular descriptors are promising tools for designing new drug formulations with improved properties. In this paper a novel approach of linear combination of descriptors screening model will be applied.

Materials and Methods Chemicals All chemicals used in this study were purchased at the highest available purity from commercial suppliers and used as received. Amides, subjected to the cocrystallization screening tests, namely urea (CAS: 57-13-6), benzamide (CAS: 55-21-0), salicylamide (CAS: 65-45-2) and ethenzamide (CAS: 938-73-8) and the following phenolic acids cocrystal formers:

2,6-dihydroxybenzoic

acid

(γ-resorcylic

acid,

CAS:

303-07-1),

3,4-

dihydroxybenzoic acid (protocatechuic acid, CAS: 99-50-3), 3,5-dihydroxybenzoic acid (α-

4 ACS Paragon Plus Environment

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

resorcylic acid, CAS: 99-10-5), ferulic acid (CAS: 1135-24-6) were obtained from SigmaAldrich. Methanol used as solvent was purchased from Avantor (Gliwice, Poland).

Experimental cocrystal screening procedure During the experimental stage of this study, solid dispersions of different compounds were prepared and analyzed in terms of new cocrystal phase formation. All measurements were performed using deposited on glass slides crystallite layers obtained via droplet evaporative crystallization (DEC) technique.32,33,41–43 As it was reported in the previous works32,33,43 this method of sample preparation can be successfully used for fast and efficient cocrystal screening. The procedure involved mixing of pure components methanolic solutions (0.3 M) in 1:1 proportion and evaporation of 20 µL droplets of these mixtures on glass microscope slides at 43°C. Crystallite deposits prepared in such way were measured using powder X-ray diffractometry (PXRD) and attenuated total reflection-Fourier-transform infrared spectroscopy (ATR-FTIR). Powder patterns were recorded by applying PW3050/60 goniometer and Empyrean XRD tube Cu LFF DK303072 (2θ angle range: 5° to 40°, width step: 0.02°). Then, the obtained raw diffractograms were cleaned from Kα2 contribution, background subtracted, smoothed and normalized. ATR-FTIR spectra were recorded using Alpha-PFT-IR instrument (Bruker Optik, Ettlingen).

Computational protocol The datasets of cocrystals and eutectics were collected based on the extensive literature studies and Cambridge Structural Database (CSD)44 survey. This collection includes solid dispersions formed by phenolic acids and their precursors (benzoic and cinnamic acid). For each compound canonical SMILES codes were retrieved from the PubChem website45 and were directly used for 1D and 2D molecular descriptors calculation using PaDEL 2.21 software.46 The values of all 1444 computed descriptors were collected and used for 5 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 6 of 39

distinguishing of binary cocrystals from simple eutectics by analyzing the absolute differences between descriptors corresponding to each pair. For this purpose 2×2 contingency statisitcs47 was applied, which relies on classification of true positives (TPs=properly classified experimentally observed pairs forming cocrystals) and true negatives (TNs=properly classified cases of simple binary eutectics or non-successful cocrystallizations). Before computing the absolute differences characterizing each pair the incomplete distributions of molecular descriptors (non-computable for all cases) as well as those comprising more than 80% zeroth values were rejected from further analysis. The remaining data were normalized and ranked in terms of potential of predictive power defined by the highest percentage of true positives TP% with the assumption of zeroth error (FP%=0). Linearly correlated sets (Pearsons’ correlation coefficient R2>=0.6) were rejected, which narrowed further analysis to only 13 orthogonal descriptors. Based on these variables the linear combinations were defined and their coefficients were optimized during regression analysis. Hence, the computational problem was reduced to weighting of descriptors contributions for maximization of TP%. Hence, for each size of the problem the regression function (RF) was defined as follows

∑   = ∑    |∆  | 

(1)

where |∆  | = |  ℎ −   | and where DkI denotes the value of k-th descriptor characterizing i-th compound being either phenolic acid (PhA) or other excipients. Among the latter there are many active pharmaceutical ingredients (API). Next, threshold value defining false positives percentage equal zero (FP%=0) was estimated within 58 eutectics (EU) subset and TP% was computed for remaining 168 cocrystals (CC). The values of regression parameters were optimized in each model for maximizing the value of TP%. Hence, the applied procedure leads to identification of the threshold level guaranteeing the highest TP%

6 ACS Paragon Plus Environment

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

values with as low as possible false positives in the eutectics subset. In the case of models with lower dimensions than 13 the values of corresponding regression coefficients were set to zero. This led to a total of 8178 models resulting from all combinations of 13 parameters examined in this study. The best models were analyzed in details using internal validation procedure. This step, aiming for verification sensitivity of each model on accidental data collection, relied on random exclusion of 10%, 20% and 30% of cases and re-optimization of suitable regression parameters included in the particular model. This procedure was repeated 100 times for each validation step and for every model. The error of each coefficient was calculated based on Student’s t-test approach at 99.9% confidence level.48 As a result, the set of 300 solutions of regression function was used for formulation of parameters and their uncertainties for every model. Finally, the model of the most predictive power of distinguishing cocrystals from simple binary eutectics was identified based on maximum value of TP%-FP% provided 95% confidence of prediction. Hence, models with FP% higher than 5% were excluded as too noisy. This means that in each proposed model the probability of inclusion of false positive cases is less than 5%. Details characterizing the top 20 solutions are provided in Supporting Information in Table S1.

Results and discussion Both experimental and theoretical screening of cocrystallization propensities of phenolic acids with variety of compounds was performed. In the first step the models were generated and validated according to the three-step procedure described in the methodology part. Final part presents the potential range of applications for practical purposes of selecting the most promising candidates for cocrystallization. It is worth mentioning that typical screening procedures are usually performed for seeking suitable coformers for given API. Here the reverse strategy is utilized in which for range of coformers the potential drug-like

7 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 8 of 39

candidates are searched. This approach seems to be more promising since for a homogeneous set of molecules common molecular descriptors are expected to be found. Of course this narrows the range of potential applications but still seems to be very effective for practical screening of new solid forms of drugs.

Dataset characteristics Theoretical model of cocrystal screening was developed utilizing dataset comprising binary solid dispersions formed by phenolic acids with different compounds. We were able to collect 216 systems found in the literature augmented with ten new mixtures prepared and characterized in this study. The detailed characteristic of the whole training set is provided in supplementary as Table S1. It is worth emphasizing that the phenolic acids are quite often used excipients for various classes of APIs. For examples they were used in formulations with non-steroidal anti-inflammatory drugs (e.g. meloxicam49), iron chelators (e.g. deferiprone50), cardiovascular agents (e.g. caffeine51), muscle relaxants (e.g. chlorzoxazone52), anti-microbial agents (e.g. sulfamethazine53), anti-asthmatic agents (e.g. theophylline34) and anti-depressants (e.g. agomelatine54). According to our literature survey, among all mixtures, 58 systems are immiscible in the solid state. This subclass comprises simple eutectics of all considered types of phenolic acids including their precursors, namely benzoic acid26,55–67 and cinnamic acid57,68,69 and their hydroxylated analogues.26,62–66,70–77 Although in simple eutectics counterpart the non-substituted benzoic acids are dominant there are, however, some interesting exclusions. Indeed examples of solid state associates formed by these coformers include many pharmaceutically relevant complexes, such as carbamazepine-cinnamic acid molecular complex78 and benzoic acid cocrystals with many various APIs and drug-like compounds

such

as

aromatic

amides

(Isonicotinamide,79

4-(1H-pyrazol-1-

ylmethyl)benzamide80), amines (1,2-diaminobenzene,81 1,2-diaminobenzene81) diazepams

8 ACS Paragon Plus Environment

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

(nevirapine82),

methylxanthines

(caffeine,83

etofylline84)

and

other

heterocycles

(trimethroprim,84 flucytosine85). The abundant hydroxybenzoic acids subclass contains cocrystals

of

4-hydroxybenzoic

acid,17,86–90

2-hydroxybenzoic

acid,4,27,90–104

3-

hydroxybenzoic acid,105–108 3,5-dihydroxybenzoic acid,107–109 2,4-dihydroxybenzoic acid,51–53 2,5-dihydroxybenzoic acid,110–112 2,6-dihydroxybenzoic acid,113–115 3,4-dihydroxybenzoic acid,33,63,88 4-hydroxy-3-methoxybenzoic acid116–118 and 3,4,5-trihydroxybenzoic acid.119–122 Hydroxycinnamic acids subclass contains binary mixtures formed by 4-hydroxycinnamic acid, ferulic, caffeic and sinapic acids.3,123–125 The dataset used for the model formulation contains also three hydroxynaphthoic acids (1-hydroxy-2-naphthoic acid, 2-hydroxy-1naphthoic acid and 6-hydroxy-2-naphthoic acid) solid dispersions.76,126–130 Although these compounds are not naturally occurring nutraceuticals, they can be assigned to phenolic acids because of their structure. Due to the presence of two condensed benzene rings, these compounds seem to be less attractive as pharmaceutical cocrystal formers than the analogous hydroxybenzoic acids. Nevertheless, hydroxynaphthoic acids were cocrystalized with various APIs, namely sulfamethazine,53 caffeine,127,129 theophylline,34,131 gabapentin,21 meloxicam,49 nicotinamide (vitamin B3),126 pyrazinamide118 and stanozolol.128

Theoretical crystals screening The binomial classification has been already successfully applied30 for crystals screening using one parameter model based on the values of mixing enthalpy. The main difference of the approach documented in this paper is related to the alternate formulation of the target function utilizing the multivariable equation. The maximization of TP% with varying regression parameters allows, in principle, for much more precise tuning of the target function and more accurate separation of the cases belonging to two considered subsets. The whole procedure requires relative values of some descriptors, which are suitable for adequate

9 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 39

classification of pairs. Hence, the main paradigm standing behind the procedure discussed here is the assumption that there are molecular descriptors being able to capture structural, energetic or electronic diversities of components in formulation of either physical mixtures or intermolecular complex in the solid state. It is worth emphasizing that this intuitively appealing assumption may however, not be so exact as it seems from the first glance. For example, lidocaine is able to form 1:1 molecular complex with L-menthol resulting in stable cocrystals in the molten mixture.132,133 However, liquids of the components in unimolar proportions do not systematically recrystallize even by the thaw-melt method.134 This is an important note since DSC experiments provide different thermograms depending on whether equimolar stoichiometric compound of L-menthol and lidocaine is present in solution or not. Hence, in this case cocrystallization requires doping with cocrystals, or else a simple binary eutectic mixture is obtained instead. Fortunately such strange cases are quite rare although the interplay between kinetic and thermodynamic factors must be always considered, if a given pair of compounds is to be classified into cocrystals or eutectics subsets. Keeping this in mind, the training set was carefully inspected and authors believe that each of 216 considered pairs is properly assigned to the either subset. As mentioned in the methodology part, the molecular descriptor selected for the purpose of classification into cocrystallizing and non-crystallizing pairs comes from computable PaDEL descriptors set. It happened that in the cases of considered populations of binary eutectics mixtures and cocrystals the following types of molecular information were found to be important, namely: •

WTPT-5 weighted path, which is the summed path lengths starting from nitrogen atoms135



naaN is one of eccentric connectivity indices,136–138 counting E-State of :NH atom-type

10 ACS Paragon Plus Environment

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



SHsOH, SsssN characterize atom type electrotopological state index values for hydrogen atoms of OH groups and >N- atoms (sp3 nitrogen atoms) respectively136–138



nHeteroRing belongs to ring count subset characterizing the number of hererocyclic rings



ATSC1i, AATSC1v descriptors defined using autocorrelation functions40 using averaged centered Broto-Moreau autocorrelation weighted by first ionization potential or van der Waals volumes, respectively



MATS8m is Moran autocorrelation with lag 8 weighted by mass40



GATS1i is Geary autocorrelation - lag 1 / weighted by first ionization potential40



MLFER_A associated with molecular linear free energy relations139 summarizing solute hydrogen bond acidity



maxHBint6 is atom type electrotopological state Maximum E-State descriptors of strength for potential Hydrogen Bonds of path length 6136–138



MDEN-22 parameters quantifying the molecular distance edges140



C1SP2 counting number of sp2 carbon atom types

It is worth emphasizing that distributions of these values are orthogonal and using all of them in the regression analysis justified. All possible combinations of descriptors, starting from selection of two and ending on inclusion of all 13, resulted in as high as 8178 regression equations. For each mode of size from 2 up to 13 the values of coefficients were optimized. The results of this first step of model seeking are documented in Fig. 1. It represents the distribution of obtained TP% values plotted against dimension of the regression function for all computed models.

11 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 39

Fig. 1. The distribution of predictive powers of all possible 8178 models for distinguishing cocrystals (TP) from simple binary eutectics (TN) formulated at the first stage of the procedure (threshold for FP%=0).

From results presented in Fig. 1 it can be inferred that there are several models worth considering, which include only few parameters. Indeed, the most surprising is that the simplest two-parameter model has quite efficient power of distinguishing CC from EU. It comprises atom counts descriptors interrelated via the following formula:  = 42.021.44 ∙ |∆naaN| # 9.680.33 ∙ |∆SHsOH|

(2)

The predicative power of the simplest model is provided in Fig. 2. Even such a crude representation of the components’ structure is quite efficient in assessments of crystallization probability. About 69.7% of TP are already embraced by the model with only 1.7% of the false positive cases. Furthermore, such simplified two-parameter quantification has quite predictive power as it was documented in left panel of Fig. 2. The classification of EU and CC according to distribution of regression function values is sufficient for practical purposes. It is

12 ACS Paragon Plus Environment

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

worth mentioning that one of our cases of unsuccessful crystallization is properly classified and among nine newly synthesized cocrystals, 6 are properly assigned. Another surprising factor comes from comparison of this simple model with approach using Hmix as criterions of crystallization. In such approach29–34 the affinity of potential cocrystals formers is quantified in terms of enthalpy of mixing of hypothetical super-cooled liquids under ambient conditions. Although this model was applied with some success29,30,32,33 it seems to be less adequate for phenolic acids characteristics. Indeed, as it is documented in Fig. 2 much steeper plot of FP% cases is observed for -Hmix distribution compared to the one resulting from two-parameter regression. The same behavior can be attributed to corresponding TP% plots.

Fig. 2. Predicative power of the simples two-parameters regression model. Right panel documents the distribution of true positive cases (cocrystals) and false negative cases (simple binary eutectic mixtures) included in the regressions solution of the highest predictive power. Additionally the plot of Matthews correlation coefficient (MCC) was provided. Left panel shows results of screening for new cocrystals based on regression A. Open circles denote data from training set, while red circles stand for results of our experiments.

13 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 39

The above presented two-parameter model does not belong to the top ten solutions obtained from the calculations. They are characterized in Supporting Information in Table S2. It has been found that the best model (denoted as C) which outperforms most other both in terms of maximizing TP% and minimizing FP% happens to be rich in eleven descriptors. The detailed regression formula has the following form:

 = −0.700.02 ∙ |∆WTPT- 5| # 43.750.05 ∙ |∆naaN| # 8.200.01 ∙ |∆SHsOH| # 13.070.01 ∙ |∆nHeteroRing| # 1.260.01 ∙ |∆ATSC1i| # 1.660.01 ∙ |∆MLFER_A| # 10.410.01 ∙ |∆GATS1i| # 3.730.00 ∙ |∆maxHBint6| − 16.380.01 ∙ |∆SsssN| # 31.470.03 ∙ |∆MDEN- 22| − 15.730.02 ∙ |∆C1SP2| (3)

This model has much higher precision and can properly account for about 80% positive cases (%TP=79.9%, Table S2). Interestingly, model denoted as H (Table S2) consists of only six parameters and it is slightly less accurate than the best model (%TP=79.4). Therefore this classification system, represented by the following equation also deserves attention due to its simplicity.

 = 43.112.43 ∙ |∆WTPT- 5| # 72.672.63 ∙ |∆naaN| # 20.411.70 ∙ |∆SHsOH| − 31.481.68 ∙ |∆MLFER_A| # 21.221.07 ∙ |∆MATS8m| # 42.993.93 ∙ |∆MDEN- 22|

(4)

Since relative values computed for each pair were normalized for given molecular descriptor, the values of regression coefficients directly inform about the significance of each parameter. In the case of all models with the highest predicting power, the most significant positive contribution to cocrystallization comes from the number of nitrogen atoms included in the

14 ACS Paragon Plus Environment

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

aromatic ring, naaN. This is not surprising since carboxylic and phenolic groups can be actively involved in intermolecular hydrogen bonding with such proton acceptor centers. This conclusion is further strengthened by contribution coming from nHeteroRing descriptor, which essentially brings the same kind of information. Another quite intuitive factor positively affecting ability of cocrystallization in the case of the best model represented by eq. 3 is included in the SHsOH descriptor. However, in this case the high contributions can be attributed rather to phenolic acids than to other coformers. Since SHsOH descriptor characterizing electrotopological state of hydrogen atoms from OH groups is frequently used in phenolic compounds QSAR studies its significance is probably related to the hydrogen bonding ability of ligands to the receptor.141–143 Similar contribution and also positive comes from Geary autocorrelation descriptor weighted by first ionization potential, GATS1i. The other positive influence comes from maxHBint6. There were negatives contributions to the cocrystallization propensities of the phenolic acids, the strongest of which was the number of sp2 carbons atoms quantized by the C1sp2 descriptor. This suggests that the presence of apolar regions is not beneficial for cocrystallization of PhA. Also closed structures of nitrogen atoms in sp3 configuration make it difficult for PhA to form intermolecular complexes. Other negative contributions of minor importance can be attributed to distributions of Broto-Moreau autocorrelation function weighted by van der Waals volumes, ATSC1v. Among ten systems tested experimentally for the purpose of this work this model properly assigns to EU subset the only case of unsuccessful cocrystallization. Two of nine newly synthesized cocrystals are misclassified, but this is within accuracy of the model. It predicts exclusion of about 30% of positive cases what is in full agreement with performed experiments. In summarizing it is worth emphasizing that the best regression model selected among many possible ones leads not only to efficient classification of pairs into cocrystals and simple binary eutectics (Fig. 3) but also provides very clear structural rationale of this effectiveness.

15 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 16 of 39

This shows that parameters of the proposed models are closely related to the intermolecular interactions formation potential, which is a similar concept to Hansen solubility parameters analysis.27,28 Hence restricting the interest to homogeneous classes of cocrystals formers seems to be beneficial despite the fact that the obtained model is probably not suited in general for screening purposes. However, it is directly applicable to finding such active pharmaceutical ingredient, for which synthesis of new cocrystals is expected to be very probable.

Fig. 3. Predictive power of the best model found among all 2036 studied. Notation is the same as in Fig. 2. Detailed values are provided in Supporting Information.

Application of the regression model The real strength of any QSPR model comes from its predictive potential and practical applications. Thus, full screening was performed for selected APIs using eq. 1. In the list of considered systems enumerated in Table S1 one can find several approved drugs. After selection of 48 of them, all possible combinations of pairs with 21 phenolic acids were generated and values of regression function were computed using the best model provided by 16 ACS Paragon Plus Environment

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

eq. 1. The resulting distributions were graphically collected in Fig. 4. It is evident that phenolic acids are supposed to be very good cocrystals formers of many drugs. Indeed, according to our prediction each of 21 phenolic acids can be used for cocrystallization with the following drugs: Voriconazole, Fluconazole, Adenine, Sulfamethazine, Epoxiconazole, Trimethroprim, Didanosine, Phenazine, Pyrazinamide, Theophylline, Caffeine, Nevirapine, Paliperidone, Meloxicam, Danazol, Etofylline, Stanozolol, Pentoxifylline, Piroxicam, Lornoxicam, Metronidazole, Iloperidone, Etodolac, Tadalafil, Imazamox, Imazethapyr, Isoniazid or Pyraclostrobin. On the contrary, none of PhA is suitable for cocrystallziation with curcumin and only few can be chosen for synthesis of molecular complexes in the solid state with Carbamazepine, Indomethacin, Progesterone or Pregnenolone. The full list of predictions made using eg.1 is provided in Supporting Information in Table S3 and S4.

Fig. 4. Distribution of the regression function values computed using eg. 1 for the set of potential cocrystals of 21 phenolic acids with 46 API. Red line represents threshold above which the propensity of cocrystallization is claimed to by high.

Experimental cocrystal screening In our previous works,32,33,43 cocrystallization abilities of phenolic acids with several popular drugs and drug-like compounds belonging to the class of amides were examined using fast 17 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 18 of 39

and efficient DEC method and theoretical COSMO-RS based calculations. These studies and many other reports22,62,92,110 showed that amides exhibit significant affinity to the phenolic acids. However, there are some exceptions. For instance benzamide cocrystalizes with salicylic acid, but on the other hand cannot form multicomponent crystals with structurally similar benzoic acid and 4-hydroxybenzoic acid.26 At this stage of the study new mixtures that have not been reported in the literature were examined. As a result of this, 9 new cocrystals were found. Noteworthy, only in the case of the benzamide-ferulic acid system was the cocrystal phase not detected. The results of experimental cocrystal screening were summarized in the Table 1. In Fig. 5 exemplary PXRD patterns and FTIR spectra of ethenzamide-ferulic acid were presented. The data for the other compounds in this study are summarized in the Supporting Information (Fig. S1-S9). The new cocrystal phase can be detected by the appearance of new diffraction peaks on the PXRD pattern, which cannot be assigned to the pure coformers and by the absorption bands shifts corresponding to the new hydrogen bonds formation on FTIR-ATR spectra. In the case of the ethenzamide-ferulic acid system, the cocrystal can be identified by the appearance of intense diffraction peak at 2θ=27.1° and two less abundant peaks at 2θ=11.3° and 13.6° (Fig. 5). Cocrystal formation can be also easily confirmed by the characteristic absorption bands shifts. As one can see from Fig. 5, new intermolecular interactions between ethenzamide and ferulic acid are manifested by the blue shifts of N-H absorption stretching modes, ν(NH) at 3169 and 3368 cm-1 on pure ethenzamide FTIR-ATR spectra, leading to the appearance of three peaks at 3186, 3289 and 3427 cm-1. Similarly a characteristic ν(NH) band shifts can be observed for benzamide-3,4dihydroxybenzoic acid (Fig. S3), nicotinamide-2,6-dihydroxybenzoic (Fig. S6), nicotinamide3,4-dihydroxybenzoic acid (Fig. S7) and nicotinamide-3,5-dihydroxybenzoic acid (Fig S8) solid dispersions. On the other hand, in the case of several other positive examples, due to the overlapping of absorption bands coming from the pure components, it is difficult to assign

18 ACS Paragon Plus Environment

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

unambiguously absorption bands to the amide ν(NH) stretching modes or to ν(OH) of phenolic acids. Nevertheless, in all these cases new absorption bands appeared on the spectra of the mixture which confirms new molecular complex formation. As one can see from Table. 1, new crystal phase can be observed in almost all cases except the benzamide-ferulic acid solid dispersion. In the case of this system all peaks corresponding to the mixture are overlapping with pure components diffraction signals. Another confirmation of immiscibility in the solid state comes from FTIR-ATR spectra, since no absorption bands shifts can be observed (Fig. S4).

Table 1. The results of experimental cocrystal screening, “+” denotes positive cases (cocrystal formation), while “-“ denotes negative cases. amide

phenolic acid

screening characteristic score

urea

benzamide

new crystal phase 2θ signals [°]

3,4-dihydroxybenzoic acid +

10.8, 17.2, 27.1

ferulic acid

8.8, 19.6, 24.4, 28.0

+

3,4-dihydroxybenzoic acid +

23.2, 27.9

ferulic acid

-

-

salicylamide

ferulic acid

+

16.1, 18.0, 19.6, 23.1, 25.9, 27.2

ethenzamide

ferulic acid

+

11.3, 13.6, 27.1

nicotinamide 2,6-dihydroxybenzoic acid +

26.0, 27.8

3,4-dihydroxybenzoic acid +

8.2, 16.5

3,5-dihydroxybenzoic acid +

9.6, 14.0, 14.6, 16.8, 17.3

ferulic acid

19.2, 25.3, 26.1, 27.0

+

19 ACS Paragon Plus Environment

Crystal Growth & Design

mixture ferulic acid

Intensity

ethenzamide

5.0

10.0

15.0

20.0

25.0

30.0

35.0

40.0

2θ [° ] 1.0

normalized transmitance

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 39

0.8

0.6

0.4

mixture

0.2

ferulic acid ethenzamide 0.0 2400

2600

2800

3000

3200

3400

3600

wave number [cm-1]

Fig. 5. PXRD pattern and FTIR-ATR spectra of ethenzamide-ferulic acid system.

Conclusions In this paper novel screening approach involving simple molecular indices for assessment of phenolic acids cocrystal potential was developed. The selection of these particular molecular compounds was guided by their practical applicability since phenolic acids containing solid dispersions have been extensively studied. Proper selection of pairs of

20 ACS Paragon Plus Environment

Page 21 of 39 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

chemical compounds for successful cocrystallization is not a trivial task. Hence, many theoretical models of cocrystals thermodynamic properties prediction involving combination of structural features and descriptors have been proposed.35–39 This paper presents new strategy for theoretical cocrystals screening which utilizes 1D and 2D molecular descriptors. This method offers several unique and advantageous features. First of all, this approach enables for finding adequate molecular descriptors properly assigning both similarities and diversions of structural, electronic and energetic contributions for homogenous set of cocrystals. The crucial step is addressed to proper definition of the target function enabling for maximization of involved true positive cases and at the same time minimization of number of false positives included in the prediction. Here it is done by optimization of weighting coefficients of relative values of selected molecular descriptors. As a result, a series of equations were generated with very promising predictive power. The real strength of the proposed approach is clearly documented by finding very simple two-parameter regression formula, which is only slightly less precise compared to the best model defined using eleven molecular descriptors. The proposed and validated procedure is very simple and does not require complex computations since all necessary information about the structure is decoded from canonical SMILES string. Since the focus was restricted exclusively to phenolic acids as cocrystals formers proposed model is probably not adequate for other sets of compounds. However, it is general in the sense that broad range of active pharmaceutical ingredients can be tested for cocrystallization with considered here group of coformers. The extension to other classes of compounds seems to be straightforward and worth further exploration.

SUPPORTING INFORMATION

21 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 22 of 39

Supporting information contains phenolic acids cocrystals and eutectics database, classification models details, theoretical cocrystal screening results and collection of PXRD patterns and FTIR-ATR spectra recorded for newly synthesized dispersions.

ASSOCIATED CONTENT AUTHOR INFORMATION Corresponding Author *

E-mail: [email protected]

ACKNOWLEDGMENT This study was supported by the Nicolaus Copernicus University in Toruń (grant MN2/WF/2017). REFERENCES (1)

Chattopadhyay, B.; Jacobs, L.; Panini, P.; Salzmann, I.; Resel, R.; Geerts, Y. Accessing Phase-Pure and Stable Acetaminophen Polymorphs by Thermal Gradient Crystallization. Cryst. Growth Des. 2018, 18 (3), 1272–1277.

(2)

Szilágyi, B.; Nagy, Z. K. Population Balance Modeling and Optimization of an Integrated Batch Crystallizer–Wet Mill System for Crystal Size Distribution Control.

Cryst. Growth Des. 2018, 18 (3), 1415–1424. (3)

Chadha, R.; Bhalla, Y. Cocrystallization in Nutraceuticals. In Recrystallization in

Materials Processing; InTech, 2015; pp 35–50. (4)

Zaworotko, M.; Clarke, H.; Kapildev, A.; Kavuru, P.; Shytle, R. D.; Pujari, T.; Marshall, L.; Ong, T. T. Nutraceutical Co-Crystal Compositions, U.S. Patent

22 ACS Paragon Plus Environment

Page 23 of 39 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

20100204204A1, August 12, 2010. (5)

Laranjinha, J.; Vieira, O.; Almeida, L.; Madeira, V. Inhibition of metmyoglobin/H2O2Dependent Low Density Lipoprotein Lipid Peroxidation by Naturally Occurring Phenolic Acids. Biochem. Pharmacol. 1996, 51 (4), 395–402.

(6)

Yeh, K. L.; Lee, T. Intensified Crystallization Processes for 1:1 Drug-Drug Cocrystals of Sulfathiazole-Theophylline, and Sulfathiazole-Sulfanilamide. Cryst. Growth Des. 2018, 18 (3), 1339–1349.

(7)

Sathya, P.; Vidyalaksmi, Y.; Pugazhendhi, S.; Gopalakrishnan, R. Benzotriazole pHydroxybenzoic Acid: Physicochemical and Biological Evaluation of an Organic Cocrystal. Mater. Res. Innov. 2016, 8917, 1–7.

(8)

Putra, O. D.; Umeda, D.; Nugraha, Y. P.; Nango, K.; Yonemochi, E.; Uekusa, H. Simultaneous Improvement of Epalrestat Photostability and Solubility via Cocrystallization: A Case Study. Cryst. Growth Des. 2018, 18 (1), 373–379.

(9)

Badhani, B.; Sharma, N.; Kakkar, R. Gallic Acid: A Versatile Antioxidant with Promising Therapeutic and Industrial Applications. RSC Adv. 2015, 5 (35), 27540– 27557.

(10)

Silva, F. A. M.; Borges, F.; Guimarães, C.; Lima, J. L. F. C.; Matos, C.; Reis, S. Phenolic Acids and Derivatives: Studies on the Relationship among Structure, Radical Scavenging Activity, and Physicochemical Parameters. J. Agric. Food Chem. 2000, 48 (6), 2122–2126.

(11)

Khuwijitjaru, P.; Plernjit, J.; Suaylam, B.; Samuhaseneetoo, S.; Pongsawatmanit, R.; Adachi, S. Degradation Kinetics of Some Phenolic Compounds in Subcritical Water and Radical Scavenging Activity of Their Degradation Products. Can. J. Chem. Eng. 2014, 92 (5), 810–815.

(12)

Zhao, H.; Zhang, H.; Yang, S. Phenolic Compounds and Its Antioxidant Activities in

23 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 39

Ethanolic Extracts from Seven Cultivars of Chinese Jujube. Food Sci. Hum. Wellness 2014, 3 (3–4), 183–190. (13)

Nagarani, G.; Abirami, A.; Siddhuraju, P. Food Prospects and Nutraceutical Attributes of Momordica Species: A Potential Tropical Bioresources – A Review. Food Sci. Hum.

Wellness 2014, 3 (3–4), 117–126. (14)

Natella, F.; Nardini, M.; Di Felice, M.; Scaccini, C. Benzoic and Cinnamic Acid Derivatives as Antioxidants: Structure−Activity Relation. J. Agric. Food Chem. 1999,

47 (4), 1453–1459. (15)

Alam, M. A.; Subhan, N.; Hossain, H.; Hossain, M.; Reza, H. M.; Rahman, M. M.; Ullah, M. O. Hydroxycinnamic Acid Derivatives: A Potential Class of Natural Compounds for the Management of Lipid Metabolism and Obesity. Nutr. Metab. 2016,

13 (1), 27. (16)

Sreekanth, B. R.; Vishweshwar, P.; Vyas, K. Supramolecular Synthon Polymorphism in 2 : 1 Co-Crystal of 4-Hydroxybenzoic Acid and 2,3,5,6-Tetramethylpyrazine. Chem.

Commun. 2007, No. 23, 2375. (17)

Mukherjee, A.; Desiraju, G. R. Synthon Polymorphism and Pseudopolymorphism in Co-Crystals. The 4,4′-bipyridine–4-Hydroxybenzoic Acid Structural Landscape. Chem.

Commun. 2011, 47 (14), 4090. (18)

Swapna, B.; Maddileti, D.; Nangia, A. Cocrystals of the Tuberculosis Drug Isoniazid: Polymorphism, Isostructurality, and Stability. Cryst. Growth Des. 2014, 14 (11), 5991– 6005.

(19)

Zhang, T.; Yang, Y.; Zhao, X.; Jia, J.; Su, H.; He, H.; Gu, J.; Zhu, G. Dissolution and Pharmacokinetic Properties of Two Paliperidone Cocrystals with 4-Hydroxybenzoic and 4-Aminobenzoic Acid. CrystEngComm 2014, 16 (33), 7667.

(20)

Brown, R. Investigating Pharmaceutical Co-Crystals as a Means to Improve the

24 ACS Paragon Plus Environment

Page 25 of 39 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

Solubility of a Drug. Thesis, University of Kansas, 2012. (21)

Sreenivas Reddy, L.; Bethune, S. J.; Kampf, J. W.; Rodríguez-Hornedo, N. Cocrystals and Salts of Gabapentin: pH Dependent Cocrystal Stability and Solubility. Cryst.

Growth Des. 2009, 9 (1), 378–385. (22)

Aitipamula, S.; Wong, A. B. H.; Chow, P. S.; Tan, R. B. H. Pharmaceutical Cocrystals of Ethenzamide: Structural, Solubility and Dissolution Studies. CrystEngComm 2012,

14 (24), 8515. (23)

Nalte, Y.; Nalte, Y. K.; Arsul, V. A.; Shep, S. G.; Sunil Bothara, D. B. Solubility Enhancement of Nevirapine by Cocrystallisation Technique. J. Pharm. Res. 2015, 99 (8), 556–561.

(24)

Ash, M.; Ash, I. Handbook of Preservatives; Synapse Information Resources: Endicott: NY, 2004.

(25)

Sanphui, P.; Kumar, S. S.; Nangia, A. Pharmaceutical Cocrystals of Niclosamide.

Cryst. Growth Des. 2012, 12 (9), 4588–4599. (26)

Seaton, C. C.; Parkin, A. Making Benzamide Cocrystals with Benzoic Acids: The Influence of Chemical Structure. Cryst. Growth Des. 2011, 11 (5), 1502–1511.

(27)

Gadade, D. D.; Pekamwar, S. S.; Lahoti, S. R.; Patni, S. D.; Sarode, M. C. Cocrystallization of Etodolac: Prediction of Cocrystallization, Synthesis, Solid State Characterization And In Vitro Drug Release. Marmara Pharm. J. 2017, 21 (1), 78–88.

(28)

Mohammad, M. A.; Alhalaweh, A.; Velaga, S. P. Hansen Solubility Parameter as a Tool to Predict Cocrystal Formation. Int. J. Pharm. 2011, 407 (1–2), 63–71.

(29)

Klamt, A. Solvent-Screening and Co-Crystal Screening for Drug Development with COSMO-RS. J. Cheminform. 2012, 4 (Suppl 1), O14.

(30)

Cysewski, P. Efficacy of Bi-Component Cocrystals and Simple Binary Eutectics Screening Using Heat of Mixing Estimated under Super Cooled Conditions. J. Mol.

25 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 26 of 39

Graph. Model. 2016, 68, 23–28. (31)

Cysewski, P.; Przybyłek, M. Selection of Effective Cocrystals Former for Dissolution Rate Improvement of Active Pharmaceutical Ingredients Based on Lipoaffinity Index.

Eur. J. Pharm. Sci. 2017, 107, 87–96. (32)

Cysewski, P.; Przybyłek, M.; Ziółkowska, D.; Mroczyńska, K. Exploring the Cocrystallization Potential of Urea and Benzamide. J. Mol. Model. 2016, 22 (5), 103.

(33)

Przybyłek, M.; Ziółkowska, D.; Mroczyńska, K.; Cysewski, P. Propensity of Salicylamide and Ethenzamide Cocrystallization with Aromatic Carboxylic Acids. Eur.

J. Pharm. Sci. 2016, 85, 132–140. (34)

Przybyłek, M.; Ziółkowska, D.; Mroczyńska, K.; Cysewski, P. Applicability of Phenolic Acids as Effective Enhancers of Cocrystal Solubility of Methylxanthines.

Cryst. Growth Des. 2017, 17 (4), 2186–2193. (35)

Wicker, J. G. P.; Crowley, L. M.; Robshaw, O.; Little, E. J.; Stokes, S. P.; Cooper, R. I.; Lawrence, S. E. Will They Co-Crystallize? CrystEngComm 2017, 19 (36), 5336– 5340.

(36)

Rama Krishna, G.; Ukrainczyk, M.; Zeglinski, J.; Rasmuson, Å. C. Prediction of Solid State Properties of Cocrystals Using Artificial Neural Network Modeling. Cryst.

Growth Des. 2018, 18 (1), 133–144. (37)

Gamidi, R. K.; Rasmuson, Å. C. Estimation of Melting Temperature of Molecular Cocrystals Using Artificial Neural Network Model. Cryst. Growth Des. 2017, 17 (1), 175–182.

(38)

Perlovich, G. L. Thermodynamic Characteristics of Cocrystal Formation and Melting Points for Rational Design of Pharmaceutical Two-Component Systems.

CrystEngComm 2015, 17 (37), 7019–7028. (39)

Wood, P. A.; Feeder, N.; Furlow, M.; Galek, P. T. A.; Groom, C. R.; Pidcock, E.

26 ACS Paragon Plus Environment

Page 27 of 39 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

Knowledge-Based Approaches to Co-Crystal Design. CrystEngComm 2014, 16 (26), 5839. (40)

Todeschini, R.; Consonni, V. Molecular Descriptors for Chemoinformatics: Vol. 1&2; Wiley-VCH, 2009.

(41)

Cysewski, P.; Przybyłek, M.; Miernik, T.; Kobierski, M.; Ziółkowska, D. On the Origin of Surfaces-Dependent Growth of Benzoic Acid Crystal Inferred through the Droplet Evaporation Method. Struct. Chem. 2015, 26 (3), 705–712.

(42)

Przybyłek, M.; Cysewski, P.; Pawelec, M.; Ziółkowska, D.; Kobierski, M. On the Origin of Surface Imposed Anisotropic Growth of Salicylic and Acetylsalicylic Acids Crystals during Droplet Evaporation. J. Mol. Model. 2015, 21 (3), 49.

(43)

Przybyłek, M.; Ziółkowska, D.; Kobierski, M.; Mroczyńska, K.; Cysewski, P. Utilization of Oriented Crystal Growth for Screening of Aromatic Carboxylic Acids Cocrystallization with Urea. J. Cryst. Growth 2016, 433, 128–138.

(44)

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 (2), 171– 179.

(45)

Kim, S.; Thiessen, P. A.; Bolton, E. E.; Chen, J.; Fu, G.; Gindulyte, A.; Han, L.; He, J.; He, S.; Shoemaker, B. A.; Wang, J.; Yu, B.; Zhang, J.; Bryant, S. H. PubChem Substance and Compound Databases. Nucleic Acids Res. 2016, 44 (D1), D1202– D1213.

(46)

Yap, C. W. PaDEL-Descriptor: An Open Source Software to Calculate Molecular Descriptors and Fingerprints. J. Comput. Chem. 2011, 32 (7), 1466–1474.

(47)

Matthews, B. W. Comparison of the Predicted and Observed Secondary Structure of T4 Phage Lysozyme. BBA - Protein Struct. 1975, 405 (2), 442–451.

(48)

Ellison, S. L. R.; Barwick, V. J.; Farrant, T. J. D. Practical Statistics for the Analytical

27 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 28 of 39

Scientist; Royal Society of Chemistry: Cambridge, 2009. (49)

Cheney, M. L.; Weyna, D. R.; Shan, N.; Hanna, M.; Wojtas, L.; Zaworotko, M. J. Supramolecular Architectures of Meloxicam Carboxylic Acid Cocrystals, a Crystal Engineering Case Study. Cryst. Growth Des. 2010, 10 (10), 4401–4413.

(50)

Zhang, X.; Tian, Y.; Jia, J.; Zhang, T.; Zhu, G. Synthesis, Characterization and Dissolution of Three Pharmaceutical Cocrystals Based on Deferiprone. J. Mol. Struct. 2016, 1108, 560–566.

(51)

Bučar, D. K.; Henry, R. F.; Lou, X.; Duerst, R. W.; MacGillivray, L. R.; Zhang, G. G. Z. Cocrystals of Caffeine and Hydroxybenzoic Acids Composed of Multiple Supramolecular Heterosynthons: Screening via Solution-Mediated Phase Transformation and Structural Characterization. Cryst. Growth Des. 2009, 9 (4), 1932– 1943.

(52)

Childs, S. L.; Hardcastle, K. I. Cocrystals of Chlorzoxazone with Carboxylic Acids.

CrystEngComm 2007, 9 (5), 364. (53)

Ghosh, S.; Bag, P. P.; Reddy, C. M. Co-Crystals of Sulfamethazine with Some Carboxylic Acids and Amides: Co-Former Assisted Tautomerism in an Active Pharmaceutical Ingredient and Hydrogen Bond Competition Study. Cryst. Growth Des. 2011, 11 (8), 3489–3503.

(54)

Lahamage, S. R.; Darekar, A. B.; Zalte, A. G.; Saudagar, R. B. Co-Crystal Approach for Improving Poorly Soluble Drug Agomelatine by Using Acid Coformer for Tablets.

Int. J. Pharm. Pharm. Res. 2016, 7 (71), 19–34. (55)

Singh, N. B.; Das, S. S.; Singh, N. P.; Agrawal, T. Computer Simulation, Thermodynamic and Microstructural Studies of Benzamide-Benzoic Acid Eutectic System. J. Cryst. Growth 2008, 310 (11), 2878–2884.

(56)

Sharma, B. L.; Kant, R.; Sharma, R.; Tandon, S. Deviations of Binary Organic Eutectic

28 ACS Paragon Plus Environment

Page 29 of 39 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

Melt Systems. Mater. Chem. Phys. 2003, 82 (1), 216–224. (57)

Rai, U. S.; Shekhar, H. Chemistry of Organic Eutectics: Phenanthrene — Benzoic Acid and Phenanthrene — Cinnamic Acid Systems. Cryst. Res. Technol. 1990, 25 (7), 771– 779.

(58)

Postnikov, V. A. Phase Diagram of the Eutectic Benzoic Acid-Naphthalene System in the Temperature Range of 300–400 K. Russ. J. Phys. Chem. A 2014, 88 (8), 1336– 1339.

(59)

Reschke, T.; Zherikova, K. V; Verevkin, S. P.; Held, C. Benzoic Acid and Chlorobenzoic Acids: Thermodynamic Study of the Pure Compounds and Binary Mixtures With Water. J. Pharm. Sci. 2016, 105 (3), 1050–1058.

(60)

Kumar, R.; Kumar, R.; Dixit, A. Thermal Phase Diagram of Acetamide-Benzoic Acid and Benzoic Acid-Phthalimide Binary Systems for Solar Thermal Applications. In AIP

Conference Proceedings; 2016; Vol. 1728, p 20687. (61)

Wang, T. C.; Deng, K. C. Solid-Liquid Equilibria for Three Binary Mixtures of Benzoic Acid with Heptanedioic Acid, 3-Methylpentanedioic Acid, and 2,3Dimethylbutanedioic Acid. J. Chem. Eng. Data 2012, 57 (1), 78–81.

(62)

Manin, A. N.; Voronin, A. P.; Drozd, K. V.; Manin, N. G.; Bauer-Brandl, A.; Perlovich, G. L. Cocrystal Screening of Hydroxybenzamides with Benzoic Acid Derivatives: A Comparative Study of Thermal and Solution-Based Methods. Eur. J.

Pharm. Sci. 2014, 65, 56–64. (63)

Kaur, R.; Gautam, R.; Cherukuvada, S.; Guru Row, T. N. Do Carboximide-Carboxylic Acid Combinations Form Co-Crystals? The Role of Hydroxyl Substitution on the Formation of Co-Crystals and Eutectics. IUCrJ 2015, 2 (3), 341–351.

(64)

Kojima, T.; Tsutsumi, S.; Yamamoto, K.; Ikeda, Y.; Moriwaki, T. High-Throughput Cocrystal Slurry Screening by Use of in Situ Raman Microscopy and Multi-Well Plate.

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

Page 30 of 39

Int. J. Pharm. 2010, 399 (1–2), 52–59. (65)

Sangster, J. Phase Diagrams and Thermodynamic Properties of Binary and Ternary Systems Based on Nitroaromatic Compounds. J. Phys. Chem. Ref. Data 1997, 26 (2), 351–502.

(66)

Sharma, S.; Kumar, A.; Singh, V. Dependence of Eutectic - Composition on Densities of Components of Binary Eutectic Mixtures of Solid Carboxylic Acids. J. Appl. Chem. 2012, 3 (1), 18–22.

(67)

Rai, U. S.; Rai, R. N. Physical Chemistry of Organic Eutectics. J. Therm. Anal.

Calorim. 1998, 53 (3), 883–893. (68)

Sharma, B. L.; Jamwal, R.; Kant, R. Thermodynamic and Lamella Models Relationship for the Eutectic System Benzoic Acid - Cinnamic Acid. Cryst. Res. Technol. 2004, 39 (5), 454–464.

(69)

Singh, N. B.; Kumar, P. Solidification Behavior of the Cinnamic Acid-P -Nitrophenol Eutectic System. J. Chem. Eng. Data 1986, 31 (4), 406–408.

(70)

Cui, D. X.; Sugibayashi, K.; Morimoto, Y.; Li, F. L. Study on Transdermal Absorption of Borneol-Salicylic Acid Eutectic Mixture. Acta Pharm. Sin. 1989, 24 (7), 556–559.

(71)

Meltzer, V.; Pincu, E. A DSC Study for Binary Mixture of 2-Chlorobenzoic Acid with Salicylic Acid. Rev. Roum. Chim. 2009, 54 (5), 333–338.

(72)

Liao, X.; Gautam, M.; Grill, A.; Zhu, H. J. Effect of Position Isomerism on the Formation and Physicochemical Properties of Pharmaceutical Co-Crystals. J. Pharm.

Sci. 2010, 99 (1), 246–254. (73)

Mayer, M. M.; Howell, W. J.; Tomasko, D. L.; Eckert, C. A.; Mayer, M. M.; Howell, W. J.; Eckert, C. A. Solid-Liquid Equilibria in the Systems Thianthrene + Phenanthrene, Salicylic Acid + Phenanthrene, and 3-Hydroxybenzoic Acid + Phenanthrene. J. Chem. Eng. Data 1990, 35 (4), 446–449.

30 ACS Paragon Plus Environment

Page 31 of 39 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

(74)

Zuń, M.; Wojciechowska, K.; Dwornicka, D.; Poleszak, E. The Influence of the Eutectic Mixtures: Salicylic Acid - Menthol and Benzocaine - Menthol on Physical Properties of the Creams with Fluconazole. Curr. Issues Pharm. Med. Sci. 2013, 26 (4), 457–460.

(75)

Sathisaran, I.; Dalvi, S. V. Crystal Engineering of Curcumin with Salicylic Acid and Hydroxyquinol as Coformers. Cryst. Growth Des. 2017, 17 (7), 3974–3988.

(76)

Friscić, T.; Lancaster, R. W.; Fábián, L.; Karamertzanis, P. G. Tunable Recognition of the Steroid Alpha-Face by Adjacent Pi-Electron Density. Proc. Natl. Acad. Sci. U. S. A. 2010, 107 (30), 13216–13221.

(77)

Goud, N. R.; Suresh, K.; Sanphui, P.; Nangia, A. Fast Dissolving Eutectic Compositions of Curcumin. Int. J. Pharm. 2012, 439 (1–2), 63–72.

(78)

Keramatnia, F.; Shayanfar, A.; Jouyban, A. Thermodynamic Solubility Profile of Carbamazepine-Cinnamic Acid Cocrystal at Different pH. J. Pharm. Sci. 2015, 104 (8), 2559–2565.

(79)

Seaton, C. C.; Parkin, A.; Wilson, C. C.; Blagden, N. Controlling the Formation of Benzoic Acid: Isonicotinamide Molecular Complexes. Cryst. Growth Des. 2009, 9 (1), 47–56.

(80)

Aakeröy, C. B.; Desper, J.; Scott, B. M. T. Balancing Supramolecular Reagents for Reliable Formation of Co-Crystals. Chem. Commun. (Camb). 2006, 0 (13), 1445–1447.

(81)

Sangster, J. Phase Diagrams and Thermodynamic Properties of Binary Organic Systems Based on 1,2-, 1,3-, 1,4-Diaminobenzene or Benzidine. J. Phys. Chem. Ref.

Data 1994, 23 (2), 295–338. (82)

Gujar, P. P.; Kokil, A. A.; Karekar, P. S.; Gurav, Y. A.; Yadav, A. V. Improvement in Physicochemical Properties of Nevirapine. Int. J. Pharm. Technol. IJPT 2013, 4 (4), 4831–4842.

31 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

(83)

Page 32 of 39

Bučar, D.-K.; Day, G. M.; Halasz, I.; Zhang, G. G. Z.; Sander, J. R. G.; Reid, D. G.; MacGillivray, L. R.; Duer, M. J.; Jones, W. The Curious Case of (Caffeine)·(benzoic Acid): How Heteronuclear Seeding Allowed the Formation of an Elusive Cocrystal.

Chem. Sci. 2013, 4 (12), 4417. (84)

Sangster, J. Phase Diagrams and Thermodynamic Properties of Binary Systems of Drugs. J. Phys. Chem. Ref. Data 1999, 28 (4), 889–930.

(85)

Da Silva, C. C. P.; Pepino, R. de O.; De Melo, C. C.; Tenorio, J. C.; Ellena, J. Controlled Synthesis of New 5-Fluorocytosine Cocrystals Based on the P K a Rule.

Cryst. Growth Des. 2014, 14 (9), 4383–4393. (86)

Kumar, S. S. Polymorphism and Solubility of Selected Active Pharmaceutical Ingredients. PhD Thesis, University of Hyderabad, 2014.

(87)

Nauha, E. Crystalline Forms of Selected Agrochemical Actives: Design and Synthesis of Cocrystals. PhD Thesis, University of Jyväskylä, 2012.

(88)

Lou, M.; Mao, S. H.; Luo, Y. H.; Zhao, P.; Sun, B. W. Synthesis, Co-Crystal Structure and Characterization of Pyrazinamide with m-Hydroxybenzoic Acid, pHydroxybenzoic Acid and 3,4-Dihydroxy Benzolic Acid. Res. Chem. Intermed. 2015,

41 (5), 2939–2951. (89)

Wang, W.; Liu, B. W.; Liu, J.; Ren, R. 4-Hydroxy-Benzoic Acid-1H-Imidazole (1/1).

Acta Crystallogr. Sect. E Struct. Reports Online 2009, 65 (6), o1205–o1205. (90)

Zhou, Z.; Chan, H. M.; Sung, H. H.-Y.; Tong, H. H. Y.; Zheng, Y. Identification of New Cocrystal Systems with Stoichiometric Diversity of Salicylic Acid Using Thermal Methods. Pharm. Res. 2016, 33 (4), 1030–1039.

(91)

Lee, K.-S.; Kim, K.-J.; Ulrich, J. Formation of Salicylic Acid/4,4′-Dipyridyl Cocrystals Based on the Ternary Phase Diagram. Chem. Eng. Technol. 2015, 38 (6), 1073–1080.

(92)

Manin, A. N.; Voronin, A. P.; Shishkina, A. V.; Vener, M. V.; Churakov, A. V.;

32 ACS Paragon Plus Environment

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

Perlovich, G. L. Influence of Secondary Interactions on the Structure, Sublimation Thermodynamics, and Solubility of Salicylate:4-Hydroxybenzamide Cocrystals. Combined Experimental and Theoretical Study. J. Phys. Chem. B 2015, 119 (33), 10466–10477. (93)

Lee, K.; Kim, K. J.; Ulrich, J. N-H⋯O, O-H⋯O Hydrogen Bonded Supramolecular Formation in the Cocrystal of Salicylic Acid with N-Containing Bases. Cryst. Res.

Technol. 2016, 51 (3), 197–206. (94)

Hsu, I. N.; Gellert, R. W. The 1:1 Crystal Complex of Salicylic Acid and Urea. J.

Crystallogr. Spectrosc. Res. 1983, 13 (1), 43–48. (95)

Alatas, F.; Soewandhi, S. N.; Sasongko, L.; Ismunandar; Uekusa, H. Cocrystal Formation between Didanosine and Two Aromatic Acids. Int. J. Pharm. Pharm. Sci. 2013, 5, 275–280.

(96)

Kastelic, J.; Lah, N.; Kikelj, D.; Leban, I. A 1:1 Cocrystal of Fluconazole with Salicylic Acid. Acta Crystallogr. Sect. C Cryst. Struct. Commun. 2011, 67 (9), o370– o372.

(97)

Hathwar, V. R.; Pal, R.; Guru Row, T. N. Charge Density Analysis of Crystals of Nicotinamide with Salicylic Acid and Oxalic Acid: An Insight into the Salt to Cocrystal Continuum. Cryst. Growth Des. 2010, 10 (8), 3306–3310.

(98)

Byres, M.; Cox, P. J.; Kay, G.; Nixon, E. Supramolecular Structures of Six AdenineCarboxylic Acid Complexes. CrystEngComm 2009, 11 (1), 135–142.

(99)

Nanjwade, V. K.; Manvi, F. V.; Ali, M. S.; Nanjwade, B. K. Development and Evaluation of Prulifloxacin Tablet. Int. J. Drug Formul. Res. 2011, 2 (6), 302–314.

(100) Zhang, S.; Chen, H.; Rasmuson, Å. C. Thermodynamics and Crystallization of the Theophylline – Salicylic Acid Cocrystal. CrystEngComm 2013, 0 (22), 1–3. (101) Vinesha, V.; Sevukarajan, M.; Rajalakshmi, R.; Chowdary, G. T.; Haritha, K.;

33 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 34 of 39

Vidyanikethan, S.; Dt, C.; Pradesh, A. Enhancement of Solubility of Tadalafil By Cocrystal Approach. Int. Res. J. Pharm 2013, 4 (4), 218–223. (102) Emsley, J.; Reza, N. M.; Kuroda, R. Hydrogen Bonding of Urea-Salicylic Acid, U·SA.

J. Crystallogr. Spectrosc. Res. 1986, 16 (1), 57–69. (103) Lipert, M. P.; Roy, L.; Childs, S. L.; Rodríguez-Hornedo, N. Cocrystal Solubilization in Biorelevant Media and Its Prediction from Drug Solubilization. J. Pharm. Sci. 2015,

104 (12), 4153–4163. (104) Patel, D. M.; Shah, H. R.; Patel, R. J.; Patel, C. N. Preparation and Characterization of Lornoxicam Co-Crystals. World J. Pharm. Pharm. Sci. 2014, 3 (6), 713–732. (105) Du, Y.; Zhang, H.; Xue, J.; Fang, H.; Zhang, Q.; Xia, Y.; Li, Y.; Hong, Z. Raman and Terahertz Spectroscopical Investigation of Cocrystal Formation Process of Piracetam and 3-Hydroxybenzoic Acid. Spectrochim. Acta - Part A Mol. Biomol. Spectrosc. 2015,

139, 488–494. (106) Lemmerer, A.; Adsmond, D. A.; Esterhuysen, C.; Bernstein, J. Polymorphic CoCrystals from Polymorphic Co-Crystal Formers: Competition between Carboxylic Acid⋯pyridine and Phenol⋯pyridine Hydrogen Bonds. Cryst. Growth Des. 2013, 13 (9), 3935–3952. (107) Mashhadi, S. M. A.; Yunus, U.; Bhatti, M. H.; Tahir, M. N. Isoniazid Cocrystals with Anti-Oxidant Hydroxy Benzoic Acids. J. Mol. Struct. 2014, 1076, 446–452. (108) Zhang, T.; Yang, Y.; Wang, H.; Sun, F.; Zhao, X.; Jia, J.; Liu, J.; Guo, W.; Cui, X.; Gu, J.; Zhu, G. Using Dissolution and Pharmacokinetics Studies of Crystal Form to Optimize the Original Iloperidone. Cryst. Growth Des. 2013, 13 (12), 5261–5266. (109) Luo, T. J. M.; Palmore, G. T. R. Engineering Crystalline Architecture with Supramolecular Tapes: Studies on Secondary Donor-Acceptor Interactions in Cocrystals of the Cyclic Dipeptide of Glycine. Cryst. Growth Des. 2002, 2 (5), 337–

34 ACS Paragon Plus Environment

Page 35 of 39 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

350. (110) Aitipamula, S.; Chow, P. S.; Tan, R. B. H. Trimorphs of a Pharmaceutical Cocrystal Involving Two Active Pharmaceutical Ingredients: Potential Relevance to Combination Drugs. CrystEngComm 2009, 11 (9), 1823. (111) Shattock, T. R. Crystal Engineering of Co-Crystals and Their Relevance to Pharmaceutical Forms. PhD Thesis, University of South Florida, 2007. (112) Mukherjee, A.; Desiraju, G. R. Combinatorial Exploration of the Structural Landscape of Acid-Pyridine Cocrystals. Cryst. Growth Des. 2014, 14 (3), 1375–1385. (113) Back, K. The Crystallisation of Conformationally Flexible Molecules, PhD Thesis, The University of Manchester, 2011. (114) Gdaniec, M.; Połoński, T. Phenazine-2,6-Dihydroxybenzoic Acid (1/1). Acta

Crystallogr. Sect. E Struct. Reports Online 2007, 63 (2), o707–o709. (115) Cohen, D. E.; Benedict, J. B.; Morlan, B.; Chiu, D. T.; Kahr, B. Dyeing Polymorphs: The MALDI Host 2,5-Dihydroxybenzoic Acid. Cryst. Growth Des. 2007, 7 (3), 492– 495. (116) Childs, S. L.; Kandi, P.; Lingireddy, S. R. Formulation of a Danazol Cocrystal with Controlled Supersaturation Plays an Essential Role in Improving Bioavailability. Mol.

Pharm. 2013, 10 (8), 3112–3127. (117) Jacobs, A.; Amombo Noa, F. M. Co-Crystals and Co-Crystal Hydrates of Vanillic Acid. CrystEngComm 2015, 17 (1), 98–106. (118) Adalder, T. K.; Sankolli, R.; Dastidar, P. Homo- or Heterosynthon? A Crystallographic Study on a Series of New Cocrystals Derived from Pyrazinecarboxamide and Various Carboxylic Acids Equipped with Additional Hydrogen Bonding Sites. Cryst. Growth

Des. 2012, 12 (5), 2533–2542. (119) Childs, S. L. Metronidazole Cocrystals and Imipramine Cocrystals, U.S. Patent

35 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 36 of 39

US8163790B2, April 24, 2012. (120) Sekhon, B. Nutraceutical Cocrystals: An Overview. RGUHS J. Pharm. Sci. 2012, 2 (2), 16–25. (121) Mukherjee, S. Crystal Engineering of Pharmaceutical Co-Crystals. Thesis, University of South Florida, 2011. (122) Dabir, T. O.; Gaikar, V. G.; Jayaraman, S.; Mukherjee, S. Thermodynamic Modeling Studies of Aqueous Solubility of Caffeine, Gallic Acid and Their Cocrystal in the Temperature Range of 303 K–363 K. Fluid Phase Equilib. 2018, 456, 65–76. (123) Sinha, A. S.; Rao Khandavilli, U. B.; O’Connor, E. L.; Deadman, B. J.; Maguire, A. R.; Lawrence, S. E. Novel Co-Crystals of the Nutraceutical Sinapic Acid. CrystEngComm 2015, 17 (26), 4832–4841. (124) Tan, Z.; Zhu, E.; Luo, L.; Lin, Z.; Yan, R. 3-(4-Hy-Droxy-3-Meth-Oxy-PhenYl)acrylic Acid-2,3,5,6-Tetra-Methyl-Pyrazine (2/1). Acta Crystallogr. Sect. E. Struct.

Rep. Online 2011, 67 (Pt 2), o424. (125) Bevill, M. J.; Vlahova, P. I.; Smit, J. P. Polymorphic Cocrystals of Nutraceutical Compound P -Coumaric Acid with Nicotinamide: Characterization, Relative SolidState Stability, and Conversion to Alternate Stoichiometries. Cryst. Growth Des. 2014,

14 (3), 1438–1448. (126) Lou, B.; Hu, S. Different Hydrogen-Bonded Interactions in the Cocrystals of Nicotinamide with Two Aromatic Acids. J. Chem. Crystallogr. 2011, 41 (11), 1663– 1668. (127) Bučar, D. K.; Henry, R. F.; Lou, X.; Duerst, R. W.; Borchardt, T. B.; MacGillivray, L. R.; Zhang, G. G. Z. Co-Crystals of Caffeine and Hydroxy-2-Naphthoic Acids: Unusual Formation of the Carboxylic Acid Dimer in the Presence of a Heterosynthon. Mol.

Pharm. 2007, 4 (3), 339–346.

36 ACS Paragon Plus Environment

Page 37 of 39 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

(128) Tilborg, A.; Springuel, G.; Norberg, B.; Wouters, J.; Leyssens, T. How Cocrystallization Affects Solid-State Tautomerism: Stanozolol Case Study. Cryst.

Growth Des. 2014, 14 (7), 3408–3422. (129) Bučar, D. K.; Henry, R. F.; Duerst, R. W.; Lou, X.; MacGillivray, L. R.; Zhang, G. G. Z. A 1:1 Cocrystal of Caffeine and 2-Hydroxy-1-Naphthoic Acid Obtained via a Slurry Screening Method. J. Chem. Crystallogr. 2010, 40 (11), 933–939. (130) Zong, Y.; Shao, H.; Pang, Y.; Wang, D.; Liu, K.; Wang, L. Multicomponent Hydrogen-Bonding Organic Solids Constructed from 6-Hydroxy-2-Naphthoic Acid and N-Heterocycles: Synthesis, Structural Characterization and Synthon Discussion. J.

Mol. Struct. 2016, 1115, 187–198. (131) Padrela, L.; Rodrigues, M. A.; Tiago, J.; Velaga, S. P.; Matos, H. A.; Azevedo, E. G. de. Tuning Physicochemical Properties of Theophylline by Cocrystallization Using the Supercritical Fluid Enhanced Atomization Technique. J. Supercrit. Fluids 2014, 86, 129–136. (132) Wright, F. E. The Crystallization of Menthol. J. Am. Chem. Soc. 1917, 39 (8), 1515– 1524. (133) Kang, L.; Jun, H. W.; McCall, J. W. Physicochemical Studies of Lidocaine-Menthol Binary Systems for Enhanced Membrane Transport. Int. J. Pharm. 2000, 206 (1–2), 35–42. (134) Corvis, Y.; Négrier, P.; Lazerges, M.; Massip, S.; Léger, J. M.; Espeau, P. Lidocaine/ L -Menthol Binary System: Cocrystallization versus Solid-State Immiscibility. J. Phys.

Chem. B 2010, 114 (16), 5420–5426. (135) Randic, M. On Molecular Identification Numbers. J. Chem. Inf. Comput. Sci. 1984, 24 (3), 164–175. (136) Gramatica, P.; Corradi, M.; Consonni, V. Modelling and Prediction of Soil Sorption

37 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 38 of 39

Coefficients of Non-Ionic Organic Pesticides by Molecular Descriptors. Chemosphere 2000, 41 (5), 763–777. (137) Liu, R.; Sun, H.; So, S.-S. Development of Quantitative Structure−Property Relationship Models for Early ADME Evaluation in Drug Discovery. 2. Blood-Brain Barrier Penetration. J. Chem. Inf. Comput. Sci. 2001, 41 (6), 1623–1632. (138) Hall, L. H.; Kier, L. B. Electrotopological State Indices for Atom Types: A Novel Combination of Electronic, Topological, and Valence State Information. J. Chem. Inf.

Comput. Sci. 1995, 35 (6), 1039–1045. (139) Platts, J. A.; Butina, D.; Abraham, M. H.; Hersey, A. Estimation of Molecular Linear Free Energy Relation Descriptors Using a Group Contribution Approach. J. Chem. Inf.

Comput. Sci. 1999, 39 (5), 835–845. (140) Liu, S.; Cao, C.; Li, Z. Approach to Estimation and Prediction for Normal Boiling Point (NBP) of Alkanes Based on a Novel Molecular Distance-Edge (MDE) Vector, λ.

J. Chem. Inf. Comput. Sci. 1998, 38 (3), 387–394. (141) Jiang, D. X.; Li, Y.; Li, J.; Wang, G. X. Prediction of the Aquatic Toxicity of Phenols to Tetrahymena Pyriformis from Molecular Descriptors. Int. J. Environ. Res. 2011, 5 (4), 923–938. (142) John, A. P.; Uzairu, A.; Shallangwa, G. A. Insilico Design of Highly Potent AntiSalmonella Typhi Drug Candidates from Schiff Bases. Alban. Med. J. 2016, 4, 17–28. (143) Basic, J.; Kalinić, M.; Ivkovic, B.; Eric, S.; Milenkovic, M.; Vladimirov, S.; Vujić, Z. Synthesis, QSAR Analysis and Mechanism of Antybacterial Activity of Simple 2′Hydroxy Chalcones. Dig. J. Nanomater. Biostructures 2014, 9 (4), 1537–1546.

38 ACS Paragon Plus Environment

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

For Table of Contents Use Only Distinguishing of cocrystals from simple eutectic mixtures: phenolic acids as potential pharmaceutical coformers Chair and Department of Physical Chemistry, Faculty of Pharmacy, Collegium Medicum of Bydgoszcz, Nicolaus Copernicus University in Toruń, Kurpińskiego 5, 85-950 Bydgoszcz, Poland

Synopsis: The model comprising 1D and 2D molecular descriptors was proposed and validated for phenolic acids cocrystal screening. This approach based on the methodology involving optimization of regression coefficients for maximization of true positives percentage leads to very valuable clues for predicting cocrystallization with probability about 80%.

39 ACS Paragon Plus Environment