Let $f(z)=\sum_{k=0}^\infty a_k(z-z_0)^k$ with radius of convergence $R$ then $f$ is analytic on the open disk around $z_0$ with radius $R$.
What I was thinking about is an approach based on this sketch:
I want to find a new power series based on a new disk with center $\widetilde{z_0}$ and variable $z$. Therefore I am using $z-z_0 = (\widetilde{z_0}-z_0)+(z-\widetilde{z_0})$ which yields
$$ f(z) = \sum_{k=0}^\infty a_k\left( (\widetilde{z_0}-z_0)+(z-\widetilde{z_0}) \right)^k = \sum_{k=0}^\infty a_k\sum_{n=0}^k\binom{k}{n}(\widetilde{z_0}-z_0)^n(z-\widetilde{z_0})^{k-n}. $$
I am looking for power series of the form $\sum_{k=0}^\infty b_k(z-\widetilde{z_0})^k$ but haven't been able to come up with the next step. What is your suggestion based on the method of using the binomial formula (explicitly no derivatives of power series etc. if possible)?

Here are two variants. We start with proceeding with your approach.
Comment:
In (2) we change the index summation $n \rightarrow k-n$
In (3) we change the order of summation according to (1)
Another variant could be based upon Taylor series expansion.