We present a systematic theoretical study of bismuth selenide electronic structure in the presence of spin-orbit interaction. It is confirmed that the inversion of Bi pz and Se pz states in the bulk band structure is induced by spin-orbit coupling. The energy gap in the bulk electronic structure is closed by pz states associated mainly with Bi atoms from the first quintuple layer. We elucidate the influence of Bi substitution in the fifth atomic layer on the local density of electronic states in the top quintuple layer. Finally, we describe how it develops into the additional p states in the surface electronic structure, which produce the triangular protrusion observed in Scanning Tunnelling Microscopy measurements and characteristic narrow features in Scanning Tunnelling Spectroscopy results.