Uncovering molecular details of urea crystal growth in the presence of additives.

Controlling the shape of crystals is of great practical relevance in fields like pharmacology and fine chemistry. Here we examine the paradigmatic case of urea which is known to crystallize from water with a needle-like morphology. To prevent this undesired effect, inhibitors that selectively favor or discourage the growth of specific crystal faces can be used. In urea the most relevant faces are the {001} and the {110} which are known to grow fast and slow, respectively. The relevant growth speed difference between these two crystal faces is responsible for the needle-like structure of crystals grown in water solution. To prevent this effect, additives are used to slow down the growth of one face relative to another, thus controlling the shape of the crystal. We study the growth of fast {001} and slow {110} faces in water solution and the effect of shape controlling inhibitors like biuret. Extensive sampling through molecular dynamics simulations provides a microscopic picture of the growth mechanism and of the role of the additives. We find a continuous growth mechanism on the {001} face, while the slow growing {110} face evolves through a birth and spread process, in which the rate-determining step is the formation on the surface of a two-dimensional crystalline nucleus. On the {001} face, growth inhibitors like biuret compete with urea for the adsorption on surface lattice sites; on the {110} face instead additives cannot interact specifically with surface sites and play a marginal sterical hindrance of the crystal growth. The free energies of adsorption of additives and urea are evaluated with advanced simulation methods (well-tempered metadynamics) allowing a microscopic understanding of the selective effect of additives. Based on this case study, general principles for the understanding of the anisotropic growth of molecular crystals from solutions are laid out. Our work is a step toward a rational development of novel shape-affecting additives.

[1]  K. Roberts,et al.  Application of Bravais-Friedel-Donnay-Harker, attachment energy and Ising models to predicting and understanding the morphology of molecular crystals , 1991 .

[2]  Alessandro Laio,et al.  “Site Binding” of Ca2+ Ions to Polyacrylates in Water: A Molecular Dynamics Study of Coiling and Aggregation , 2007 .

[3]  Michele Parrinello,et al.  Freezing of a Lennard-Jones fluid: from nucleation to spinodal regime. , 2006, Physical review letters.

[4]  B. Roux,et al.  Energetics of ion conduction through the gramicidin channel , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[5]  L. Addadi,et al.  Molecular Recognition at Crystal Interfaces , 1991, Science.

[6]  Massimiliano Bonomi,et al.  PLUMED: A portable plugin for free-energy calculations with molecular dynamics , 2009, Comput. Phys. Commun..

[7]  S. Kuyucak,et al.  Response to “Comment on ‘Free energy simulations of single and double ion occupancy in gramicidin A’ ” [J. Chem. Phys.128, 227101 (2008)] , 2008 .

[8]  William L. Jorgensen,et al.  Do denaturants interact with aromatic hydrocarbons in water , 1993 .

[9]  B. Roux,et al.  Comment on "Free energy simulations of single and double ion occupancy in gramicidin A" [J. Chem. Phys. 126, 105103 (2007)]. , 2008, The Journal of chemical physics.

[10]  J. Kirkwood.  Statistical Mechanics of Fluid Mixtures , 1935 .

[11]  Francis S. Romanski,et al.  Atomistic simulation study of surfactant and polymer interactions on the surface of a fenofibrate crystal. , 2011, European journal of pharmaceutical sciences : official journal of the European Federation for Pharmaceutical Sciences.

[12]  S. Piana,et al.  Simulating micrometre-scale crystal growth from solution , 2005, Nature.

[13]  Marco Mazzotti,et al.  Monitoring the particle size and shape in the crystallization of paracetamol from water , 2010 .

[14]  Marco Mazzotti,et al.  Slowing the Growth Rate of Ibuprofen Crystals Using the Polymeric Additive Pluronic F127 , 2011 .

[15]  M. Parrinello,et al.  Polymorphic transitions in single crystals: A new molecular dynamics method , 1981 .

[16]  Haiping Fang,et al.  Urea-induced drying of hydrophobic nanotubes: comparison of different urea models. , 2011, The journal of physical chemistry. B.

[17]  Michele Parrinello,et al.  Effect of urea on the β-hairpin conformational ensemble and protein denaturation mechanism. , 2011, Journal of the American Chemical Society.

[18]  S. C. Parker,et al.  Growth modification of seeded calcite using carboxylic acids: atomistic simulations. , 2010, Journal of colloid and interface science.

[19]  P. Kollman,et al.  A well-behaved electrostatic potential-based method using charge restraints for deriving atomic char , 1993 .

[20]  Kwok Chow,et al.  Engineering of pharmaceutical materials: an industrial perspective. , 2008, Journal of pharmaceutical sciences.

[21]  Fujio Izumi,et al.  VESTA: a three-dimensional visualization system for electronic and structural analysis , 2008 .

[22]  C. S. Strom,et al.  On the Prediction of Crystal Morphology. I. The Hartman-Perdok Theory Revisited , 1998 .

[23]  M. Parrinello,et al.  Binding of calcium and carbonate to polyacrylates. , 2009, Journal of Physical Chemistry B.

[24]  W. Yuan,et al.  Crystallization of zinc lactate in presence of malic acid , 2010 .

[25]  E. Vlieg,et al.  Towards an atomic-scale understanding of crystal growth in solution. , 2007, Faraday discussions.

[26]  M. Parrinello,et al.  Well-tempered metadynamics: a smoothly converging and tunable free-energy method. , 2008, Physical review letters.

[27]  François Puel,et al.  Industrial batch crystallization of a plate-like organic product. In situ monitoring and 2D-CSD modelling: Part 1: Experimental study , 2007 .

[28]  In-Line Analysis of Impurity Effects on Crystallisation , 2005 .

[29]  Erik E. Santiso,et al.  A general set of order parameters for molecular crystals. , 2011, The Journal of chemical physics.

[30]  Carsten Kutzner,et al.  GROMACS 4:  Algorithms for Highly Efficient, Load-Balanced, and Scalable Molecular Simulation. , 2008, Journal of chemical theory and computation.

[31]  M. Parrinello,et al.  A molecular dynamics study of the early stages of calcium carbonate growth. , 2009, The journal of physical chemistry. B.

[32]  R. Tan,et al.  Molecular simulation study of the effect of various additives on salbutamol sulfate crystal habit. , 2011, Molecular pharmaceutics.

[33]  S. Piana,et al.  Understanding the barriers to crystal growth: dynamical simulation of the dissolution and growth of urea from aqueous solution. , 2005, Journal of the American Chemical Society.

[34]  J. Anwar,et al.  Mode of action and design rules for additives that modulate crystal nucleation. , 2009, Angewandte Chemie.

[35]  Linear models for prediction of ibuprofen crystal morphology based on hydrogen bonding propensities , 2009 .

[36]  K. Roberts,et al.  Molecular and Solid-State Modeling of the Crystal Purity and Morphology of ε-Caprolactam in the Presence of Synthesis Impurities and the Imino-Tautomeric Species Caprolactim , 2003 .

[37]  J. Anwar,et al.  Uncovering molecular processes in crystal nucleation and growth by using molecular simulation. , 2011, Angewandte Chemie.

[38]  R. K. McMullan,et al.  The crystal structure and molecular thermal motion of urea at 12, 60 and 123 K from neutron diffraction , 1984 .

[39]  Lorna J. Smith,et al.  Computer Simulation of Urea−Water Mixtures: A Test of Force Field Parameters for Use in Biomolecular Simulation , 2004 .

[40]  M. Parrinello,et al.  Isothermal-isobaric molecular dynamics using stochastic velocity rescaling. , 2009, The Journal of chemical physics.

[41]  Stéphane Veesler,et al.  The Evolution of Crystal Shape During Dissolution: Predictions and Experiments , 2008 .

[42]  J. Ulrich,et al.  Predicting the effect of impurity adsorption on crystal morphology , 2007 .

[43]  Narayan Variankaval,et al.  From Form to Function: Crystallization of Active Pharmaceutical Ingredients , 2008 .

[44]  A. Laio,et al.  Escaping free-energy minima , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[45]  Jacopo Tomasi,et al.  Evaluation of Solvent Effects in Isotropic and Anisotropic Dielectrics and in Ionic Solutions with a Unified Integral Equation Method: Theoretical Bases, Computational Implementation, and Numerical Applications , 1997 .

[46]  M. Parrinello,et al.  Asprich peptides are occluded in calcite and permanently disorder biomineral crystals. , 2010, Journal of the American Chemical Society.

[47]  P. Bennema,et al.  The attachment energy as a habit controlling factor: I. Theoretical considerations , 1980 .

[48]  Junmei Wang,et al.  Development and testing of a general amber force field , 2004, J. Comput. Chem..

[49]  Åke C. Rasmuson,et al.  Influence of Additives on Nucleation of Vanillin: Experiments and Introductory Molecular Simulations , 2004 .

[50]  P. Kollman,et al.  A Second Generation Force Field for the Simulation of Proteins, Nucleic Acids, and Organic Molecules , 1995 .