Crystal Engineering of Naturally Occurring Seselin to Obtain a Co

Sciences, University of Karachi, Karachi-75270, Pakistan. ‡. Department of Biochemistry, Faculty of Science, King Abdulaziz University, Jeddah-21412...
0 downloads 0 Views 1MB Size
Subscriber access provided by - Access paid by the | UCSB Libraries

Article

Crystal Engineering of Naturally Occurring Seselin to Obtain a Co-crystal with Enhanced Antileishmanial Activity, Hirshfeld Surface Analysis and Computational Insight Ejaz Hussain, Rajesh Kumar, M. Iqbal Choudhary, and Sammer Yousuf Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.8b00602 • Publication Date (Web): 05 Jul 2018 Downloaded from http://pubs.acs.org on July 9, 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 24 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

Crystal Engineering of Naturally Occurring Seselin to Obtain Co-crystal with Enhanced Antileishmanial Activity, Hirshfeld Surface Analysis and Computational Insight

Ejaz Hussain †, Rajesh Kumar †, M.Iqbal Choudhary †‡, SammerYousuf †*.



H.E.J. Research Institute of Chemistry, International Center for Chemical and Biological

Sciences, University of Karachi, Karachi-75270, Pakistan. ‡

Department of Biochemistry, Faculty of Science, King Abdulaziz University, Jeddah-21412,

Saudi Arabia.

*Correspondence author: H. E. J. Research Institute of Chemistry, International Center for Chemical and Biological Sciences, University of Karachi, Karachi-75270, Pakistan.

Tel: +92-21-34824924-5

Fax: +92-21-4819018-9

E-mail address: [email protected]

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

ABSTRACT Seselin (8,8-dimethyl-2H,8Hbenzo[1,2-b:3,4-b']dipyran-2-one) is a naturally occurring pyrano coumarin, isolated from Seseli diffusum (Roxb. ex Sm.) Sant. & Wagh.We have synthesized cocrystal of seselin with thiourea by using liquid-assisted grinding and solution methods (via refluxing). The purity and homogeneity of seselin-thiourea (1:1) co-crystal was confirmed by analysis based on single-crystal X-ray diffraction technique. A deep insight into the geometry of co-crystal demonstrated that 1:1 co-crystal stoichiometry is sustained by N-H…O hydrogen bonding between amine (-NH2) groups of thiourea and carbonyl group of seselin. Synthesized co-crystal exhibited good anti-leishmanial activity invitro (IC50 = 13.2 + 1.2 µg/mL) in comparison to seselin (IC50 = 29.4 + 1.0 µg/mL) and thiourea (IC50 = 86.06+0.6µg/mL), against the tested standards, pentamidine (IC50 = 5.09 + 0.09 µg/mL) and amphotericin B (IC50 = 0.12 + 0.15µg/mL). Geometrical values of computationally optimized structure were calculated by density functional theory (DFT) at B3LYP/6-31(d,p) to obtain HOMOs and LUMOs of cocrystal, and found to be consistent with XRD data. Hence, this study further highlights the significance of co-crystallization towards improving the efficacy of medicinal and industrial materials. Furthermore, the consistency of computational results with crystallographic values

2 ACS Paragon Plus Environment

Page 2 of 24

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

ensured the computational modeling of co-crystal is equally important to predict the stable geometries, and next-generation applications of co-crystals. Key Words: Co-crystal, Seselin, Sesseli diffusium, Leishmaniasis, Single-crystal X-ray diffraction, Density functional theory (DFT), Hirshfeld surface analysis.

1. INTRODUCTION 2. Co-crystals represent a class of multi-component compounds having a significant effect on physiochemical and pharmaceutical properties1. The non-covalent interactions, predominantly hydrogen bonding2, halogen bonding3,4or heterosynthones5 (in supramolecularco-crystals) trigger the assembly of the components in a fixed stoichiometry at ambient conditions. Crystal engineering approaches via co-crystallization are mean to alter the physiochemical properties of active pharmaceutical ingredients (APIs) to enhance their stability, efficacy, and bioavailability6,7. In fact, co-crystallization is used extensively in industries to associate APIs with excipient (s) for designing effective drugs dosage forms. Their application ranges from drug designing to the analysis of active pharmaceutical ingredients (APIs)8. Co-crystals have significant difference in physiochemical and biological properties from their individual components. Indeed, the purpose of co-crystallization is to obtain superior physiochemical and biological properties from the existing APIs, while the original properties of APIs should sustaine9. In this context, numerous pharmaceutical key drugs were crystallized with suitable coformers to fine-tune the desired properties10,11. Many crystal engineering approaches as analytical techniques have been introduced for the development of co-crystals12. An antitumor pro-drug temozolmide (TZM) that has broad spectrum anti-neoplastic activity was co-crystallized

with

amide

coformers

(nicotineamide,

p-hydroxybenzamide,

pyrazinamide, saccharine, and caffeine) to obtain stable co-crystals in alkaline medium13. Nalidixic acid is a quinolone antibiotic used in the treatment of urinary tract infection. The co-crystals of nalidixic acid were developed to improve the bioavailability due to replacement of intermolecular drug-drug interaction in polymorph with strong drug3 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

coformer hydrogen bonds in co-crystals14. The limited in-vivo efficacy and poor bioavailability of quercetin was improved by co-crystallization with caffeine and isonicotinamide as conformers15. Along with their pharmaceutical importance, co-crystals are equally important for sensing applications. Co-crystal of 10-methylphenthiazine and 1,3-dinitrobenzene was adequately utilized for the optical sensing of TNT based explosives16. Micro- and nano-crystalline co-crystals exhibit single-crystal-to-singlecrystal chemical reactivity that is promising for sensors, and high-density data storage devices17. Currently co-crystallization is now a preferred technique to improve dissolution, solubility, and bioavailability of less water-soluble drugs18-20. Carbamazepine co-crystals with nicotinamide (NCT), malonic acid (MLN), glutaric acid (GTA), saccharin (SAC), anhydrous oxalic acid (OXA), succinic acid (SUC), and salicylic acid (SLC) exhibit

152 times greater solubility than the stable carbamazepine dihydride form21. Hence, the co-crystallization has been an adequate subject for the researchers from many past years because of their use in many industries which includes pharmaceuticals, paper, textile, photographic, chemical processing, and electronics22. There are many commercially important compounds that can be co-crystallized with suitable coformers to fine-tune the stable, and superior properties for next-generation applications. The US-FDA approval for co-crystals to be considered as drug intermediates has opened new vistas, both in pharmaceutical and academic research sector. The article entitled, “Reflection paper” by the European medicines agency has discussed in detail the regularity implementation of the co-crystals and clearly mention that any co-crystal will not be consider as a new drug if its safety and efficacy found to be equivalent to the active pharmaceutical ingredient of the considered co-crystal23. The parasitic infection leishmaniasis is transmitted through bite of infected sand fly (Phlebotomus or Lutzomyia). The infection is responsible for more than 2 million new cases, and 70,000 deaths annually along with 350 million people at risk in various regions of world, predominantly in 89 countries of American, Asian, and African subcontinents, and Mediterranean region24. The severity of infection depends upon the type and the penetrating ability of the parasite ranging from facial and cutaneous to subcutaneous and the most severe visceral leishmaniasis25. According to a worldwide survey of the WHO to identify the world-wide regions with distribution and major burdens of various types of 4 ACS Paragon Plus Environment

Page 4 of 24

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

leishmaniasis in 102 countries, Pakistan has identified among those having high burden of both cutaneous and viral leishmaniasis (World Health Organization: Weekly Epidemiological Record (WER). 2016). Antimonial drugs are the most effective and remain the only primary treatment in various effected regions of world, followed by amphotericin B, and pentamidine as second line of treatment 26. However, the severe side effects and development of drug resistance due to complex immune response of the causative agents open new vistas for the search of alternatives. That has been resulted in identification and approval of miltefosine as the only oral treatment for leishmaniasis27. The low efficacy, high price and increased drug resistance developed in last 10 to 20 years against amphotericin B, pentamidine not only increasing the economic burden but also demand a vigorous search towards more effective treatment from available resources including locally available natural products, and through repurposing of existing drugs. Seselin

(8,8-dimethyl-2H,8H

benzo[1,2-b:3,4-b']dipyran-2-one)

is

a

natural

pyranocoumarin, isolated from Seseli diffusem(Roxb. ex Sm.) Sant. & Wagh.Seselin is widely important due to its biological activities, including anti-oxidant and antispasmodic properties28. The photo-responsive coumarin-based multifunctional nanoparticles were used for the controlled delivery of anticancer drug chlorambucil29. Improvement of coumarin derivatives via co-crystallization might open new avenues for next generation applications. For instance,coumarin derivatives bear luminescent properties such as coumarin-3-carboxylic acid were modified by co-crystallization with different conformers30. Theoretical insight to co-crystals by means of computational modeling enables to predict the geometry,stability, andbioavailability31. The stability of 1,3,5trinitrobenzene assembled by halogen bonding were evaluated computationally, and cocrystal forms exhibited strong stability than itsindividualcomponents32. The binding energies even associated with weak interaction were evaluated for some of coumarin derivatives33. In fact,computational predictions helped in the targeted synthesis of cocrystals for next-generation applications. By keeping above facts in mind, co-crystals of naturally occurring seselin with thiourea (Figure 1) was synthesized sucessfully. The main objectives were to assimilate the smaller molecules (thiourea in our case) in heterocyclic solid (seselin) that results in 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 24

alteration the crystalline packing, to modify the biological properties, and to study the function of non-covalent forces that are responsible for these modifications. The 1:1 stoichiometry of co-crystal was characterized based on single-crystal X-ray diffraction analysis. The in vitro antileishmanial activity against L. major (causative agent of cutaneous leishmaniasis) promastegoates was evaluated and results showed improvement in co-crystal form as compared to seselin and thiourea. Therefore, as per approved guidelines of US-FDA the synthesized co-crystal is a potential lead for future studies against cutaneous leishmaniasis. The stable geometrical parameters were evaluated by density functional theory (DFT). The computed non-covalent interactions were in good agreement with those obtained by diffraction values. Hirshfeld surface analysisis, a powerful technique to determine intermolecular interactions qualitatively and quantitatively, was also successfully employed to quantitatively analyze the role of various interactions towards stability of arrangements in crystal lattice33,34.

Figure1. Structure of seselin-thiourea co-crystal in 1:1 stoichiometry. 3. EXPERIMENTAL SECTION 3.1.

Materials:

Commercially available solvents and chemicals were utilized without any further purification. Seselin-(8,8-dimethyl-2H,8H-benzo[1,2-b:3,4-b']dipyran-2-one)

was

isolated

from

Seseli

diffusum(Roxb. ex Sm.) Sant. & Wagh,as described by Ahmed and co-workers28.Thiourea was purchased from Merck, Germany (Index No: 612-082-0000). A small sized pestle, and mortar was used for grinding. 3.2.

Synthesis of seselin : thiourea co-crystal:

6 ACS Paragon Plus Environment

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

Seselin (8,8-dimethyl-2H,8H-benzo[1,2-b:3,4-b']dipyran-2-one)was successfully co-crystallized with thiourea by using both solution reflux and grindingmethods, in the stoichiometric ratio (1:1). After dissolving the grinded mixture in methanol, left on slow evaporation at ambient temperature and pressure to obtain homogenous crystals. a. Solution Crystallization: 68.8 mg (0.3014 m. mol) seselin, and 24.5 mg (0.32 m. mol) thiourea in 1:1 molar ratio were allowed to dissolve in 5 mL methanol solution and allowed to reflux at 70°C for 3hrs. the solution kept for evaporation at RT for 2 days in order to obtain fine and shiny co-crystals. b. Grinding Method: 68.8 mg (0.3014 m. mol) seselin, and 24.5 mg (0.32 m. mol) thiourea in 1:1 molar ratio were carried out through mechanical neat grinding for 1hrin a mortar pestal, followed by methanol(2 mL) assisted grinding for 30 minutes. Resulting mixture was then dissolved in 10 mL methanol and kept on slow evaporation for 2 days to obtain suitable co-crystals for single-crystal X-ray diffraction analysis.

3.3.

Single Crystal X-ray Diffraction:

Using a compound microscope, suitable crystalof synthesized co-crystal(0.81 x 0.26 x 0.07 mm) wasselected andthen mounted on a copper pin by mean of epoxy glue for data collection. Single crystal X-ray diffraction data was collected on Bruker Smart APEX II, fitted with CCD detector diffractometer,installed with Bruker AXS kryoflex open flow cryostate(Mo Kα (0.71073 Å) radiations).Data integration, and reductionwere performed by SAINT program. The structure of compound was solved by utilizing direct method, and Fourier transformation techniques, respectively by using SHELXL97 program.Structureswere refined by using full-matrix leastsquare calculation on F2. All non-hydrogen atoms were refined anisotropic thermal parameters and were refined at geometrically idealized positions. SHELXL and PLATON were used for final refinement and various statistics. Using the ORTEP program, 3D view of molecules,and inter-molecular interaction with crystal packing were produced. 2.4. Hirshfeld surface analysis: Hirshfeld and their functions were determined by incorporating crystal explorer35, using standard settings from the calculated crystal structure parameters.Furthermore, 2D fingerprint plot was generated for visualizing percentage contribution of each contact. dnorm, di, andde surface 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

generated on the surface scale of -0.404__1.424 Å, 0.849__2.552 Å, and 0.848__2.597Å for dnorm, di, andde, respectively. 2.5.

Biological Screening of Co-Crystal and Seselin:

In-vitro anti-leishmanial activity against major promastigotes was evaluated by using the same procedure as used by36. 2.6.

Computational Modelling:

The geometry optimization of co-crystal was performed by density functional theory (DFT) at B3LYP/6-31(d,p) by using Gaussian 09W at super computer37. 3. RESULTS AND DISCUSIONS 3.5.

Synthesis of Co-crystal

Despite the variation of methods, our targeted co-crystal was successfully prepared by adopting two simple approaches. Neat mechanical grinding followed by methanol assisted grinding and refluxing of equiv-molar amounts of components i.e. seselin and thiourea (1:1). Both approaches were ended up to left methanolic solution of co-crystal for slow evaporation to obtain suitable crystals for X-ray diffraction analysis. The synthesis scheme for co-crystal has been illustrated in (Figure 2).

Figure2. Schematic representation of co-crystal synthesis.

8 ACS Paragon Plus Environment

Page 8 of 24

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

3.6.

Crystal Structure Analysis

Figure 3. ORTEP drawing of co-crystal of seselin, and thiourea [(C14H12O3) (CH4N2S)] drawn at 50% probability level. Table-1: The summary of crystal data, data collection and structure refinement parameters for synthesized co-crystal.

Crystal Parameters

Co-crystal

Empirical formula

C15H16N2O3S

Formula weight

304.36

Temperature

273(2) K

Wavelength

0.71073 Å

Crystal system

Monoclinic,

Space group

P2(1)/c

Unit cell dimensions

a = 7.1579(16) Å, b = 8.1507(19) Å, c = 24.808(6) Å, β = 96.128(5)°

Volume

1439.1(6)Å3

Z

4

Calculated density

1.405 mg/m3

Absorption coefficient

0.237 mm-1

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

F(000)

640

Crystal size

0.816 x 0.266 x 0.078mm

Theta range for data collection

1.651 to 25.50°

Limiting indices

h = -88, k = -99, l = -2330

Reflections Collected/ Unique

8282/2594

(Rint)

0.045

Completeness to theta = 25.49°

96.9 %

Refinement Method

Full-matrix least-squares on F2

Data / restraints / parameters

2594 / 0 / 206

Final R indices [I > 2s(I)]

R1 = 0.0388, wR2 = 0.0893

R indices (all data) Goodness of fit on F2

R1 = 0.0527, wR2 = 0.0975 1.040

Largest difference. Peak and hole.

0.297 and -0.208 e.Å-3

Structural Features: The structure of the co-crystal found to be composed of a pyranocoumarin seselin and thiourea moietyadhere in crystal lattice viaN1---H1...O3, and N2---H3...O3 noncovalent interaction(Figure. 4). The dihedral angle between coumarin ring (O1-O2/C1-C9)and six member nonaromatic pyran ring (C7/O2/C8-C11) was found to be 9.48(7)°with maximum deviation of 0.260(2)° for C-8. The bond length of 1.698(2) Å was found for S1-C14. Carbonyl group C1-O1 has bond length of 1.215(3) Å. The torsional angle of 168.31(16)° and 178.15(19)° were found for C7-O3-C12-C14 and O1-C1 -C2-C3.The summary of crystal data, data collection, and structure refinement parameters are summarized in Table-1. Crystal Packing: An overview of crystal packing of co-crystal revealed that incorporated thiourea moiety in the lattice of co crystal interacted with carbonyl functionality of seselin molecules via protons of amino group (N-H…O) without bringing any configurational changes. In the crystal lattice, the co-crystal moieties arranged in two-dimensional pattern that further stabilized by π-π interaction {(Cg1···Cg1)(O2-C1-C2-C3-C4-C9)} having minimum centroidcentroid distance of 3.691Å. Each thiourea moiety found to be connected to two neighboring thiourea molecules via N1-H1B…S1, and N2-H2B…S1 interactions to form S8 graph set repeating motif that extended to from chains running along b-axis (Figure 4). The onedimensional chains of thiourea moieties further linked with neighboring seselin molecules to form two dimensional network via N1-H1A…O1, and N2-H2B…O1 interactions, with hydrogen 10 ACS Paragon Plus Environment

Page 10 of 24

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

bond distances of 2.25(2) Å and 2.21(2) Å, respectively, summarized in Table2. The existence of glide

plane

symmetry

can

be

clearly

observed

in

the

Figure

4.

Figure 4. Crystal packing diagram of co-crystals [(C14H12O3) (CH4N2S)]

Table-2: Hydrogen Bond Distances in Å. D

H

A

N1

H1A H1B

N1

D-H

H...A

D...A

D-H…A

O1a

0.85(2)

2.25(2)

3.016(3)

150(2)

b

0.89(2)

3.310(2)

3.326(2)

169.7(19)

0.87(2)

2.21(2)

3.009(3)

153(2)

0.87(2)

2.55(2)

3.407(2)

168.7(19)

S1

N2

H2A

O1

N2

H2B

S1c

a

Symmetry codes #a: -1-x,-y,1-z #b: -x,-1/2+y,1/2-z #c: -x,1/2+y,1/2-z

3.7.

Antileishmanial Activity In vitro:

In vitro antileishmanial activities of the pure seselin, thiourea and its co-crystal were investigated against the L. major promastigoates. A significant enhancement in biological activity of co-crystal (IC50 = 13.2 + 1.2 µg/mL) in comparison to that of seselin (IC50 = 29.4

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 24

+ 1.0 µg/mL) and thiourea (IC50 = 86.06+0.6µg/mL) against the tested standard pentamidine (IC50 = 5.09+0.09µg/mL) and amphotericin B (IC50 = 0.12 +0.15 µg/mL) was observed. Significant increase in anti-leishmanial activity clearly demonstrate the synergic effect of both seselin and thiourea moieties. However, detailed in vivo and mechanistic studies are required to study the mode of action. The results of the biological assays are summarized in the Table-3. .

Table-3: Biological activities of co-crystals. Sample

Leishmanicidal activity (IC50(µ µg/ml))

Seselin

29.4+1.0

Thiourea

86.06 + 0.6

Co-crystal (1)

13.2+1.2

Standard drugs

3.8.

Amphotericin B

0.29+0.05

Pentamidine

5.09+0.09

Computational Insights:

Molecular structure of co-crystal was evaluated by using density functional theory (DFT) at B3LYP/6-31 G(d,p) level. The stable molecular geometry of co crystal obtained by computational modeling is illustrated in Figure 5. The hydrogen bonding data for the components of co-crystal was evaluated theoretically. The theoretical studies supported two hydrogen bonding (N1—H1A…O1, and N2—H2A…O2) responsible for the stabilization of the co-crystal geometry with bond distance of 2.28 Å and 2.29(2) Å, respectively. Theoretical

12 ACS Paragon Plus Environment

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

hydrogen bond values were found consistent with experimental values of XRD analysis (Table 4).

Figure 5. Molecular geometry of co-crystal obtained by computational modeling. Table 5: Hydrogen bonding data obtained by computational simulation D

H

A

H...A (Theoretical)

H...A (Exp)

D-H…A Theoretical

D-H…A (Exp)

N1

H1A

O1

2.28

2.25(2)

150

151

N2

H2A

O1

2.29

2.21(2)

154.24

154

The HOMO and LUMO molecular orbitals were depicted by density functional theory (DFT) at B3LYP/6-31(d,p) are presented in Figure 6. Highest occupied molecular orbitals (HOMOs) were localized on thiourea moiety that act as hydrogen bond donor (N1-H1A, N2-H2B). While LUMOs extended along the seselin skeleton acts as hydrogen bond acceptor i.e. O1as illustrated in Figure 6.

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

Figure 6.HOMO (left) and LUMO (right) of co-crystal depicted by DFT. Hirshfeld Surface Analysis:

Figure7. Hirshfeld surface generated on different parameters de, and di for co-crystal

14 ACS Paragon Plus Environment

Page 14 of 24

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

Figure8. Hirshfeld surface generated on dnorm parameters for co-crystal.

Figure9. Interaction of co crystal with the generated Hirshfeld surface in co-crystal.

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

Figure10. (a) Packing diagram of Hirshfeld surface with neighboring molecules (b) close packing diagram of Hirshfeld surface for co-crystal. The comparison of surface mapped on the values of dnorm with the di and de reveals the changing of contacts formed with the Hirshfeld surface (Figure 7). Highlighted area over di and de surfaces clearly shows the contacts; bright red spots vary due to the interacting atoms from the outside and inside of the surfaces. Hirshfeld surface mapped on dnorm value of compound, shows the area highlighted with bright red spots, where close contacts were formed34. Figure9 shows the red spots indicating the close contacts with Hirshfeld surface due to two hydrogen bonds O---H, and S---H, in which N2---H2B…O1 interaction is the strongest.

Figure11. (a) Shape index (b) Curvedness mapped over Hirshfeld surface for co-crystal. 16 ACS Paragon Plus Environment

Page 16 of 24

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

Two additional properties mapped over Hirshfeld surface reveals the probability of interaction with the neighboring molecules and planar stacking. Blue and red adjacent triangles indicate the π--π interactions, However, in our case highlighted area shows that these triangles are not clear due to the weak π--π interactions.Encircled area at one side indicates the π-hole. Curvedness surface indicates the location of surface depending upon the probability of interaction with neighboring molecules34. Packing diagram of Hirshfeld surface reveals about the packing pattern, Figure 10(a) with the neighboring molecules while Figure 10(b) shows the closed packing, depending upon the clearance of the surface packing only one-sided interactions have been generated. Figure12.1,and12.2 show the spikes are of close contacts over the Hirshfeld surface, and the highlighted blue color area, and length of spikes are due to the strength of contacts present in crystal packing. FPs indicates that most of the contributions over the Hirshfeld surface are due to H---H contact (45.2%)as the molecular surface is comprised of a sea of H atoms. Large number of H-atoms over surface supported the C-H contact to be the second one having 14.8%. Beside these interactions hydrogen bond contacts, on which the interactions are totally dependent are SH and O-H having percentages of 14.3%, and 12%, respectively38. O-H is the main contact on which crystal packing is extended further but it is contributing less than S-H. Possible reason is that S atom is contributing with Hirshfeld surface from two different sides, while O atom has only one site. Among all contacts that were contributing in crystal packing, other interaction are O-C, N-H, C-C, O-O, and O-S having 4.7%, 3.8%, 3.4%, 1.2%, and 0.6%, respectively.

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 24

__

Figure12.1 . 2D fingerprint plots according to dnorm value (-0.404 1.424 Å) and Hirshfeld surface view about these interactions, having (a), (b), (b), (d) and (e) for total, H-H, C-H, S-H and O-H interactions in co-crystal.

__

Figure12.2 . 2D fingerprint plots according to dnorm value (-0.404 1.424 Å) and Hirshfeld surface view about these interactions, having (f), (g), (h), (i) and (j) for total, O-C, N-H, C-C, OO, and O-S interactions in co-crystal.

18 ACS Paragon Plus Environment

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

Electrostatic potential:

Figure13.Electrostatic potential mapped over hirshfeld surface for co-crystal.

Electrostatic property (Figure13), mapped over Hirshfeld surface39, reveals the positive and negative end of the molecules, involved in interaction with neighboring molecules. Diagram has been extended to the neighboring molecules to elaborate the electrostatic potential. Red color represents negative potential while the blue represents the positive potential. The highlighted areas in circles represent the negative potential, which are because of sulfur atom, making S---H contact with the positive part. The rectangular highlighted area is due to oxygen of carbonyl moiety bonded viaN---H…O contact, here H atom in blue color showing positive potential.

4. CONCLUSION In conclusion co-crystal of naturally occurring seselin with thiourea has synthesized successfully. The antileishmanial activity of synsthesized co-crystal was found to be several folds more than the seselin and thiourea individually. The obtained promising results have clearly demonstrated the synergic effects of the carbonyl and amine functionalities of two molecules and promoted the ability of synthesized co-crystal as a new lead for further studies against cutaneous leishmaniasis. Furthermore, via crystallographic studies we were able to study the crystal parameters, configurations and packing of molecules in crystal lattice. The non-covalent forces involved in molecular geometry were also evaluated. The computational studies indicated the stable geometries as well as HOMO and LUMO of co-crystal. Hirshfeld surface analysis further 19 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

elaborate the qualitative and quantitative contributions of inter-molecular interactions towards crystal packing, and in addition to this electrostatic potential study showed the potential sites of the molecule for further structural modifications. The presented study further emphasizes the need of co-crystallization of the natural compounds with suitable conformers to enhance their medicinal properties.

5. CONFLICT OF INTEREST The authors has declared no competing conflict of financial interest.

REFERENCES 1. Sanphui, P.; Goud, N. R.; Khandavilli, U. R.; Nangia, A., Fast dissolving curcumin cocrystals. Cryst. Growth Des.2011,11, 4135-4145. 2. Sanphui, P.; Devi, V. K.; Clara, D.; Malviya, N.; Ganguly, S.; Desiraju, G. R., Cocrystals of hydrochlorothiazide: solubility and diffusion/permeability enhancements through drug– coformer interactions. Mol. Pharmaceutics2015,12, 1615-1622. 3. Goroff, N. S.; Curtis, S. M.; Webb, J. A.; Fowler, F. W.; Lauher, J. W., Designed cocrystals based on the pyridine− iodoalkyne halogen bond. Org. Lett.2005,7, 1891-1893. 4. Baldrighi, M.; Cavallo, G.; Chierotti, M. R.; Gobetto, R.; Metrangolo, P.; Pilati, T.; Resnati, G.; Terraneo, G., Halogen bonding and pharmaceutical cocrystals: the case of a widely used preservative. Mol. Pharmaceutics2013,10, 1760-1772. 5. Goud, N. R.; Babu, N. J.; Nangia, A., Sulfonamide− Pyridine-N-oxide Cocrystals. Cryst. Growth Des.2011,11, 1930-1939. 6. Sanphui, P.; Kumar, S. S.; Nangia, A., Pharmaceutical cocrystals of niclosamide. Cryst. Growth Des.2012,12, 4588-4599. 7. Wang, L.; Tan, B.; Zhang, H.; Deng, Z., Pharmaceutical cocrystals of diflunisal with nicotinamide or isonicotinamide. Org. Process Res. Dev.2013,17, 1413-1418. 8. Schultheiss, N.; Newman, A., Pharmaceutical cocrystals and their physicochemical properties. Cryst. Growth Des.2009,9, 2950-2967. 9. Surov, A. O.; Voronin, A. P.; Manin, A. N.; Manin, N. G.; Kuzmina, L. G.; Churakov, A. V.; Perlovich, G. L., Pharmaceutical cocrystals of diflunisal and diclofenac with theophylline. Mol. Pharmaceutics2014,11, 3707-3715. 10. Chen, J.-M.; Li, S.; Lu, T.-B., Pharmaceutical cocrystals of ribavirin with reduced release rates. Cryst. Growth Des.2014,14, 6399-6408. 11. E. Castro, R. A.; Ribeiro, J. D.; Maria, T. M.; Ramos Silva, M.; Yuste-Vivas, C.; Canotilho, J.; Eusébio, M. E. S., Naproxen cocrystals with pyridinecarboxamide isomers. Cryst. Growth Des.2011,11, 5396-5404.

20 ACS Paragon Plus Environment

Page 20 of 24

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

12. Charron, D. M.; Ajito, K.; Kim, J.-Y.; Ueno, Y., Chemical mapping of pharmaceutical cocrystals using terahertz spectroscopic imaging. Anal. Chem.2013,85, 1980-1984. 13. Sanphui, P.; Babu, N. J.; Nangia, A., Temozolomide cocrystals with carboxamide coformers. Cryst. Growth Des.2013,13, 2208-2219. 14. Gangavaram, S.; Raghavender, S.; Sanphui, P.; Pal, S.; Manjunatha, S. G.; Nambiar, S.; Nangia, A., Polymorphs and cocrystals of nalidixic acid. Cryst. Growth Des.2012,12, 4963-4971. 15. Smith, A. J.; Kavuru, P.; Wojtas, L.; Zaworotko, M. J.; Shytle, R. D., Cocrystals of quercetin with improved solubility and oral bioavailability. Mol. Pharmaceutics2011,8, 18671876. 16. McNeil, S. K.; Kelley, S. P.; Beg, C.; Cook, H.; Rogers, R. D.; Nikles, D. E., Cocrystals of 10-Methylphenthiazine and 1, 3-Dinitrobenzene: Implications for the Optical Sensing of TNTBased Explosives. ACS Appl. Mater. Interfaces2013,5, 7647-7653. 17. Bucar, D.-K.; MacGillivray, L. R., Preparation and reactivity of nanocrystalline cocrystals formed via sonocrystallization. J. Am. Chem. Soc.2007,129, 32-33. 18. Yan, Y.; Chen, J.-M.; Geng, N.; Lu, T.-B., Improving the solubility of agomelatine via cocrystals. Cryst. Growth Des.2012,12, 2226-2233. 19. Alhalaweh, A.; Sokolowski, A.; Rodríguez-Hornedo, N. r.; Velaga, S. P., Solubility behavior and solution chemistry of indomethacin cocrystals in organic solvents. Cryst. Growth Des.2011,11, 3923-3929. 20. Alhalaweh, A.; Roy, L.; Rodríguez-Hornedo, N.; Velaga, S. P., pH-dependent solubility of indomethacin–saccharin and carbamazepine–saccharin cocrystals in aqueous media. Mol. Pharmaceutics2012,9, 2605-2612. 21. Good, D. J.; Rodríguez-Hornedo, N. r., Solubility advantage of pharmaceutical cocrystals. Cryst. Growth Des.2009,9, 2252-2264. 22. Stahly, G. P., A survey of cocrystals reported prior to 2000. Cryst. Growth Des.2009,9, 4212-4229. 23. Gadade, D. D.; Pekamwar, S. S., Pharmaceutical Cocrystals: regulatory and strategic aspects, design and development.Adv. Pharm. Bull.2016,6, 479. 24. Torres-Guerrero, E.; Quintanilla-Cedillo, M. R.; Ruiz-Esmenjaud, J.; Arenas, R., Leishmaniasis: a review. F1000Research 2017,6. 25. Andrade-Narváez, F. J.; Vargas-González, A.; Canto-Lara, S. B.; Damián-Centeno, A. G., Clinical picture of cutaneous leishmaniases due to Leishmania (Leishmania) mexicana in the Yucatan peninsula, Mexico. Mem. Inst. Oswaldo Cruz2001,96, 163-167. 26. Haldar, A. K.; Sen, P.; Roy, S., Use of antimony in the treatment of leishmaniasis: current status and future directions. Mol. Biol. Int.2011,2011. 27. Ponte-Sucre, A.; Gamarro, F.; Dujardin, J.-C.; Barrett, M. P.; López-Vélez, R.; GarcíaHernández, R.; Pountain, A. W.; Mwenechanya, R.; Papadopoulou, B., Drug resistance and treatment failure in leishmaniasis: A 21st century challenge. PLoS Neglected Trop. Dis.2017,11 , e0006052. 28. Abbaskhan, A.; Choudhary, M. I.; Ghayur, M. N.; Parween, Z.; Shaheen, F.; Gilani, A. u. H.; Maruyama, T.; Iqbal, K.; Tsuda, Y., Biological Activities of Indian Celery, Seseli diffusum (Roxb. ex Sm.) Sant. & Wagh. Phytother. Res.2012,26, 783-786.

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

29. Karthik, S.; Puvvada, N.; Kumar, B. P.; Rajput, S.; Pathak, A.; Mandal, M.; Singh, N. P., Photoresponsive coumarin-tethered multifunctional magnetic nanoparticles for release of anticancer drug. ACS Appl. Mater. Interfaces2013,5, 5232-5238. 30. Yan, D.; Delori, A.; Lloyd, G. O.; Patel, B.; Friščić, T.; Day, G. M.; Bučar, D.-K.; Jones, W.; Lu, J.; Wei, M., Modification of luminescent properties of a coumarin derivative by formation of multi-component crystals. CrystEngComm 2012,14, 5121-5123. 31. Song, K.-p.; Zhang, S.-h.; Shi, W.-j., Theoretical insights into the stabilities, detonation performance, and electrostatic potentials of cocrystals containing α-or β-HMX and TATB, FOX7, NTO, or DMF in various molar ratios. J. Mol. Model.2016,22, 249. 32. Bennion, J. C.; Vogt, L.; Tuckerman, M. E.; Matzger, A. J., Isostructural cocrystals of 1, 3, 5-trinitrobenzene assembled by halogen bonding. Cryst. Growth Des.2016,16, 4688-4693. 33. Seth, S. K.; Maity, G. C.; Kar, T., Structural elucidation, Hirshfeld surface analysis and quantum mechanical study of para-nitro benzylidene methyl arjunolate. J. Mol. Struct.2011,1000, 120-126. 34. Spackman, M. A.; Jayatilaka, D., Hirshfeld surface analysis. CrystEngComm 2009,11, 19-32. 35. Wolff, S.; Grimwood, D.; McKinnon, J.; Turner, M.; Jayatilaka, D.; Spackman, M., Crystal explorer. The University of Western Australia Perth, Australia: 2012. 36. Choudhary, M. I.; Yousuf, S.; Shah, S. A. A.; Ahmed, S., Biotransformation of physalin H and leishmanicidal activity of its transformed products. Chem. Pharm. Bull.2006,54, 927-930. 37. Mohamed, S.; Tocher, D. A.; Price, S. L., Computational prediction of salt and cocrystal structures—Does a proton position matter? Int. J. Pharm.2011,418, 187-198. 38. McKinnon, J. J.; Jayatilaka, D.; Spackman, M. A., Towards quantitative analysis of intermolecular interactions with Hirshfeld surfaces. Chem. Commun.2007, (37), 3814-3816. 39. Spackman, M. A.; McKinnon, J. J.; Jayatilaka, D., Electrostatic potentials mapped on Hirshfeld surfaces provide direct insight into intermolecular interactions in crystals. CrystEngComm 2008,10, 377-388. 40. SMART, A. I., Data Collection Software, version 2.1. Bruker AXS Inc.: Madison, WI 2005. 41. Bruker, S.; SMART, B. A., Inc., Madison, Wisconsin, USA, 2001;(b) GM Sheldrick. Acta Crystallogr., Sect. A: Found. Crystallogr 2008,64, 112. 42. Sheldrick, G. M., Crystal structure refinement with SHELXL. Acta Crystallogr., Sect. C: Struct. Chem.2015,71, 3-8. 43. Sheldrick, G., SHELXL-97; University of Göttingen: Göttingen, Germany, 1997. There is no corresponding record for this reference 2010. 44. Farrugia, L. J., ORTEP-3 for Windows-a version of ORTEP-III with a Graphical User Interface (GUI).J. Appl. Crystallogr.1997,30, 565-565.

22 ACS Paragon Plus Environment

Page 22 of 24

Page 23 of 24 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 Content Use Only

Crystal Engineering of Naturally Occurring Seselin to Obtain Co-crystal with Enhanced Antileishmanial Activity, Hirshfeld Surface Analysis and Computational Insight Ejaz Hussain †, Rajesh Kumar †, M. Iqbal Choudhary †‡, SammerYousuf †*. †

H.E.J. Research Institute of Chemistry, International Center for Chemical and Biological Sciences, University of Karachi, Karachi-75270, Pakistan. ‡

Department of Biochemistry, Faculty of Science, King Abdulaziz University, Jeddah-21412, Saudi Arabia.

*Correspondence author: H. E. J. Research Institute of Chemistry, International Center for Chemical and Biological Sciences, University of Karachi, Karachi-75270, Pakistan. Tel: +92-21-34824924-5 Fax: +92-21-4819018-9 E-mail address: [email protected]

Current study demonstrated the synthesis, structural studies and anti-leishmanial activity evaluation of the co-crystal of naturally occurring seslin with thiourea. Synthesized co-crystal showed improved anti-leishmanial activity against promastegoates of L. major responsible for cutaneous leishmaniasis. DFT was employed to obtain HOMOs and LUMOs of co-crystal followed by Hirshfeld surface analysis to quantitively analyze various intermolecular interactions involved in crystal packing.

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

24 ACS Paragon Plus Environment

Page 24 of 24