A flexible docking scheme to explore the binding selectivity of PDZ domains

Z. Nevin Gerek, Sefika Ozkan

Research output: Contribution to journalArticle

41 Citations (Scopus)

Abstract

Modeling of protein binding site flexibility in molecular docking is still a challenging problem due to the large conformational space that needs sampling. Here, we propose a flexible receptor docking scheme: A dihedral restrained replica exchange molecular dynamics (REMD), where we incorporate the normal modes obtained by the Elastic Network Model (ENM) as dihedral restraints to speed up the search towards correct binding site conformations. To our knowledge, this is the first approach that uses ENM modes to bias REMD simulations towards binding induced fluctuations in docking studies. In our docking scheme, we first obtain the deformed structures of the unbound protein as initial conformations by moving along the binding fluctuation mode, and perform REMD using the ENM modes as dihedral restraints. Then, we generate an ensemble of multiple receptor conformations (MRCs) by clustering the lowest replica trajectory. Using ROSETTALIGAND, we dock ligands to the clustered conformations to predict the binding pose and affinity. We apply this method to postsynaptic density-95/Dlg/ZO-1 (PDZ) domains; whose dynamics govern their binding specificity. Our approach produces the lowest energy bound complexes with an average ligand root mean square deviation of 0.36 Å. We further test our method on (i) homologs and (ii) mutant structures of PDZ where mutations alter the binding selectivity. In both cases, our approach succeeds to predict the correct pose and the affinity of binding peptides. Overall, with this approach, we generate an ensemble of MRCs that leads to predict the binding poses and specificities of a protein complex accurately. Published by Wiley-Blackwell.

Original languageEnglish (US)
Pages (from-to)914-928
Number of pages15
JournalProtein Science
Volume19
Issue number5
DOIs
StatePublished - May 2010

Fingerprint

Post-Synaptic Density
Molecular Dynamics Simulation
Conformations
Molecular dynamics
Binding Sites
Ligands
Protein Binding
Cluster Analysis
Proteins
Docks
Peptides
Mutation
Trajectories
Sampling
Computer simulation

Keywords

  • Docking
  • PDZ domain
  • Protein flexibility
  • RosettaLigand
  • Selectivity

ASJC Scopus subject areas

  • Biochemistry
  • Molecular Biology

Cite this

A flexible docking scheme to explore the binding selectivity of PDZ domains. / Nevin Gerek, Z.; Ozkan, Sefika.

In: Protein Science, Vol. 19, No. 5, 05.2010, p. 914-928.

Research output: Contribution to journalArticle

@article{58612682a4e64734b6f06f4a7b2f2bb5,
title = "A flexible docking scheme to explore the binding selectivity of PDZ domains",
abstract = "Modeling of protein binding site flexibility in molecular docking is still a challenging problem due to the large conformational space that needs sampling. Here, we propose a flexible receptor docking scheme: A dihedral restrained replica exchange molecular dynamics (REMD), where we incorporate the normal modes obtained by the Elastic Network Model (ENM) as dihedral restraints to speed up the search towards correct binding site conformations. To our knowledge, this is the first approach that uses ENM modes to bias REMD simulations towards binding induced fluctuations in docking studies. In our docking scheme, we first obtain the deformed structures of the unbound protein as initial conformations by moving along the binding fluctuation mode, and perform REMD using the ENM modes as dihedral restraints. Then, we generate an ensemble of multiple receptor conformations (MRCs) by clustering the lowest replica trajectory. Using ROSETTALIGAND, we dock ligands to the clustered conformations to predict the binding pose and affinity. We apply this method to postsynaptic density-95/Dlg/ZO-1 (PDZ) domains; whose dynamics govern their binding specificity. Our approach produces the lowest energy bound complexes with an average ligand root mean square deviation of 0.36 {\AA}. We further test our method on (i) homologs and (ii) mutant structures of PDZ where mutations alter the binding selectivity. In both cases, our approach succeeds to predict the correct pose and the affinity of binding peptides. Overall, with this approach, we generate an ensemble of MRCs that leads to predict the binding poses and specificities of a protein complex accurately. Published by Wiley-Blackwell.",
keywords = "Docking, PDZ domain, Protein flexibility, RosettaLigand, Selectivity",
author = "{Nevin Gerek}, Z. and Sefika Ozkan",
year = "2010",
month = "5",
doi = "10.1002/pro.366",
language = "English (US)",
volume = "19",
pages = "914--928",
journal = "Protein Science",
issn = "0961-8368",
publisher = "Cold Spring Harbor Laboratory Press",
number = "5",

}

TY - JOUR

T1 - A flexible docking scheme to explore the binding selectivity of PDZ domains

AU - Nevin Gerek, Z.

AU - Ozkan, Sefika

PY - 2010/5

Y1 - 2010/5

N2 - Modeling of protein binding site flexibility in molecular docking is still a challenging problem due to the large conformational space that needs sampling. Here, we propose a flexible receptor docking scheme: A dihedral restrained replica exchange molecular dynamics (REMD), where we incorporate the normal modes obtained by the Elastic Network Model (ENM) as dihedral restraints to speed up the search towards correct binding site conformations. To our knowledge, this is the first approach that uses ENM modes to bias REMD simulations towards binding induced fluctuations in docking studies. In our docking scheme, we first obtain the deformed structures of the unbound protein as initial conformations by moving along the binding fluctuation mode, and perform REMD using the ENM modes as dihedral restraints. Then, we generate an ensemble of multiple receptor conformations (MRCs) by clustering the lowest replica trajectory. Using ROSETTALIGAND, we dock ligands to the clustered conformations to predict the binding pose and affinity. We apply this method to postsynaptic density-95/Dlg/ZO-1 (PDZ) domains; whose dynamics govern their binding specificity. Our approach produces the lowest energy bound complexes with an average ligand root mean square deviation of 0.36 Å. We further test our method on (i) homologs and (ii) mutant structures of PDZ where mutations alter the binding selectivity. In both cases, our approach succeeds to predict the correct pose and the affinity of binding peptides. Overall, with this approach, we generate an ensemble of MRCs that leads to predict the binding poses and specificities of a protein complex accurately. Published by Wiley-Blackwell.

AB - Modeling of protein binding site flexibility in molecular docking is still a challenging problem due to the large conformational space that needs sampling. Here, we propose a flexible receptor docking scheme: A dihedral restrained replica exchange molecular dynamics (REMD), where we incorporate the normal modes obtained by the Elastic Network Model (ENM) as dihedral restraints to speed up the search towards correct binding site conformations. To our knowledge, this is the first approach that uses ENM modes to bias REMD simulations towards binding induced fluctuations in docking studies. In our docking scheme, we first obtain the deformed structures of the unbound protein as initial conformations by moving along the binding fluctuation mode, and perform REMD using the ENM modes as dihedral restraints. Then, we generate an ensemble of multiple receptor conformations (MRCs) by clustering the lowest replica trajectory. Using ROSETTALIGAND, we dock ligands to the clustered conformations to predict the binding pose and affinity. We apply this method to postsynaptic density-95/Dlg/ZO-1 (PDZ) domains; whose dynamics govern their binding specificity. Our approach produces the lowest energy bound complexes with an average ligand root mean square deviation of 0.36 Å. We further test our method on (i) homologs and (ii) mutant structures of PDZ where mutations alter the binding selectivity. In both cases, our approach succeeds to predict the correct pose and the affinity of binding peptides. Overall, with this approach, we generate an ensemble of MRCs that leads to predict the binding poses and specificities of a protein complex accurately. Published by Wiley-Blackwell.

KW - Docking

KW - PDZ domain

KW - Protein flexibility

KW - RosettaLigand

KW - Selectivity

UR - http://www.scopus.com/inward/record.url?scp=77951588413&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=77951588413&partnerID=8YFLogxK

U2 - 10.1002/pro.366

DO - 10.1002/pro.366

M3 - Article

C2 - 20196074

AN - SCOPUS:77951588413

VL - 19

SP - 914

EP - 928

JO - Protein Science

JF - Protein Science

SN - 0961-8368

IS - 5

ER -