This thesis focus on the analysis of spin-lattice NMRD relaxation profilesmeasured in various complex systems such as proteins, zeolites and ionicliquids. Proton, deuterium and fluoride T1-NMRD relaxation profiles wereobtained from a fast-field cycling (FFC) instrument. It is found that alsopossible to obtain NMRD profiles from the molecular dynamics (MD)simulation trajectories. NMRD Profiles were analyzed by using differentrelaxation models, such as the Solomon-Bloembergen-Morgan (SBM) theoryand the Stochastic Liouville (SL) theory. Paper I described the hydration of protein PrxV obtained from a MDsimulation, and compared with the picture emerges from an analysis byusing a generally accepted relaxation model [appendix C]. The result showsthat the information from NMRD analysis is an averaged picture of watermolecules with similar relaxation times; and the MD simulations containsinformation of all types of interested water molecules with differentresidence times. In paper II NMRD profiles have been used to characterize the hydration ofthe oxygen-evolving complex in state S1 of photosystem II. NMRDexperiments were performed on both intact protein samples and Mndepletedsamples, and characteristic dispersion difference were foundbetween 0.03 MHz to 1 MHz; approximately. Both the SBM theory and theSL theory have been used to explain this dispersion difference, and it isfound that this is due to a paramagnetic enhancement of 1-2 water moleculesnearby ~10 Å from the spin center of the Mn4CaO5 cluster. The result showsthe reorientation of the molecular cluster is in μs time interval. Whencompare these two theories, the SL theory presented a better interpretationbecause parameters obtained from the SBM theory shows they didn’t fulfilthe presupposed perturbation criterion (the Kubo term). Paper III deals with the water dynamics in the restricted/confined spaces inthe zeolite samples (H-ZSM-5 and NH4-ZSM-5) and obtained by proton anddeuterium spin-lattice NMRD profiles. The results show that the spin-latticeNMRD can be used to characterize various zeolites. The temperature has aweak effect on the relaxation rate R1, but the change of different counter ionsmay change the hydration and the translational diffusion pores and givedifferent R1. Proton and fluoride NMRD profiles and MD simulations were both used tostudy the dynamics of BMIM[PF6] in paper IV. Results indicate the reorientation of the molecules are in the ns time regime, and the effectivecorrelation time obtained from 1H and 19F are the same. From the MDsimulation it is found the reorientation of [PF6]- ions is much faster (in ps)compare with BMIM+ ion which moves in the ns time range. With previous results, the FFC NMRD profiles are indeed very informativetools to study the molecular dynamics of complex systems. The MDsimulation can be used as a complementary method to obtain detailedinformation. By combine these two methods, it provide a more colorfulpicture in the study of protein hydration and liquid molecular dynamics.