Here we present the results of our unbiased searches of glycine polymorphs obtained using the genetic algorithms search implemented in MGAC, modified genetic algorithm for crystals, coupled with the local optimization and energy evaluation provided by Quantum Espresso. We demonstrate that it is possible to predict the crystal structures of a biomedical molecule using solely first principles calculations. We were able to find all the ambient pressure stable glycine polymorphs, which are found in the same energetic ordering as observed experimentally and the agreement between the experimental and predicted structures is of such accuracy that the two are visually almost indistinguishable. (C) 2015 Elsevier B.V. All rights reserved.