I present the recent results obtained using the updated version of MG-MAMPOSSt, a code that constrains modified gravity (MG) models viable at cosmological scales using determination of galaxy cluster mass profiles with kinematics and lensing analyses. I will discuss limitations and future developments of this method in view of upcoming imaging and spectroscopic surveys, as well as the possibilities of including X-ray data to break degeneracy among model parameters. Finally I will show preliminary results about the constraints that can be obtained on the inner slope of dark matter profiles when adding the velocity dispersion of the Brightest Central Galaxy (BCG) in the dataset of MG-MAMPOSSt.