We use (2)H NMR spectroscopy to investigate the rotational motion of glycerol molecules in matrices provided by the connective tissue proteins elastin and collagen. Analyzing spin-lattice relaxation, line-shape properties, and stimulated-echo decays, we determine the rates and geometries of the motion as a function of temperature and composition. It is found that embedding glycerol in an elastin matrix leads to a mild slowdown of glycerol reorientation at low temperatures and glycerol concentrations, while the effect vanishes at ambient temperatures or high solvent content. Furthermore, it is observed that the nonexponential character of the rotational correlation functions is much more prominent in the elastin matrix than in the bulk liquid. Results from spin-lattice relaxation and line shape measurements indicate that, in the mixed systems, the strong nonexponentiality is in large part due to the existence of distributions of correlation times, which are broader on the long-time flank and, hence, more symmetric than in the neat system. Stimulated-echo analysis of slow glycerol dynamics reveals that, when elastin is added, the mechanism for the reorientation crosses over from small-angle jump dynamics to large-angle jump dynamics and the geometry of the motion changes from isotropic to anisotropic. The results are discussed against the background of present and previous findings for glycerol and water dynamics in various protein matrices and compared with observations for other dynamically highly asymmetric mixtures so as to ascertain in which way the viscous freezing of a fast component in the matrix of a slow component differs from the glassy slowdown in neat supercooled liquids.