We present the result of the search for the rare B meson decay of B+ -> l(+)nu e gamma with l = e, mu. For the search the full data set recorded by the Belle experiment of 711 fb(-1) integrated luminosity near the Upsilon(4S) resonance is used. Signal candidates are reconstructed for photon energies E-gamma larger than 1 GeV using a novel multivariate tagging algorithm. The novel algorithm fully reconstructs the second B meson produced in the collision using hadronic modes and was specifically trained to recognize the signal signature in combination with hadronic tag-side B meson decays. This approach greatly enhances the performance. Background processes that can mimic this signature, mainly charmless semileptonic decays and continuum processes, are suppressed using multivariate methods. The number of signal candidates is determined by analyzing the missing mass squared distribution as inferred from the signal side particles and the kinematic properties of the tag-side B meson. No significant excess over the background-only hypothesis is observed and upper limits on the partial branching fraction Delta B with E-gamma > 1 GeV individually for electron and muon final states as well as for the average branching fraction of both lepton final states are reported. We find a Bayesian upper limit of Delta B(B+ -> l(+)nu el) < 3.0 x 10(-6) at 90% CL and also report an upper limit on the first inverse moment of the light-cone distribution amplitude of the B meson of lambda(B) > 0.24 GeV at 90% CL.