There is demand for safer and bio-based solvents, brought on by legislation and sustainability objectives. The prediction of physical properties is highly desirable to help design new molecules. Here we present an in silico approach to obtain calculated Kamlet–Abboud–Taft solvatochromic parameters using virtual experiments. The tautomerisation equilibrium of methyl acetoacetate and dimedone was calculated in different solvents with COSMO-RS theory and converted into estimates of solvent dipolarity and hydrogen bond accepting ability, respectively. Hydrogen bond donating ability was calculated as a function of the electron deficient surface area on protic solvents. These polarity descriptors correlate with rate constants and equilibria, and so ability of calculated Kamlet–Abboud–Taft solvatochromic parameters to recreate experimental free energy relationships was tested with sixteen case studies taken from the literature. The accuracy of the calculated parameters was also satisfactory for solvent selection, as demonstrated with a 1,4-addition reaction and a multicomponent heterocycle synthesis.