Porphyrins are a class of photosensitizers used in photodynamic therapy (PDT). Understanding the interaction of porphyrins with membrane cells components is important in order to improve this therapy. Many analytical methods can be used for this aim. High performance liquid chromatography (HPLC) was used for the separation of porphyrins on RP and HILIC stationary phases as well as on a biomimetic membrane IAM phase. Twenty-six tetraphenyl porphyrins (TPP) were successfully separated on an IAM column, a C18 Gravity RP column, a C8 Gravity RP column, a PolarTec RP column and a HILIC column. Stationary phases were chosen as the most appropriate to cover the study of different types of interactions. Elution was performed with a 45 min linear gradient. Obtained gradient retention times were converted to gradient chromatography hydrophobicity index (CHI) and to an apparent retention factor (kapp). The partition coefficients (logP) of the 26 compounds were measured in a 2-octanol/PBS system and estimated in silico. Correlation between kapp values was studied. Moreover, a multivariate analysis was performed to explain columns relationships. Obtained results show that porphyrins are separated mainly according to hydrophobic interactions that are relative to their structure (sugar number and the disposition around the porphyrin macrocycle).