Let $L:=-\Delta +V$ be a Schrödinger operator on ${ℝ}^n$ with $n\ge 3$ and $V\ge 0$ satisfying ${\Delta }^{-1}V\in L^{\infty }\left({ℝ}^n\right)$. Assume that $\varphi :{ℝ}^n\times [0,\infty )\to [0,\infty )$ is a function such that $\varphi (x,·)$ is an Orlicz function, $\varphi (·,t)\in {𝔸}_{\infty }\left({ℝ}^n\right)$ (the class of uniformly Muckenhoupt weights). Let $w$ be an $L$-harmonic function on ${ℝ}^n$ with $0<C_1\le w\le C_2$, where $C_1$ and $C_2$ are positive constants. In this article, the author proves that the mapping $H_{\varphi ,L}\left({ℝ}^n\right)\ni f\mapsto wf\in H_{\varphi }\left({ℝ}^n\right)$ is an isomorphism from the Musielak-Orlicz-Hardy space associated with $L$, $H_{\varphi ,L}\left({ℝ}^n\right)$, to the Musielak-Orlicz-Hardy space $H_{\varphi }\left({ℝ}^n\right)$ under some assumptions on $\varphi$. As applications, the author further obtains the...