I came here before asking myself the same question and, as I figured it out, I'd like to answer just for the record.
Respecting the former answer, with plain and usual commutators (e.g. see this question) in the Schrödinger picture with locality and causality, if $X(\xi)=\sum_i\left(\alpha_i(\xi){U}+b_i(\xi){V}\right)$, the condition $[X(\xi),X(\xi^\prime)]=0$ (which comes from definition of the usual commutator) implies $[U,V]=0$ (yes, explicitly you will get something non-zero, but this must be equated to the zero that comes from the definition of the commutator of the operator with itself) and thus also $[\partial_\xi{X}(\xi),X(\xi^\prime)]=0$.
Now, in general for any field $\phi$ in the Heisenberg picture,
$$[\partial_i\phi(x),\phi(y)]=\partial_i[\phi(x),\phi(y)]+[\partial_i,\phi(y)]\phi(x)\tag{1}$$ where $x$ and $y$ are 4-vectors and $i$ stands for $x^i$ derivatives in the spatial components $i=1,2,3$. Now, $[\partial_i,\phi(y)]\phi(x)=0$ as all derivatives are in $x^i$, and thus, simply
$$[\partial_i\phi(x),\phi(y)]=\partial_i[\phi(x),\phi(y)]=0\tag{2}$$
as $[\phi(x),\phi(y)]=0$ for spacelike separation and particularly for equal times $x^0=y^0$, which as far as I've seen are usually assumed in the Heisenberg picture unless said otherwise. As user Pisanty commented in the former answer, at least in the Schrödinger picture, scalars always commute with themselves and with their derivatives. The case of the time component is different depending on spacetime separation; if you want to compute $[\dot\phi,\phi]$ directly (at equal times) I guess you should either know the Fourier decomposition of $\dot\phi$ or compute $\dot\phi=i[H,\phi]$ using the Hamiltonian $H$, at least you get from $[\partial_0,\phi(y)]\phi(x)$ that $[\dot\phi(x),\phi(y)]$ should be an even operator function in space.
And so, if you know that $[A(x),B(y)]=C(x,y)$, then $[\partial_xA(x),B(y)]=\partial_xC(x,y)$.
Let me exemplify this with the equations of motion of the real scalar theory (you should be able to extend it for your complex scalar one):
$$\mathcal{L}=\frac{1}{2}\partial_\mu\phi\partial^\mu\phi-\frac{1}{2}m^2\phi^2\,\Longrightarrow\,H=\frac{1}{2}\int{d^3x}\left(\pi^2+(\nabla\phi)^2+m^2\phi^2\right)\tag{3}$$
Then you get
\begin{align}\dot\phi(x)=i[H,\phi(x)]&=\frac{i}{2}\int{d^3y}[\pi^2(y)+(\nabla\phi)^2(y)+m^2\phi^2(y),\phi(x)]\\&=\frac{i}{2}\int{d^3y}[\pi^2(y),\phi(x)]=\pi(x)\tag{4}\end{align}
where I used $[(\nabla\phi)^2(y),\phi(x)]=0$ because of the reasons above. From here you also get in particular for this theory $[\dot\phi,\phi]=[\pi,\phi]=-i\delta(\vec{x}-\vec{y})$. And for $\pi$,
\begin{align}\dot\pi(x)=i[H,\pi(x)]&=\frac{i}{2}\int{d^3y}[\pi^2(y)+(\nabla\phi)^2(y)+m^2\phi^2(y),\pi(x)]\\&=\frac{i}{2}\int{d^3y}[(\nabla\phi)^2(y)+m^2\phi^2(y),\pi(x)]\tag{5}\end{align}
from where, using the above statements,
\begin{align}\int{d^3y}[(\nabla\phi)^2(y),\pi(x)]&=\int{d^3y}\left(\nabla\phi\cdot[\nabla\phi,\pi]+[\nabla\phi,\pi]\cdot\nabla\phi\right)\\
&=2\int{d^3y}(\nabla\phi\cdot\nabla[\phi,\pi])\\
&=2i\int{d^3y}\left(\nabla\phi\cdot\nabla\delta(\vec{x}-\vec{y})\right)\\
&=2i\nabla\phi\cdot\vec\delta(\vec{x}-\vec{y})\bigg|_{\pm\infty}-2i\int{d^3y}\nabla^2\phi\,\delta(\vec{x}-\vec{y})\\
&=-2i\nabla^2\phi(x)\tag{6}\end{align}
where I used the commutation of the Euclidean dot product and good old integration by parts, I also stopped writing the variable dependence but just remember that $\nabla$ acts on $y$'s only. So this way one gets $\dot\pi=\nabla^2\phi-m^2\phi^2$ which using $(4)$ gives the Klein-Gordon equation for $\phi$ as expected.
