Plotting the sequence 1, 4, 22, 139, ... as a scatter diagram, one can fit a 3rd degree polynomial curve with equation y = 14x^3 - 34.5x^2 + 23.5x + 1 (x = 0, 1, 2, 3). Using this graph as an estimator, one can guess that the number of shapes will be in the region of 439 if four are removed. Therefore the guestimate for 500 is not far out!
T25 appears to be 9 (see attached image) and I agree that T26 is 1. The roots indicate the distance between the two cubelet-middle points in the grid (x = y = z = 1).
Thus, would the sequence end with ..., 9, 1, 0 (the 0 indicating no shape at all)?
Will determine T24 a bit later!