Thalidomide is a chiral molecule, which was sold in the 1960s as a sedative in its (S,R)-racemic form. The tragedy was that the (S)-isomer was tetragenic, and only the (R) enantiomer acts as a sedative. What was not appreciated at the time is that interconversion of the (S)- and (R) forms takes place quite quickly in aqueous media. Nowadays, quantum modelling can provide good in-silico estimates of the (free) energy barriers for such processes, which in this case is a simple keto-enol tautomerism. In a recently published article[1], just such a simulation is reported. By involving two explicit water molecules in the transition state, an (~enthalpic) barrier of 27.7 kcal/mol was obtained. The simulation was conducted just with two water molecules acting as solvent, and without any additional continuum solvation applied. So I thought I would re-evaluate this result by computing it at the ωB97XD/6-311G(d,p)/SCRF=water level (a triple-ζ basis set rather than the double-ζ used before[1]), and employing a dispersion-corrected DFT method rather than B3LYP.
References
- C. Tian, P. Xiu, Y. Meng, W. Zhao, Z. Wang, and R. Zhou, "Enantiomerization Mechanism of Thalidomide and the Role of Water and Hydroxide Ions", Chemistry - A European Journal, vol. 18, pp. 14305-14313, 2012. http://dx.doi.org/10.1002/chem.201202651