Determining the Crystal Structures of Peptide Analogs of Boronic Acid

6 days ago - wavelength (Å). 0.399978(2). 0.399963(3) θ range for data processing (deg). 0.502−22.498. 0.502−22.498. Rp; Rwp. 0.0442; 0.0643. 0...
0 downloads 3 Views 1MB Size
Subscriber access provided by UNIVERSITY OF THE SUNSHINE COAST

Article

Determining the crystal structures of peptide analogs of boronic acid in the absence of single crystals: Intricate motifs of ixazomib citrate revealed by XRPD guided by ss-NMR Michal Husak, Alexandr Jegorov, Jan Rohlicek, Andrew N. Fitch, Jiri Czernek, Libor Kobera, and Ji#í Brus Cryst. Growth Des., Just Accepted Manuscript • DOI: 10.1021/acs.cgd.8b00402 • Publication Date (Web): 30 Apr 2018 Downloaded from http://pubs.acs.org on May 2, 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 36 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

Cover Page

Determining the crystal structures of peptide analogs of boronic acid in the absence of single crystals: Intricate motifs of ixazomib citrate revealed by XRPD guided by ss-NMR Michal Hušák,a Alexandr Jegorov,b Jan Rohlíček,c Andrew Fitch,d Jiří Czernek,e Libor Kobera,e and Jiří Brus*e a)

Department of Solid State Chemistry, University of Chemistry and Technology, Technicka 5, 166 28 Prague 6, Czech Republic. b) Teva Pharmaceuticals CR, s.r.o, Branisovska 31, 370 05 Ceske Budejovice, Czech Republic. c) Department of Structural Analysis, Institute of Physics of the Czech Academy of Sciences, Na Slovance 2, Praha 8, 182 21, Czech Republic. d) ESRF- The European Synchrotron, CS40220, F-38043 Grenoble Cedex 9, France. e) Department of NMR Spectroscopy, Institute of Macromolecular Chemistry CAS, Heyrovsky sq. 2, 162 06 Prague 6, Czech Republic. ABSTRACT: Uncertainties in the structure determination of peptide analogs of boronic acid, exacerbated by the many coordination modes of boron, represent an obstacle in understanding their role in living organisms and thus also in developing the next generation of anticancer drugs. For that reason, we present here a general experimentalcomputational strategy allowing structure determination of complex boronic acid derivatives with extensive conformational variability. We demonstrate successful solution of the crystal structures of two non-solvated polymorphs of ixazomib citrate directly from synchrotron powder diffraction data, which is challenging because the two molecules in the asymmetric unit cell exhibit 32 degrees of conformation freedom push the limits of current solution procedures. We used a novel two-step Rietveld refinement based on DFT-D restraints to improve information quality derived from powder diffraction data to be comparable with that of single-crystal solutions. NMR crystallography was applied to verify the crystal structures, and the high potential value of using 11B NMR parameters toward the solution of unknown structures was demonstrated. Evolution of 11B-11B double-quantum coherences allows probing of interatomic distances up to 7 Å. Overall, we present an integrated approach that applies several techniques in conjunction to provide otherwise unavailable structural information.

For Table of Contents Use Only

Synopsis: General experimental-computational strategy for exploring the assembly of peptidic boronic acid derivatives in microcrystalline state is demonstrated. By combining the simulated annealing and two-step Rietveld refinement of XRPD data with the comparative analysis of NMR parameters the crystal structures of two non-solvated polymorphs of ixazomib citrate, each consisting of two symmetry-independent molecules in unit cell with 32 free parameters, were solved and refined.

Corresponding Author J.B. *Email: [email protected] http://www.imc.cas.cz/en/umch/o_nmr.htm

ACS Paragon Plus Environment

1

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 36

Determining the crystal structures of peptide analogs of boronic acid in the absence of single crystals: Intricate motifs of ixazomib citrate revealed by XRPD guided by ss-NMR Michal Hušák,a Alexandr Jegorov,b Jan Rohlíček,c Andrew Fitch,d Jiří Czernek,e Libor Kobera,e and Jiří Brus*e a)

Department of Solid State Chemistry, University of Chemistry and Technology, Technicka 5, 166 28 Prague 6, Czech Republic.

b)

c)

Teva Pharmaceuticals CR, s.r.o, Branisovska 31, 370 05 Ceske Budejovice, Czech Republic.

Department of Structural Analysis, Institute of Physics of the Czech Academy of Sciences, Na Slovance 2, Praha 8, 182 21, Czech Republic. d)

e)

ESRF- The European Synchrotron, CS40220, F-38043 Grenoble Cedex 9, France.

Department of NMR Spectroscopy, Institute of Macromolecular Chemistry CAS, Heyrovsky sq. 2, 162 06 Prague 6, Czech Republic.

KEYWORDS: XRPD • NMR crystallography • peptide and boronic acid analogs • ixazomib citrate.

ABSTRACT: Uncertainties in the structure determination of peptide analogs of boronic acid, exacerbated by the many coordination modes of boron, represent an obstacle in understanding

ACS Paragon Plus Environment

2

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

their role in living organisms and thus also in developing the next generation of anticancer drugs. For that reason, we present here a general experimental-computational strategy allowing structure determination of complex boronic acid derivatives with extensive conformational variability. We demonstrate successful solution of the crystal structures of two non-solvated polymorphs of ixazomib citrate directly from synchrotron powder diffraction data, which is challenging because the two molecules in the asymmetric unit cell exhibit 32 degrees of conformation freedom push the limits of current solution procedures. We used a novel two-step Rietveld refinement based on DFT-D restraints to improve information quality derived from powder diffraction data to be comparable with that of single-crystal solutions. NMR crystallography was applied to verify the crystal structures, and the high potential value of using 11

B NMR parameters toward the solution of unknown structures was demonstrated. Evolution of

11

B-11B double-quantum coherences allows probing of interatomic distances up to 7 Å. Overall,

we present an integrated approach that applies several techniques in conjunction to provide otherwise unavailable structural information.

Introduction Boron-containing compounds have long been recognized as potentially active pharmaceutical ingredients.1-3 As recent investigations have resulted in the discovery of many promising pharmaceuticals exhibiting anticancer and antibacterial activity such as bortezomib, MLN978, CEP-18770, and GSK2251052, research on peptidic derivatives with boronic acid fragments is rapidly gaining in intensity.4

ACS Paragon Plus Environment

3

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 36

Ixazomib citrate, an active pharmaceutical ingredient (API) consisting of a peptidic boronic acid derivative and citric acid, is a prodrug for the treatment of multiple myeloma (a type of white blood cell cancer). Ixazomib, the active part, is formed by hydrolysis of boronic ester under physiological conditions. However, due to the many coordination possibilities for boron atoms, the structures of ixazomib and its complexes are still under extensive investigation. Several possible coordination schemes for citric acid have been proposed, including those featuring five- or six-membered rings (Scheme 1).5 However, the coordination possibilities of boron are not entirely exhausted by these options; which coordination takes part in each crystalline form of ixazomib citrate is unclear, as is whether the coordination modes are equivalent when several symmetry-independent molecules are present in the crystal unit.

Scheme 1. Expected Structures of Ixazomib Citrate Complexes (C20H23BCl2N2O8).

In this regard, the recently derived single-crystal structures of two solvates indicated that the more preferred arrangement involves a 5-membered ring.6 In practice, however, most crystalline forms do not form single crystals of sufficient quality for crystal structure determination, and only attempting to solve the structure from powder diffraction data remains an option. However, prediction of the complicated solid-state structures of the peptide analogs of boronic acid, which result from a complex pathway involving reversible covalent and non-covalent interactions, is particularly challenging. The future design of new drugs based on boronic acid compounds thus requires new tools for their precise structural characterization.

ACS Paragon Plus Environment

4

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

Hence, the purpose of our study is to demonstrate an extended experimental-computational strategy for exploring the structure of complex peptide and boronic acid analogs with many degrees of conformational freedom. By combining advanced analysis of powder diffraction data with comparative analysis of NMR parameters and high-level quantum chemical calculations, complete crystal structures of two non-solvated polymorphs of ixazomib citrate (Form II and Form X), each consisting of two symmetry-independent molecules in a unit cell with 32 free parameters, were successfully determined.

Results and Discussion Structure solution and refinement from powder data. In the absence of suitable single crystals allowing for complete X-ray diffraction (XRD) analysis, we focused our attention on the advanced analysis of powder diffraction data (Figure 1). The data were preprocessed by synchrotron beamline-specific software7 as described in experimental section. The reflection positions were determined via DASH software,8 and indexing of all peaks was done via DICVOL06.9 DASH automatic space group determination indicated a P21/c space group for both forms of ixazomib citrate (Form II and Form X), which is inconsistent with the chiral character of the molecule. Detailed investigation of the diffraction record for Form II showed several weak reflections violating the P21/c group extinction, proving that the space group is P21 with two molecules in the asymmetric unit (very probably exhibiting pseudo-centrosymmetric behavior). The diffraction record of Form X did not violate the P21/c group extinction in any way—the correct P21 space group was confirmed during refinement. Both Forms II and X have two molecules of ixazomib citrate in their asymmetric unit.

ACS Paragon Plus Environment

5

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 36

X-ray synchrotron powder records a)

Form II

b) Form X

Figure 1. Comparison of the Measured and Calculated Powder Diffractograms for Ixazomib Form II and Form X (a and b, respectively). For Full-size Diffractograms, See Supporting Information, Figure S1.

Correct selection of an initial model is the key step in the structure determination from XRPD data. In the case of ixazomib citrate it was necessary to decide between the two possible structures consisting either of a five-membered ring or the six-membered arrangement (see Scheme 1). Therefore, to select the most suitable initial model for the subsequent structure solution the model of ixazomib citrate complex was taken from the known crystal structure of the ixazomib citrate ethanol solvate.6 However, the in vacuo RI-MP2 optimized models (IXAModel-4 and Model-5, see discussion below) showing the lowest relative potential energies and excellent agreement with experimental 11B NMR isotropic chemical shifts can be also utilized,

ACS Paragon Plus Environment

6

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

because in all cases the structures featuring five-membered rings were preferred. Consequently, in the absence of single crystal data, the ab initio optimized models featuring lowest potential energy represent quite suitable initial approximations. Since the second axial chiral center on boron is created during the crystallization, predicting whether the structure contains R,R, R,S enantiomers (as in the ethanol solvate), R,R, R,R enantiomers or the R,S, R,S combination (where the first chiral center corresponds to ixazomib and the second is generated on citric acid as a consequence of its coordination to boronic acid) was not possible. Testing all three possible combinations was required. Each molecule of ixazomib citrate has 16 DOF (degrees of freedom). Two molecules of ixazomib citrate in the asymmetric unit give 32 DOF, which is a complex problem whose solution is on the edge of current powder-solution computational methodology. The three possible combinations of conformations were each used 1000 times as a starting point for DASH simulated annealing runs, and 2.0e+7 steps were used per run. The computation was executed on 16 CPU cores in parallel, and the required equivalent single CPU execution time was approximately 1 year for each form. Finding the correct torsion angles was accelerated by statistical bias based on similar molecules’ torsion angles as found in the Cambridge Structural Database (CSD). Three identical solutions were found for Form II, all for the R,R, R,S conformation combination with the lowest cost functions. Fourteen identical solutions were found for Form X, all for the R,R, R,S conformations.

ACS Paragon Plus Environment

7

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 36

Table 1. Crystal Data for Form II and Form X, Both Restrained by the DFT-Calculated Molecular Geometry. Form II

Form X

Formula

C20H23BCl2N2O9

M.w.

517.119

Crystal system

Monoclinic

Monoclinic

Space group, Z

P21 (No. 4), 4

P21 (No. 4), 4

a, Å

16.0169(2)

16.1521(2)

b, Å

18.26999(15)

17.03683(12

c, Å

8.52827(6)

8.51817(6)

β, (°)

104.9730(7)

94.1479(8)

V, Å3

2410.88(4)

2337.90(4)

T(K)

293

293

Absorption coefficient (mm-1)

0.08

0.08

Wavelength (Å)

0.399978(2)

0.399963(3)

θ range for data processing (°)

0.502 - 22.498

0.502 - 22.498

Rp; Rwp

0.0442; 0.0643

0.0361, 0.0535

Goof

6.06

4.23

Parameters refined

254

246

Number of restraints and constraints

180, 266

180, 266

March-Dollase correction and direction

1.107(1 0 0)

1.045(1 0 0)

Software used for structure solution, DASH, DICVOL06 DASH, DICVOL06, indexation, interpretation.

refinement

and Jana2006, 2005

MCE Jana2006,

MCE

2005

The final refinement of the structure was done in Jana2006 software.10 Both structural models were checked and corrected by placing atoms via the difference Fourier maps in MCE software.11 Positions of all non-hydrogen atoms were refined together using one isotropic

ACS Paragon Plus Environment

8

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

displacement parameter shared by carbon, nitrogen, oxygen and boron atomic types; one harmonic displacement parameter shared by all chlorine atoms; unit cell parameters; background (36 points in the case of Form II and 28 points in the case of Form X); zero shift and the MarchDollase parameter. Hydrogen atoms were kept in their positions calculated from the geometry. In the first step of refinement, the bond-length and bond-angle restraints based on the ixazomib citrate isopropyl alcohol solvate6 were used. The resulting model was checked by DFT calculations (described below). In the final refinement, bond-length and bond-angle restraints based on the DFT-calculated geometry were used. All restraints were kept at the same initial weight (w = 1/σ2) with expected deviations of σ = 0.001 for bond-lengths and σ = 0.01 for bondangles from the average. For details of refinement, see Table 1. The accuracy of the X-ray powder diffraction (XRPD) structures of both Forms II and X was first tested by periodic DFT-D geometry optimization using CASTEP12 with the GGA PBE functional, TS dispersion correction, "Fine" settings, and an energy cut-off of 520 eV. "Fine" setting corresponds to mapping k-points with 0.07 Å-1 in reciprocal space - the Monkhorst-Pack parameters are derived to give the specified separation between neighboring grid points. The optimization was done in a two-step manner as suggested by van de Streek and Neumann.13 The geometry of the structure was first optimized, followed by the cell parameters. The maximum number of optimization steps was limited to 100. The convergence to predefined criteria is summarized in Supporting Information Table S1. After the cell optimization, the "Fine" convergence criteria were satisfied for Form X. However, for Form II, we were able to fully satisfy only the "Coarse" criteria for maximal displacement. Better results for Form II were unobtainable, so we used the data satisfying the "Coarse" criteria for maximum displacement, "Medium" criteria for maximal force and "Fine"

ACS Paragon Plus Environment

9

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 36

criteria for the rest. (For further details defining “Coarse”, “Medium” and “Fine” criteria see Supporting Information S2). The issues with DFT calculation convergence can be related to a very flat energy minimum. The flat minimum hypothesis corresponds well with the occurrence of two similar polymorphs. To compare the original and DFT-D-optimized geometries, rootmean-square (r.m.s.) Cartesian displacement (RMSCD) as previously defined13 was calculated in Crystal CMP software.14 The results are summarized in Supporting Information Table S2.

Ixazomib citrate Form II

Ixazomib citrate Form X

Figure 2. DFT-Optimized Refined XRPD Crystal Structures of Ixazomib Citrate Forms II and X.

ACS Paragon Plus Environment

10

Page 11 of 36 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 both Forms II and X, the RMSCD values calculated without H atoms (cell fixed as well as cell optimized) are less than 0.25 Å, which is typical for correctly solved structures. The RMSCD values calculated with H atoms are ca. 0.30 Å, which indicates slightly incorrect hydrogen placement based on only single-crystal restraints. The introduction of DFT-D-based restraints decreased both the r.m.s. deviation (r.m.s.d.) and the R factor. The refined crystal structures of both Forms II and X (Figure 2) are deposited in the Cambridge Crystallographic Data Centre (CCDC) with provisional numbers CCDC 1822510 and 1822511.

Independent validation of the crystal structures. As discussed above, determining the crystal structures of complex systems consisting of multiple molecules with unknown conformations directly from powder XRD data represents a great challenge. Independent verification of the refined XRPD structures is thus of particular interest. In this regard, NMR crystallography displays a remarkable potential because of the synergy of experimental and computational approaches.15-25 In general, solid-state NMR data can be used in two ways to complement crystal structure refinement. One approach is based on the systematic comparison of the NMR parameters calculated using periodic DFT methods on the refined XRPD structure with the corresponding experimentally determined values. Agreement between the calculated and experimental values then validates the crystal structure. The second, more straightforward route utilizes measurements of dipole-dipole interactions because they provide direct information on specific internuclear distances, molecular conformations and bonding arrangements.15-25 Specifically, we first focused on the analysis of

11

B NMR isotropic chemical shifts and

11

B

quadrupolar coupling constants to validate the refined XRPD structures of ixazomib citrate

ACS Paragon Plus Environment

11

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 36

Forms II and X. In this manner, we investigated the coordination geometry of the boron centers. Subsequently, we measured

11

B-11B and

1

H-1H dipolar interactions to probe long-range

interatomic distances and intermolecular contacts. Finally, we compared the DFT-calculated and experimentally determined 1H and

13

C NMR isotropic chemical shifts to determine the overall

consistency of the global architecture of both polymorphs of ixazomib citrate derived from XRPD data.

Determination of coordination geometry of boron centers. NMR parameters for the

11

B

nucleus such as the isotropic chemical shift (δiso) and the quadrupolar coupling constant (QC) represent receptive probes of the coordination geometries around boron atoms.26-31 Although these parameters cover relatively narrow ranges (0 - 40 ppm and 0 - 4 MHz, respectively), the resonance frequencies and line shapes made discriminating a variety of boron species possible, as illustrated in Figure 3.

ACS Paragon Plus Environment

12

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

11

B MAS NMR Spectra

symmetric BIV with 6-/6- rings deformed BIV OH 5-/6- rings B deformed BIV B B 5-/5- rings -

-

-

symmetric BIII

B-

HO B OH

HO

deformed BIII R B OH

HO

35

Figure 3.

11

30

25

20

15

10

5

0

-5

ppm

B MAS NMR Lineshapes of Model Compounds Simulated using Experimental

Parameters Found in the Literature:30,32 (2-methylpropyl)boronic acid (δiso = 32.7 ppm, QC = 3.1 MHz and η = 0.44; red line); boric acid (δiso = 19.2 ppm, QC = 2.1 MHz and η = 0.17, blue line); bortezomib Form I (δiso = 10.5 ppm, QC = 1.65 MHz and η = 0.82, purple line); the tetrahedral site in borax (δiso = 2.8 ppm, QC = 0.97 MHz and η = 0.90, green line) and the 6-membered cyclic borate diesters (δiso = 0.4 ppm, QC = 0.99 MHz and η = 0.70, violet line). For direct comparison, the

11

B MAS NMR spectrum of ixazomib citrate Form II simulated from the

experimentally determined parameters (δiso = 12.4 ppm, QC = 1.66 MHz and η = 0.38, black line) is also included.

For ixazomib citrate, the 2D

11

B triple-quantum (TQ)/MAS NMR spectra (Figure 4a)

demonstrate that the structural sensitivity of

11

B resonances is high enough to distinguish

different polymorphs with similar crystal structures. Moreover, in accord with quantum chemical calculations discussed in detail below, the obtained quadrupolar coupling constants (QC = 1.671.75 MHz) and isotropic chemical shifts (δiso = 12.4-13.5 ppm) (Tables 2 and 3) indicate that

ACS Paragon Plus Environment

13

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 36

boron atoms in both polymorphs adopt 4-fold coordination with significantly perturbed tetrahedral symmetry.

Table 2. Experimental

11

B NMR Isotropic Chemical Shifts and

11

B Quadrupolar Coupling

Constants (Qcc) of Boron Species in Ixazomib Citrate Forms II and X. Experimental 11B NMR parameters

δ(11B), ppm

Qc, MHz

η

Site 1

Site 2

Site 1

Site 2

Site 1

Site 2

F-II

12.45

12.45

1.67

1.67

0.39

0.39

F-X

13.11

13.54

1.75

1.74

0.45

0.48

Table 3. Calculated

11

B NMR Isotropic Chemical Shifts and

11

B Quadrupolar Coupling

Constants (Qc) of Boron Species in Ixazomib Citrate Forms II and X. The calculated chemical shifts were estimated from the GIPAW-PBE computed chemical shielding using the conversion factor.22 Calculated 11B NMR parameters

δ(11B), ppm

Qc, MHz

η

Site 1

Site 2

Site 1

Site 2

Site 1

Site 2

F-II

14.07

14.02

1.96

2.00

0.37

0.38

F-X

15.23

15.75

2.05

2.06

0.45

0.48

ACS Paragon Plus Environment

14

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

11

2D B TQ/MAS NMR Ixazomib citrate Form X

a) Ixazomib citrate Form II ppm

ppm

14

14

15

15

16

16 14

12

10

8

6

14

ppm

12

10

8

6

ppm

13

C CP/MAS NMR Ixazomib citrate Form X

b) Ixazomib citrate Form II C-OH (citric acid)

150

CH, CH3 (IXA)

100

50

C-OH (citric acid)

CH, CH3 (IXA)

150

ppm

100

50

ppm

1

H CRAMPS NMR

c) Ixazomib citrate Form II

Ixazomib citrate Form X COOH(2)

COOH(2) COOH(1)

COOH(1)

14

Figure 4. 2D

11

B TQ/MAS,

12

13

10

8

6

4

2

ppm

14

12

10.0 ppm

10

8

6

4

2

ppm

C CP/MAS and 1H CRAMPS NMR Spectra of Ixazomib Citrate

Forms II and X (Red and Blue Lines, Respectively). Spectral regions containing the C-OH and CH3 resonances of citric acid and ixazomib molecules are expanded in boxes above the fullrange spectra. When considering sources of deformation of tetrahedral symmetry, we first analyzed

11

B

isotropic chemical shifts. Literature data show that tetra-coordinated boron atoms shared by two 6-membered rings in spirocyclic borate diesters exhibit very low isotropic chemical shifts ranging from 4 to -2 ppm30 and that boron atoms shared by a 6-membered ring and a 5membered ring have chemical shifts of 9-11 ppm.32-34 Boron atoms shared by two 5-membered rings resonate at higher frequencies (12-14 ppm).27

ACS Paragon Plus Environment

15

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 36

As exactly the same trends are reproduced by high-level quantum chemical calculations of 11

B NMR parameters,30-35 we optimized the geometries of a set of different conformations of

ixazomib citrate complexes by the RI-MP2/TZVP method and calculated their isotropic chemical shifts (Figure 5). We found that the boron atoms shared by two 5-membered rings in ixazomib citrate complexes resonate at ca. 13±0.5 ppm, depending on the global conformation. In contrast, the presence of a 6-membered ring considerably decreases the calculated isotropic chemical shifts, to ca. 10±0.5 ppm. Boron atoms in ixazomib citrate Forms II and X are therefore incorporated into spirocyclic borate fragments consisting of two 5-membered rings.

IXA-Model-1 5-/5-membered

IXA-Model-2 6-/5-membered

δι(11B) = 12.9 ppm

δι(11B) = 11.0 ppm

δι(11B) = 12.7 ppm

δι(11B) = 13.2 ppm

δι(11B) = 14.2 ppm

∆E = 19.9 kJ/mol

∆E = 11.8 kJ/mol

∆E = 16.0 kJ/mol

∆E = 0.0 kJ/mol

∆E = 8.9 kJ/mol

IXA-Model-3 5-/5-membered

IXA-Model-4 5-/5-membered

IXA-Model-5 5-/5-membered

Figure 5. Typical Examples of RI-MP2/TZVP Optimized Geometries of Ixazomib Citrate Complexes with the 5-/5- and 5-/6-Membered Cyclic Borate Fragments. The corresponding values of calculated

11

B NMR isotropic chemical shifts (δi(11B)) and the relative energies with

respect to the most stable conformer (∆E) are also included.

In the subsequent step, we calculated absolute total energies, which allowed us to estimate relative stabilities for each conformer as demonstrated in Figure 5. Notably the most stable structure, IXA-Model-4, is fairly similar to the geometry of one of the crystallographically

ACS Paragon Plus Environment

16

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

independent ixazomib citrate complexes clipped out of the unit cell of the experimentally determined structure described in the file IXA Form II (the r.m.s.d. of aligned positions is 1.664 Å), whereas the second most stable structure, IXA-Model-5, is similar to the geometry of one such complex from IXA Form X, with an r.m.s.d. of aligned positions of 1.379 Å. It thus appears that in this particular case the most stable configurations consisting of two 5-membered rings are favored. Moreover, the structures IXA-Model-4 and IXA-Model-5 are also clearly suitable models for structure refinement of XRPD data.

Measurements of long-range intermolecular distances.

11

B nuclei, owing to their high

gyromagnetic ratio and high natural isotopic abundance, are particularly suited to probe intermolecular dipolar contacts. As demonstrated previously, the evolution of

11

B-11B double-

quantum (DQ) coherences allow determining 11B…11B distances up to ca. 4 Å.32 In the case of ixazomib citrate, the build-up of 11B-11B DQ coherence shows that the maximum signal intensity is reached at ca. tm = 5.6 ms for both polymorphic forms (Figure 6). By applying the previously derived relation,28 r = 0.23tm0.38, where r represents the 11B…11B interatomic distance and tm is the recoupling time for reaching maximum signal intensity, the obtained tm values indicate that the shortest

11

B…11B interatomic distances are typically 6.3-6.5 Å. Although this estimation is

approximate, the obtained distances correspond well with the values of 6.2-6.4 Å derived from the refined XRPD crystal structures (Figure 2).

ACS Paragon Plus Environment

17

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

11

Page 18 of 36

11

B- B DQ coherence build-ups

I, %

2.5 Å Borax

IXA Form-II

3.6 Å H3BO3

IXA Form-X

t, µs 0.38

r, Å

r = 0.23tm

IXA Form-X IXA Form-II

6.4-6.5 Å 3.6 Å 2.5 Å 1.7 Å

tmax, µs

Figure 6. 11B-11B DQC Build-Ups Recorded for Ixazomib Form II (Red Line) and IXA Form X (Black Line). DQ coherence build-ups for short

11

B…11B pairs (3.6 Å and 2.5 Å) are also

displayed for comparison (upper panel). The relation between the recoupling time for maximum DQ coherence intensity and the interatomic 11B…11B distance is displayed in the bottom panel.

Furthermore, when looking for additional markers of the crystal structures, hydrogen atoms located at molecular peripheries represent sensitive receptors of intermolecular contacts, which can be traced in 1H-1H correlation spectra.36 For ixazomib citrate, 1D 1H DUMBO NMR spectra reveal the existence of two types of COOH protons in citric acid residues (Figure 4c). In the corresponding 2D 1H-1H single-quantum (SQ)/DQ MAS NMR correlation spectra (Figure 7), the high-frequency carboxyl protons (COOH(1)) resonating at 12.3 and 10.7 ppm for Forms II and X, respectively, exhibit clear autocorrelation signals, indicating the presence of strongly hydrogen-bonded COOH(1)…(1)HOOC pairs. In contrast, the low-frequency carboxyl protons (COOH(2)) show clear dipolar contact with the aromatic protons, thus revealing secondary

ACS Paragon Plus Environment

18

Page 19 of 36 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 bonding of citric acid residues with amide C=O groups attached to aromatic rings. These experimental findings are in excellent accord with the proposed crystal structures (Figure 8). 1

1

H- H SQ/DQ DUMBO MAS NMR

Ixazomib citrate Form II

Ixazomib citrate Form X

CH3

Ar(2) COOH(1)

COOH(2)

B-CH

Ar(1) NH

ppm

COOH(2) COOH(1) NH

ppm

-5

Ar(2) Ar(1)

B-CH

CH3

CH2

N-CH

0

0

5

COOH(1)…CH3

5

10

COOH(1)…CH3

10

15

COOH(2)…Ar(1)

15

ssb

20

20

COOH(2)…Ar(1) 25

COOH(1)…HOOC(1)

25

15 14 13 12 11 10

9

8

7

6

5

4

3

2

1

1

0

-1

COOH(1)…HOOC(1) 15 14 13 12 11 10

ppm

9

8

7

6

5

4

3

2

1

ppm

CH3

CH-OH

CH-OH

COOH(2)

-1

13

H- C HETCOR MAS NMR

Ar(1)

COOH(1)

0

Ar(1)

COOH(1)

ppm

ppm

0

0

CH3

2

2 4

CH-OH

4

=C-H(Ar1)

6

6

8

=C-H(Ar2)

8

10

COOH(2) 12

COOH(2) 10

COOH(1)

14 180

160

140

120

100

80

60

40

20

ppm

COOH(1)

12 180

160

140

120

100

80

60

40

20

ppm

Figure 7. 1H-1H DQ/SQ DUMBO MAS NMR and 1H-13C HETCOR MAS NMR Spectra of Ixazomib Citrate Form II and X.

ACS Paragon Plus Environment

19

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

Ixazomib citrate Form II

Page 20 of 36

Ixazomib citrate Form X

Weak H-bonding COOH(2)…O=C(Ar)

Strong H-bonding COOH(1)…(1)HOOC

Figure 8. Hydrogen-Bonding Motifs in Ixazomib Citrate Forms II and X.

Verifying the global molecular arrangement. As shown recently, solid-state NMR spectra also represent fingerprints from which complete crystallographic data can be extracted.37-39 First, solid-state NMR spectra, via the site-specific structural receptivity, allow for identification of different polymorphic forms and determination of the number of symmetry independent molecules in the unit cell. In our particular case, the signal splitting observed in

13

C CP/MAS

NMR spectra (Figure 4b; and Supporting Information S1, Figure S2) shows that the unit cells of both polymorphic modifications consist of two inequivalent ixazomib molecules and two inequivalent molecules of citric acid. The 11B TQ/MAS NMR spectra (Figure 4a) then reveal that the nonequivalence of boron centers is relatively small for Form II, whereas much more substantial differences in the coordination geometries exist in ixazomib citrate Form X. These findings thus clearly support the generated structural models. Second, systematic comparison of experimentally determined 1H NMR isotropic chemical shifts with such shifts calculated with DFT from the corresponding crystal structure models

ACS Paragon Plus Environment

20

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

allows selection of the most suitable candidates.37-39 Agreement between theory and experiment was analyzed using linear regression when the quantitative measure was the r.m.s.d. From the literature data, the correct structure is clearly characterized by the smallest 1H r.m.s.d., less than 0.5 ppm. Although 13C chemical shifts were found to be less sensitive, the 13C r.m.s.d. indicating the correct crystal structure was always smaller than 2.5 ppm. Therefore, we used the same approach to validate refined XRPD crystal structures in the next step. First, the 1H and

13

C isotropic chemical shifts were calculated for both refined XRPD

models (Forms II and X) and then compared with the corresponding experimental data extracted from 2D 1H-13C HETCOR and 1H-1H DQ/SQ DUMBO MAS NMR spectra (Figure 7). Due to the splitting of

13

C NMR resonances caused by the presence of two symmetry-inequivalent

molecules in the unit cells, the signal assignment was made on the basis of an analysis of 1H-13C and 1H-1H correlation signals with the aid of automated prediction of 1H-13C correlation patterns.40 The obtained similarity parameters are summarized in Table 4, while further details are given in Supporting Information SI5.

Table 4. Comparison of Experimental and Calculated 1H and

13

C NMR Parameters for the

Refined XRPD Crystal Structures of Ixazomib Citrate Forms II and X. NMR similarity parameters Form II

Form X

r.m.s.d. (13C, COOH excl.), ppm

1.18

1.85

r.m.s.d. (1H, OH/NH excl.), ppm

0.47

0.40

r.m.s.d. (13C, all), ppm

1.99

2.55

r.m.s.d. (1H, all), ppm

0.65

0.93

covariance for 1H-13C pairs, ppm2

0.24

0.11

ACS Paragon Plus Environment

21

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 36

While determining the precise location of labile hydrogen atoms is difficult, even from synchrotron powder diffraction data, we excluded these labile species (COOH, OH and NH) from the comparative analysis. We thus obtained similarity parameters that clearly fall below the threshold limits, confirming that the refined XRPD backbone conformations and global molecular packing of ixazomib citrate data are consistent with the experimentally determined solid-state NMR data. However, when the labile species COOH, OH and NH are included in the comparative analysis, the r.m.s.d.s increased, reflecting certain inconsistencies in the locations of hydrogen-bonded protons. Whereas refined XRPD structures predict very strong hydrogen bonding between citric acids (COOH(1)…(1)HOOC), which is typically characterized by high isotropic chemical shifts (reaching ca. 16 ppm; 15.5 and 14.5 ppm were obtained from DFT calculations for Forms II and X, respectively), the corresponding experimental values are considerably lower (12.3 and 10.7 ppm), indicating weaker hydrogen bonding and longer intermolecular distances in these molecular fragments or the presence of thermal motions. Analysis of 1H NMR spectra thus provides additional input which may possibly improve structure determination from XRPD data.

Conclusions Boron-containing compounds have long been recognized as potentially active pharmaceutical ingredients. Uncertainties in the structure determination of these peptide analogs, however, represent an obstacle in developing the next generation of anticancer drugs. Determination of the atomic-resolution structures of these multicomponent solids in the absence of single crystals remains a challenge. When focusing on generally improving structure determination from X-ray powder diffraction data, we considered several aspects of current experimental and

ACS Paragon Plus Environment

22

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

computational approaches. As a result of this analysis, we formulated and successfully applied an extended experimental-computational strategy to explore the crystal structures of complex peptide analogs of boronic acid derivatives, even those exhibiting extensive conformational variability. By combining an advanced analysis of powder diffraction data with a comparative analysis of NMR parameters and high-level quantum chemical calculations, we solved the complete crystal structures of two previously undescribed non-solvated polymorphs of ixazomib citrate. As both polymorphic forms consisted of two molecules in the asymmetric unit cell and thus exhibited 32 degrees of conformational freedom, solving the structure from powder diffraction data represented a challenging problem on the edge of the capabilities of current methodology. Consequently, a novel two-step Rietveld refinement based on DFT-D restraints was used to improve information quality derived from powder diffraction data. NMR crystallography was applied to verify the refined crystal structures. In this regard, the high potential value of 11B NMR parameters was demonstrated. Whereas the evolution of 11B-11B DQ coherences allows probing long-range molecular arrangements with interatomic distances up to ca. 7 Å, a comparative analysis of NMR parameters that were experimentally determined and calculated for the refined crystal structures gives strong support for the accuracy of the XRPDdetermined structures. Overall, an integrated approach that applies several techniques in conjunction, thus providing otherwise unavailable structural information, was developed and successfully demonstrated, resulting in an efficient tool for further advancement of highly active anticancer drugs based on peptidic analogs of boronic acid.

ACS Paragon Plus Environment

23

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 36

Experimental Section Materials. Ixazomib citrate polymorphs Teva Form II and Teva Form X were obtained from Teva Pharmaceuticals CR (Czech Republic). The applied labeling of the investigated polymorphs strictly follows the terminology of Teva Pharmaceutical Industries Ltd. used in the process of polymorphs screening. Solid-state NMR spectroscopy. Solid-state NMR spectra were measured at 11.7 T using a Bruker Avance III HD 500 US/WB NMR spectrometer (Karlsruhe, Germany, 2013). The following techniques were applied: i) one-dimensional (1D)

11

B and 13C MAS, CP/MAS, and

CPPI/MAS NMR experiments;41 ii) two-dimensional (2D) experiments;42 iii) triple-quantum (TQ) experiments (2D

11

1

H-13C FSLG HETCOR

B TQ/MAS NMR);43 iv)

11

B-11B

double-quantum (DQ) experiments with the BR212 recoupling sequence;44,45 and v) 2D DQ/SQ 1

H-1H MAS NMR experiment46 with SPC5 DQ recoupling47 and DUMBO homodecoupling.48

Frictional heating49,50 of the spinning samples was offset by active cooling, and the temperature calibration was performed with Pb(NO3)2. For all experimental details, see Supporting Information S6. X-ray synchrotron powder data measurement. Powdered samples were placed in a thin-walled 1/1.5 mm capillary (Hilgenberg low-absorption borosilicate glass No. 50). The powder diffraction measurement was made at the ESRF synchrotron radiation source in Grenoble. Data were collected on the high-resolution powder X-ray diffraction beamline ID22 at wavelengths of 0.399978(2) Å (Form II) and 0.399963(3) Å (Form X) and were calibrated against NIST 640c standard silicon. The sample was spun on the axis of the diffractometer to optimize the powder averaging. Measurements were made at room temperature in continuous-scanning mode, from -5 to 35 degrees for 2theta (Form II) and from -5 to 28 degrees for 2theta (Form X) at 10 degrees

ACS Paragon Plus Environment

24

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

per minute, sampling the diffraction pattern every 3 ms (0.0005 degrees). As a precaution against radiation damage, the sample was translated axially by 1.35 mm after each scan to expose fresh sample to the beam. Twenty-six scans were normalized, summed and rebinned into steps of 0.002 degrees, taking account of the exact offsets and efficiencies of the nine channels of detection and the decay in the storage ring current during the acquisition via an incident-beam monitor.7 Computational methods for isolated molecules. The in vacuo potential energy surface (PES) of ixazomib citrate was scanned in the part of the conformational space relevant for the coordination of the boron moiety. Thus, an extensive search for the minima of the PES of ixazomib citrate was performed using the DFT-based B3LYP/6-311G** method by employing numerous starting structures (not shown). This search yielded five minimal geometries that were subject to geometrical optimization at the resolution-of-the-identity (RI) second-order Møller– Plesset (MP2) level of quantum chemical theory. The resulting RI-MP2/TZVP structures were used to approximate their complete basis set (CBS) total energies (see Supporting Information SI3 for details). The CBS-extrapolated MP2 energy was then combined with the zero-point vibrational energy obtained for the corresponding B3LYP/6-311G** minimum to arrive at estimates of the absolute total energy, E. The Gaussian 09 and TURBOMOLE program packages were used for the DFT and MP2 calculations, respectively. Periodic computational methods. The obtained refined XRPD crystal structure geometries for both Forms II and X were the input for the NMR chemical shielding calculations, which were carried out by combining the PBE functional with the gauge-including projector augmented wave (GIPAW) method51,52 as implemented in the CASTEP-NMR module.12 In all calculations (the geometrical optimizations, followed by the calculations of the NMR parameters), the

ACS Paragon Plus Environment

25

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 36

corresponding “Fine” level of the CASTEP settings was adopted. Additional computational details are provided in Supporting Information SI4.

ASSOCIATED CONTENT Supporting Information. This material is available free of charge via the Internet at http://pubs.acs.org. Full-size diffractograms; Computational methods; Details of calculations of the total energy in vacuo; Details of DFT calculations of isotropic chemical shifts; The HETCOR spectra peakassignments aided with DFT calculations; Solid-State NMR Spectroscopy CIF file for XRPD refined crystal structure of Ixazomib Citrate Forms II and X (CSD Code CCDC 1822510 and 1822511). Accession Codes CCDC 1822510 and 1822511 contain the supplementary crystallographic data for this paper. These data can be obtained free of charge via www.ccdc.cam.ac.uk/data_request/cif, or by emailing [email protected], or by contacting The Cambridge Crystallographic Data Centre, 12, Union Road, Cambridge CB2 1EZ, UK; fax: +44 1223 336033.

AUTHOR INFORMATION Corresponding Author *Email: [email protected]

ACS Paragon Plus Environment

26

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

Author Contributions The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript. Notes The authors declare no competing financial interests.

ACKNOWLEDGMENT The authors thank the Czech Science Foundation (grant nos. GA16-10035S and GA1604109S) and the Ministry of Education, Youth and Sports of the Czech Republic within the National Sustainability Program I (NPU I), Project LO1507 POLYMAT, for their financial support. Computational resources were partially provided under program LM2010005, and through the Centre CERIT Scientific Cloud, part of the Operational Program Research and Development for Innovations, reg. no. CZ.1.05/3.2.00/08.0144.

REFERENCES (1)

Bryan, J.; Kantarjian, H.; Jabbour, E. Decitabine for the treatment of acute myeloid

leukemia. Expert Opin. Orphan Drugs 2013, 1, 661–673. (2)

Gnyszka, A.; Jastrzebski, Z.; Flis, S. DNA Methyltransferase Inhibitors and Their

Emerging Role in Epigenetic Therapy of Cancer. Anticancer Res. 2013, 33, 2989–2996.

ACS Paragon Plus Environment

27

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

(3)

Page 28 of 36

Singh, V.; Sharma, P. ; Capalash, N. DNA Methyltransferase-1 Inhibitors as Epigenetic

Therapy for Cancer. Curr. Cancer Drug Targets 2013, 13, 379–399. (4)

Rogstad, D. K.; Herring, J. L.; Theruvathu, J. A.; Burdzy, A.; Perry, C. C.; Neidigh, J.

W.; Sowers, L. C. Chemical Decomposition of 5-Aza-2 '-deoxycytidine (Decitabine): Kinetic Analyses and Identification of Products by NMR, HPLC, and Mass Spectrometry. Chem. Res. Toxicol. 2009, 22, 1194–1204. (5)

Elliott, E. L.; Ferdous, A. J.; Kaufman, M. J.; Komar Lay, S. A.; Mazaik, D. L.;

Mccubbin, Q. J.; Nguyen, P. M.; Palaniappan, V. Boronate Ester Compounds And Pharmaceutical Compositions Thereof. US 8,859,504 B2, October 14, 2014. (6)

Skořepová, E.; Čerňa, I.; Vlasáková, R.; Zvoníček, V.; Tkadlecová, M.; Dušek, M.

Spirocyclic character of ixazomib citrate revealed by comprehensive XRD, NMR and DFT study. J. Mol. Struct. 2017, 1148, 22–27. (7)

Wright, J.; Vaughan, G.; Fitch. A. Merging data from a multi-detector continuous

scanning powder diffraction system. IUCr Computing Commission Newsletter, No. 1, 2003, pp. 92–96. (8)

David, W. I. F.; Shankland, K.; van de Streek, J.; Pidcock, E. ; Motherwell, W. D. S.;

Cole, J. C. DASH: a program for crystal structure determination from powder diffraction data. J. Appl. Crystallogr. 2006, 39, 910–915. (9)

Boultif, A.; Louër, D. Powder pattern indexing with the dichotomy method. J. Appl.

Crystallogr. 2004, 37, 724–731.

ACS Paragon Plus Environment

28

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

(10) Petříček, V.; Dušek, M.; Palatinus, L. Crystallographic Computing System JANA2006: General features Z. Kristallogr. 2014, 229, 345–352. (11) Rohlíček, J.; Hušák, M. MCE2005 - a new version of a program for fast interactive visualization of electron and similar density maps optimized for small molecules. J. Appl. Crystallogr. 2007, 40, 600–601. (12) Clark, S. J.; Segall, M. D.; Pickard, C. J.; Hasnip, P. J.; Probert, M. I. J.; Refson, K.; Payne, M. C. First principles methods using CASTEP. Zeitschrift für Kristallographie Crystalline Materials 2005, 220, 567–570. (13) van de Streek, J.; Neumann, M. A. Validation of molecular crystal structures from powder diffraction data with dispersion-corrected density functional theory (DFT-D). Acta Crystallogr., Sect. B: Struct. Sci., Cryst. Eng. Mater. 2014, 70, 1020–1032. (14) Rohlíček, J.; Skořepová, E.; Babor, M.; Čejka, J. CrystalCMP: an easy-to-use tool for fast comparison of molecular packing. J. Appl. Crystallogr. 2016, 49, 2172–2183. (15) Hughes, C.E.; Reddy, G.N.M.; Masiero, S.; Brown, S. P.; Williams, P. A.; and Harris, K.D.M. Determination of a complex crystal structure in the absence of single crystals: analysis of powder X-ray diffraction data, guided by solid-state NMR and periodic DFT calculations, reveals a new 20-deoxyguanosine structural motif. Chem. Sci. 2017, 8, 3971-3979. (16) NMR crystallography; Harris, R.K., Wasylishen, R.E., Duer, M.J., Eds.; John Wiley & Sons, Ltd. 2009.

ACS Paragon Plus Environment

29

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 36

(17) Brus, J.; Jegorov, A. Through-bonds and through-space solid-state NMR correlations at natural isotopic abundance: Signal assignment and structural study of simvastatin. J. Phys. Chem. A 2004, 108 (18), 3955-3964. (18) Aguilar-Granda, A.; Perez-Estrada, S.; Roa, A. E.; Rodriguez-Hernandez, J.; HernandezOrtega, S.; Rodriguez, M.; Rodriguez-Molina, B. Synthesis of a Carbazole-[pi]-carbazole Molecular Rotor with Fast Solid State Intramolecular Dynamics and Crystallization-Induced Emission. Cryst. Growth Des. 2016, 16 (6), 3435-3442. (19) Geppi, M.; Mollica, G.; Borsacchi, S.; Veracini, C.A. Solid-state NMR studies of pharmaceutical systems. Appl. Spectrosc. Rev. 2008, 43, 202-302. (20) Putra, O. D.; Yoshida, T.; Umeda, D.; Higashi, K.; Uekusa, H.; and Yonemochi, E. Crystal Structure Determination of Dimenhydrinate after More than 60 Years: Solving SaltCocrystal Ambiguity via Solid-State Characterizations and Solubility Study. Cryst Growth Des 2016, 16 (9), 5223−5229. (21) Costantino, F; Sassi, P.; Geppi, M.; Taddei, M. Multitechnique Experimental Insight on an Unusual

Crystal-to-Crystal

High

Temperature

Solid

State

Reaction

in

Zirconium

Carboxypyridinephosphonates: From One-Dimensional Chains to Two-Dimensional Hybrid Layers Through HF Elimination. Cryst. Growth Des. 2012, 12, 5462-5470. (22) Reddy, G. N. M.; Marsh, A.; Davis, J. T.; Masiero, S.; Brown, S. P. Interplay of Noncovalent Interactions in Ribbon-like Guanosine Self-Assembly: An NMR Crystallography Study. Cryst. Growth Des. 2015, 15 (12), 5945-5954.

ACS Paragon Plus Environment

30

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

(23) Bryce, D.L.; Taulelle, F., NMR Crystallography, Acta Crystallographica Section C: Structural Chemistry 2017, 73, 126-127. (24) Bryce, D.L., NMR crystallography: Structure and properties of materials from solid-state nuclear magnetic resonance observables, IUCrJ 2017, 4, 350-359. (25) Beran, G.J.O., Hartman, J.D., Heit, Y.N. Predicting Molecular Crystal Properties from First Principles: Finite-Temperature Thermochemistry to NMR Crystallography, Accounts of Chemical Research 2017, 49, 2501-2508. (26) Perras, F. A.; Bryce, D. L. Measuring dipolar and J coupling between quadrupolar nuclei using double-rotation NMR. J. Chem. Phys. 2013, 138, 174202. (27) Peters, G. M.; Skala, L. P.; Plank, T. N.; Hyman, B. J.; Reddy, G. N. M.; Marsh, A.; Brown, S. P.; Davis, J. T. A G(4)center dot K+ Hydrogel Stabilized by an Anion. J. Am. Chem. Soc. 2014, 136, 12596–12599. (28) Feng, N.; Wang, Q.; Zheng, A.; Zhang, Z.; Fan, J.; Liu, S.-B.; Amoureux, J.-P.; Deng, F. Understanding the High Photocatalytic Activity of (B, Ag)-Codoped TiO2 under Solar-Light Irradiation with XPS, Solid-State NMR, and DFT Calculations. J. Am. Chem. Soc. 2013, 135, 1607–1616. (29) Vogt, F. G.; Williams, G. R.; Copley, R. C. B. Solid-State NMR Analysis of a BoronContaining Pharmaceutical Hydrochloride Salt. J. Pharm. Sci. 2013, 102, 3705–3716. (30) Kobera, L.; Czernek, J.; Strečková, M.; Urbanova, M.; Abbrent, S.; Brus, J. Structure and Distribution of Cross-Links in Boron-Modified Phenol-Formaldehyde Resins Designed for Soft

ACS Paragon Plus Environment

31

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 32 of 36

Magnetic Composites: A Multiple-Quantum B-11-B-11 MAS NMR Correlation Spectroscopy Study. Macromolecules 2015, 48, 4874–4881. (31) Angel Wong, Y.-T.; Bryce, D.L. Recent Advances in

11

B Solid-State Nuclear Magnetic

Resonance Spectroscopy of Crystalline Solids, Annual Reports on NMR Spectroscopy 2018, 93, 145-212. (32) Brus, J.; Czernek, J.; Kobera, L.; Urbanová, M.; Jegorov, A. An efficient 2D 11B–11B solid-state NMR spectroscopy strategy for monitoring covalent self-assembly of boronic acidderived com-pounds: the transformation and unique architecture of bortezomib molecules in the solid state. Phys. Chem. Chem. Phys. 2017, 19, 487-495. (33) Matejicek, P.; Brus, J.; Jigounov, A.; Plestil, J.; Uchman, M.; Prochazka, K.; Gradzielski, M. On the structure of polymeric composite of metallacarborane with poly(ethylene oxide). Macromolecules 2011, 44 (10), 3847-3855. (34) Brus, J.; Zhigunov, A.; Czernek, J.; Kobera, L.; Uchman, M.; Matejicek, P. Control over the self-assembly and dynamics of metallacarborane nanorotors by the nature of the polymer matrix: A solid-state NMR study. Macromolecules 2014, 47 (18), 6343-6354. (35) Czernek, J.; Brus, J. On the predictions of the B-11 solid state NMR parameters. Chem. Phys. Lett. 2016, 655-656, 66–70. (36) Brown, S. P. Applications of high-resolution H-1 solid-state NMR. Solid State Nucl. Magn. Reson. 2012, 41, 1–27.

ACS Paragon Plus Environment

32

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

(37) Baias, M.; Dumez, J. N.; Svensson, P. H.; Schantz, S.; Day, G. M.; Emsley, L. De Novo Determination of the Crystal Structure of a Large Drug Molecule by Crystal Structure Prediction-Based Powder NMR Crystallography. J. Am. Chem. Soc. 2013, 135, 17501-17507. (38) Baias, M.; Widdifield, C. M.; Dumez, J. N.; Thompson, H. P. G.; Cooper, T. G.; Salager, E.; Bassil, S.; Stein, R. S.; Lesage, A.; Day, G. M.; Emsley, L. Powder crystallography of pharmaceutical materials by combined crystal structure prediction and solid-state 1H NMR spectroscopy. Phys. Chem. Chem. Phys. 2013, 15, 8069-8080. (39) Brus, J.; Czernek, J.; Kobera, L.; Urbanová, M.; Abbrent, S.; Husak, M. Predicting the Crystal Structure of Decitabine by Powder NMR Crystallography: Influence of Long-Range Molecular Packing Symmetry on NMR Parameters. Cryst. Growth Des. 2016, 16, 7102–7111. (40) Czernek, J.; Brus, J. The covariance of the differences between experimental and theoretical chemical shifts as an aid for assigning two-dimensional heteronuclear correlation solid-state NMR spectra. Chem. Phys. Lett. 2014, 608, 334-339. (41) Wu, X. L.; Burns, S. T.; Zilm, K. W. Spectral Editing in CPMAS NMR. Generating Subspectra Based on Proton Multiplicities. J. Magn. Reson. Ser. A 1994, 111, 29-36. (42) vanRossum, B. J.; Forster, H.; deGroot, H. J. M. High-Field and High-Speed CPMAS13C NMR Heteronuclear Dipolar-Correlation Spectroscopy of Solids with FrequencySwitched Lee–Goldburg Homonuclear Decoupling. J Magn. Reson. 1997, 124, 516-519. (43) Amoureux, J.-P.; Fernandez, C.; Steuernagel, S. Z filtering in MQMAS NMR. J. Magn. Reson., Ser. A 1996, 123, 116–118.

ACS Paragon Plus Environment

33

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 36

(44) Edén, M.; Zhou, D.; Yu, J. Improved double-quantum NMR correlation spectroscopy of dipolar-coupled quadrupolar spins. Chem. Phys. Lett. 2006, 431, 397–403. (45) Wang, Q.; Hu, B.; Lafon, O.; Trébosc, J.; Deng, F.; Amoureux, J. P. Double-quantum homonuclear NMR correlation spectroscopy of quadrupolar nuclei subjected to magic-angle spinning and high magnetic field. J. Magn. Reson. 2009, 200, 251–260. (46) Brown, S. P.; Lesage, A.; Elena, B.; Emsley, L. Probing Proton−Proton Proximities in the Solid State:  High-Resolution Two-Dimensional 1H−1H Double-Quantum CRAMPS NMR Spectroscopy. J. Am. Chem. Soc. 2004, 126, 13230-13231. (47) Hohwy, M.; Rienstra, C. M.; Jaroniec, C. P.; Griffin, R. G. Fivefold symmetric homonuclear dipolar recoupling in rotating solids: Application to double quantum spectroscopy. J. Chem. Phys. 1999, 110, 7983-7992. (48) Sakellariou, D.; Lesage, A.; Emsley, L. Proton−Proton Constraints in Powdered Solids from 1H−1H−1H and 1H−1H−13C Three-Dimensional NMR Chemical Shift Correlation Spectroscopy. J. Am. Chem. Soc. 2001, 123, 5604-5605. (49) Langer, B.; Schnell, L.; Spiess, H. W.; Grimmer, A. R. Temperature Calibration under Ultrafast MAS Conditions. J. Magn. Reson. 1999, 138, 182-186. (50) Brus, J. Heating of samples induced by fast magic-angle spinning. Solid State Nucl. Magn. Reson. 2000, 16, 151-160. (51) Pickard, C. J.; Mauri, F. All-electron magnetic response with pseudopotentials: NMR chemical shifts. Phys. Rev. B 2001, 63, 245101.

ACS Paragon Plus Environment

34

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

Crystal Growth & Design

(52) Yates, J. R.; Pickard, C. J.; Mauri, F. Calculation of NMR chemical shifts for extended systems using ultrasoft pseudopotentials. Phys. Rev. B 2007, 76, 024401.

ACS Paragon Plus Environment

35

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 36

For Table of Contents Use Only

Determining the crystal structures of peptide analogs of boronic acid in the absence of single crystals: Intricate motifs of ixazomib citrate revealed by XRPD guided by ss-NMR

Michal Hušák, Alexandr Jegorov, Jan Rohlíček, Andrew Fitch, Jiří Czernek, Libor Kobera, and Jiří Brus*

General experimental-computational strategy for exploring the assembly of peptidic boronic acid derivatives in microcrystalline state is demonstrated. By combining the simulated annealing and two-step Rietveld refinement of XRPD data with the comparative analysis of NMR parameters the crystal structures of two non-solvated polymorphs of ixazomib citrate, each consisting of two symmetry-independent molecules in unit cell with 32 free parameters, were solved and refined.

ACS Paragon Plus Environment

36