The general change of Numeraire formula gives the following Radon-Nikodym derivative:
$$ \frac{dN_2}{dN_1}(t)|\mathcal{F}_{t_0}=\frac{N_1(t_0)N_2(t)}{N_1(t)N_2(t_0)} $$
I am able to derive this Radon-Nikodym for specific examples, such as changing from the risk-neutral measure $Q$ to the T-Forward Measure associated with a zero-coupon bond $P(t_0,t)$: in this case, we have under $Q$:
$$ \frac{S_0}{N_Q(t_0)=1}=\mathbb{E}^Q\left[\frac{S_t}{N_Q(t)=e^{rt}}|\mathcal{F}_{t_0}\right] $$
So that:
$$ (i) S_0 = \mathbb{E}^Q\left[S_t\frac{N_Q(t_0)=1}{N_Q(t)=e^{rt}}|\mathcal{F}_{t_0}\right] $$
Under the T-forward Bond numeraire:
$$ \frac{S_0}{N_{P}(t_0)=P(t_0,t)}=\mathbb{E}^{P_t}\left[\frac{S_t}{N_P(t)=1}|\mathcal{F}_{t_0}\right]$$
So that:
$$(ii) S_0 = P(t_0,t)\mathbb{E}^{P_t}\left[\frac{S_t}{N_P(t)=1}|\mathcal{F}_{t_0}\right]$$
Equating (i) to (ii) we get:
$$\mathbb{E}^Q\left[S_t\frac{N_Q(t_0)}{N_Q(t)}|\mathcal{F}_{t_0}\right]=N_P(t_0)\mathbb{E}^{P_t}\left[\frac{S_t}{N_P(t)}|\mathcal{F}_{t_0}\right]$$
Since $N_P(t)$ at time $t$ is by definition constant (equal to one), it is easy to take it out of the expectation and group all the Numeraire terms on the LHS, so that:
$$ \mathbb{E}^Q\left[S_t\frac{N_Q(t_0)N_P(t)}{N_Q(t)N_P(t_0)}|\mathcal{F}_{t_0}\right]=\mathbb{E}^{P_t}\left[S_t|\mathcal{F}_{t_0}\right] $$
And the result follows be inspection.
Note: in general, the numeraire $N_2(t)$ would not be a constant at time $t$, as is the case for the numeraire associated with the T-forward maturing bond. So it would not be possible to take $N_2(t)$ out of the expectation $\mathbb{E}_{t_0}^{N_2}[]$ as in the case above. It would therefore not be so straight forward to group all the numeraire terms and deduce the Radon-Nikodym derivative by inspection.
Question: How can the change of Numeraire Radon-Nikodym formula be derived or proved in the general case? (not thinking about specific numeraires as in the case above).