Adsorption is a fundamental step in catalysis on a solid surface, and therefore its understanding is important to explaining its behavior. This work investigated the adsorption of various small molecules, including H2, N2, CO, CO2, CH4, NH3, H2O, H2S, DMSO, alkanes, alkenes, alkynes, aromatic compounds, alcohols, aldehydes, ketones, nitriles, carboxylic acids, amides, and amines, on the anatase (101) and rutile (110) surfaces of TiO2, using periodic density functional theory calculations and statistical methods. Adsorption energies were computed at the same level of theory to obtain a clean and consistent data set. A linear relationship was observed between the adsorption energies of these molecules and their highest occupied molecular orbital (HOMO) levels, whereas no obvious correlation was evident for the lowest unoccupied molecular orbital (LUMO) levels. Improved correlations between the adsorption energies and the HOMO levels were generated by dividing these molecules into two subgroups: hydrocarbons and heteroatom-containing compounds. Interactions between frontier molecular orbitals and the surfaces were considered, to gain a better understanding of the significant correlations that were identified. The data show that these relationships can be primarily ascribed to interactions between the HOMO of the small molecule and conduction state of the TiO2 surface. Statistical analysis using machine learning demonstrated that the HOMO and the dipole moment are the first and second most important properties, respectively, in terms of rationalizing and predicting the adsorption energies.