Modeling Olanzapine Solution Growth Morphologies

Modeling Olanzapine Solution Growth. Morphologies. Yuanyuan Sun. 1,2. , Carl J. Tilbury. 2. , Susan M. Reutzel-Edens. 3. , Rajni M. Bhardwaj. 3. , Jin...
1 downloads 0 Views 1MB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

Article

Modeling Olanzapine Solution Growth Morphologies yuanyuan sun, Carl J. Tilbury, Susan M. Reutzel-Edens, Rajni M. Bhardwaj, Jinjin Li, and Michael F. Doherty Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.7b01389 • Publication Date (Web): 12 Dec 2017 Downloaded from http://pubs.acs.org on December 14, 2017

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

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

Page 1 of 28 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

Modeling Olanzapine Solution Growth Morphologies Yuanyuan Sun1,2, Carl J. Tilbury2, Susan M. Reutzel-Edens3, Rajni M. Bhardwaj3, Jinjin Li1,*, & Michael F. Doherty 2,* 1

National Key Laboratory of Science and Technology on Micro/Nano Fabrication, Department of

Micro/Nano Electronics, Shanghai Jiao Tong University, Shanghai 200240, China. 2

Department of Chemical Engineering, University of California, Santa Barbara, CA 93106, United States.

3

Small Molecule Design & Development, Eli Lilly and Company, Indianapolis, Indiana 46285, United

States.

Correspondence

and

requests

for

materials

should

be

addressed

to

J.L.

([email protected]) or to M. F. D. (mfd@ ucsb.edu)

Abstract The ability to predict crystal growth habits is an important component of drug design, enabling a targeted sweep of optimal growth conditions that confer desirable properties. This article presents an investigation into the shape of olanzapine crystals grown from various solvents, exemplifying how mechanistic models of spiral growth can be applied to small molecule therapeutics. Olanzapine is recognized as the most effective treatment for schizophrenia, but a mechanistic treatment of the underlying crystal growth has yet to be established. We model spiral growth of olanzapine form I from five solvents (acetone, ethyl acetate, toluene, methyl isobutyl ketone and n-butyl acetate), considering a dimeric growth unit and periodic bond chains

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 28

consisting of inter-dimer bonds. The centrosymmetric dimers are stabilized by multiple C-H… contacts. The {1 0 0} face family dominates the predicted crystal habits, in agreement with our experiments; this morphology stems from the in-plane hydrogen bonds that are exposed on the {1 0 0} surface. The close agreement between predicted morphologies and experimental determinations lends support to the hypothesis that olanzapine grows as a dimer from these solvents. This mechanistic treatment can be readily applied to other compounds.

2 ACS Paragon Plus Environment

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

Introduction Schizophrenia is a severe chronic psychotic disorder that causes significant functional impairment and a severe reduction in the quality of life.1-3 Both pharmacologic and psychological treatments are employed to address these symptoms, and the former is typically viewed as more effective, with antipsychotic medications as the mainstay.4,5 After much exploration, olanzapine has emerged as the therapeutic of choice for treating schizophrenia, due to its demonstrated efficacy.6 Although olanzapine represents the state-of-the-art, its growth mechanism and crystallization behavior are not fully understood, in part due to the complexity of its molecular structure. Olanzapine

(2-methyl-4-(4-methyl-l-piperazinyl)-10H-thieno-[2,3b][1,5]benzo-

diazepine) can form polymorphs, hydrates and solvates, resulting in a complicated manufacturing process and diverse solid-state forms.7

The building block of all known

olanzapine crystals is a molecular dimer stabilized by multiple weak C–H…π contacts.8-10 There are three fused rings (diazepine, phenyl and thiophene) and one additional piperazine ring in the olanzapine molecule as shown in Fig. 1, with the middle seven membered diazepine ring adopting a distorted boat conformation.11 The dispersion-bound dimer is rather compact, its packing leaves exposed H-bonding groups (donor NH and acceptor N) and void spaces which can accommodate many different small organic solvent molecules, leading to solvate solid forms. Olanzapine crystallizes in three polymorphic anhydrates (I, II, III), three polymorphic hydrates, a higher hydrate and many pure and mixed solvates.12 Among the anhydrous forms, form I is monotropically the most stable and can be obtained by direct crystallization from suitable anhydrous solvents that are less amenable to solvate formation, such as ethyl acetate, acetone, toluene, n-butyl acetate and methyl isobutyl ketone. Metastable forms II and III have 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 28

been concomitantly produced by desolvation of solvates, as well as by sublimation; pure samples of form II and III have not been produced to date. The solubility of olanzapine in water is very low, which is why an organic solvent is used in the manufacturing process. Since form I is most stable, it is the preferred pharmaceutical solid form. The crystallization process influences many fundamental properties of the resultant drug, including chemical purity and composition, polymorphic state, crystal size and shape distribution, etc. Size and shape distributions in turn influence many solid properties, such as end-use efficacy, e.g., bioavailability for pharmaceuticals,13 reactivity for catalysts,14 flowability, wettability,15-17 and adhesion.18 As a consequence of the impact of crystal morphology on a material’s properties, control of crystallization is essential for effective downstream blending and compaction of drug particles. In order to scientifically engineer the shape of crystalline solids, significant efforts have been undertaken in recent decades to understand and effectively model the underlying crystal growth mechanisms.19-24,26-28 For the controlled production of pharmaceutical compounds from solution, the 2D nucleation and spiral growth mechanisms are most relevant. For organic molecular crystals, these models provide useful insight to design materials and guide experiments. Our group at UCSB has contributed to the development of a general mechanistic model to predict the crystal morphologies of organic molecules (for both centrosymmetric and non-centrosymmetric growth units) grown from different environments (including vapor or solutions), via the surface growth mechanisms of spirals and twodimensional nucleation.19-25 These models capture the effect of factors such as temperature, solvent and supersaturation by accounting for the fundamental surface growth physics and interfacial chemistry. These models have furthermore been successfully applied by Koo and coworkers for many crystalline explosives.26-28 The Small Molecule Design and Development 4 ACS Paragon Plus Environment

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

team at Eli Lilly and Co. has extensive experience crystallizing olanzapine from solution, where olanzapine form I is grown commercially by cooling crystallization from ethyl acetate.29 Our two teams have joined in this collaboration to crystallize olanzapine from several solvents, to predict its morphology from those same solvents, and to compare the results. Although the technological processes for olanzapine are relatively mature, from the perspective of fundamental understanding, fewer reports exist on experimental morphologies or the growth mechanism; correspondingly, to our knowledge there are no crystal habit predictions of olanzapine grown from solution. Crystal growth happens at the interface between the crystal surface and the solution or vapor. During the process of crystallization, the building blocks can be a single molecule (i.e., monomer) or several molecules stabilized by e.g. intermolecular bonds (a molecular cluster, typically existing as a dimer). The propensity to crystallize as a monomer or dimer can depend on both polymorph and solution environment. For example, as evidenced by solution diffusion data and surface diffraction, the crystallization of α-glycine from aqueous solution proceeds via incorporation of cyclic, hydrogen-bonded dimers,30-31 while the γ-form in acid or base solutions crystallizes as a monomer. For benzoic acid and tetrolic acid, FTIR spectroscopy also supports the conclusion that nucleation of a dimer is favored kinetically as a result of the preexistence of dimers in solution.32 The aforementioned dimers are hydrogenbonded with strong intermolecular (intra-dimer) interactions stabilizing the molecular pair and rendering it as the most favorable lattice building block. According to theoretical calculations, the binding energy of the olanzapine dimer is 8.68 kcal/mol,8 which indicates a similarly strong intermolecular interaction, although it is not formed from hydrogen bonds. Furthermore, PIXEL calculations indicate that this dispersive-interactionstabilized dimer could be the growth unit in olanzapine crystallization.7 Based on the strong 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 28

intermolecular binding interaction, we consider this dimeric growth unit in modeling both the vapor crystallization and solution crystallization of olanzapine form I from acetone, ethyl acetate, toluene, methyl isobutyl ketone and n-butyl acetate. Due to its complicated molecular structure and ability to form various solvates and hydrates, it presents an interesting, industrially relevant system to study and benchmark the predictive utility of these mechanistic models. Of particular interest is accounting for dimeric growth, which we expect given the strength of the intra-dimer interactions and propensity for dimers to exist in solution. Thus, the application of mechanistic models and periodic bond chains (PBCs) consisting of dimeric growth units is a principal focus of this article. The following sections describe the experimental morphologies and the morphological predictions using the spiral growth model for dimeric olanzapine form I.

Experimental Morphologies Form I single crystals were grown by slow evaporation from a short list of solvents that have yet to form solvates with olanzapine:7 ethyl acetate, n-butyl acetate, acetone (dry), methyl isobutyl ketone, and toluene. Olanzapine solutions were freshly prepared, filtered through 0.45 micron Millex PTFE filters, then stored in covered petri dishes in a dark refrigerator at 5 °C (to minimize chemical decomposition) to allow for slow solvent evaporation. Depending on the solvent, single crystals grew in hours to days, in some cases to mm dimensions, from the solutions. The morphology of the form I crystals was captured, in the mother liquors when possible, by polarized light microscopy using an Olympus SZH Zoom stereo microscope, equipped with an Olympus SZH-D.F. Plan 0.75X objective and interfaced with a PaxCAM2 digital camera and PC. In all cases, the crystals of form I were rhombic and the {1 0 0} faces dominant. 6 ACS Paragon Plus Environment

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

X-ray powder diffraction (XRPD) data were collected for samples (5-10 mg) placed on Si based low background sample holders, using a Bruker D8 Advance reflection diffractometer, equipped with a CuKa source (λ=1.54056 Å) and a Lynxeye detector, and operating at 40 kV and 40 mA, with a 0.2 mm divergence slit and an auto-changer. Data were collected in the range 240° with a 0.02° 2θ step size and 0.2 s step-1 count time. Experimental data were processed using MDI Jade 2010 (vs. 3.5.7). The XRPD patterns for crystals grown from these five solvents are reported in Fig. 2, together with the reference pattern for form I. The similarity of these patterns to the reference indicates that form I was grown from each solvent.

Predicted Morphologies Olanzapine form I crystallizes in the space group P21/c with four molecules in the unit cell. The cell parameters are a = 10.38 Å, b = 14.83 Å, c = 10.56 Å, β = 100.62° (UNOGIN01).9 We first calculated the intermolecular interactions within the lattice, using ADDICT (Advanced Design and Development of Industrial Crystallization Technology)22 and applying the AMBER force field33,34 as implemented by our group at UC Santa Barbara. The partial charges for each atom were obtained using GAUSSIAN0335 with the restrained electrostatic potential (RESP) model.36 The centrosymmetric olanzapine dimers with mass centers located at (0.5, 0.5, 0.5) and (0.5, 0, 1) were selected for the unit cell, as shown in Fig. 3.9 We then calculated interactions between olanzapine dimers in the solid state in order to determine periodic bond chains representing the strongest directional repeated interactions. The calculated lattice energy for dimer growth units is -26.73 kcal/mol (per mol of monomers), while it is -35.28 kcal/mol for a monomeric growth unit. The difference between these two lattice energies provides the binding energy for olanzapine dimers: -8.54 kcal/mol. Wawrzycka-Gorczyca et al. calculated the intra7 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 28

dimer binding energies as -8.68 kcal/mol, at MP2 level of theory;8 the close agreement indicates that the applied AMBER force field is accurate enough for the present system. In order to predict the crystal shape for olanzapine, normalized step velocities on each Fface in the solid state were calculated and used to form expressions for spiral rotation times. The rotation times were used to calculate relative growth rates between various face families and determine the steady-state crystal shape via the Frank-Chernov condition (see eq. S13-S14).37,38 F-faces are crystal faces with two or more periodic bond chains. To determine periodic bond chains, we first calculated the inter-dimer solid-state interactions. Calculating aggregate interdimer interactions requires summation of four intermolecular interactions (each molecule in dimer 1 interacts with each molecule in dimer 2). For solution growth, the interfacial energy is calculated to account for the presence of solvent molecules. This treatment considers both pure solute and solvent cohesive energy penalties, in addition to the favorable contribution from the work of adhesion at the interface.19,39 The internal crystalline and solvent bonds must be broken to form a surface, which explains the cohesive energy penalty, and new crystal-solvent bonds form and deliver the adhesive energy reward (see eqs. S16-S18). The pure surface energies are divided into dispersive and acid-base components in order to quantify the degree to which matching interactions can exist at the surface. This division uses the force field components for the solid-state interactions and solubility parameters for the solvent-side interactions (see eq. S15, and values of solubility parameters in Table 1).40 As shown in eq. S17-S18, the interfacial energy between solid and solvent is calculated as the cohesive energy of each phase minus the interfacial work of adhesion. After obtaining bond energies in solution, periodic bond chains (PBCs) constructed from the strong repeating inter-dimer bonds were generated, and F-faces with two or more PBCs were found. We only include bond chains for which | | ≥ 8 ACS Paragon Plus Environment

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

0.3 / (≈ ⁄2, ≈ 0.5  per bond), because edges with lower bond energy do not form identifiable steps on the crystal surface. There are seven F-face families with two or more PBCs:{1 0 0}, {1 1 1}, {0 1 1}, {1 1 0}, {1 2 0}, {0 2 0}, and {1 0 2}. Among these F-face families, {1 0 0}, {1 1 1} and {0 1 1} are the three slowest growing. The PBC bond energies for these faces are listed in Table 2 (the force field components for the PBC bonds are reported in Table S1), and the bonding structures for them are shown in Fig. 4. The corresponding results for the other faces can be seen in Fig. S1. As shown in Table 2, compared with the base solid-state bond energies (which correspond to vapor growth), the values for most PBCs are reduced for each solvent, indicating the stabilizing influence of solvent on crystal growth. The work of adhesion can be large if either significant dispersive interaction matching between the solid and solvent exists across the interface, or if there is potential for favorable crystal-solvent hydrogen bonding. From Table 2, it can also be seen that the [011] edge family ([011], [011], [01 1], [011]) typically represents the most stable edges (strongest kinkdirection interactions) under all solvents. The strength of the [011] edge family stems from the solid-state hydrogen-bonding interactions that exist between dimers (this distribution of hydrogen bonds is illustrated in Fig. 5 for the {100} face). Table 2 also indicates which edges remain on the steady-state spiral shapes predicted for each face, under each solvent. Eqs. S1-S2 were used to obtain kink densities and step velocities resulting from these solid-state interactions that form PBCs of centrosymmetric olanzapine dimers; these simple expressions are applicable due to the symmetric interaction environments present. Relative growth rates of F-faces are then calculated using eqs. S10-S14. Kink energies, kink densities, critical lengths, and face growth rates (relative to the {100} family) for the {1 0 0}, {1 1 1} and {0 1 1} faces are presented in Table 3, with intermediate data and corresponding results for the other faces reported in the SI. 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 28

In vapor and all solvents studied, the {1 0 0} face family represents the slowest growing faces with greatest expression on the predicted steady-state morphologies. The surrounding faces in vapor growth are predicted to be the {0 1 1} family, generating a regular platelet crystal habit (shown in Fig. 6). For solution growth, however, the surrounding faces are instead predicted to be the {1 1 1} and {0 1 1} face families, except in toluene; this still produces rhombic crystals with the {1 0 0} faces dominant (see Fig. 7). For toluene, the only existing surrounding face family is {1 1 1}, due to different solvent effects. The predicted spiral shapes for vapor growth are shown in the SI. For all studied growth conditions, the remaining F-face families ({1 1 0}, {1 2 0}, {0 2 0} and {1 0 2 }) are predicted to grow out of the steady-state shape. These morphologically unimportant faces have high relative growth rates because of the generally lower values of kink energies within the slice, which produce high kink densities and low critical lengths, both of which lead to relatively short spiral rotation times, and thus larger normal growth rates. According to a recent article,41 crystals of olanzapine form I grown from ethyl acetate solution were indexed. The experimental crystal shape was observed to be rhombic with dominant {1 0 0} faces, and surrounding faces belonging to the {1 1 1} and {0 1 1} families. Therefore, our predictions agree well with the experiments. From the spiral growth rate expression there are three natural groupings of mechanistic parameters that primarily depend on solvent, crystal face, and supersaturation, respectively (see eq. S22). All quantities can be readily calculated except for  (a solvent dependent parameter). To compare predicted growth rates between solvents we report  ⁄  for the dominant face family (see Table 4). The detachment work anisotropy (i.e., differences in ∆ between solvents, or equivalently, differences in solubility between solvents) is the primary cause of differences

10 ACS Paragon Plus Environment

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

between these predicted normalized growth rates. Table 4 indicates that on the basis of values of  ⁄  , the (100) face grows fastest in toluene and slowest in acetone. Absolute growth rates will depend on the individual values of  in each solvent (note,  is related to the exponential of the free energy barrier for attachment of growth units at kink sites). Thus, we also report relative  attachment rate constant values (  , arbitrarily assigning  in toluene as unity) that would

make the (100) face grow equally fast in each solvent. In acetone, for example,  would need to be ~ 106 times larger than  for toluene in order for both solvents to produce similar absolute growth rates for the (100) face. Until such time as  can be routinely calculated, Table 4 is the best estimate of relative growth rates for olanzapine crystals in different solvents. To test the accuracy of our predicted steady-state crystal habits, we grew olanzapine form I from five different solvents—acetone, ethyl acetate, toluene, methyl isobutyl ketone and nbutyl acetate by slow evaporation. These morphologies are shown in Fig. 7, alongside the theoretical predictions under spiral growth. For each solvent, the model predictions align well with the experimentally observed rhombic platelets.

Conclusion The successful prediction of olanzapine form I crystal habits grown from solution lends support to the hypothesis of a dimeric growth unit and operation of a spiral growth mechanism. This strongly bonded dimer is stabilized by multiple C—H…π contacts. The {100} face family dominates predicted shapes for all solvents. For solvents with larger polar and hydrogen bonding solubility parameters, the surrounding face families belong to {0 1 1} and {1 1 1}, while for nonpolar toluene, there is only surrounding face family {1 1 1}. Solvent modified bond energies

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 28

are highest in acetone and smallest in toluene, suggesting the rate of crystallization may be faster for toluene (though anisotropic desolvation kinetics could still reverse that trend). The {100} face family exposes solid-state hydrogen-bonds, offering a mechanistic explanation for its morphological dominance. Ultimately, mechanistic modeling offers a powerful tool for engineering pharmaceutical compounds, providing insight into potential shape modification through fundamental theoretical understanding.

Supporting Information Mechanistic growth model and parameters used to predict the crystal shapes in vapor growth and all five solvents, bonding structures for all the F-faces, predicted spirals in vapor and additional illustrations of the crystal shape grown from all five solvents.

Acknowledgements The authors are grateful for the financial support provided by Eli Lilly and Company, Pfizer, Novartis, the National Natural Science Foundation of China (No. 51672176) and the Intergovernmental International Scientific and Technological Cooperation of Shanghai (No. 17520710200). We thank an anonymous reviewer for very helpful suggestions.

References 1. Nordstroem, A. L.; Talbot, D., Bernasconi, C.; Berardo, C. G.; Lalonde, J. Burden of Illness of People with Persistent Symptoms of Schizophrenia: A Multinational Crosssectional Study. Int J Soc Psychiatry. 2017, 63, 139-150. 2. Csernansky, J. G.; Schuchart, E. K. Relapse and Rehospitalisation Rates in Patients with 12 ACS Paragon Plus Environment

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

Schizophrenia: Effects of Second Generation Antipsychotics. CNS Drugs, 2002, 16, 473484. 3. Insel T. R. Rethinking schizophrenia. Nature, 2010, 468,187-93. 4. Biedermann, F.; Fleischhacker, W. W. Emerging Drugs for Schizophrenia. Expert Opin. Emerg. Drugs 2011, 16, 271-282. 5. Jarskog, L. F.; Miyamoto, S.; Lieberman, J. A. Schizophrenia: New Pathological Insights and Therapies. Annu. Rev. Med. 2007, 58, 49-61. 6. Attard, A.; Taylor, D. M. Comparative Effectiveness of Atypical Antipsychotics in Schizophrenia CNS Drugs 2012, 26, 491-508. 7. Bhardwaj, R. M.; Price, L. S.; Price, S. L.; Reutzel-Edens, S. M.; Miller, G. J.; Oswald, I. D. H.; Johnston, B.; Florence, A. J. Exploring the Experimental and Computed Crystal Energy Landscape of Olanzapine. Cryst. Growth Des. 2013, 13, 1602-1617. 8. Wawrzycka-Gorczyca, I.; Borowski, P.; Osypiuk-Tomasik, J.; Mazur, L.; Koziol, A. E. Crystal Structure of Olanzapine and its Solvates. Part 3. Two and Three-component Solvates with Water, Ethanol, Butan-2-ol and Dichloromethane J. Mol. Struct. 2007, 830, 188-197. 9. Wawrzycka-Gorczyca, I.; Koziol A. E.; Glice, M.; Cybulski, J. Polymorphic Form II of 2methyl-4-(4-methyl-1-piperazinyl)-10H-thieno-[2,3-b][1,5]-benzodiazepine. Acta Cryst. 2004, E60, o66-o68. 10. Wawrzycka-Gorczyca, I.; Mazur, L.; Koziol, A. E. 2-Methyl-4-(4-methyl-1-piperazinyl)10H-thieno-[2,3-b][1,5]-benzodiazepine Methanol Solvate Acta Cryst. 2004, E60, o69−o71. 11. Hendrickson, J. B. Molecular Geometry. VII. Modes of Interconversion in the Medium 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 28

Rings J. Am. Chem. Soc. 1967, 89, 7047-7061. 12. Reutzel-Edens, S. M.; Bush, J. K.; Magee, P. A.;, Stephenson, G. A.; Byrn, S. R. Anhydrates and Hydrates of Olanzapine: Crystallization, Solid-State Characterization, and Structural Relationships. Cryst. Growth Des. 2003, 3, 897-907. 13. Yin, J.-C.; Zhou, J.-S.; Sun, J.; Qiu, Y.; Wei, D.-Z.; Shen, Y.-L. Study of the crystal shape and its influence on the anti-tumor activity of tumor necrosis factor-related apoptosisinducing ligand (Apo2L/TRAIL). Cryst. Res. Technol. 2008, 43, 888-893. 14. Yang, H. G.; Sun, C. H.; Qiao, S. Z.; Zou, J.; Liu, G.; Smith, S. C.; Cheng, H. M.; Lu, G. Q. Anatase TiO2 single crystals with a large percentage of reactive facets. Nature 2008, 453, 638-641. 15. Prestidge, C. A.; Tsatouhas, G. Wettability studies of morphine sulfate powders. Int. J. Pharm. 2000, 198, 201-212. 16. Heng, J.; Williams, D. Wettability of Paracetamol Polymorphic Forms I and II. Langmuir 2006, 22, 6905-6909. 17. Heng, J.; Bismarck, A.; Lee, A.; Wilson, K.; Williams, D. Anisotropic Surface Energetics and Wettability of Macroscopic Form I Paracetamol Crystals. Langmuir 2006, 22, 27602769. 18. Riebel, U.; Kofler, V;, Lffler, F. Shape Characterization of Crystals and Crystal Agglomerates. Part. Part. Syst. Charact. 1991, 8, 48-54. 19. Snyder, R. C.; Doherty, M. F. Predicting Crystal Growth by Spiral Motion. Proc. R. Soc. London, Ser. A. 2009, 465, 1145-1171. 20. Kuvadia, Z. B.; Doherty, M. F. Spiral Growth Model for Faceted Crystals of Noncentrosymmetric Organic Molecules Grown from Solution. Cryst. Growth Des. 2011, 11, 14 ACS Paragon Plus Environment

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

2780-2802. 21. Lovette, M. A.; Doherty, M. F. Predictive Modeling of Supersaturation dependent Crystal Shapes. Cryst. Growth Des. 2012, 12, 656-669. 22. Li, J. J.; Tilbury, C. J.; Kim, S.H.; Doherty, M. F. A Design Aid for Crystal Growth Engineering. Prog. Mater Sci. 2016, 82, 1-38. 23. Tilbury, C. J.; Doherty, M. F. Modeling Layered Crystal Growth at Increasing Supersaturation by Connecting Growth Regimes. AIChE J. 2017, 63, 1338-1352. 24. Tilbury, C. J.; Joswiak, M. N.; Peters, B.; Doherty, M. F. Modeling Step Velocities and Edge Surface Structures during Growth of Non-Centrosymmetric Crystals. Cryst. Growth Des. 2017, 17, 2066-2080. 25. Li, J. J.; Abramov, Y. A.; Doherty, M. F. New Tricks of the Trade for Crystal Structure Refinement. ACS Cent. Sci. DOI: 10.1021/acscentsci.7b00130. 26. Shim, H.-M.; Koo, K.-K. Crystal Morphology Prediction of Hexahydro-1,3,5-trinitro1,3,5-triazine by the Spiral Growth Model. Cryst. Growth Des. 2014, 14, 1802−1810. 27. Shim, H.-M.; Koo, K.-K. Prediction of Growth Habit of b-cyclotetramethylenetetranitramine Crystals by the First-principles Models. Cryst. Growth Des. 2015, 15, 3983-3991. 28. Shim, H.-M.; Kim, H.-S.; Koo, K.-K. Molecular Modeling on Supersaturationdependent Growth Habit of 1,1-diamino-2,2-dinitroethylene. Cryst. Growth Des. 2015, 15, 18331842. 29. Bunnell, C. A.; Hendriksen, B. A.; Larsen, S. D. Olanzapine Polymorph Crystal Form. US Patent 5736541, 1998. 30. Myerson, A. S.; Lo, P. Y. Cluster Formation and Diffusion in Supersaturated Binary and 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 28

Ternary Amino Acid Solutions. J. Cryst. Growth 1991, 110, 26-33. 31. Gidalevitz, D.; Feidenhans’l, R.; Matlis, S.; Smilgies, D.-M.; Christensen, M. J.; Leiserowitz, L. Monitoring In Situ Growth and Dissolution of Molecular Crystals: Towards Determination of the Growth Units. Angew. Chem., Int. Ed. Engl. 1997, 36, 955959. 32. Davey, R. J.; Dent, G.; Mughal, R. K.; Parveen, S. Concerning the Relationship between Structural and Growth Synthons in Crystal Nucleation:  Solution and Crystal Chemistry of Carboxylic Acids As Revealed through IR Spectroscopy. Cryst. Growth. Des. 2006, 6, 1788-1796. 33. Case, D. A. AMBER 10; University of California: San Francisco, 2008. 34. Cornell, W. D., Cieplak, P., Bayly, C. I., Gould, I. R., Merz, K. M., Ferguson, D. M., Spellmeyer, D. C., Fox, T., Caldwell, J. W., Kollman, P. A. J. Am. Chem. Soc. 1995, 117, 5179-5197. 35. Frisch, M. J. Gaussian 03, Revision C.02. 2003. 36. Bayly, C. I.; Cieplak, P.; Cornell, W.; Kollman, P. A. A Well-Behaved Electrostatic Potential Based Method Using Charge Restraints for Deriving Atomic Charges: The RESP Model. J. Phys. Chem. 1993, 97, 10269-10280. 37. Chernov, A. A. Crystal Growth Forms and Their Kinetic Stability. Sov. Phys. Cryst. 1963, 7, 728-730. 38. Frank, F. C. In: Growth and Perfection of Crystals; Doremus, R. H.; Roberts, B. W.; Turnbull, D. Eds.; Wiley: New York, 1958. 39. Tilbury, C. J., Green, D. A., Marshall, W. J., Doherty, M. F. Predicting the Effect of Solvent on the Crystal Habit of Small Organic Molecules. Cryst. Growth Des. 2016, 16, 16 ACS Paragon Plus Environment

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

2590-2604. 40. Barton, A. F. M. Solubility Parameters. Chem. Rev. 1975, 75, 731-753. 41. Warzecha, M; Guo, R.; Bhardwaj, R. M.; Reutzel-Edens, S. M.; Price, S. L.; Lamprou, D. A.; Florence, A. J. Direct observation of templated two-step nucleation mechanism during olanzapine hydrate formation. Cryst. Growth Des. DOI: 10.1021/acs.cgd.7b01060

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 28

Figure 1. Molecular structure of olanzapine.

18 ACS Paragon Plus Environment

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

Crystal Growth & Design

Figure 2. XRPD data for olanzapine form I grown from toluene, n-butyl acetate, methyl isobutyl ketone, ethyl acetate, and acetone. (Form I-reference means the calculated XRPD pattern using the form I crystal structure (UNOGIN03) taken from the CSD)

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

Page 20 of 28

Figure 3. The olanzapine unit cell containing two centrosymmetric dimers.

Figure 4. Bonding structures of different F-faces for olanzapine dimers (a) {1 0 0} face, (b) {1 11} face, (c) {0 1 1} face (with graphic axes shown in the picture). 20 ACS Paragon Plus Environment

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

Crystal Growth & Design

Figure 5. The distribution of hydrogen bonds between olanzapine dimers within the {1 0 0} face.

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 28

Figure 6. Crystal habit prediction for sublimation growth of olanzapine form I. ({1 0 0} face family in yellow; and {0 1 1} face family in pink)

22 ACS Paragon Plus Environment

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

Crystal Growth & Design

Figure 7. Comparison between predicted and experimentally observed morphologies of olanzapine form I grown from (a) acetone, (b) ethyl acetate, (c) toluene, (d) methyl isobutyl ketone and (e) n-butyl acetate. ({100} face family in yellow, {1 11} face family in blue, and {0 1 1} face family in pink)

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 28

Table 1. Solvent information: molar volume (V" , cc/mol), dispersion solubility parameter (δ) , (cal/ mol)-/. ), polar solubility parameter ( δ/0 , (cal/mol)-/. ), hydrogen bonding solubility parameter (

δ1 , (cal/mol)-/.

),

scale

factor

(f),

and

surface

energies

(

2303 , erg⁄cm.). Solvent

V"

δ)

δ789

δ1

f

2303

Acetone

74.0

7.6

5.1

3.4

0.10

20.39

Methyl isobutyl ketone

125.8

7.5

3.0

2.0

0.10

17.52

Ethyl acetate

98.5

7.7

2.6

3.5

0.12

22.74

n-butyl acetate

132.5

7.7

1.8

3.1

0.12

21.90

Toluene

106.8

8.8

0.7

1.0

0.153

28.49

24 ACS Paragon Plus Environment

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

Crystal Growth & Design

Table 2. Energetic interactions for olanzapine, organized by face and edge. Solid-state (vapor growth) interaction strengths are listed, in addition to solvent-modified interactions. Bond energies refer to interdimer interactions and each PBC edge is responsible for two opposing edges of the resulting spiral surface structures (i.e., positive and negative directions). Bold typeface indicates the edges that remain in the steady-state spiral surface structures. Total PBC bond energy (kcal/mol) Face

{1 0 0}

{1 1 1}

{0 1 1}

Edge Vapor

Acetone

Ethyl acetate

Toluene

n-butyl acetate

Methyl isobutyl ketone

[001]

-5.76

-3.42

-2.34

-1.41

-1.85

-1.41

[011]

-14.60

-5.50

-4.61

-4.70

-4.37

-5.00

[01 1]

-14.60

-5.50

-4.61

-4.70

-4.37

-5.00

[101]

-3.02

-2.81

-2.23

-2.08

-1.94

-1.35

[211]

-2.80

-1.42

-0.92

-0.39

-0.70

-0.58

[011]

-14.60

-5.50

-4.61

-4.70

-4.37

-5.00

[100]

-7.70

-3.04

-1.92

-0.20

-1.47

-1.75

[211]

-2.80

-1.42

-0.92

-0.39

-0.70

-0.58

[011]

-14.60

-5.50

-4.61

-4.70

-4.37

-5.00

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 28

Table 3. Predicted relative growth rates of olanzapine form I grown from ethyl acetate, acetone, toluene, n-butyl

acetate

and

methyl

isobutyl

ketone.

(h

is

the

step

height

in

Å;

Φ< is the kink energy in kcal/mol; ρ is the kink density; lc is the critical length in Å; R is growth rate of

each face relative to the {100} family under the same solvent) Faces

{1 0 0}

{1 1 1}

{0 1 1}

h

10.21

7.08

8.50

Edges

Vapor

[001]

[011]

[01 1]

[101]

[211]

[011]

[100]

[211]

[011]

Φ