The phenyl ring dynamics of [2H5]Phe4-labeled Leu5- and Met5-enkephalin molecules in crystals grown from four solvents were examined using solid state 2H NMR spectroscopy. 2H NMR powder patterns clearly indicated the presence of 180° flip motions about the Cβ−Cγ bond axis of the phenyl rings. Frequencies of the 180° flip motions were estimated to be 5.0 × 103, 3.0 × 104, and 2.4 × 106 Hz for Leu5-enkephalin crystallized from H2O, methanol/H2O, and N,N-dimethylformamide (DMF)/H2O, respectively and 1.0 × 104 Hz for Met5-enkephalin crystallized from ethanol/H2O at ambient temperature. The difference of the frequencies for the motion was attributed to the manner of their moleculer packing in the crystals as determined by X-ray diffraction. Because the correlation times determined from the 2H spin−lattice relaxation times (T1D values) were much shorter than those of the 180° flip motions, it was shown that the phenyl rings of these four crystals have small-amplitude librations. Therefore, we concluded that the T1D values were dominated by the librations, even for the ring deuterium. These motions became slower at lowered temperatures and caused the change of the peak intensities and increased quadrupole splittings which were observed in each 2H NMR spectra. Isotropic sharp signals due to naturally abundant solvent molecules were observed at the center of the 2H NMR spectra. The stepwise loss of the signal intensity was interpreted in terms of differential temperatures of freezing of motions of both bound water or organic solvent and the mixed solvent as the temperature was lowered, consistent with the buildup of solvent peaks in the 13C CP-MAS NMR spectra. It is suggested that there are a number of bound mobile solvent molecules in the crystals and the freezing of the solvents causes considerable changes in the conformations and dynamics of enkephalin molecules.