(74j) On the Interplay between Conformational Complexity, Solution Structure, and Polymorphism in Succinic Acid Nucleation from Solution.
Obtaining an insight into the interplay between conformational transitions, the structure of the liquid parent phase, and the nucleation process leading to the appearance of a specific polymorph, is thus key in order to improve our overall understanding of nucleation of complex materials from a liquid solution.
In this work we investigate the relation between conformational transitions of succinic acid (1,4-Butanedioic acid), its solution structure, and the nucleation of two conformational polymorphs. The first polymorph is the notorious Î² form, stable at ambient temperature, in which succinic acid is in a planar conformation. The second is the recently discovered Î³ polymorph, in which succinic acid appears in a non-planar conformation. The Î³ polymorph has been serendipitously discovered during a co-crystallization experiment with Leu-Leu dipeptide, and was produced in concomitance with the Î² form. Despite CSP studies have shown that the Î³ form is indeed stable and thermodynamically plausible, extensive attempts at reproducing Î³ were not successful, classifying this system as a case in which disappearing polymorphism is observed. It is not uncommon in literature to encounter this phenomenon for molecules with conformational freedom.
We begin our work by analysing the stability at finite temperature and pressure of the bulk Î² and Î³ phases. As the stability of polymorphs changes at finite-temperature and pressure, molecular dynamics (MD) simulations are performed on Î² and Î³ crystals at 300K. These simulations show that while Î² is stable and the results are in reasonable agreement with experiments and CSP predictions, Î³ displays an incipient deformation of the cell shape. To further investigate such a structural relaxation, we perform metadynamics (MetaD) and well-tempered metadynamics (WTMetaD) simulations on Î³. As a result, we observe that the fluctuation of an angle of the supercell induces firstly the formation of local defects accompanied by conformational transitions of individual succinic acid molecules, then local melting, which ultimately leads to the irreversible melting of the entire cell.
We move then to analyse the conformational behaviour of succinic acid in solution. We initially carry out a metadynamics simulation of a single molecule in aqueous solution, biasing three characteristic dihedral angles. From the resulting free energy surface, we identify 10 distinct conformers, including the folded isomer found in Î³ and the planar one found in Î². A Markov State Model (MSM), informed by 18 100ns-long unbiased simulations allowed us to evaluate the equilibrium distribution of the conformers as well as the relaxation time of their population. The equilibrium distribution obtained from the MSM is in agreement with that obtained from metadynamics. Most notably we observe that Î³ conformer is dominant in solution (~70% probability), while the Î² planar conformer has a probability lower than 10%. The relaxation time for the conformer population is however fast, with a characteristic time of the order of 200 ps, hinting that in this case the conformational rearrangement does not represent a rate determining step in the nucleation process.
By repeating this analysis for increasing concentrations of succinic acid in water, we note that both the equilibrium probability and the relaxation time are affected by the transient formation of spatial domains locally dense in solute molecules.
To further investigate the interplay between conformational transitions and the nucleation process we study the homogeneous nucleation of both Î² and Î³ succinic acid crystals from solution through WTMetaD simulations. We expect such simulations to lead predominantly to the formation of Î² crystals, as it is well-known that this polymorph forms spontaneously even at low supersaturations; however, uncovering the mechanism of nucleation can shed a light upon possible role of Î³ as a metastable intermediate according to Ostwaldâs rule of stages.
In conclusion, we present a simulation work motivated by the discovery of a new elusive polymorph of succinic acid, Î³. This disappearing polymorph is compatible with the CSP potential energy landscape, which identify it as plausible; MD simulations point towards the interpretation of Î³ as a labile structure, possibly prone to melting as seen with metadynamics simulations. In addition, the conformer found in Î³ is the most abundant in solution. Nucleation studies will possibly uncover the mechanism followed and thus lead to a better understanding of the role of Î³ in this scenario.
 A. J. Cruz-Cabeza and J. Bernstein, âConformational Polymorphism,â Chem. Rev., vol. 114, no. 4, pp. 2170â2191, Feb. 2014.
 S. L. Price, âWhy donât we find more polymorphs?,â Acta Crystallogr. Sect. B Struct. Sci. Cryst. Eng. Mater., vol. 69, no. 4, pp. 313â328, Aug. 2013.
 J. D. Dunitz and J. Bernstein, âDisappearing Polymorphs,â Acc. Chem. Res., vol. 28, no. 4, pp. 193â200, Apr. 1995.
 A. Laio and M. Parrinello, âEscaping free-energy minima,â Proc. Natl. Acad. Sci., vol. 99, no. 20, pp. 12562â12566, Oct. 2002.
 A. Barducci, G. Bussi, and M. Parrinello, âWell-tempered metadynamics: A smoothly converging and tunable free-energy method,â Phys. Rev. Lett., vol. 100, no. 2, pp. 1â4, 2008.
 P. Lucaioli, E. Nauha, I. Gimondi, L. S. Price, R. Guo, L. Iuzzolino, I. Singh, M. Salvalaglio, S. L. Price, N. Blagden, submitted, 2018