Microplastics (MPs) are ubiquitously distributed in freshwater systems and they can determine the environmental fate of organic pollutants (OPs) via sorption interaction. However, the diverse physicochemical properties of MPs and the wide range of OP species make a deeper understanding of sorption mechanisms challenging. Traditional isotherm-based sorption models are limited in their universality since they normally only consider the nature and characteristics of either sorbents or sorbates individually. Therefore, only specific equilibrium concentrations or specific sorption isotherms can be used to predict sorption. To systematically evaluate and predict OP sorption under the influence of both MPs and OPs properties, we collected 475 sorption data from peer-reviewed publications and developed a poly-parameter-linear-free-energy-relationship-embedded machine learning method to analyze the collected sorption datasets. Models of different algorithms were compared, and the genetic algorithm and support vector machine hybrid model displayed the best prediction performance (R of 0.93 and root-mean-square-error of 0.07). Finally, comparison results of three feature importance analysis tools (forward step wise method, Shapley method, and global sensitivity analysis) showed that chemical properties of MPs, excess molar refraction, and hydrogen-bonding interaction of OPs contribute the most to sorption, reflecting the dominant sorption mechanisms of hydrophobic partitioning, hydrogen bond formation, and π-π interaction, respectively. This study presents a novel sorbate-sorbent-based ML model with a wide applicability to expand our capacity in understanding the complicated process and mechanism of OP sorption on MPs.