How to find $\mathrm{\sum_{n=1}^\infty E_{n}(n)= \lim_{k\to\infty}\int_1^\infty \frac{e^{kt+t}t^k-e^t}{e^{kt+2t}t^{k+1}-e^{kt+t}t^k}dt=.26929…}$? shows a nice result from your integral which will be generalized here using geometric series and gamma type functions. The result is simple, but has the restriction of x<Ω=W(1) which is the Omega constant:
$$\int \frac{dx}{1-xe^x}=-\int \sum_{n=0}^\infty \left(xe^x\right)^ndx=-\sum_{n=0}^\infty \int x^n e^{xn}dx$$
This is a close definition of the gamma function and exponential integral function mentioned in the bolded link:
$$\sum_{n=0}^\infty \int x^n e^{xn}dx=C-\sum_{n=0}^\infty \frac{ (e^{-n x} (e^x x)^n (-n x)^{-n} Γ(n + 1, -n x))}{n}\mathop=^{x\in\Bbb R} C-\sum_{n=0}^\infty x^{n+1} \text E_{-n}(-n x)=C+\sum_{n=0}^\infty \frac{(-1)^nΓ(1 + n, -n x)}{ n^{n+1}} $$
If there was some function, like a hypergeometric function to give a closed form, please post it. Please correct me and give me feedback!