This paper compares the performance of two recently developed algorithms and methods for peak alignment of first-order NMR data of complex biological samples. The NMR spectra of such samples exhibit variations in peak position and peak shape due to variations in the sample matrix and to instrumental instabilities. The first method comprises an alignment of spectral segments with linear interpolation and shift correction to accommodate correspondence between a target and a test spectrum by a beam search or genetic algorithm. The second method is based on peak picking and needle vector representation of the NMR data with subsequent breadth-first search to establish shift corrections between the target and the test spectrum. The two proposed peak alignment methods and their respective merits are discussed for a real metabonomics application. Both alignment methods have been shown to enhance the interpretability of the resulting multivariate models, thereby increasing the prospect of detecting and following the onset of subtle biological changes reflected in the NMR data.