Computational Study of Conformation Dynamics and Allostery in WW Protein Domains

Proteins continually and naturally incur evolutionary selection through mutagenesis that optimizes their fitness, which is primarily determined by their function. It is known that allosteric regulation alters a protein's conformational dynamics leading to functional changes. We have computationally introduced a mutation at a predicted regulatory site of a short, 46 residue-long, protein interaction module composed of a WW domain and corresponding polyproline ligand (PDB id: 1k9r). The dynamic flexibility index (DFI) was computed for the binding site of the wild type and mutant WW domains to quantify the mutations effect on the rigidity of the binding pocket. DFI is used as a metric to quantify the resilience of a given position to perturbation along the chain. Using steered molecular dynamics (SMD), we also measure the effect of the point mutation on allosteric regulation by approximating the binding free energy of the system calculated using Jarzynski's Equality. Calculation of the DFI shows that the overall flexibility of the protein complex increases as a result of the distal point mutation. Total change in DFI percentile of the binding site showed a 0.067 increase suggesting an allosteric, loss of function mutation. Furthermore, we see that the change in the binding free energy is greater for that of the mutated complex supporting the idea that an increase in flexibility is correlated to a decrease in proteinlig and binding affinity. We show that sequence mutation of an allosteric site affects the mechanical stability and functionality of the binding pocket.