w = Total[
QuantumState /@ {"10000", "01000", "00100", "00010", "00001"}];
s = QuantumPartialTranspose[w, {2}];
s["DensityMatrix"] // MatrixForm
QuantumPartialTranspose[state,{j}]
with j the qubit with respect to which you do the partial transpose
Note the density matrix of w is a 2^5x2^5 matrix.