Genotyping of the phenylalanine hydroxylase (PAH) gene offers a new tool for characterizing patients with phenylketonuria (PKU), refining the diagnosis and aiding in the prediction of the clinical outcome and in the implementation of a more adequate treatment. The primary goal of this work was the detailed study of the different allele combinations and the metabolic phenotypes in Spanish PKU patients in order to understand better the clinical heterogeneity of PAH deficiency in our population. The results show that the disease phenotype is a consequence of a combination of mutations at the PAH locus and this observation is valid throughout the spectrum of clinical and biochemical varieties found in Spanish PKU patients. A stronger correlation was found between the predicted residual activity, when known from previous in vitro studies of the mutant proteins, and the Phe tolerance than between the predicted residual activity and the inverse of Phe levels at diagnosis. The observed genotype-phenotype correlations and the available data on the in vitro residual activity of the mutant proteins has enabled the estimation of the severity of most of the mutations found in Spain. This study includes relevant data for clinicians and pediatricians adding to the present knowledge which relates allelic PAH genotypes to biological phenotypes.