I prefer to work in four-vector format: (and use [itex]m=m_0[/itex] 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 [itex]g_{\mu\nu} \sim diag(1,-1,-1,-1)[/itex] indicating the metric measures proper time (rather than proper distance).
The action then is the path integral:
[tex]S = \int_{x(\tau)}[(mc)ds + e/c A_\mu dx^\mu] = \int L d\tau[/tex]
with [itex]L=mc\sqrt{\dot{x}^\mu \dot{x}_\mu} + e/c A_\mu \dot{x}^\mu[/itex].
Note that [itex]\tau[/itex] is a path parameter not necessarily proper time.
The canonical momentum-energy four vector components are:
[tex]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[/tex]
The Euler-Lagrange equations are then:[itex]\dot{P}_\mu = (e/c)\dot{x}^\nu A_{\nu,\mu}[/itex] or...
[tex]\dot{p}_\mu = (e/c) \dot{x}^\nu [A_{\nu,\mu}-A_{\mu,\nu}] = (e/c)\dot{x}^\nu F_{\mu\nu}[/tex]
OK Now I will attempt toapply Noether's theorem. We need to write the variation of the action in the form:[itex]\delta S = \int Q_\mu \delta\dot{x}^\mu d\tau[/itex] and then [itex]\dot{Q}_\mu = 0[/itex] along paths satisfying the dynamics.
[itex]\delta S = \int d\tau [P_\mu \delta\dot{x}^\mu + e/c\dot{x}^\nu A_{\nu,\mu} \delta x^\mu][/itex]
By expanding [itex]e/c \frac{d}{d\tau}(x^\nu A_{\nu,\mu} dx^\mu)[/itex] and substitution with integration by parts (zeroing variation on the boundary) this yields:
[tex]\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[/tex]
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:
[itex]Q_\mu = P_\mu -( e/c )x^\nu A_{\nu,\mu}[/itex] are conserved.
Now the fact that the four potential is linear gives us:[itex]A_\nu = W_{\nu\mu}x^\mu[/itex] with the [itex]W[/itex]'s constant and that the e-m field tensor components are:
[tex]F_{\nu\mu} = A_{\mu,\nu}-A_{\nu,\mu} = W_{\mu\nu}-W_{\nu\mu}\equiv 2W_{[\mu\nu]}[/tex]
The bracketed index indicates the anti-symmetric component of W: [itex]W_{\mu\nu} = W_{[\mu\nu]}+W_{(\mu\nu)}[/itex] 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. [itex]\Lambda = (1/2)W_{(\mu\nu)}x^\mu x^\nu[/itex], [itex]A_\mu \to A_\mu - \partial_\mu \Lambda[/itex].
This means that [itex]A_{\nu,\mu} = -A_{\mu,\nu}[/itex] and thus that [itex]x^\nu A_{\nu,\mu} = - x^\nu A_{\mu,\nu} = -A_\mu[/itex].
So the conserved quantities are:
[tex]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}[/tex]
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.)