Solving boundary value problems with boundary element methods requires specific Green functions suited to the boundary conditions of the problem. Using vector algebra, one often needs to use a Green function for the Helmholtz equation whereas it is a solution of the perturbed Dirac equation that is required for solving electromagnetic problems using Clifford algebra. A wealth of different Green functions of the Helmholtz equation are already documented in the literature. However, perturbed Dirac equation is only solved for the generic case and only its fundamental solution is reported. In this paper, we present a simple framework to use known Green functions of Helmholtz equation to construct the corresponding Green functions of perturbed Dirac equation which are essential in finding the appropriate kernels for integral equations of electromagnetic problems. The procedure is further demonstrated in a few examples.