Differential evolution Markov Chain Monte Carlo for fitting transit timing variations with N-body integrations