I. INTRODUCTION
Levomefolic acid is a metabolite of folic acid (Vitamin B9) and is a major active form of folate found in foods and in the blood circulation. It is transported across membranes, including the blood–brain barrier, and plays an essential role in DNA and protein synthesis. Levomefolate is approved as a food additive and is designated as a GRAS (generally recognized as safe) compound. It is available commercially as a crystalline Ca salt (trade name Metafolin), which has the required stability for use as a supplement (https://www.drugbank.ca/salts/DBSALY001276). The IUPAC name (CAS Registry number 151533-22-1 for the anhydrous salt) is (2S)-2-[[4-[[(6S)-2-amino-5-methyl-4-oxo-3,6,7,8-tetrahydropteridin-6-yl]methylamino]benzoyl]amino]pentanedioate calcium trihydrate. A two-dimensional molecular diagram is shown in Figure 1.
The 2D molecular structure of calcium L-5-methyltetrahydrofolate.
Stable crystalline salts of 5-methyltetrahydrofolic acid are disclosed and claimed in US Patent 6,441,168 (Müller et al., Reference Müller, Moser and Egger2002; Eprova AG). The patent includes claims for “a water of crystallization of at least one equivalent per equivalent of 5-methyltetrahydrofolic acid” and “≥3 equivalents of water”. Commercial samples are generally described as the trihydrate, although the pentahydrate is also available commercially. Powder diffraction data for additional (6S)-5-methyltetrahydrofolate calcium salts are reported in Chinese Patent CN 104530051 A (Wang and Cheng, Reference Wang and Cheng2015; Beijing Jinkang Hexin Pharmaceutical Technology Co.). A connectivity search in the Cambridge Structural Database (Groom et al., Reference Groom, Bruno, Lightfoot and Ward2016) for derivatives of levomefolic acid yielded no hits. Name and formula searches in the Powder Diffraction File (Gates-Rector and Blanton, Reference Gates-Rector and Blanton2019) yielded no hits.
II. EXPERIMENTAL
Calcium L-5-methyltetrahydrofolate trihydrate (Lot #WN02-130227) was supplied by Virtus Pharmaceuticals. A laboratory pattern (measured on a Bruker D2 Phaser using Cu Kα radiation; Kα 1 = 1.540593 Å, Kα 2 = 1.544451 Å) could be indexed on a primitive monoclinic unit cell with a = 13.314(3), b = 26.868(6), c = 8.187(2) Å, β = 94.33(2)°, V = 2563.6(8) Å3, and Z = 4. This cell predicted a peak at 3.3° 2θ, which was confirmed by changing the configuration of the diffractometer. Attempts to solve the structure using this cell were unsuccessful.
The same pale yellow 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.413691(2) Å from 0.5° to 50° 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 twelve silicon crystal analyzers that allow 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 synchrotron pattern was indexed with difficulty using DICVOL14 (Louër and Boultif, Reference Louër and Boultif2014) on a primitive monoclinic unit cell having a = 6.9575, b = 6.5372, c = 53.7699 Å, β = 92.320°, V = 2444.33 Å3, and Z = 4. Analysis of the systematic absences using EXPO2014 (Altomare et al., Reference Altomare, Cuocci, Giacovazzo, Moliterni, Rizzi, Corriero and Falcicchio2013) suggested the space group P21. Attempts to solve the structure using this cell yielded some plausible structures, but they all contained some molecular overlap. Analysis of this cell using PLATON (Spek, Reference Spek2009, Reference Spek2020) showed that it was not the conventional monoclinic cell, which has a = 7.202, b = 6.548, c = 53.837 Å, and β = 90.5 4°. The fact that the β angle was close to 90° suggested that we explore orthorhombic unit cells. EXPO2014 yielded a better Le Bail fit using an orthorhombic cell, and suggested the space group P 2 1 2 1 2 1, which has a 4-fold general position. A reduced cell search in the Cambridge Structural Database (Groom et al., Reference Groom, Bruno, Lightfoot and Ward2016) yielded 36 hits, but no folate derivatives.
The levomefolic acid molecule was downloaded from PubChem (Kim et al., Reference Kim, Chen, Cheng, Gindulyte, He, He, Li, Shoemaker, Thiessen, Yu, Zaslavsky, Zhang and Bolton2023) as Conformer3D_CID_135398561.mol2. In this molecule, both C14 and C25 have the S stereochemistry. The structure was difficult to solve, using several different strategies and programs. The successful solution came by using FOX (Favre-Nicolin and Černý, Reference Favre-Nicolin and Černý2002). The maximum sinθ/λ used was 0.25 Å−1. The anion (converted to a Fenske-Hall Z-matrix using OpenBabel (O'Boyle et al., Reference O'Boyle, Banck, James, Morley, Vandermeersch and Hutchison2011)), a Ca atom, and 3 O atoms (for the water molecules) were used as fragments, along with 001 preferred orientation (the extreme anisotropy of the unit cell makes preferred orientation likely) and a bump penalty (which increases the cost factor when pairs of atoms become closer than specified distances). In the lowest-cost solution, two of the O atoms were too close to each other, so one was removed from the model. The third O was added in a void, detected by Mercury (Macrae et al., Reference Macrae, Sovago, Cottrell, Galek, McCabe, Pidcock, Platings, Shields, Stevens, Towler and Wood2020). The default mode of FOX includes enough molecular flexibility to invert the chirality of a carbon atom occasionally, as was the case here. In the lowest-cost model, C25 was R, so H47 was removed and reintroduced to the other side of this C atom using Materials Studio (Dassault, 2021). This model was subjected to a molecular mechanics optimization (fixed unit cell) using the Forcite module of Materials Studio. This optimized model was the start of a DFT optimization using CRYSTAL14 (Dovesi et al., Reference Dovesi, Erba, Orlando, Zicovich-Wilson, Civalleri, Maschio, Rérat, Casassa, Baima, Salustro and Kirtman2018). The basis sets for the H, C, N, and O atoms were those of Gatti et al. (Reference Gatti, Saunders and Roetti1994), and the basis set for Ca was that of Peintinger et al. (Reference Peintinger, Vilela Oliveira and Bredow2013). The calculation was run on eight 2.1 GHz Xeon cores (each with 6 GB RAM) of a 304-core Dell Linux cluster at IIT, using 8 k-points and the B3LYP functional, and took ~150 h.
Rietveld refinement was carried out using GSAS-II (Toby and Von Dreele, Reference Toby and Von Dreele2013). Only the 0.7–18.0° portion of the diffraction pattern was included in the refinement (d min = 1.322 Å). All non-H bond distances and angles in the folate anion 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) of the molecule. The Mogul average and standard deviation for each quantity were used as the restraint parameters. The Ca–O distances were restrained manually at 2.48(10) Å. The restraints contributed 6.87% to the overall χ 2. The hydrogen atoms were included in calculated positions, which were recalculated during the refinement using Materials Studio. The U iso were grouped by chemical similarity, and the U iso of each H atom was constrained to be 1.3× that of the heavy atom to which it is attached. The background was modeled using a 6-term shifted Chebyshev polynomial, along with one peak at 4.56° to model the scattering from the Kapton capillary and an amorphous component. The peak for the Kapton capillary generally occurs ~5.2°.
The final refinement of 136 variables using 17 303 observations and 91 restraints yielded the residuals R wp = 0.1011 and GOF = 1.77. The largest peak (0.06 Å from Ca34) and hole (1.82 Å from C21) in the difference Fourier map were 0.17 and −0.12(3) eÅ−3, respectively. The largest errors in the fit (Figure 2) are in the shapes of some of the peaks. The data did not support refinement of the generalized microstrain model, so a uniaxial model (with 001) as the unique axis was used.
The Rietveld plot for the refinement of calcium L-5-methyltetrahydrofolate trihydrate. The blue crosses represent the observed data points, and the green line is the calculated pattern. The cyan curve is the error plot, and the red line is the background curve. The vertical scale has been multiplied by a factor of 5× for 2θ > 3.0°.
Additionally, a dihydrate model (without a coordinated water molecule, and derived from an optimization using VASP (Kresse and Furthmüller, Reference Kresse and Furthmüller1996)) was refined. This model yielded poorer residuals: R wp = 0.1252 and GOF = 2.13.
III. RESULTS AND DISCUSSION
This synchrotron power pattern of calcium L-5-methyltetrahydrofolate trihydrate is in excellent agreement with that reported by Müller et al. (Reference Müller, Moser and Egger2002) (Figure 3). The agreement is good enough to conclude that the patterns represent the same material, known as Type I. We should note that the patent pattern does not include the strongest peak of the pattern, because the pattern was not measured to low-enough an angle. This is not a problem for characterizing the phase for patent purposes, but is a useful cautionary tale.
Comparison of the synchrotron pattern of calcium L-5-methyltetrahydrofolate trihydrate (magenta) to that of Form I reported by Müller et al. (Reference Müller, Moser and Egger2002; black). Note that the patent pattern does not include the lowest-angle (and strongest) peak of the pattern. The literature pattern (measured using Cu Kα radiation) was digitized using UN-SCAN-IT (Silk Scientific, 2013) and converted to the synchrotron wavelength of 0.458963(2) Å using JADE Pro (MDI, 2022). Image generated using JADE Pro (MDI, 2022).
The refined atom coordinates of L-5-methyltetrahydrofolate trihydrate Type I and the coordinates from the DFT optimization are reported in the CIFs deposited with ICDD. The root-mean-square (rms) Cartesian displacement of the non-hydrogen atoms in the Rietveld-refined and DFT-optimized structures is 0.744 Å (Figure 4), outside the normal range for correct structures (van de Streek and Neumann, Reference van de Streek and Neumann2014). The absolute difference in the Ca position in the two structures is 1.025 Å. The major differences are in the orientation of the C21–C29 phenyl ring, the orientations of the side chains, and the absolute position of the molecule in the unit cell. In the population of samples of this compound, this is an exceptionally crystalline one (hence the progress described here), but the peaks are relatively broad and the sample contains an amorphous component. As we will see later, the sample exhibits significant texture. All of these “features” mean that we can expect poorer accuracy of the refined structure than usual. Perhaps this should be considered a “proposed” structure for this compound; it is certainly better than no structure. The discussion below concentrates on the DFT-optimized structure, as we believe it is more reliable. The asymmetric unit (with atom numbering) is illustrated in Figure 5, and the crystal structure is presented in Figure 6.
Comparison of the Rietveld-refined (red) and VASP-optimized (blue) structures of calcium L-5-methyltetrahydrofolate trihydrate. The rms Cartesian displacement is 0.744 Å. Image generated using Mercury (Macrae et al., Reference Macrae, Sovago, Cottrell, Galek, McCabe, Pidcock, Platings, Shields, Stevens, Towler and Wood2020).
The asymmetric unit of calcium L-5-methyltetrahydrofolate trihydrate, with the atom numbering. The atoms are represented by 50% probability spheroids/ellipsoids. Image generated using Mercury (Macrae et al., Reference Macrae, Sovago, Cottrell, Galek, McCabe, Pidcock, Platings, Shields, Stevens, Towler and Wood2020).
The crystal structure of calcium L-5-methyltetrahydrofolate trihydrate, viewed down the b-axis. Image generated using Diamond (Crystal Impact, 2022).
The structure is characterized by alternating hydrophobic and hydrophilic layers along the c-axis (Figure 6). The Ca cations are 7-coordinate (one water molecule and six carboxylate oxygen atoms), and share edges to form chains along the b-axis (Figure 7). The bond valence sum of the Ca is 1.89. As expected, the carboxylate groups bind to the Ca. The group C32–O3–O4 bridges two Ca34, while the group C33–O5–O6 chelates to one Ca34, and bridges two others. One water molecule (O35) is coordinated to the Ca, while the other two (O36 and O37) are zeolitic. Additional small voids can be located by decreasing the probe radius in Mercury to 1.0 Å (from the default value of 1.2 Å) (Figure 8). The presence of additional water molecules in some samples can thus be rationalized easily. Attempts to refine the occupancies of water molecules placed in these voids resulted in values insignificantly different from zero, consistent with the formulation as a trihydrate.
The chains of edge-sharing CaO7 coordination polyhedra, viewed down the c-axis. Image generated using Diamond (Crystal Impact, 2022).
Potential additional voids in the crystal structure of calcium L-5-methyltetrahydrofolate trihydrate, obtained by decreasing the probe radius to 1.0 Å. Image generated using Mercury (Macrae et al., Reference Macrae, Sovago, Cottrell, Galek, McCabe, Pidcock, Platings, Shields, Stevens, Towler and Wood2020).
Thermogravimetric analysis (TGA) indicated that the sample was a hexahydrate. The TGA, however, was measured 5 years after the sample had been acquired and the synchrotron pattern measured, so it is uncertain how relevant that water content is to the specimen used to measure the data. At the time of this writing, the diffraction pattern of the sample differed from the synchrotron pattern, indicating that the sample changed over time. Comparison of the DFT-optimized trihydrate and dihydrate structures (Figure 9) reveals significant differences. The water content can apparently affect the structure significantly, perhaps resulting in the poorer-than-usual agreement of the refined and optimized structures.
Comparison of the DFT-optimized crystal structure of calcium L-5-methyltetrahydrofolate trihydrate (blue) to that of the dihydrate (green). Image generated using Mercury (Macrae et al., Reference Macrae, Sovago, Cottrell, Galek, McCabe, Pidcock, Platings, Shields, Stevens, Towler and Wood2020).
Almost all of the bond distances, 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 O57–C28–C30 angle of 109.3° (average = 105.8(10)°, Z-score = 3.6) is flagged as unusual. This angle lies slightly outside a narrow distribution of a few similar angles. The torsion angle C9–O54–C11–C12 is flagged as unusual; this represents the linking of the two portion of the molecule, and it is likely that crystal packing forces influence the molecular conformation. The torsion angles involving rotation around the C37–C38 bond (such as C33–C37–C38–C39) are flagged as unusual. This angle lies near a minor gauche population of mainly trans angles.
Quantum chemical geometry optimization of the isolated anion (DFT/B3LYP/6-31G*/water) using Spartan ‘18 (Wavefunction, 2020) indicated that the observed conformation of the cation is 36.5 kcal/mol higher in energy than the local minimum conformation. The conformational differences are spread throughout the molecule. The global minimum-energy conformation of the anion (MMFF force field) curls up on itself so that the ring systems are roughly parallel, and intramolecular N–H⋯O hydrogen bonds form between the amino group N13 and carboxylate oxygen atoms. Intermolecular interactions are thus important in determining the solid-state conformation.
Analysis of the contributions to the total crystal energy of the structure using the Forcite module of Materials Studio (Dassault Systèmes, 2021) suggests that angle distortion terms dominate the intramolecular deformation energy, but that bond and torsion distortion terms are also significant. The intermolecular energy is dominated by electrostatic attractions, which in this force field analysis also include hydrogen bonds. The hydrogen bonds are better analyzed using the results of the DFT calculation.
Hydrogen bonds are prominent in the crystal structure (Table I). Each of the water molecules acts as a donor in two hydrogen bonds. The coordinated water molecule O35 makes two strong intermolecular O–H⋯O hydrogen bonds to carboxyl and carbonyl groups. The zeolitic water molecules O36 and O37 form weaker hydrogen bonds, to carbonyl O atoms, ring N atoms, and aromatic C atoms. Several N–H⋯O/N hydrogen bonds, as well as C–H⋯O hydrogen bonds, also contribute to the lattice energy. The energies of the O–H⋯O hydrogen bonds were calculated using the correlation of Rammohan and Kaduk (Reference Rammohan and Kaduk2018), and the energies of the N–H⋯O hydrogen bonds were calculated using the correlation of Wheatley and Kaduk (Reference Wheatley and Kaduk2019).
Hydrogen bonds (CRYSTAL14) in calcium L-5-methyltetrahydrofolate trihydrate Type I.
a Intramolecular.
The Bravais–Friedel–Donnay–Harker (Bravais, Reference Bravais1866; Friedel, Reference Friedel1907; Donnay and Harker, Reference Donnay and Harker1937) morphology suggests that we might expect platy morphology for calcium L-5-methyltetrahydrofolate trihydrate, with {002} as the principal faces. A sixth-order spherical harmonic model was included in the refinement. The refined texture index was 1.563(10), indicating that preferred orientation was significant in this rotated capillary specimen. We should expect that preferred orientation would be significant in most specimens of this material, especially on a Bragg-Brentano diffractometer.
IV. DEPOSITED DATA
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
The 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. We thank Lynn Ribaud and Saul Lapidus for their assistance in the data collection, Andrey Rogachev for the use of computing resources at IIT, and Rhett Daniels of Virtus Pharmaceuticals for the sample.
CONFLICTS OF INTEREST
The authors have no conflicts of interest to declare.
Open access
