Whether you use second quantization formalism or whether you are even talking about classical or quantum systems current is defined via the continuity equation for some quantity, $\hat{O}$,
$$\frac{\partial \hat{O}}{\partial t} + \nabla \cdot \mathbf{\hat{J}} = 0,$$
where I have used hat to denote we are talking about quantum mechanical observables.
The question is can we find a pair of observables for which the above equation holds.
For integrable models, such as the Hubbard model, Heisenberg spin chain model, free fermions the answer is yes. We can identify local conserved charges for which the above equation holds.
Now, in the Heisenberg picture we have,
$$\frac{d}{dt}\hat{O}(t)=\frac{i}{\hbar}[H,\hat{O}(t)]$$
So if you have some Hamiltonian and some corresponding local conserved charge you compute it's commutator with the Hamiltonian and use that to find the current operator.
For instance, for the Hubbard model,
$$ H = -t \sum_{\langle i,j \rangle,\sigma}( c^{\dagger}_{i,\sigma} c^{}_{j,\sigma}+ h.c.) + U \sum_{i=1}^{N} n_{i\uparrow} n_{i\downarrow}$$
the number density current at site $i$, $n_i=c^\dagger_i c_i$ can be easily found by a discretized version of the continuity equation,
$$\frac{i}{\hbar}[H,n_i(t)]=- (J_{i+1}-J_{i}) = -t(i c^\dagger_{i+1} c_i+ h.c.)+t(i c^\dagger_{i} c_{i-1}+ h.c.),$$
which allows us to identify,
$$J_i=-t(-i c^\dagger_{i} c_{i-1}+ h.c.) $$
as the particle number current density at site $i$.
Please note that the total (sum for all sites) particle number density commutes with the Hamiltonian, assuming periodic boundary conditions. This is very important for integrability.
For more examples (e.g., energy current in the Hubbard model), see this paper:
http://arxiv.org/abs/cond-mat/9611007
