A systematic study has been performed of the proton affinity of a large family of agonists and antagonists of the alpha1-adrenoceptor at the B3LYP/6-31G* level of theory. After a conformational search, all the N atoms were considered as protonation sites and protonation energy values were determined. The inclusion of solvation by means of the Onsager model yielded stabilization in the proton affinity values obtained. In addition, a good correlation was found between the proton affinity values corresponding to the first protonation in gas phase of some of the compounds and their corresponding experimental affinity constants K(i) for the alpha1A adrenergic receptor.