You are trying to do this in a needlessly complicated and muddled way. It is better to do it the simple easy way (see below) and then realize the relation between electric potential and the hypothetical work of electrostatic field on unit positive charge that undergoes infinite motion. Also, I recommend deprecating that non-standard underscore notation for vectors, the formulae using it are hard to read.

So here is the simple derivation; first, it is important to realize that the relation

$$ \mathbf E = - \nabla V \tag{*} $$ between electric field and electric potential isn't a generally valid relation, but it is only valid for the Coulomb electric field $\mathbf E_C$ of all charges (in general, this is only a part of the total electric field). The Coulomb field is defined as integral over all space:

$$ \mathbf E_C(\mathbf x) = \int_{space} K\rho(\mathbf x') \frac{\mathbf x - \mathbf x'}{|\mathbf x - \mathbf x'|^3}\,d^3\mathbf x' $$ where $\rho$ is density of electric charge.

The electric potential in (*) is one of infinite number of possible potentials that are related to the Coulomb potential of all charges in space via arbitrary shift in its value $V_0$: $$ V(\mathbf x) = V_0 + \int_{space} K\rho(\mathbf x') \frac{1}{|\mathbf x - \mathbf x'|}\,d^3\mathbf x'. $$

Now, using these two definitions, it is easy to verify that the two quantities obey the relation (*). The nabla operator applied to $V$ means that derivatives of $V$ with respect to components of vector $\mathbf x$ are to be taken.

$$ -\nabla V (\mathbf x) = \left [ -\frac{\partial V}{\partial x_1}, -\frac{\partial V}{\partial x_2}, -\frac{\partial V}{\partial x_3} \right ] $$ You can easily show that $-\partial V/\partial x_1$ is the same as $\mathbf E_{C,1}$ and so on for components 2,3.

You are trying to do this in a needlessly complicated and muddled way (work integral with a variable limit). It isn't wrong but it more complicated and confusing. 

It is better to do it the simple easy way (see below) and then realize the relation between electric potential and the work integral with variables limit (hypothetical work of electrostatic field on unit positive charge that undergoes infinite motion). Also, I recommend deprecating that non-standard underscore notation for vectors, the formulae using it are hard to read.

So here is the simple derivation; first, it is important to realize that the relation

$$ \mathbf E = - \nabla V \tag{*} $$ between electric field and electric potential isn't a generally valid relation, but it is only valid for the Coulomb electric field $\mathbf E_C$ of all charges (in general, this is only a part of the total electric field). The Coulomb field is defined as integral over all space:

$$ \mathbf E_C(\mathbf x) = \int_{space} K\rho(\mathbf x') \frac{\mathbf x - \mathbf x'}{|\mathbf x - \mathbf x'|^3}\,d^3\mathbf x' $$ where $\rho$ is density of electric charge.

The electric potential in (*) is one of infinite number of possible potentials that are related to the Coulomb potential of all charges in space via arbitrary shift in its value $V_0$: $$ V(\mathbf x) = V_0 + \int_{space} K\rho(\mathbf x') \frac{1}{|\mathbf x - \mathbf x'|}\,d^3\mathbf x'. $$

Now, using these two definitions, it is easy to verify that the two quantities obey the relation (*). The nabla operator applied to $V$ means that derivatives of $V$ with respect to components of vector $\mathbf x$ are to be taken.

$$ -\nabla V (\mathbf x) = \left [ -\frac{\partial V}{\partial x_1}, -\frac{\partial V}{\partial x_2}, -\frac{\partial V}{\partial x_3} \right ] $$ You can easily show that $-\partial V/\partial x_1$ is the same as $\mathbf E_{C,1}$ and so on for components 2,3.