1H NMR spectra of biofluids provides a wealth of biochemical information on the metabolic status of an organism. Through the application of pattern recognition and classification algorithms, the data have been shown to provide information on disease diagnosis and the beneficial and adverse effects of potential therapeutics. Here, a novel approach is described for identifying subsets of spectral patterns in databases of NMR spectra, and it is shown that the intensities of these spectral patterns can be related to the onset and recovery from a toxic lesion in both a time-related and dose-related fashion. These patterns form a new type of combination biomarker for the biological effect under study. The approach is illustrated with a study of liver toxicity in rats using NMR spectra of urine following administration of a model hepatotoxin hydrazine.