Hostname: page-component-745bb68f8f-5r2nc Total loading time: 0 Render date: 2025-02-06T08:57:45.635Z Has data issue: false hasContentIssue false

Ab initio structure determination of 3,4-diaminopyridin-1-ium dihydrogen phosphate

Published online by Cambridge University Press:  05 March 2012

A. Le Bail*
Affiliation:
Laboratoire des Oxydes et Fluorures, CNRS UMR 6010, Université du Maine, avenue O. Messiaen, 72085 Le Mans Cedex 9, France
L’. Smrčok
Affiliation:
Institute of Inorganic Chemistry, Slovak Academy of Sciences, Dúbravská cesta 9, SK-845 36 Bratislava, Slovak Republic
*
a)Electronic mail: armel.le_bail@univ-lemans.fr
Rights & Permissions [Opens in a new window]

Abstract

The structure of 3,4-diaminopyridin-1-ium dihydrogen phosphate, [C5H3(NH)(NH2)2]+ (H2PO4), is solved from conventional X-ray powder diffraction data in direct space (monoclinic unit cell with a = 16.0725(9) Å, b = 7.7301(3) Å, c = 14.6189(9) Å, β = 96.869(1)°, V = 1803.2(2) Å3, Z = 8, and space group I2/c), and optimized by energy minimization in the solid state. In the crystal structure of the title compound, dihydrogenphosphate tetrahedra are linked by strong hydrogen O-H…O bonds forming chains running parallel to the b-axis. Antiparallelly π–π stacked DAP cations form molecular columns in the spaces between the chains. Although the dominant interaction of the molecules with their surroundings is electrostatic, their bonding are further enhanced by N-H…O and C-H…O hydrogen bonds.

Type
Technical Articles
Copyright
Copyright © Cambridge University Press 2011

INTRODUCTION

The 3,4-diaminopyridine (3,4-DAP) molecule, whose crystal structure was recently established (Rubin-Preminger and Englert, Reference Rubin-Preminger and Englert2007), can be used pharmaceutically (amifampridine) for the treatment of various diseases including the Lambert-Eaton myasthenic syndrome (LEMS) (McEvoy et al., Reference McEvoy, Windebank, Daube and Low1989). Tartrate and phosphate salts of 3,4-DAP were shown to have very high stability in storage conditions complying with international recommendations contrarily to 3,4-DAP (Guyon et al., Reference Guyon, Pradeau, Le Hoang and Houri2002). The 3,4-DAP phosphate (3,4-DAPP) formulation has obtained the orphan medicinal product status both in the European Union and in the Unites States (Quartel et al., Reference Quartel, Turbeville and Lounsbury2010).

The aim of the present work is to provide an answer for the 3,4-DAPP to the usual questions for pharmaceutical compounds: Is a single phase structure established by crystallographic investigation or are there polymorphs in the sample produced? Is there formation of a salt or cocrystallization of the API (active pharmaceutical ingredient) with H3PO4? In spite of the absence of a suitable single crystal, solving the structure of 3,4-DAPP is realized from powder diffraction data by using methodologies which continuously have demonstrated their efficiency during repeated blind tests (Le Bail et al., Reference Le Bail, Cranswick, Adil, Altomare, Avdeev, Cerny, Cuocci, Giacovazzo, Halasz, Lapidus, Louwen, Moliterni, Palatinus, Rizzi, Schilder, Stephens, Stone and van Mechelen2009).

STRUCTURE DETERMINATION, RIETVELD REFINEMENTS, AND DFT STRUCTURE OPTIMIZATION

The powder diffraction pattern, recorded from a sample as received from the SERATEC Company (www.serateclab.com), could be indexed in a monoclinic cell by the MCMAILLE software (Le Bail, Reference Le Bail2004). The FoMs (figures of merit) for the pattern selected for this study were found to be M 20 = 24.7 (de Wolf, Reference de Wolf1968), F 20 = 55.4 (0.0049, 73) (Smith and Snyder, Reference Smith and Snyder1979), McM 20 = 362.0 (Le Bail, Reference Le Bail, Clearfield, Reibenspies and Bhuvanesh2008), and Z = 8. Using Le Bail fitting (Le Bail, Reference Le Bail2005) for evaluating the reliability of the indexing and for intensities extraction, the I2/c space group was determined unambiguously. The structure solution was done in the direct space using the ESPOIR software (Le Bail, Reference Le Bail2001) moving in the unit cell, by a Monte Carlo process, a 3,4-diaminopyridine molecule and a PO4 tetrahedron as rigid bodies (for a total of 12 degrees of freedom). This led to R = 14% on the first 100 “|F obs|” extracted at low diffracting angles (<42°2θ). Refinements were undertaken by the Rietveld (Reference Rietveld1969) FULLPROF software (Rodriguez-Carvajal, Reference Rodriguez-Carvajal1993). The hydrogen atom positions were suggested by the SHELXL software (Sheldrick, Reference Sheldrick2008) in the hypothesis of cocrystallization. By using soft restraints on N-C and C-C interatomic distances, the 3,4-diaminopyridine molecule was maintained in the shape similar to that established by Rubin-Preminger and Englert (Reference Rubin-Preminger and Englert2007). The O and P atoms of the H3PO4 group were refined freely, but all H atoms were fixed at their positions suggested by the procedure implemented in the SHELXL program.

Refined atomic coordinates were further optimized by energy minimization in the solid state using the VASP program (Kresse and Hafner, Reference Kresse and Hafner1993; Kresse and Furthmüller, Reference Kresse and Furthmüller1996a, Reference Kresse and Furthmüller1996b). The electron exchange-correlation interaction was described in the generalized gradient approximation (Perdew and Wang, Reference Perdew and Wang1992). Plane waves formed the basis set, and calculations were performed using the projector-augmented wave method (Blöchl, Reference Blöchl1994; Kresse and Joubert, Reference Kresse and Joubert1999) and atomic pseudopotentials

TABLE I. Experimental and Rietveld refinement details for the 3,4-diaminopyridin-1-ium dihydrogen phosphate.

(Kresse and Hafner, Reference Kresse and Hafner1994). The energy cut-off controlling the accuracy of the calculation was set to 500 eV, representing an extended basis set and consequently highly accurate calculations. The positions of all atoms were optimized by means of the conjugated gradient method in the Γ-point of the Brillouin zone with the unit-cell parameters fixed. In the course of the calculation, one of the hydrogen atoms from the H3PO4 group jumped to the molecule ring N atom forming

Figure 1. (Color online) Rietveld plot for 3,4-diaminopyridin-1-ium dihydrogen phosphate, [C5H3(NH)(NH2)2]+ (H2PO4). Calculated diffraction pattern as solid line, observed as dots, difference curve in the lower part, and positions of reflections as vertical bars.

Figure 2. (Color online) Atom numbering scheme for the moieties of 3,4-diaminopyridin-1-ium dihydrogen phosphate.

pyridinium cation. This new conformation was then confirmed by the independent calculations starting with a different initial geometry. The final formula of the title compound hence is [C5H3(NH)(NH2)2]+ (H2PO4). The same behavior conducting to a 3,4-pyridinium salt is observed when the phosphoric acid is replaced by succinic (Fun and Balasubramini, Reference Fun and Balasubramani2009), squaric (Koleva et al., Reference Koleva, Tsanev, Kolev, Mayer-Figge and Sheldrick2007), or tartaric (Koleva et al., Reference Koleva, Kolev, Tsanev, Kotov, Mayer-Figge, Seidel and Sheldrick2008) acids.

The basic crystallographic data, experimental, and Rietveld-refinement details are summarized in Table I. The Rietveld fit, atoms numbering, and crystal packing are depicted in Figures 13, respectively. The atomic coordinates obtained by both methods are listed in Table II, and the essential bond distances, bond angles are detailed in Table III. Hydrogen bonds geometry is summarized in Table IV, and X-ray powder diffraction data are given in Table V. Both sets of atomic coordinates were deposited (CIF) at the Crystallography Open Database (www.crystallography.net) (Gražulis et al., Reference Gražulis, Chateigner, Downs, Yokochi, Quirós, Lutterotti, Manakova, Butkus, Moeck and Le Bail2009) and in the supplementary date.

DISCUSSION

According to the optimized atomic coordinates, the molecule of 3,4-DAP+ is essentially planar with the maximum

TABLE II. Fractional atomic coordinates and isotropic displacement parameters (Å2) for the 3,4-diaminopyridin-1-ium dihydrogen phosphate. The coordinates obtained by energy minimization in the solid state are listed in the bottom lines. For H atoms, the first line gives the coordinates as suggested by SHELXL in the cocrystal hypothesis.

Figure 3. (Color online) Crystal packing of the 3,4-diaminopyridin-1-ium dihydrogen phosphate viewed along the b-axis. Dashed lines indicate hydrogen bonds.

deviation from the least-squares plane fit to all nonhydrogen atoms being only 0.039 Å (due to N3). The optimized bond distances and bond angles are in good agreement with the values derived from the single crystal data (Rubin-Preminger and Englert, Reference Rubin-Preminger and Englert2007), though it is N2 which is significantly displaced by 0.131 Å from the plane in the 3,4-DAP. The phosphate group in 3,4-DAPP is a slightly deformed tetrahedra with a minimum spread of the O-P-O bond angles

TABLE III. Selected bond distances (Å) and angles (°) for the 3,4-diaminopyridin-1-ium dihydrogen phosphate obtained by the PLATON program (Spek, Reference Spek2003). The values derived from energy minimization in the solid state are listed in the bottom lines.

(Table III). As expected, the P-O bond distances fall into two groups, and the protonated P-O vertices show the expected lengthening of the P=O bonds (Blessing, Reference Blessing1988; Mahmoudkhani and Langer, Reference Mahmoudkhani and Langer2002; Demir et al., Reference Demir, Yilmaz and Harrison2003; Smrčok et al., Reference Smrčok, Havlíček, Kaman and Rundlőf2009; Balamurugan et al., Reference Balamurugan, Jagan and Sivakumar2010; Kaman et al., Reference Kaman, Smrčok, Císařová and Havlíček2011; Marouani et al., Reference Marouani, Al-Deyab and Rzaigui2011). The structure of the 3,4-DAPP consists of dihydrogen phosphate tetrahedra linked by moderately strong (Jeffrey, Reference Jeffrey1997) O-H…O hydrogen bonds [Table IV and Figure 4(a)] involving all oxygen

TABLE IV. Hydrogen bond geometry of the 3,4-diaminopyridin-1-ium dihydrogen phosphate from the energy optimized coordinates.

a 1/2 − x, 1/2 − y, 1/2 − z.

b 1/2 − x, −1/2 − y, 1/2 − z.

c x, −y, 1/2 + z.

dx, −y, 1 − z.

a −1/2 + x, 1/2 − y, z.

TABLE V. X-ray diffraction data of 3,4-diaminopyridin-1-ium dihydrogen phosphate (Cu Kα1).

atoms to the chains running along the b-axis. Antiparallelly, π–π stacked 3,4-DAP+ cations form molecular columns in the spaces between the chains. The molecules are mutually slipped, but the perpendicular distances of the molecules are only ∼3.31 Å. Although the dominant interaction of the molecules with their surroundings is electrostatic, their bonding is further enhanced by the formation of N-H…O and C-H…O hydrogen bonds [Table IV and Figure 4(b)]. As a result, the architecture of hydrogen bonds comprise oxygen atoms, which are either acceptors or both donors and acceptors. Specifically, while the O1 and O2 atoms play only the role of acceptors of relatively strong O-H…O and N-H…O hydrogen bonds, the O3, except for playing a donor role in the phosphate chain, also receives a pair of a strong N-H and a weak C-H bonds. Finally, the donor O4 atom is also acceptor of a weak C-H…O bond. The N3…O1(O2) separations are comparable to the estimated mean values (∼2.85 Å) of the distribution of N…O contact distances extracted from CCSD database (Kumara Swamy et al., Reference Kumara Swamy, Kumaraswamy and Kommana2001). On the contrary, the hydrogen bonds involving the

Figure 4. (Color online) Details of hydrogen bonding in the structure of the 3,4-diaminopyridin-1-ium dihydrogen phosphate: (a) O-H…O bonds in the chains of dihydrogen phosphate moieties; (b) hydrogen bonds formed by the molecules of 3,4-diaminopyridin-1-ium cation. Symmetry codes are given at the bottom of Table IV.

N2 atom are remarkably longer and hence weaker. The lengths of two C-H…O are comparable, but the bond angle in C5-H3…O4 is close to linear, and the bond is hence probably somewhat stronger.

CONCLUSION

In the absence of a suitable single crystal, modern powder diffraction methodologies have the power sometimes to reveal crystal structures, including those of important pharmaceutical compounds. However, peak overlapping leads to unfavorable data/parameter ratios, and the crystal structures are often quite inaccurate without using restraints and constraints and even when using them. Moreover, not locating clearly H atoms may lead to difficulties to give the correct names to molecular compounds (in the cocrystallization hypothesis, the name would have been 3,4-diaminopyridine phosphoric acid). The presence of a single phase (no polymorph) in the sample studied is thus confirmed. Optimizing the atomic coordinates by quantum chemical calculations looks to be a complementary approach that will probably be soon generalized in such cases.

ACKNOWLEDGEMENT

One of the authors (L.S.) wishes to express his thanks to Slovak Grant Agency VEGA for a partial financial support of this study under the contract No. 2/150/09.

References

Balamurugan, P., Jagan, R., and Sivakumar, K. (2010). “Dihydrogen phosphate mediated supramolecular frameworks in 2- and 4-chloroanilinium dihydrogen phosphate salts,” Acta Crystallogr., Sect. C: Cryst. Struct. Commun. 66, o109o113.10.1107/S0108270110001940/gd3320sup1.cifGoogle Scholar
Blessing, R. H. (1988). “New analysis of the neutron diffraction data for anhydrous orthophosphoric acid and the structure of H3PO4 molecules in crystals,” Acta Crystallogr. Sect. B: Struct. Sci. 44, 334340.10.1107/S0108768188001429Google Scholar
Blöchl, P. E. (1994). “Projector augmented-wave method.” Phys. Rev. B 50, 1795317979.10.1103/PhysRevB.50.17953Google Scholar
Demir, S., Yilmaz, V. T., and Harrison, W. T. A. (2003). “2-(Hydrogenmethyl) pyridinium dihydrogenphosphate,” Acta Crystallogr., Sect. C: Cryst. Struct. Commun. 59, o378o380.10.1107/S0108270103011077CrossRefGoogle Scholar
de Wolf, P. M. (1968). “A simplified criterion for the reliability of a powder pattern indexing,” J. Appl. Crystallogr. 1, 108113.10.1107/S002188986800508XCrossRefGoogle Scholar
Fun, H.-K. and Balasubramani, K. (2009). “3,4-Diaminopyridinium hydrogen succinate,” Acta Crystallogr., Sect. E: Struct. Rep. Online 65, o1531o1532.10.1107/S1600536809021205/sj2630sup1.cifGoogle Scholar
Gražulis, S., Chateigner, D., Downs, R. T., Yokochi, A. F. T., Quirós, M.,Lutterotti, L., Manakova, E., Butkus, J., Moeck, P., and Le Bail, A. (2009). “Crystallography Open Database - an open-access collection of crystal structures,” J. Appl. Crystallogr. 42, 726729.10.1107/S0021889809016690CrossRefGoogle ScholarPubMed
Guyon, F., Pradeau, D., Le Hoang, M. D., and Houri, J.-J.(2002). International Patent No. WO2002062760 (15 August).Google Scholar
Jeffrey, G. A. (1997). An Introduction to Hydrogen Bonding (Oxford University Press, New York), p. 12.Google Scholar
Kaman, O., Smrčok, L’., Císařová, I., and Havlíček, D. (2011). “Dihydrogen phosphate and hydrogen sulphate of 1,4-dimethyl-1,4-diazabicyclo[2.2.2]octane-1,4-diium: Crystal structures, hydrogen bonding and infrared spectra,” J. Chem. Crystallogr. 41, 1539154610.1007/s10870-011-0137-0CrossRefGoogle Scholar
Koleva, B., Tsanev, T., Kolev, T., Mayer-Figge, H., and Sheldrick, W. S. (2007). “3,4-diaminopyridinium hydrogen squarate,” Acta Crystallogr., Sect. E: Struct. Rep. Online 63, o3556.10.1107/S1600536807031170/bt2412sup1.cifGoogle Scholar
Koleva, B., Kolev, T., Tsanev, T., Kotov, S., Mayer-Figge, H., Seidel, R. W., and Sheldrick, W. S. (2008). “Spectroscopic and structural elucidation of 3,4-diaminopyridine and its hydrogentartarate salt: Crystal structure of 3,4-diaminopyridinium hydrogentartarate dihydrate,” J. Mol. Struct. 881, 146155.10.1016/j.molstruc.2007.09.006Google Scholar
Kresse, G. and Furthmüller, J. (1996a). “Efficient iterative scheme for ab initio total energy calculations using a plane-wave basis set,” Phys. Rev. B, 54, 1116911186.10.1103/PhysRevB.54.11169Google Scholar
Kresse, G. and Furthmüller, J. (1996b). “Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set,” Comput. Mater. Sci. 6, 1550.10.1016/0927-0256(96)00008-0CrossRefGoogle Scholar
Kresse, G. and Hafner, J. (1993). “Ab initio molecular dynamics for open-shell transition metals,” Phys. Rev. B 48, 1311513118.10.1103/PhysRevB.48.13115Google Scholar
Kresse, G. and Hafner, J. (1994). “Norm-conserving and ultrasoft pseudopotentials for first-row and transition elements,” J. Phys.: Condens. Matter 6, 82458527.10.1088/0953-8984/6/40/015Google Scholar
Kresse, G. and Joubert, J. (1999). “From ultrasoft potentials to the projector augmented wave method,” Phys. Rev. B 59, 17581775.10.1103/PhysRevB.59.1758CrossRefGoogle Scholar
Kumara Swamy, K. C., Kumaraswamy, S., and Kommana, P. (2001). “Very strong C-H…O, N-H…O, and O-H…O hydrogen bonds involving a cyclic phosphate,” J. Am. Chem. Soc. 123, 1264212643.10.1021/ja010713xGoogle Scholar
Le Bail, A. (2001). “ESPOIR: A program for solving structures by Monte Carlo from powder diffraction data,” Mater. Sci. Forum 378–381, 6570.10.4028/www.scientific.net/MSF.378-381.65Google Scholar
Le Bail, A. (2004). “Monte Carlo indexing with MCMAILLE,” Powder Diffr. 19, 249254.10.1154/1.1763152Google Scholar
Le Bail, A. (2005). “Whole powder pattern decomposition methods and applications—A retrospection,” Powder Diffr. 20, 316326.10.1154/1.2135315CrossRefGoogle Scholar
Le Bail, A. (2008). “Structure solution,” in Principles and Applications of Powder Diffraction, edited by Clearfield, A., Reibenspies, J., and Bhuvanesh, N. (Wiley, New York), pp. 261309.Google Scholar
Le Bail, A., Cranswick, L. M. D, Adil, K., Altomare, A., Avdeev, M., Cerny, R., Cuocci, C., Giacovazzo, C., Halasz, I., Lapidus, S. H., Louwen, J. N., Moliterni, A., Palatinus, L., Rizzi, R., Schilder, E. C., Stephens, P. W., Stone, K. H., and van Mechelen, J. (2009). “Third structure determination by powder diffractometry round robin (SDPDRR-3),” Powder Diffr. 24, 254262.10.1154/1.3200881Google Scholar
Mahmoudkhani, A. H. and Langer, V. (2002) “Phenylphosphonic acid as a building block for two-dimensional hydrogen-bonded supramolecular array,” J. Mol. Struct. 609, 97108.10.1016/S0022-2860(01)00954-1Google Scholar
Marouani, H., Al-Deyab, S. S., and Rzaigui, M. (2011). “2-Aminopyrimidinium dihydrogen phosphate monohydrate,” Acta Crystallogr., Sect. E: Struct. Rep. Online 67, o970o971.10.1107/S1600536811010658/pv2399sup1.cifGoogle Scholar
McEvoy, K. M., Windebank, A. J., Daube, J. R., and Low, P. A. (1989). “2,4-Diaminopyridine in the treatment of Lambert-Eaton mysathenic syndrome,” N. Engl. J. Med. 321, 15671571.10.1056/NEJM198912073212303CrossRefGoogle Scholar
Perdew, J. P. and Wang, Y. (1992). “Accurate and simple analytic representation of the electron-gas correlation energy,” Phys. Rev. B 45, 1324413249.10.1103/PhysRevB.45.13244Google Scholar
Quartel, A., Turbeville, S., and Lounsbury, D. (2010). “Current therapy for Lambert-Eaton myasthenic syndrome: Development of 3,4-diaminopyridine phosphate salt as first-line symptomatic treatment,” Curr. Med. Res. Opin. 26, 13631375.10.1185/03007991003745209Google Scholar
Rietveld, H. M. (1969). “A profile refinement method for nuclear and magnetic structures,” J. Appl. Crystallogr. 2, 6571.10.1107/S0021889869006558Google Scholar
Rodriguez-Carvajal, J. (1993). “Recent advances in magnetic-structure determination by neutron powder diffraction,” Physica B, 192, 5569.10.1016/0921-4526(93)90108-IGoogle Scholar
Rubin-Preminger, J. M. and Englert, U. (2007). “3,4-Diaminopyridine,” Acta Crystallogr., Sect. E: Struct. Rep. Online 64, o757o758.10.1107/S1600536807001444/bt2244sup1.cifGoogle Scholar
Sheldrick, G. (2008). “A short history of SHELX,” Acta Crystallogr., Sect. A: Found. Crystallogr. 64, 112122.10.1107/S0108767307043930CrossRefGoogle Scholar
Smith, G. S. and Snyder, R. L. (1979). “F N: A criterion for rating powder diffraction patterns and evaluating the reliability of powder-pattern indexing,” J. Appl. Crystallogr. 12, 6065.10.1107/S002188987901178XCrossRefGoogle Scholar
Smrčok, Ľ., Havlíček, D., Kaman, O., and Rundlőf, H. (2009). “1,4-Diazabicyclo[2.2.2]octane-1,4-diium dihydrogenphosphate monohydrate from X-ray and neutron data,” Z. Kristallogr. 224, 174178.10.1524/zkri.2009.1127CrossRefGoogle Scholar
Spek, A. L. (2003). “Single crystal structure validation with the program PLATON,” J. Appl. Crystallogr. 36, 713.10.1107/S0021889802022112Google Scholar
See supplementary material at http://dx.doi.org/10.1154/1.3660160 E-PODIE2-26-011104 for the crystallographic information file (CIF).Google Scholar
Figure 0

TABLE I. Experimental and Rietveld refinement details for the 3,4-diaminopyridin-1-ium dihydrogen phosphate.

Figure 1

Figure 1. (Color online) Rietveld plot for 3,4-diaminopyridin-1-ium dihydrogen phosphate, [C5H3(NH)(NH2)2]+ (H2PO4). Calculated diffraction pattern as solid line, observed as dots, difference curve in the lower part, and positions of reflections as vertical bars.

Figure 2

Figure 2. (Color online) Atom numbering scheme for the moieties of 3,4-diaminopyridin-1-ium dihydrogen phosphate.

Figure 3

TABLE II. Fractional atomic coordinates and isotropic displacement parameters (Å2) for the 3,4-diaminopyridin-1-ium dihydrogen phosphate. The coordinates obtained by energy minimization in the solid state are listed in the bottom lines. For H atoms, the first line gives the coordinates as suggested by SHELXL in the cocrystal hypothesis.

Figure 4

Figure 3. (Color online) Crystal packing of the 3,4-diaminopyridin-1-ium dihydrogen phosphate viewed along the b-axis. Dashed lines indicate hydrogen bonds.

Figure 5

TABLE III. Selected bond distances (Å) and angles (°) for the 3,4-diaminopyridin-1-ium dihydrogen phosphate obtained by the PLATON program (Spek, 2003). The values derived from energy minimization in the solid state are listed in the bottom lines.

Figure 6

TABLE IV. Hydrogen bond geometry of the 3,4-diaminopyridin-1-ium dihydrogen phosphate from the energy optimized coordinates.

Figure 7

TABLE V. X-ray diffraction data of 3,4-diaminopyridin-1-ium dihydrogen phosphate (Cu Kα1).

Figure 8

Figure 4. (Color online) Details of hydrogen bonding in the structure of the 3,4-diaminopyridin-1-ium dihydrogen phosphate: (a) O-H…O bonds in the chains of dihydrogen phosphate moieties; (b) hydrogen bonds formed by the molecules of 3,4-diaminopyridin-1-ium cation. Symmetry codes are given at the bottom of Table IV.