Configuring Q1 as voltage follower will not give 1kV.
As for the original circuit, the first thing is that stability will depend on Q1 stage gain and overall closed loop gain, i.e. the ratios R1/R2 and R5/R4. So set them to something realistic (100~200x or so, I guess) if you want your prototype to be realistic. Well, this may screw up DC bias of the prototype, so if it can't work, at least try to make the ratio (R1/R2)/(R5/R4) be realistic - I think near unity would be OK. This ratio goes into loop gain, so lower is more stable. But you need to mind things like Q1 gate voltage and available output swing of U1.
Then I don't know, perhaps try 10kΩ between the control voltage source and U1, plus a fairly substantial (dunno, 100nF or more?) local feedback capacitor around U1 (OUT to IN-). This is where playing with SPICE AC analysis could make sense. The big unknown is Q1 - it will have flat gain and close to 180° phase up to some frequency, then start to fall off and lag 270°. You want your overall loop gain to fall below unity by the time it happens, and the direct feedback around U1 should provide for that. I think.
An alternative option for high voltage is the old "boosted opamp" circuit which takes signals from the opamp's supply pins. Maybe slightly more complex than the above, provided that the above would work.