Euler-Lagrange simulations of gas-solid flows in unbounded domains have been performed to study sub-grid modeling of the filtered drag force for non-cohesive and cohesive particles. The filtered drag forces under various microstructures and flow conditions were analyzed in terms of various sub-grid quantities: the sub-grid drift velocity, which stems from the sub-grid correlation between the local fluid velocity and the local particle volume fraction, and the scalar variance of solid volume fraction, which is a measure to identify the degree of local inhomogeneity of volume fraction within a filter volume. The results show that the drift velocity and the scalar variance exert systematic effects on the filtered drag force. Effects of particle and domain sizes, gravitational accelerations, and mass loadings on the filtered drag are also studied, and it is shown that these effects can be captured by both sub-grid quantities. Additionally, the effect of cohesion force through the van der Waals interaction on the filtered drag force is investigated, and it is found that there is no significant difference on the dependence of the filtered drag coefficient of cohesive and non-cohesive particles on the sub-grid drift velocity or the scalar variance of solid volume fraction. The assessment of predictabilities of sub-grid quantities was performed by correlation coefficient analyses in a priori manner, and it is found that the drift velocity is superior. However, the drift velocity is not available in "coarse-grid" simulations and a specific closure is needed. A dynamic scale-similarity approach was used to model drift velocity but the predictability of that model is not entirely satisfactory. It is concluded that one must develop a more elaborate model for estimating the drift velocity in "coarse-grid" simulations.