Depending on the chemical scaffold, the bioactive species could reflect the interplay between ionization and tautomerism, often complicated by the possibility to populate different conformational states in the case of flexible ligands. In this context, theoretical methods can be valuable to discern the role of these factors, as shown here for βenamino esters of 1,2,3,6-tetrahydroazocino fused ring systems, some of which had proven to be suitable scaffolds for designing novel acetylcholinesterase inhibitors. The compounds investigated herein form two clusters with distinctive experimental pK a values (i.e., α,β-diesters and β-esters ranging within 6.1-7.3 and 8.2-9.0 pK a intervals, respectively), which implies a drastic difference in the most populated species at physiological conditions. While chemoinformatic tools did not provide a consistent description of the actual pK a values, the theoretical analysis performed for the protonated and neutral species of these compounds revealed a marked change in the tautomeric preference of the tetrahydroazocine moiety upon (de)protonation. Excellent agreement between calculated and experimental pK a values was found when the tautomeric preference of protonated and neutral species was considered. Overall, this study highlights the potential use of high-level computational methods to disclose the mutual influence between ionization, tautomerism and conformational preferences in multifunctional (bio)organic compounds.