The three aromatic amino acid hydroxylases (phenylalanine, tyrosine, and tryptophan hydroxylase) and nitric oxide synthase (NOS) all utilize (6R)-l-erythro-5,6,7,8-tetrahydrobiopterin (BH(4)) as cofactor. The pterin binding site in the three hydroxylases is well conserved and different from the binding site in NOS. The structures of phenylalanine hydroxylase (PAH) and of NOS in complex with BH(4) are still the only crystal structures available for the reduced cofactor-enzyme complexes. We have studied the enzyme-bound and free conformations of BH(4) by NMR spectroscopy and molecular docking into the active site of the three hydroxylases, using endothelial NOS as a comparative probe. We have found that the dihydroxypropyl side chain of BH(4) adopts different conformations depending on which hydroxylase it interacts with. All the bound conformations are different from that of BH(4) free in solution at neutral pH. The different bound conformations appear to result from specific interactions with nonconserved amino acids at the BH(4) binding sites of the hydroxylases, notably the stretch 248-251 (numeration in PAH) and the residue corresponding to Ala322 in PAH, i.e., Ser in TH and Ala in TPH. On the basis of analysis of molecular interaction fields, we discuss the selectivity determinants for each hydroxylase and explain the high-affinity inhibitory effect of 7-tetrahydrobiopterin specifically for PAH.