H-1 hyperfine (hf) coupling constants of semiquinone radical anions of 1,4-naphthoquinone, 2-methyl-1,4-naphthoquinone, and 2-methyl-3-phytyl-1,4-naphthoquinone in frozen alcoholic solutions were measured using pulse Q-band electron nuclear double resonance spectroscopy. The resolved signals of the quinone protons as well as from hydrogen bond and solvent shell protons were analyzed and assigned. Both in-plane and out-of-plane hydrogen bonding with respect to the pi-plane of the radical is observed. Interactions with nonexchangeable protons from the surrounding matrix are detected and assigned to solvent protons above and below the quinone plane. Density functional theory was used to calculate spin Hamiltonian parameters of the radical anions. Solvent molecules of the first solvent shell that provide hydrogen bonds to the quinones were included in the geometry optimization. The conductor-like screening model was employed to introduce additional effects of the solvent cage. From a comparison of the experimental and calculated hf tensors it is concluded that four solvent molecules are coordinated via hydrogen bonds to the quinone oxygens. For all radicals very good agreement between experimental and calculated data is observed. The influence of different substituents on the spin density distribution and hydrogen bond geometries is discussed.