Could you try removing step = pm.Metropolis() and just call trace = pm.sample(1000)?
step = pm.Metropolis()
trace = pm.sample(1000)