We study the cosmological dynamics of a recently proposed infrared modification of the Einstein equations, based on the introduction of a nonlocal term constructed with m2gμν□-1R, where m is a mass parameter. The theory generates automatically a dynamical dark energy component, that can reproduce the observed value of the dark energy density without introducing a cosmological constant. Fixing m so to reproduce the observed value ΩDE≃0.68, and writing wDE(a)= w0+(1-a)wa, the model provides a neat prediction for the equation of state parameters of dark energy, w0≃-1.042 and wa≃-0.020, and more generally provides a pure prediction for wDE as a function of redshift. We show that, because of some freedom in the definition of □-1, one can extend the construction so to define a more general family of nonlocal models. However, in a first approximation this turns out to be equivalent to adding an explicit cosmological constant term on top of the dynamical dark energy component. This leads to an extended model with two parameters, ΩΛ and m. Even in this extension the EOS parameter w0 is always on the phantom side, in the range -1.33 ≲w0≤-1, and there is a prediction for the relation between w0 and wa.