Variational optimization of the 2DM: approaching three-index accuracy using extended cluster constraints