Replacement of Oxygen by Sulfur in Small Organic Molecules. 3. Theoretical Studies on the Tautomeric Equilibria of the 2OH and 4OH-Substituted Oxazole and Thiazole and the 3OH and 4OH-Substituted Isoxazole and Isothiazole in the Isolated State and in Solution

This follow-up paper completes the author’s investigations to explore the in-solution structural preferences and relative free energies of all OH-substituted oxazole, thiazole, isoxazole, and isothiazole systems. The polarizable continuum dielectric solvent method calculations in the integral-equation formalism (IEF-PCM) were performed at the DFT/B97D/aug-cc-pv(q+(d))z level for the stable neutral tautomers with geometries optimized in dichloromethane and aqueous solution. With the exception of the predictions for the predominant tautomers of the 3OH isoxazole and isothiazole, the results of the IEF-PCM calculations for identifying the most stable tautomer of the given species in the two selected solvents agreed with those from experimental investigations. The calculations predict that the hydroxy proton, with the exception for the 4OH isoxazole and 4OH isothiazole, moves preferentially to the ring nitrogen or to a ring carbon atom in parallel with the development of a C=O group. The remaining, low-fraction OH tautomers will not be observable in the equilibrium compositions. Relative solvation free energies obtained by the free energy perturbation method implemented in Monte Carlo simulations are in moderate accord with the IEF-PCM results, but consideration of the ΔG(solv)/MC values in calculating ΔG(s)(tot) maintains the tautomeric preferences. It was revealed from the Monte Carlo solution structure analyses that the S atom is not a hydrogen-bond acceptor in any OH-substituted thiazole or isothiazole, and the OH-substituted isoxazole and oxazole ring oxygens may act as a weak hydrogen-bond acceptor at most. The molecules form 1.0−3.4 solute−water hydrogen bonds in generally unexplored numbers at some specific solute sites. Nonetheless, hydrogen-bond formation is favorable with the NH, C=O and OH groups.