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.

Combining Enhanced Sampling with Experiment Directed Simulation of the GYG peptide

Enhanced Conformational Sampling using Replica Exchange with Collective-Variable Tempering

Replica Exchange with Solute Tempering for Protein Conformational Sampling

Enhanced Sampling in Molecular Dynamics Using Metadynamics, Replica-Exchange, and Temperature-Acceleration

Enhanced sampling methods for molecular dynamics simulations

Enhanced sampling techniques in molecular dynamics simulations of biological systems

Boosting the Conformational Sampling by Combining Replica Exchange with Solute Tempering and Well-Sliced Metadynamics

Enhanced Sampling Method in Molecular Simulations

Simulated Solute Tempering 2

Parallel Implementation of the Replica Exchange Molecular Dynamics Algorithm on Blue Gene/L

Enhanced Sampling with Sub-optimal Collective Variables: Reconciling Accuracy and Convergence Speed

Replica exchange molecular dynamics method for protein folding simulation.

Replica Exchange Molecular Dynamics: A Practical Application Protocol with Solutions to Common Problems and a Peptide Aggregation and Self-Assembly Example.

Efficient Sampling over Rough Energy Landscapes with High Barriers: A Combination of Metadynamics with Integrated Tempering Sampling.

Enhanced Sampling Method with Coarse Graining of Conformational Space.

Enhancing torsional sampling using fully adaptive simulated tempering

Enhanced sampling techniques in biomolecular simulations

Enhanced Sampling Simulations of RNA-peptide Binding using Deep Learning Collective Variables

An interoperable implementation of collective‐variable based enhanced sampling methods in extended phase space within the OpenMM package

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

Implementing Dimer Metadynamics using GROMACS