In order to describe self-bound systems, one needs a nonperturbative approach. We discuss the relativistic bound state equations of QED and QCD formulated in Basis Light-Front Quantization. In this approach, the light-front direction is discretized, and two-dimensional harmonic oscillator basis functions are used for the transverse direction. At present, the Fockspace in our calculations is limited to the minimal Fock sector, but the extension to an arbitrary number of (anti-)fermions and gauge bosons is in principle straightforward. We present initial results for the energies and distribution functions of two-body bound states obtained within this approach and discuss convergence issues.