We use heavy baryon chiral perturbation theory to evaluate the two-photon exchange corrections to the low-energy elastic lepton-proton scattering at next-to-leading order accuracy, i.e., O(α, M −1 ), including a non-zero lepton mass. We consider the elastic proton intermediate state in the two-photon exchange invoking soft photon approximation. The infrared singular contributions are projected out using dimensional regularization. The resulting infrared singularity-free two-photon exchange contribution is in good numerical agreement with existing predictions based on standard diagrammatic soft photon approximation evaluations.