Phenylalanine ammonia-lyases (PALs) deaminate L-phenylalanine to trans-cinnamic acid and ammonium and have widespread application in chemo-enzymatic synthesis, agriculture, and medicine. In particular, the PAL from Anabaena variabilis has garnered significant attention as the active ingredient in Pegvaliase®, the only FDA-approved drug treating classical phenylketonuria (PKU). Although an extensive body of literature exists on structure, substrate-specificity, and catalytic mechanism, protein-wide sequence determinants of function remain unknown, which limits the ability to rationally engineer these enzymes. Previously, we developed a high-throughput screen (HTS) for PAL, and here, we leverage it to create a detailed sequence-function landscape of PAL by performing deep mutational scanning (DMS). Our method revealed 79 hotspots that affected a positive change in enzyme fitness, many of which have not been reported previously. Using fitness values and structure-function analysis, we picked a subset of residues for comprehensive single- and multi-site saturation mutagenesis to improve the catalytic activity of PAL and identified combinations of mutations that led to improvement in reaction kinetics in cell-free and cellular contexts. To understand the mechanistic role of the most beneficial mutations, we performed QM/MM and MD and observed that different mutants confer improved catalytic activity via different mechanisms, including stabilizing first transition and intermediate states and improving substrate diffusion into the active site, and decreased product inhibition. Thus, this work provides a comprehensive sequence-function relationship for PAL, identifies positions that improve PAL activity when mutated and assess their mechanisms of action.