The retention of a solute in RP chromatography is a very complex process which depends on many factors. Therefore, the study of the influence of a mobile phase modifier concentration on the retention in different reversed phase chromatographic systems is very important for understanding the rules governing retention and mechanisms of substance separation in a chromatographic process. Composition changes and the nature of mobile phases enable tuning of the separated analytes' retention over a wide range of retention parameters and optimization of the chromatographic process as well. Optimization of the chromatographic process can be achieved by several different methods; one of them is the so-called interpretative strategy. The key approach adopted in this strategy is the implementation of adequate retention models that couple the retention of solute with the composition of a mixed mobile phase. The use of chemically bonded stationary phases composed of partially non-bonded silica matrix and organic ligands bonded to its surface in everyday chromatography practice leads to questions of the correct definition of the retention model and the dominant retention mechanism in such chromatographic systems. The retention model for an accurate prediction of retention factor as a function of modifier concentration and the heterogeneity of the adsorbent surface should be taken into consideration. In this work the influence of mobile-phase composition on the retention of sixteen model substances such as phenols, quinolines, and anilines used as test analytes in different RP-TLC systems with CN-, NH2-, and Diol-silica polar bonded stationary phases has been studied. The aim of this study is to compare the performance of three valuable retention models assumed as the partition, adsorption/partition, and adsorption mechanism of retention. All the models were verified for different RP-TLC systems by three statistical criteria. The results of investigations presented in this work demonstrate that the best agreement between the experimental and calculated Rf values was obtained by the use of new-generation retention models, which assume heterogeneity of adsorbent surface. The results reported here show that heterogeneity of the adsorbent surface may be important in analysis of the elution process in liquid chromatography. Consideration of the goodness of fit for the experimental data to the examined retention models is in conformity with the adsorption mechanism of retention on all polar bonded stationary phases in most eluent systems for most investigated compounds.