The results of electronic structure calculations based on density functional theory (DFT) for protic ionic liquids (PILs) consisting of triethanolammonium cation paired with anion of different sulfonic acids are reported. The influence of the anion nature on the structure and interactions in the ion pairs that are formed in these PILs is discussed in detail. Multiple H-bonding interactions exist between the protons in the NH/OH groups of the cation and different oxygen atoms of the acid anion in the ion pairs. The quantum theory of "atoms in molecules" has been used to estimate the individual contributions of each hydrogen bond to the stability of the ion pair. The hydrogen-bonding interactions in the ion pairs vary in their strength ranging from weak to moderately strong. In addition to these hydrogen bonds, there are other dispersion and electrostatic-dominant interactions that play an important role in the overall stability of PILs and their physicochemical properties. Aided by results from our previous DFT studies of triethanolammonium class of PILs with inorganic anions, these new data allow us to gain an improved understanding of the structure-property relationships in the studied ionic liquids. Close to linear correlation, in particular, has been found between the melting points and the binding energies of the cation and anion in the ion pairs.