Find hamilton cycle in a directed graph reduced to sat problem
Just add the constraint that $x_{1j}=x_{nj}$ for all $j$, and make a special exception to constraint 2 so that you omit constraint 2 in the case where $i=1$ and $k=n$.

Everything works fine for a directed graph.  You just need to interpret constraint 5 appropriately.