MP2, DFT and CCSD methods with 6-311++G** and aug-cc-pvdz basis sets have been used to probe the structural changes and relative energies of E-prop-2-ynylideneamine (I), Z-prop-2-ynylideneamine (II), prop-1,2-diene-1-imine (III) and vinyl cyanide (IV). The energy near-equivalence and provenance of preference of isomers and tautomers were investigated by NBO calculations using HF and B3LYP methods with 6-311++G** and aug-cc-pvdz basis sets. All substrates have Cs symmetry. The optimized geometries were found to be mainly theoretical method dependent. All elected levels of theory have computed I/II total energy of isomerization (ΔE) of 1.707 to 3.707 kJ/mol in favour of II at 298.15 K. MP2 and CCSD methods have indicated clearly the preference of II over III; while the B3LYP functional predicted nearly similar total energies. All tested levels of theory yielded a global II/IV tautomerization total energy (ΔE) of 137.3-148.4 kJ/mol in support of IV at 298.15 K. The negative values of ΔS indicated that IV is favoured at low temperature. At high temperature, a reverse tautomerization becomes spontaneous and II is preferred. The existence of II in space was debated through the interpretation and analysis of the thermodynamic and kinetic studies of this tautomerization reaction and the presence of similar compounds in the Interstellar Medium (ISM).