Ranking protein–protein docking results using steered molecular dynamics and potential of mean force calculations

Crystallization of protein–protein complexes can often be problematic and therefore computational structural models are often relied on. Such models are often generated using protein–protein docking algorithms, where one of the main challenges is selecting which of several thousand potential predictions represents the most near‐native complex. We have developed a novel technique that involves the use of steered molecular dynamics (sMD) and umbrella sampling to identify near‐native complexes among protein–protein docking predictions. Using this technique, we have found a strong correlation between our predictions and the interface RMSD (iRMSD) in ten diverse test systems. On two of the systems, we investigated if the prediction results could be further improved using potential of mean force calculations. We demonstrated that a near‐native (<2.0 Å iRMSD) structure could be identified in the top‐1 ranked position for both systems. © 2016 Wiley Periodicals, Inc.

[1]  Renos Savva,et al.  Nucleotide mimicry in the crystal structure of the uracil-DNA glycosylase–uracil glycosylase inhibitor protein complex , 1995, Nature Structural Biology.

[2]  L. Christophorou.  Science , 2018, Emerging Dynamics: Science, Energy, Society and Values.

[3]  Mieczyslaw Torchala,et al.  The scoring of poses in protein-protein docking: current capabilities and future directions , 2013, BMC Bioinformatics.

[4]  Markus A. Lill,et al.  WATsite: Hydration site prediction program with PyMOL interface , 2014, J. Comput. Chem..

[5]  Dietrich Suck,et al.  Structural basis of yeast aminoacyl-tRNA synthetase complex formation revealed by crystal structures of two binary sub-complexes , 2006, Nucleic acids research.

[6]  Yifeng D. Yang,et al.  Multi‐LZerD: Multiple protein docking for asymmetric complexes , 2012, Proteins.

[7]  O. Dym,et al.  Thermal stabilization of the protozoan Entamoeba histolytica alcohol dehydrogenase by a single proline substitution , 2008 .

[8]  G. Kachalova,et al.  Crystallographic studies of the interaction between the ferredoxin-NADP+ reductase and ferredoxin from the cyanobacterium Anabaena: looking for the elusive ferredoxin molecule. , 2000, Acta crystallographica. Section D, Biological crystallography.

[9]  J Otlewski,et al.  The refined 2.0 Å X‐ray crystal structure of the complex formed between bovine β‐trypsin and CMTI‐I, a trypsin inhibitor from squash seeds (Cucurbita maxima) Topological similarity of the squash seed inhibitors with the carboxypeptidase A inhibitor from potatoes , 1989, FEBS letters.

[10]  W. J. Allen,et al.  Steered molecular dynamics identifies critical residues of the Nodamura virus B2 suppressor of RNAi , 2014, Journal of Molecular Modeling.

[11]  S. Harrison,et al.  Structure of SARS Coronavirus Spike Receptor-Binding Domain Complexed with Receptor , 2005, Science.

[12]  BMC Bioinformatics , 2005 .

[13]  Liao Y. Chen,et al.  Computing Protein-Protein Association Affinity with Hybrid Steered Molecular Dynamics. , 2015, Journal of chemical theory and computation.

[14]  Z. Weng,et al.  Integrating atom‐based and residue‐based scoring functions for protein–protein docking , 2011, Protein science : a publication of the Protein Society.

[15]  S. Wodak,et al.  Docking and scoring protein complexes: CAPRI 3rd Edition , 2007, Proteins.

[16]  Z. Weng,et al.  ZDOCK: An initial‐stage protein‐docking algorithm , 2003, Proteins.

[17]  R. Rosenfeld.  Nature , 2009, Otolaryngology--head and neck surgery : official journal of American Academy of Otolaryngology-Head and Neck Surgery.

[18]  Anna Maria Almerico,et al.  Molecular dynamics, dynamic site mapping, and highthroughput virtual screening on leptin and the Ob receptor as anti-obesity target , 2014, Journal of Molecular Modeling.

[19]  R J Fletterick,et al.  Compromise and accommodation in ecotin, a dimeric macromolecular inhibitor of serine proteases. , 2000, Journal of molecular biology.

[20]  Carles Pons,et al.  Present and future challenges and limitations in protein–protein docking , 2010, Proteins.

[21]  Zhiping Weng,et al.  ZRANK: Reranking protein docking predictions with an optimized energy function , 2007, Proteins.

[22]  J. Deisenhofer,et al.  A structural basis of the interactions between leucine-rich repeats and protein ligands , 1995, Nature.

[23]  T. Bhat,et al.  Bound water molecules and conformational stabilization help mediate an antigen-antibody association. , 1994, Proceedings of the National Academy of Sciences of the United States of America.

[24]  Robert Huber,et al.  Catalytic Domain Structures of MT-SP1/Matriptase, a Matrix-degrading Transmembrane Serine Proteinase* , 2002, The Journal of Biological Chemistry.

[25]  Kalle Gehring,et al.  Structural basis for ubiquitin-mediated dimerization and activation of the ubiquitin protein ligase Cbl-b. , 2007, Molecular cell.

[26]  Marc Ostermeier,et al.  Characterization of Monobody Scaffold Interactions with Ligand via Force Spectroscopy and Steered Molecular Dynamics , 2015, Scientific Reports.