Hydroxylation of the aromatic amino acids phenylalanine, tyrosine and tryptophan is carried out by a family of non-heme iron and tetrahydrobiopterin (BH4) dependent enzymes, i.e. the aromatic amino acid hydroxylases (AAHs). The reactions catalyzed by these enzymes are important for biomedicine and their mutant forms in humans are associated with phenylketonuria (phenylalanine hydroxylase), Parkinson's disease and DOPA-responsive dystonia (tyrosine hydroxylase), and possibly neuropsychiatric and gastrointestinal disorders (tryptophan hydroxylase 1 and 2). We attempt to rationalize current knowledge about substrate and inhibitor specificity based on the three-dimensional structures of the enzymes and their complexes with substrates, cofactors and inhibitors. In addition, further insights on the selectivity and affinity determinants for ligand binding in the AAHs were obtained from molecular interaction field (MIF) analysis. We applied this computational structural approach to a rational analysis of structural differences at the active sites of the enzymes, a strategy that can help in the design of novel selective ligands for each AAH.