Abstract:Abstract Peptides mediate up to 40% of known protein-protein interactions in higher eukaryotes and play an important role in cellular signaling. However, it is challenging to simulate both binding and unbinding of peptides and calculate peptide binding free energies through conventional molecular dynamics, due to long biological timescales and extremely high flexibility of the peptides. Based on the Gaussian accelerated molecular dynamics (GaMD) enhanced sampling technique, we have developed a new computational method “Pep-GaMD”, which selectively boosts essential potential energy of the peptide in order to effectively model its high flexibility. In addition, another boost potential is applied to the remaining potential energy of the entire system in a dual-boost algorithm. Pep-GaMD has been demonstrated on binding of three model peptides to the SH3 domains. Independent 1 μs dual-boost Pep-GaMD simulations have captured repetitive peptide dissociation and binding events, which enable us to calculate peptide binding thermodynamics and kinetics. The calculated binding free energies and kinetic rate constants agreed very well with available experimental data. Furthermore, the all-atom Pep-GaMD simulations have provided important insights into the mechanism of peptide binding to proteins that involves long-range electrostatic interactions and mainly conformational selection. In summary, Pep-GaMD provides a highly efficient, easy-to-use approach for unconstrained enhanced sampling and calculations of peptide binding free energies and kinetics. Significance Statement We have developed a new computational method “Pep-GaMD” for enhanced sampling of peptide-protein interactions based on the Gaussian accelerated molecular dynamics (GaMD) technique. Pep-GaMD works by selectively boosting the essential potential energy of the peptide to effectively model its high flexibility. In addition, another boost potential can be applied to the remaining potential energy of the entire system in a dual-boost algorithm. Pep-GaMD has been demonstrated on binding of three model peptides to the SH3 domains. Dual-boost Pep-GaMD has captured repetitive peptide dissociation and binding events within significantly shorter simulation time (microsecond) than conventional molecular dynamics. Compared with previous enhanced sampling methods, Pep-GaMD is easier to use and more efficient for unconstrained enhanced sampling of peptide binding and unbinding, which provides a novel physics-based approach to calculating peptide binding free energies and kinetics.

Improved Modeling of Peptide-Protein Binding Through Global Docking and Accelerated Molecular Dynamics Simulations

PepBinding: A Workflow for Predicting Peptide Binding Structures by Combining Peptide Docking and Peptide Gaussian Accelerated Molecular Dynamics Simulations

Peptide Gaussian accelerated molecular dynamics (Pep-GaMD): Enhanced sampling and free energy and kinetics calculations of peptide binding

Improving Protein-Peptide Docking Results Via Pose-Clustering and Rescoring with a Combined Knowledge-Based and MM-GBSA Scoring Function.

Molecular Dynamics Scoring of Protein–Peptide Models Derived from Coarse-Grained Docking

A Protein-Peptide Docking Program With Modeling Receptor Flexible Areas

MDockPeP: An ab‐initio protein–peptide docking server

Comparison of Two Docking Methods for Peptide-Protein Interactions.

Accurate Prediction for Protein-Peptide Binding Based on High-Temperature Molecular Dynamics Simulations

Modelling peptide–protein complexes: docking, simulations and machine learning

Peptide docking and structure-based characterization of peptide binding: from knowledge to know-how

Improved Docking of Polypeptides with Glide

Towards protein-protein docking with significant structural changes using CABS-dock

Computational modeling of cyclic peptide inhibitor–MDM2/MDMX binding through global docking and Gaussian accelerated molecular dynamics simulations

Modeling of protein-peptide interactions using the CABS-dock web server for binding site search and flexible docking

Protein-Protein Docking Using A Brownian Dynamics Simulations Approach

Hierarchical Flexible Peptide Docking by Conformer Generation and Ensemble Docking of Peptides.

Biologically enhanced sampling geometric docking and backbone flexibility treatment with multiconformational superposition

Molecular docking with conformer-dependent charges

Comprehensive Evaluation of 10 Docking Programs on a Diverse Set of Protein-Cyclic Peptide Complexes.

Docking Cyclic Peptides Formed by a Disulfide Bond Through a Hierarchical Strategy