27481 is a 4k+1 prime and so it can be expressed as a sum of two squares in one way only viz. 165^2 + 16^2. The sequence can be generated as follows:

INPUT

L=[]

for p in prime_range(40000):

  cube=p^3

  if Set(cube.digits())=={0,1,2,3,4,5,6,7,8,9}:

    L.append(p)

print(L)

print(27481^3)

OUTPUT

[5437, 6221, 7219, 8443, 10903, 11353, 15937, 17123, 18229, 19429, 20353, 20903, 20929, 21803, 21841, 21961, 22123, 22283, 22993, 23053, 23369, 23663, 24733, 25183, 25219, 25463, 26317, 26387, 26449, 27127, 27481, 28181, 28631, 28711, 28961, 29059, 29443, 29501, 30169, 31153, 31183, 32213, 32801, 33739, 33797, 33811, 33941, 34283, 35027, 35051, 35729, 35963, 36137, 36251, 36383, 36809, 36943, 37223, 37369, 37511, 37619, 37967, 38281, 38917]

20753798525641

27481 is a member of OEIS A124629: primes p such that their cubes are pandigital, here 20753798525641.

27481 forms a twin prime with 27479.

27481 can be rendered as a digit equation as follows:

(-2 +7 -4) ^ 8 = 1