Application of per-Residue Energy Decomposition to Design Peptide Inhibitors of PSD95 GK Domain

Abstract

Specific interaction between the postsynaptic density protein 95 (PSD95) and synapse-associated protein 90/postsynaptic density 95–associated protein (SAPAP) is crucial for excitatory synaptic development and plasticity. Designing inhibitors that target the guanylate kinase (GK) domain of PSD95, which is responsible for the interaction, is a promising manipulation tool for the investigation of the function of PSD95 GK and the etiology of its related psychiatric disorders. Herein, we designed new peptide inhibitors of PSD95 GK/SAPAP with higher binding affinity by using molecular dynamics simulations. First, the interactions between PSD95 GK and their reported phosphorylated and unphosphorylated peptides were explored by molecular dynamics simulations. Besides the hydrogen bonding interactions mediated by the phospho-serine (p-Ser) or corresponding phosphomimic residue Asp/Glu, the hydrophobic interactions from the other amino acids also contribute to the PSD95 GK/SAPAP interaction. As an unphosphorylated synthetic peptide with moderate binding affinity and relatively lower molecular weight, the QSF inhibitory peptide was selected for further modification. Based on per-residue energy decomposition results of the PSD95 GK/QSF complex, ten peptides were designed to enhance the binding interactions, especially the hydrophobic interactions. The top-ranked five peptides with lower binding energy were eventually synthesized. The binding affinities of the synthesized peptides were determined using fluorescence polarization (FP) assay. As expected, all peptides have higher binding affinity than the QSF peptide (Ki = 5.64 ± 0.51 μM). Among them, F10W was the most potent inhibitor (Ki = 0.75 ± 0.25 μM), suggesting that enhancement of the hydrophobic interactions is an important strategy for the design of new inhibitory peptides targeting PSD95 GK.

Simplified Derivatives of Tetrandrine as Potent and Specific P-gp Inhibitors to Reverse Multidrug Resistance in Cancer Chemotherapy

“…The binding pocket of P-gp was generated by multi-channel surface mode, which was in accordance with the binding site of Taxol as well as zosuquidar. The main protocols or parameters of docking were addressed in our previous publications. − The key amino acids among the binding interactions were analyzed by ligplot v2.2.5 and generated using PyMOL v2.52…”

Revealing the Interaction Mechanism between Mycobacterium tuberculosis GyrB and Novobiocin, SPR719 through Binding Thermodynamics and Dissociation Kinetics Analysis

“…The Molecular Mechanics-Generalized Born Surface Area (MM-GBSA) method strikes a balance between computational efficiency and accuracy [59,60] and is widely used to predict the binding affinity between proteins and ligands [61][62][63][64]. Additionally, MM-GBSA facilitates the analysis of energy contributions from each residue through binding free energy decomposition, identifying key interactions during the binding process [65,66]. In MM-GBSA, the calculation of the free energy for ligand-receptor binding is conducted as follows:…”

Identification of small molecule inhibitors of CXCR4 – an important drug target in renal fibrosis

“…The MMGBSA of the three complexes were done to nd the ligand binding a nities [42]. In Table 4, The ligand binding a nity of 4990 was better than other two ligands.…”