Thanks to Jakobian for pointing it out: some technical details require stronger hypothesis for this approach to work! It suffices to let $u\in L^1[0,1]\setminus L^{p-r}[0,1]$ for some $r>0$ small enough.
The main idea is to tweak the usual "dual" function so as to get the desired $w$.
Let $u\in L^{1}[0,1] \setminus L^{p}[0,1]$. A standard "dual" function in $L^{p'}[0,1]$ would be $v = \operatorname{sign}(u)|u|^{p-1}$ - if $u \in L^p[0,1]$. However, in this case $v$ is not in $L^{p'}[0,1]$, since
$$
\int_{0}^{1} |u(x)|^{p'(p-1)} dx = \int_{0}^{1} |u(x)|^p dx = \infty.
$$
Now, lets fix this somehow. Let $g_t = \frac{1}{|u|^t}\mathbf{1}_{\{|u|>1\}}$, and let
$$
s=\inf_{t\in [0,\infty)}\{ug_t \in L^{p}[0,1]\};$$
note that since $u\in L^{1}[0,1]\setminus L^{p}[0,1]$ we have that the infimum exists, as $t=\frac{1}{p'}$ works and $t=0$ does not.
When $u\notin L^{p-r}[0,1]$ we can deduce that $t=\frac{p-r}{p}>0$ does not work either so furthermore $s>0$.
We have that then if $t>s$, $w=v(g_t^{p-1})$ is in $L^{p'}$, as
$$
\int_{0}^{1} |v(x)(g_t(x))^{p-1}|^{p'} dx = \int_{0}^{1} |u|^{(1-t)p} \mathbf{1}_{\{|u|>1\}} dx = \int_{0}^{1} |ug_t|^p < \infty.
$$
However,
$$
\int_{0}^{1} |u(x)v(x)(g_t(x))^{p-1}| dx = \int_{0}^{1} |u(x)|^{1+(1-t)(p-1)} \mathbf{1}_{\{|u|>1\}}(x) dx = \int_{0}^{1} |u(x)|^{(1-t')p} \mathbf{1}_{\{|u|>1\}}(x) dx,
$$
where $t'=t-t/p$, as
$$
(1-t')p = (1-t)p + t = 1 + (1-t)(p-1).
$$
We now just need to take $t$ close enough to $s$ such that $t'< s$ and we get that for $w=v(g_t)^{p-1}$, $uw\notin L^{1}[0,1]$. This is possible when $s>0$ as $t(1-\frac{1}{p})$ can get arbitrarily close to $0$, but it is not clear how to proceed when $s=0$.