No it is not possible to derive the full relativistic law, because in the newtonian limit electrostatics and gravitostatics are identical, but the electrostatic potential is completely different from the gravitational potential in relativity.
The easiest generalization of Newton's nonrelativistic law:
$${dp\over dt} = e \nabla(\phi)$$
is to have a four dimensional scalar Higgs-type field $\phi$, then the relativistic equation of motion for a particle with a coupling e to the Higgs is
$$ {d\over d\tau}( (m+e\phi(x)) p_\mu) = e \partial_\mu \phi $$
In the limit of nearly constant $\phi$ (or small e) and small velocities, this reproduces the inverse square law, and the force is the gradient of the Higgs.
To get EM, you have to assume that $\phi$ is a zero component of a four-vector. This is the standard case, and I won't talk about it.
The other noteworthy case is when $\phi$ is the zero-zero component of a tensor. Then the full special relativistic equation of motion is the linearized geodesic equation
$${d \over d\tau} v_{\mu} = (\partial_\mu h_{\nu\sigma} - {1\over 2}\partial_\mu h_{\mu\nu} ) v^{\mu}v^{\nu}$$
And this is yet another generalization, linearized General Relativity, that reproduces the inverse square behavior for slow velocities, from the gradient of $h_{00}$. You can rule out possibilities 1 and 3 by the observation that there is repulsion between like charges, but this requires Maxwell's equations, and you asked to do it from the force law only.