A program SPIRIT (Simulation Program considering Incomplete Recovery of z magnetization and INEPT Transfer efficiency) has been developed to simulate three-dimensional NOESY-HSQC spectra. This program takes into account (1) different transfer efficiency during INEPT and reverse INEPT durations due to differential relaxation rates and 1J coupling constants; (2) the different effect of the sensitivity-enhancement scheme on CH, CH2 and CH3 systems; and (3) incomplete recovery of longitudinal magnetization between scans. The simulation program incorporates anisotropic tumbling mode for symmetric tops, and allows for differential external relaxation rates for protons. Some well-defined internal motions, such as the fast rotation of methyl groups, are taken into account. The simulation program also allows for input of multiple conformations and their relative populations to calculate the average relaxation matrix to account for slow internal motions. With the SPIRIT program, the sensitivity-enhanced NOESY-HSQC experiment can be used directly in the evaluation of the accuracy of structures, which can potentially be improved by direct refinement against the primary data.