Not yet reviewed by Pith; the record is open.

We compute the part of the two-loop virtual amplitude for the process $gg \to V_1 V_2 \to (l_1 \bar l'_{1}) (l_2 \bar l'_2)$, where $V_{1,2}$ are arbitrary electroweak gauge bosons, that receives contributions from loops of massless quarks. Invariant masses of electroweak bosons are allowed to be different from each other. Our result provides an important ingredient for improving the description of gluon fusion contribution to the production of four-lepton final states at the LHC.