I prefer to work in four-vector format: (and use m=m_0 avoiding reference to "relativistic mass" all together.)
You can find a nice brief clarification of the relativistic Lagrangian for choices of metric at:
http://arxiv.org/pdf/0912.0655.pdf".
I choose the convention g_{\mu\nu} \sim diag(1,-1,-1,-1) indicating the metric measures proper time (rather than proper distance).
The action then is the path integral:
S = \int_{x(\tau)}[(mc)ds + e/c A_\mu dx^\mu] = \int L d\tau
with L=mc\sqrt{\dot{x}^\mu \dot{x}_\mu} + e/c A_\mu \dot{x}^\mu.
Note that \tau is a path parameter not necessarily proper time.
The canonical momentum-energy four vector components are:
P_\mu =\partial_{\dot{x}^\mu}L= \frac{mc}{\sqrt{{ \dot{x}}_\mu \dot{x}^\mu}} \dot{x}^\mu+ e/c A_\mu= p_\mu + (e/c) A_\mu
The Euler-Lagrange equations are then:\dot{P}_\mu = (e/c)\dot{x}^\nu A_{\nu,\mu} or...
\dot{p}_\mu = (e/c) \dot{x}^\nu [A_{\nu,\mu}-A_{\mu,\nu}] = (e/c)\dot{x}^\nu F_{\mu\nu}
OK Now I will attempt toapply Noether's theorem. We need to write the variation of the action in the form:\delta S = \int Q_\mu \delta\dot{x}^\mu d\tau and then \dot{Q}_\mu = 0 along paths satisfying the dynamics.
\delta S = \int d\tau [P_\mu \delta\dot{x}^\mu + e/c\dot{x}^\nu A_{\nu,\mu} \delta x^\mu]
By expanding e/c \frac{d}{d\tau}(x^\nu A_{\nu,\mu} dx^\mu) and substitution with integration by parts (zeroing variation on the boundary) this yields:
\delta S = \int d\tau [P_\mu - e/cx^\nu A_{\nu,\mu}]\delta \dot{x}^\mu - x^\nu\dot{x}^\lambda A_{\nu,\mu\lambda}\delta x^\mu
We can take different tracks from here but let us assume that the four-potential is linear in the coordinates eliminating the second order derivatives and thus the quantities:
Q_\mu = P_\mu -( e/c )x^\nu A_{\nu,\mu} are conserved.
Now the fact that the four potential is linear gives us:A_\nu = W_{\nu\mu}x^\mu with the W's constant and that the e-m field tensor components are:
F_{\nu\mu} = A_{\mu,\nu}-A_{\nu,\mu} = W_{\mu\nu}-W_{\nu\mu}\equiv 2W_{[\mu\nu]}
The bracketed index indicates the anti-symmetric component of W: W_{\mu\nu} = W_{[\mu\nu]}+W_{(\mu\nu)} the sum of anti-symmetric and symmetric components.
Now since the symmetric components contribute nothing to the physics we can eliminate them with a gauge transformation. \Lambda = (1/2)W_{(\mu\nu)}x^\mu x^\nu, A_\mu \to A_\mu - \partial_\mu \Lambda.
This means that A_{\nu,\mu} = -A_{\mu,\nu} and thus that x^\nu A_{\nu,\mu} = - x^\nu A_{\mu,\nu} = -A_\mu.
So the conserved quantities are:
Q_\mu = P_\mu + (e/c)A_\mu = p_\mu + 2(e/c)A_\mu=p_\mu + 2(e/c)W_{[\mu\nu]}x^\nu = p_\mu + (e/c)x^\nu F_{\nu\mu}
So here is a general result for general constant electro-magnetic fields. Note that if A is constant you can just gauge it to 0 as it means zero electromagnetic fields. My earlier post was incorrect in that it is not the canonical momentum which is conserved and qA is not the "potential momentum".
I think you can generalize this to arbitrary cases and I'm sure it's written up in the texts and available online somewhere.
(Note: I'm going crosseyed trying to keep everything straight so check my work.)