Polymer composite materials require softening to reduce their glass transition temperature and improve processability. To this end, plasticizers (PLs), which are small organic molecules, are added to the polymer matrix. The miscibility of these PLs has a large impact on their effectiveness and, therefore, their interactions with the polymer matrix must be carefully considered. Many PL characteristics, including their size, topology, and flexibility, can impact their miscibility and, because of the exponentially large number of PLs, the current trial-and-error approach is very ineffective. In this work, we show that using coarse-grained molecular simulations of a small dataset of 48 PLs, it is possible to identify topological and thermodynamic descriptors that are proxy for their miscibility. Using ad-hoc molecular dynamics simulation setups that are relatively computationally inexpensive, we establish correlations between the PLs' topology, internal flexibility, thermodynamics of aggregation, and degree of miscibility, and use these descriptors to classify the molecules as miscible or immiscible. With all available data, we also construct a decision tree model, which achieves a F1 score of 0.86 ± 0.01 with repeated, stratified 5-fold cross-validation, indicating that this machine learning method can be a promising route to fully automate the screening. By evaluating the individual performance of the descriptors, we show this procedure enables a 10-fold reduction of the test space and provides the basis for the development of workflows that can efficiently screen PLs with a variety of topological features. The approach is used here to screen for apolar PLs in polyisoprene melts, but similar proxies would be valid for other polyolefins, while, in cases where polar interactions drive the miscibility, other descriptors are likely to be needed.