In this contribution, we discuss a molecular representation mode for the generation of reduced descriptions of flexible molecules in various conformations. The representations of the endothiapepsin ligands are constituted by graphs of peaks obtained through a hierarchical merging algorithm which combines the location of promolecular electron density (ED) maxima with the decomposition of the molecular structures into fragments. The representations are then aligned through the use of a Monte Carlo/Simulated Annealing procedure. The evaluation function of the alignment solutions is based on the local density values of the peaks and their inter-distances. The applications show that the alignment of a given molecule onto itself, in a different conformation, is successful when a pseudo-topological path length is considered in the evaluation function of the solution, while Cartesian distances are more adapted to the scoring for alignments of two different molecules in their co-crystallized conformation. Results are compared with the available literature data.