Let L = −∆ + V be a Schrödinger operator, where ∆ is the Laplacian on R d and the nonnegative potential V belongs to the reverse Hölder class RH q for q ≥ d. The Riesz transform associated with the operator L = −∆ + V is denoted by R = ∇(−∆ + V ) −1/2 and the dual Riesz transform is denoted by R * = (−∆ + V ) −1/2 ∇. In this paper, we first introduce some kinds of weighted Morrey spaces related to certain nonnegative potentials belonging to the reverse Hölder class RH q for q ≥ d. Then we will establish the boundedness properties of the operators R and its adjoint R * on these new spaces. Furthermore, weighted strong-type estimate and weighted endpoint estimate for the corresponding commutators [b, R] and [b, R * ] are also obtained. The classes of weights, the classes of symbol functions as well as weighted Morrey spaces discussed in this paper are larger than A p , BMO(R d ) and L p,κ (w) corresponding to the classical Riesz transforms (V ≡ 0).