379700 Predicting the Neat Polymorphs of Axitinib
In recent years, there has been significant progress in crystal structure prediction methodologies, as demonstrated by recent blind tests and other work on the prediction of the crystal structure of small pharmaceutical molecules1. In this paper, we discuss the application of one particular approach that has proved successful to the challenging case of axitinib. Axitinib is well-studied molecule developed by Pfizer, which has been shown to have 71 polymorphs2 of which 5 are neat (i.e., they contain only the axitinib molecule). The focus of the present study on neat polymorphs with one molecule in the asymmetric unit, of which 4 are known experimentally (Forms I, VI, XXV and XLI).
The approach to crystal structure prediction adopted here is based on the lattice energy minimization of a large number of trial structures at 0 K in order to generate an energy landscape, in which low-energy minima are identified based on the computational model of choice. Success is dependent upon the accuracy of the lattice energy calculations, especially in terms of electrostatics and the impact of molecular flexibility, and the ability to optimise a hundreds of thousands to millions of putative structures. Such approaches have been proven to be successful in predicting the polymorphs of compounds with one or two polymorphs in many cases3.
We use a multi-stage methodology. A model is first chosen based on analysis of the conformational flexibility of the molecule. Then 4,800,000 structures are generated and locally minimised with a computationally inexpensive model, using the CrystalPredictor4,5 algorithm. Then the most promising structures (lowest in energy) are further minimised in a two-step approach with CrystalOptimizer6,7, a program performing local minimisation, adopting a highly accurate though computationally expensive model. The axitinib molecule is found to have a richly populated energy landscape. The four known structures are predicted as lattice energy minima with the chosen model. Although Form XLI is not found to be the most stable at 0 K, contrary to expectations based on data at higher temperatures, the relative stability of other forms is predicted correctly. The relative importance of the different contributions to the energy model (intramolecular energy, electrostatics, repulsion/dispersion) is analysed and provides insights into the need for future development of predictive methodologies.
Bibliography
1. C. C. Pantelides, C. S. Adjiman, A. V. Kazantsev. “General computational algorithms for ab initio crystal structure prediction for organic molecules”. In S. Atahan-Evrenk and A. Aspuru-Guzik, editors, Prediction and Calculation of Crystal Structures, volume 345 of Topics in Current Chemistry, pages 25-58. Springer International Publishing, 2014.
2. A. M. Campeta, B. P. Chekal, Y. A. Abramov, P. A. Meenan, M. J. Henson, B. Shi, R. A. Singer, K. R. Horspool, Development of a targeted polymorph screening approach for a complex polymorphic and highly solvating API, Journal of Pharmaceutical Sciences, 99(9):3874-3886, 2010.
3. D. A. Bardwell, C. S. Adjiman, Y. A. Arnautova, E. Bartashevich, S. X. M. Boerrigter, D. E. Braun, A. J. Cruz-Cabeza, G. M. Day, R. G. Della Valle, G. R. Desiraju, B. P. van Eijck, J. C. Facelli, M. B. Ferraro, D. Grillo, M. Habgood, D. W. M. Hofmann, F. Hofmann, K. V. Jovan Jose, P. G. Karamertzanis, A. V. Kazantsev, J. Kendrick, L. N. Kuleshova, F. J. J. Leusen, A. V. Maleev, A. J. Misquitta, S. Mohamed, R. J. Needs, M. A. Neumann, D. Nikylov, A. M.Orendt, R. Pal, C. C. Pantelides, C. J. Pickard, L. S.Price, S. L. Price, H. A. Scheraga, J. van de Streek, T. S.Thakur, S. Tiwari, E. Venuti, and I. K. Zhitkov. Towards crystal structure prediction of complex organic compounds -a report on thefifth blind test. Acta Crystallographica Section B, 67(6):535-551, 2011.
4. P.G. Karamertzanis, C.C. Pantelides, “Ab initio-Crystal Structure Prediction – I. Rigid molecules”, J. Comput. Chem., 26:304-324 (2005)
5. P.G. Karamertzanis, C.C. Pantelides, “Ab initio-Crystal Structure Prediction – II. Flexible molecules”, Molecular Physics, Vol. 105 , Nos. 2-3, 273-291 (2007)
6. A.V. Kazantsev, P.G. Karamertzanis, C. S. Adjiman, C C. Pantelides, J. Chem. Theory Comput., Efficient Handling of Molecular Flexibility in Lattice Energy Minimization of Organic Crystals, 7(6):1998-2016, 2011.
See more of this Group/Topical: Separations Division