An In Silico Analysis of the Binding Modes and Binding Affinities of Small Molecule Modulators of PDZ-Peptide Interactions

Abstract

Inhibitors of PDZ-peptide interactions have important implications in a variety of biological processes including treatment of cancer and Parkinson’s disease. Even though experimental studies have reported characterization of peptidomimetic inhibitors of PDZ-peptide interactions, the binding modes for most of them have not been characterized by structural studies. In this study we have attempted to understand the structural basis of the small molecule-PDZ interactions by in silico analysis of the binding modes and binding affinities of a set of 38 small molecules with known Ki or Kd values for PDZ2 and PDZ3 domains of PSD-95 protein. These two PDZ domains show differential selectivity for these compounds despite having a high degree of sequence similarity and almost identical peptide binding pockets. Optimum binding modes for these ligands for PDZ2 and PDZ3 domains were identified by using a novel combination of semi-flexible docking and explicit solvent molecular dynamics (MD) simulations. Analysis of the binding modes revealed most of the peptidomimectic ligands which had high Ki or Kd moved away from the peptide binding pocket, while ligands with high binding affinities remained in the peptide binding pocket. The differential specificities of the PDZ2 and PDZ3 domains primarily arise from differences in the conformation of the loop connecting βB and βC strands, because this loop interacts with the N-terminal chemical moieties of the ligands. We have also computed the MM/PBSA binding free energy values for these 38 compounds with both the PDZ domains from multiple 5 ns MD trajectories on each complex i.e. a total of 228 MD trajectories of 5 ns length each. Interestingly, computational binding free energies show good agreement with experimental binding free energies with a correlation coefficient of approximately 0.6. Thus our study demonstrates that combined use of docking and MD simulations can help in identification of potent inhibitors of PDZ-peptide complexes.

Andrographolide: A potent antituberculosis compound that targets Aminoglycoside 2′-N-acetyltransferase in Mycobacterium tuberculosis

“…It was earlier reported that when flexibilities are introduced into the protein and ligand, the computed binding energy values show better correlation with experimental results [34]. Thus, the MM-BSA binding free energy value, G bind = −12.27calculated for andrographolide within the binding site of 1M4I after MD simulation revealed that 1M4I could be a probable target protein for andrographolide.…”

Rapid Detection of Avian Influenza Virus by Fluorescent Diagnostic Assay using an Epitope-Derived Peptide

“…Docking was then carried out by using PEP-FOLD3 18 and I-TASSER 20 to model the structure of peptide and antigen. PyRx of AutoDock Vina docking software was used to predict binding energies and root-mean-square deviations (RMSD) 21 , 22 . PyMOL Molecular Graphics System was used to visualize interactive residues between peptides and antigen 23 .…”

Recent developments in molecular docking technology applied in food science: a review

“…Therefore, it is evident that the bond length or angle of the noncritical part of the ligand structure can be fixed (Morris et al, 2015). Semi-flexible docking is suitable for the docking of small molecules and macromolecules, such as proteins or nucleic acids and small ligand molecules (Tiwari & Mohanty, 2013;Jiang et al, 2016). It takes into account the influence of changes on the ligand structure and can be employed in a broader range of applications.…”