The gas-phase conformational potential energy surfaces (PES) of the last, 22nd amino acid pyrrolysine and related derivatives (neutral, deprotonated, and protonated) were extensively searched for the first time. By considering all possible combinations of the single-bond rotational degrees of freedom with a semiempirical and ab initio combined computational approach, a large set of unique low-energy conformers was identified for each pyrrolysine species, and essential properties such as vibrational frequencies, dipole moments, rotational constants, and intramolecular hydrogen bonding configurations were presented and characterized. The conformational electronic energies and thermochemical properties of proton affinity/dissociation energy (PA/PDE) and gas-phase acidity/basicity (GA/GB) were determined by the density functional BHandHLYP, B3LYP, and M062X, and Moller-Plesset MP2 methods. The MP2 and DFT methods are found to predict disparate PES for neutral and protonated conformations and sufficiently different thermochemical data. The measurements of dipole moments and characteristic IR modes at low temperature as well as GA/GB are demonstrated to be feasible approaches to verify the theoretical predictions.