Decoupling of catalysis and transition state analog binding from mutations throughout a phosphatase revealed by high-throughput enzymology.

Using high-throughput microfluidic enzyme kinetics (HT-MEK), we measured over 9,000 inhibition curves detailing impacts of 1,004 single-site mutations throughout the alkaline phosphatase PafA on binding affinity for two transition state analogs (TSAs), vanadate and tungstate. As predicted by catalytic models invoking transition state complementary, mutations to active site and active-site-contacting residues had highly similar impacts on catalysis and TSA binding. Unexpectedly, most mutations to more distal residues that reduced catalysis had little or no impact on TSA binding and many even increased tungstate affinity. These disparate effects can be accounted for by a model in which distal mutations alter the enzyme's conformational landscape, increasing the occupancy of microstates that are catalytically less effective but better able to accommodate larger transition state analogs. In support of this ensemble model, glycine substitutions (rather than valine) were more likely to increase tungstate affinity (but not more likely to impact catalysis), presumably due to increased conformational flexibility that allows previously disfavored microstates to increase in occupancy. These results indicate that residues throughout an enzyme provide specificity for the transition state and discriminate against analogs that are larger only by tenths of an Ångström. Thus, engineering enzymes that rival the most powerful natural enzymes will likely require consideration of distal residues that shape the enzyme's conformational landscape and fine-tune active-site residues. Biologically, the evolution of extensive communication between the active site and remote residues to aid catalysis may have provided the foundation for allostery to make it a highly evolvable trait.

Energetically unfavorable protein angles: Exploration of a conserved dihedral angle in triosephosphate isomerase.

Over the past 3.5 billion years of evolution, enzymes have adopted a myriad of conformations to suit life on earth. However, torsional angles of proteins have settled into limited zones of energetically favorable dihedrals observed in Ramachandran plots. Areas outside said zones are believed to be disallowed to all amino acids, except glycine, due to steric hindrance. Triosephosphate isomerase (TIM), a homodimer with a catalytic rate approaching the diffusion limit, contains an active site lysine residue (K13) with dihedrals within the fourth quadrant (Φ = +51/Ψ = -143). Both the amino acid and the dihedral angles are conserved across all species of TIM and known crystal structures regardless of ligand. Only crystal structures of the engineered monomeric version (1MSS) show accepted ß-sheet dihedral values of Φ = -135/Ψ = +170 but experiments show a 1000-fold loss in activity. Based on these results, we hypothesized that adopting the unfavorable torsion angle for K13 contributes to catalysis. Using both, computational and experimental approaches, four residues that interact with K13 (N11, M14, E97, and Q64) were mutated to alanine. In silico molecular dynamics (MD) simulations were performed using 2JK2 unliganded human TIM as a starting structure. Ramachandran plots, containing K13 dihedral values reveal full or partial loss of disallowed zone angles. N11A showed no detectable catalytic activity and lost the unfavorable K13 dihedral angles across four separate force fields during simulation while all other mutants plus wild type retained activity and retained the conserved K13 dihedral angles.