I. INTRODUCTION
Cariprazine, C21H32N4O (as the hydrochloride salt C21H33N4OCl, marketed under the trade names Vraylar and Reagila, among others), is used in the treatment of a variety of psychological disorders including schizophrenia as well as manic depression or a combination of behaviors associated with bipolar disorder. Cariprazine is included in the Top 200 Small Molecule Drugs by Retail Sales in 2022 (McGrath et al., Reference McGrath, Brichacek and Njardarson2010). The systematic name of cariprazine hydrochloride (CAS Registry Number 1083076-69-0) is 3-[4-[2-[4-(2,3-dichlorophenyl)piperazin-1-yl]ethyl]cyclohexyl]-1,1-dimethylurea hydrochloride. A two-dimensional molecular diagram of the cariprazine dication is shown in Figure 1.
Cariprazine hydrochloride Form I was claimed in US Patent 7,943,621 B2 (Czibula et al., Reference Czibula, Sebok, Greiner, Domany and Csongor2011; Richter Gideon). The dihydrochloride salt was also claimed, but we are unaware of any published powder diffraction data for the dihydrochloride. Form III of cariprazine hydrochloride is claimed in US Patent 7,829,569 B2 (Liao et al., Reference Liao, Zhu and Grill2010; Forest Laboratories).
This work was carried out as part of a project (Kaduk et al., Reference Kaduk, Crowder, Zhong, Fawcett and Suchomel2014) to determine the crystal structures of large-volume commercial pharmaceuticals and include high-quality powder diffraction data for them in the Powder Diffraction File (Kabekkodu et al., Reference Kabekkodu, Dosen and Blanton2024).
II. EXPERIMENTAL
Cariprazine hydrochloride was a commercial reagent, purchased from TargetMol (Batch #159104), and was used as-received. The white powder was packed into a 1.5 mm diameter Kapton capillary and rotated during the measurement at ~50 Hz. The powder pattern was measured at 295 K at beamline 11-BM (Antao et al., Reference Antao, Hassan, Wang, Lee and Toby2008; Lee et al., Reference Lee, Shu, Ramanathan, Preissner, Wang, Beno, Von Dreele, Ribaud, Kurtz, Antao, Jiao and Toby2008; Wang et al., Reference Wang, Toby, Lee, Ribaud, Antao, Kurtz, Ramanathan, Von Dreele and Beno2008) of the Advanced Photon Source at Argonne National Laboratory using a wavelength of 0.459744(2) Å from 0.5 to 40° 2θ with a step size of 0.001° and a counting time of 0.1 s/step. The high-resolution powder diffraction data were collected using 12 silicon crystal analyzers for high angular resolution, high precision, and accurate peak positions. A mixture of silicon (NIST SRM 640c) and alumina (NIST SRM 676a) standards (ratio Al2O3:Si = 2:1 by weight) was used to calibrate the instrument and refine the monochromatic wavelength used in the experiment.
The pattern was indexed using JADE Pro (MDI, 2023) on a very high-quality (F(20) = 205, JADE figure of merit = 6) primitive monoclinic unit cell with a = 27.26552, b = 7.29133, c = 12.80730 Å, β = 99.60°, V = 2510.48 Å3, and Z = 4. The suggested space group was P2 1/n, which was confirmed by successful solution and refinement of the structure. A reduced cell search of the Cambridge Structural Database (Groom et al., Reference Groom, Bruno, Lightfoot and Ward2016) yielded seven hits but no structures of cariprazine derivatives.
The structure was solved by direct methods using EXPO2014 (Altomare et al., Reference Altomare, Cuocci, Giacovazzo, Moliterni, Rizzi, Corriero and Falcicchio2013), invoking the COVMAP option. In addition to the cariprazine molecule and the Cl ion, the solution contained another atom, which was assigned as O62 (a water molecule?). The structure also contained a small void, which was filled by another O atom (O63). The occupancy of O63 refined to 0.11, and its position moved very close to Cl1, so it was removed from the refinement. The occupancy of O62 refined to 2.1, and it moved to a reasonable position for a second Cl ion, so it was renamed Cl62. The compound is thus the dihydrochloride. H63 and H64 were added to N5 and N6, respectively, based on potential hydrogen bonds.
Rietveld refinement was carried out with GSAS-II (Toby and Von Dreele, Reference Toby and Von Dreele2013). Only the 1.9–25.0° portion of the pattern was included in the refinements (d min = 1.062 Å). All non-H-bond distances and angles were subjected to restraints based on a Mercury/Mogul Geometry Check (Bruno et al., Reference Bruno, Cole, Kessler, Luo, Motherwell, Purkis, Smith, Taylor, Cooper, Harris and Orpen2004; Sykes et al., Reference Sykes, McCabe, Allen, Battle, Bruno and Wood2011). The Mogul average and standard deviation for each quantity were used as the restraint parameters. The restraints contributed 1.8% to the final χ 2. The hydrogen atoms were included in calculated positions, which were recalculated during the refinement using Materials Studio (Dassault, 2022). The Cl atoms were refined anisotropically. The Uiso of the C, N, and O atoms were grouped by chemical similarity. The Uiso for the H atoms were fixed at 1.3× the Uiso of the heavy atoms to which they are attached. The peak profiles were described using a uniaxial microstrain model, with 010 as the unique axis.
The final refinement of 136 variables using 23,137 observations and 71 restraints yielded the residuals R wp = 0.10797 and goodness of fit = 1.68. The largest peak (0.10 Å from Cl3) and hole (0.74 Å from C13) in the difference Fourier map were 0.46(9) and −0.39(9) eÅ−3, respectively. The final Rietveld plot is shown in Figure 2. The largest features in the normalized error plot represent subtle errors in peak positions and may represent changes in the specimen during the measurement.
The crystal structure of cariprazine dihydrochloride was optimized (fixed experimental unit cell) with density functional theory techniques using VASP (Kresse and Furthmüller, Reference Kresse and Furthmüller1996) through the MedeA graphical interface (Materials Design, 2016). The calculation was carried out on 16 2.4 GHz processors (each with 4 Gb RAM) of a 64-processor HP Proliant DL580 Generation 7 Linux cluster at North Central College. The calculation used the GGA-PBE functional, a plane wave cutoff energy of 400.0 eV, and a k-point spacing of 0.5 Å−1, leading to a 2 × 2 × 1 mesh, and took ~76.7 h. Single-point density functional calculations (fixed experimental cell) and population analysis were carried out using CRYSTAL23 (Erba et al., Reference Erba, Desmarais, Casassa, Civalleri, Donà, Bush, Searle, Maschio, Daga, Cossard, Ribaldone, Ascrizzi, Marana, Flament and Kirtman2023). The basis sets for the H, C, N, and O atoms in the calculation were those of Gatti et al. (Reference Gatti, Saunders and Roetti1994), and the basis set for Cl was that of Peintinger et al. (Reference Peintinger, Vilela Oliveira and Bredow2013). The calculations were run on a 3.5 GHz PC using eight k-points and the B3LYP functional and took ~5.1 h.
III. RESULTS AND DISCUSSION
Despite being purchased as “cariprazine hydrochloride”, the sample studied here is cariprazine dihydrochloride. The asymmetric unit (Figure 3) contains one cariprazine dication and two Cl anions. The root-mean-square Cartesian displacement of the non-H atoms in the Rietveld-refined and VASP-optimized cation structures is 0.188 Å (Figure 4). The agreement is within the normal range for correct structures (van de Streek and Neumann, Reference van de Streek and Neumann2014) and confirms that the structure is correct. The remainder of this discussion will emphasize the VASP-optimized structure.
Almost all of the bond distances, bond angles, and torsion angles fall within the normal ranges indicated by a Mercury Mogul Geometry check (Macrae et al., Reference Macrae, Sovago, Cottrell, Galek, McCabe, Pidcock, Platings, Shields, Stevens, Towler and Wood2020). The torsion angles involving rotation about the C8–C9 bond are flagged as unusual. They lie on the tails of peaked 0/180° distributions and indicate that the orientation of the dimethylamino group is slightly unusual.
Quantum chemical geometry optimization of the isolated cation (DFT/B3LYP/6-31G*/water) using Spartan ‘20 (Wavefunction, 2022) indicated that the observed solid-state conformation is 8.9 kcal/mol higher in energy than a local minimum, which has an almost identical conformation. The global minimum-energy conformation (MMFF force field) is 4.2 kcal/mol lower in energy than the observed conformation and is less-linear. The difference indicates that intermolecular interactions are important in determining the observed extended conformation.
The crystal structure (Figure 5) consists of layers of cations parallel to the bc-plane. The cations stack along the b-axis. The mean plane of the cation is approximately −11,−6,3. The Mercury Aromatics Analyser indicates two strong interactions between the dichlorophenyl rings in adjacent cations. Analysis of the contributions to the total crystal energy of the structure using the Forcite module of Materials Studio (Dassault Systèmes, 2022) suggests that bond, angle, and torsion distortion terms contribute about equally to the intramolecular energy. The intermolecular energy is dominated by electrostatic attractions, which, in this force field analysis, include hydrogen bonds. The hydrogen bonds are better analyzed using the results of the DFT calculation.
Each H atom on the two protonated N atoms (N5 and N6) participates in a discrete N–H⋯Cl hydrogen bond. Cl62 acts as an acceptor in two of these bonds, while Cl1 is an acceptor in only one (Table I). The result is to link the cations and anions into columns parallel to the b-axis. One C–H⋯O hydrogen bond links cations, and several C–H⋯Cl hydrogen bonds also link cations and anions. There is one intramolecular C–H⋯Cl hydrogen bond.
a Intramolecular.
The volume enclosed by the Hirshfeld surface of cariprazine dihydrochloride (Figure 6, Hirshfeld, Reference Hirshfeld1977; Spackman et al., Reference Spackman, Turner, McKinnon, Wolff, Grimwood, Jayatilaka and Spackman2021) is 618.14 Å3, which represents 98.47% of the unit cell volume. The packing density is thus fairly typical. The only significant close contacts (red in Figure 6) involve the hydrogen bonds. The volume/non-hydrogen atom is larger than the usual at 20.9 Å3, reflecting the presence of multiple Cl atoms in the compound.
The Bravais–Friedel–Donnay–Harker (Bravais, Reference Bravais1866; Friedel, Reference Friedel1907; Donnay and Harker, Reference Donnay and Harker1937) morphology suggests that we might expect elongated morphology for cariprazine dihydrochloride, with ${\langle 010}\rangle$ as the long axis. A second-order spherical harmonic model was included in the refinement. The texture index was 1.019(0), indicating that the preferred orientation was slight in this rotated capillary specimen.
IV. DEPOSITED DATA
The powder pattern of cariprazine dihydrochloride from this synchrotron data set has been submitted to the ICDD for inclusion in the Powder Diffraction File. The Crystallographic Information Framework (CIF) files containing the results of the Rietveld refinement (including the raw data) and the DFT geometry optimization were deposited with the ICDD. The data can be requested at pdj@icdd.com.
ACKNOWLEDGEMENTS
Use of the Advanced Photon Source at Argonne National Laboratory was supported by the U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences under Contract No. DE-AC02-06CH11357. This work was partially supported by the International Centre for Diffraction Data. We thank Saul Lapidus for his assistance in the data collection.
CONFLICTS OF INTEREST
The authors have no conflicts of interest to declare.