We present recent progress in constructing a generic two-loop amplitude reduction algorithm within the computational framework HELAC. Following HELAC tree- and one-loop paradigm, we have completed the generation and validation of the two-loop amplitude-integrand. Following the well-known OPP reduction approach, we express the amplitude in terms of Feynman integrals which are further reduced to master integrals by IBP identities. We also discuss the numerical reconstruction of the amplitude in $d=4-2\epsilon$ dimensions.