Organic molecules can crystallize in multiple structures or polymorphs, yielding crystals with very different physical and mechanical properties. The prediction of the polymorphs that may appear in nature is a challenge with great potential benefits for the development of new products and processes. A multistage crystal structure prediction (CSP) methodology is applied to axitinib, a pharmaceutical molecule with significant polymorphism arising from molecular flexibility. The CSP study is focused on those polymorphs with one molecule in the asymmetric unit. The approach successfully identifies all four known polymorphs within this class, as well as a large number of other low-energy structures. The important role of conformational flexibility is highlighted. The performance of the approach is discussed in terms of both the quality of the results and various algorithmic and computational aspects, and some key priorities for further work in this area are identified. (C) 2014 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY license (hap://creativecommons.org/licenses/by/3.0/).