I will provide a brief overview on my work concerning constraints on modified gravity models obtained using galaxy cluster mass profile determinations. In particular, I will present the results of a paper in which we combined the information given by the kinematic of galaxies in cluster with the information provided by lensing analyses for 2 galaxy clusters of the CLASH-CLASH\VLT collaboration to get constraints on f(R) models. In order to discuss the applicability of the proposed method in view of future imaging and spectroscopic surveys, I will further introduce my current study of cosmological simulations, aiming at estimating and calibrating the impact of systematics.