Decoding the deactivation mechanism of R192W mutation of ZAP-70 using molecular dynamics simulations and binding free energy calculations

Context ZAP-70 (zeta-chain-associated protein of 70 kDa), serving as a critical regulator for T cell antigen receptor signaling, represents an attractive therapeutic target for autoimmunity disease. How the mechanistical mechanism of ZAP-70 to a human autoimmune syndrome-associated R192W mutation remains unclear. The results indicated that the R192W mutation of ZAP-70 clearly affected the conformational flexibility of the N-terminal ITAM-Y2P. Structural analysis unveiled that the R192W mutation of ZAP-70 caused the exposure of the N-terminal ITAM-Y2P to the solvent. MM-GBSA binding free energy calculations exhibited that the R192W mutation decreased the binding affinity of ITAM-Y2P to the ZAP-70 mutant. Residue-based free energy decomposition further revealed that the protein-peptide interaction networks involving electrostatic interactions provide significant contributions for complex formation. The energy unfavorable residues include Arg43, Arg192, Tyr240, and Lys244 from ZAP-70 and Asn301, Leu303, pY304, and pY315 from ITAM-Y2P in the R192W mutant. Our obtained results may help the understanding of the deactivation mechanism of ZAP-70 induced by the R192W mutation. Methods In the work, multiple replica molecular dynamics simulations and molecular mechanics-generalized Born surface area (MM-GBSA) method were performed to reveal the doubly phosphorylated ITAMs (ITAM-Y2P)-mediated deactivation mechanism of ZAP-70 induced by the R192W mutation.