1
|
Wang J, Koirala K, Do HN, Miao Y. PepBinding: A Workflow for Predicting Peptide Binding Structures by Combining Peptide Docking and Peptide Gaussian Accelerated Molecular Dynamics Simulations. J Phys Chem B 2024; 128:7332-7340. [PMID: 39041172 DOI: 10.1021/acs.jpcb.4c02047] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 07/24/2024]
Abstract
Predicting protein-peptide interactions is crucial for understanding peptide binding processes and designing peptide drugs. However, traditional computational modeling approaches face challenges in accurately predicting peptide-protein binding structures due to the slow dynamics and high flexibility of the peptides. Here, we introduce a new workflow termed "PepBinding" for predicting peptide binding structures, which combines peptide docking, all-atom enhanced sampling simulations using the Peptide Gaussian accelerated Molecular Dynamics (Pep-GaMD) method, and structural clustering. PepBinding has been demonstrated on seven distinct model peptides. In peptide docking using HPEPDOCK, the peptide backbone root-mean-square deviations (RMSDs) of their bound conformations relative to X-ray structures ranged from 3.8 to 16.0 Å, corresponding to the medium to inaccurate quality models according to the Critical Assessment of PRediction of Interactions (CAPRI) criteria. The Pep-GaMD simulations performed for only 200 ns significantly improved the docking models, resulting in five medium and two acceptable quality models. Therefore, PepBinding is an efficient workflow for predicting peptide binding structures and is publicly available at https://github.com/MiaoLab20/PepBinding.
Collapse
Affiliation(s)
- Jinan Wang
- Computational Medicine Program and Department of Pharmacology, University of North Carolina - Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Kushal Koirala
- Computational Medicine Program and Department of Pharmacology, University of North Carolina - Chapel Hill, Chapel Hill, North Carolina 27599, United States
- Curriculum in Bioinformatics and Computational Biology, University of North Carolina - Chapel Hill, Chapel Hill, North Carolina 27599, United States
| | - Hung N Do
- Computational Biology Program, Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, United States
| | - Yinglong Miao
- Computational Medicine Program and Department of Pharmacology, University of North Carolina - Chapel Hill, Chapel Hill, North Carolina 27599, United States
| |
Collapse
|
2
|
Bayarsaikhan B, Zsidó BZ, Börzsei R, Hetényi C. Efficient Refinement of Complex Structures of Flexible Histone Peptides Using Post-Docking Molecular Dynamics Protocols. Int J Mol Sci 2024; 25:5945. [PMID: 38892133 PMCID: PMC11172440 DOI: 10.3390/ijms25115945] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/24/2024] [Revised: 05/26/2024] [Accepted: 05/27/2024] [Indexed: 06/21/2024] Open
Abstract
Histones are keys to many epigenetic events and their complexes have therapeutic and diagnostic importance. The determination of the structures of histone complexes is fundamental in the design of new drugs. Computational molecular docking is widely used for the prediction of target-ligand complexes. Large, linear peptides like the tail regions of histones are challenging ligands for docking due to their large conformational flexibility, extensive hydration, and weak interactions with the shallow binding pockets of their reader proteins. Thus, fast docking methods often fail to produce complex structures of such peptide ligands at a level appropriate for drug design. To address this challenge, and improve the structural quality of the docked complexes, post-docking refinement has been applied using various molecular dynamics (MD) approaches. However, a final consensus has not been reached on the desired MD refinement protocol. In this present study, MD refinement strategies were systematically explored on a set of problematic complexes of histone peptide ligands with relatively large errors in their docked geometries. Six protocols were compared that differ in their MD simulation parameters. In all cases, pre-MD hydration of the complex interface regions was applied to avoid the unwanted presence of empty cavities. The best-performing protocol achieved a median of 32% improvement over the docked structures in terms of the change in root mean squared deviations from the experimental references. The influence of structural factors and explicit hydration on the performance of post-docking MD refinements are also discussed to help with their implementation in future methods and applications.
Collapse
Affiliation(s)
- Bayartsetseg Bayarsaikhan
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, H-7624 Pécs, Hungary; (B.B.); (B.Z.Z.); (R.B.)
| | - Balázs Zoltán Zsidó
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, H-7624 Pécs, Hungary; (B.B.); (B.Z.Z.); (R.B.)
| | - Rita Börzsei
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, H-7624 Pécs, Hungary; (B.B.); (B.Z.Z.); (R.B.)
| | - Csaba Hetényi
- Pharmacoinformatics Unit, Department of Pharmacology and Pharmacotherapy, Medical School, University of Pécs, Szigeti út 12, H-7624 Pécs, Hungary; (B.B.); (B.Z.Z.); (R.B.)
- National Laboratory for Drug Research and Development, Magyar tudósok krt. 2, H-1117 Budapest, Hungary
| |
Collapse
|
3
|
Chen JN, Jiang F, Wu YD. Accurate Prediction for Protein-Peptide Binding Based on High-Temperature Molecular Dynamics Simulations. J Chem Theory Comput 2022; 18:6386-6395. [PMID: 36149394 DOI: 10.1021/acs.jctc.2c00743] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
The structural characterization of protein-peptide interactions is fundamental to elucidating biological processes and designing peptide drugs. Molecular dynamics (MD) simulations are extensively used to study biomolecular systems. However, simulating the protein-peptide binding process is usually quite expensive. Based on our previous studies, herein, we propose a simple and effective method to predict the binding site and pose of the peptide simultaneously using high-temperature (high-T) MD simulations with the RSFF2C force field. Thousands of binding events (nonspecific or specific) can be sampled during microseconds of high-T MD. From density-based clustering analysis, the structures of all of the 12 complexes (nine with linear peptides and three with cyclic peptides) can be successfully predicted with root-mean-square deviation (RMSD) < 2.5 Å. By directly simulating the process of the ligand binding onto the receptor, our method approaches experimental precision for the first time, significantly surpassing previous protein-peptide docking methods in terms of accuracy.
Collapse
Affiliation(s)
- Jia-Nan Chen
- Lab of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Fan Jiang
- Lab of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China
| | - Yun-Dong Wu
- Lab of Computational Chemistry and Drug Design, State Key Laboratory of Chemical Oncogenomics, Peking University Shenzhen Graduate School, Shenzhen 518055, China.,Shenzhen Bay Laboratory, Shenzhen 518132, China.,College of Chemistry and Molecular Engineering, Peking University, Beijing 100871, China
| |
Collapse
|
4
|
Wang J, Miao Y. Peptide Gaussian accelerated molecular dynamics (Pep-GaMD): Enhanced sampling and free energy and kinetics calculations of peptide binding. J Chem Phys 2021; 153:154109. [PMID: 33092378 DOI: 10.1063/5.0021399] [Citation(s) in RCA: 66] [Impact Index Per Article: 22.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/17/2022] Open
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.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, USA
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, Kansas 66047, USA
| |
Collapse
|
5
|
Wang J, Alekseenko A, Kozakov D, Miao Y. Improved Modeling of Peptide-Protein Binding Through Global Docking and Accelerated Molecular Dynamics Simulations. Front Mol Biosci 2019; 6:112. [PMID: 31737642 PMCID: PMC6835073 DOI: 10.3389/fmolb.2019.00112] [Citation(s) in RCA: 49] [Impact Index Per Article: 9.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2019] [Accepted: 10/09/2019] [Indexed: 01/31/2023] Open
Abstract
Peptides mediate up to 40% of known protein-protein interactions in higher eukaryotes and play a key role in cellular signaling, protein trafficking, immunology, and oncology. However, it is challenging to predict peptide-protein binding with conventional computational modeling approaches, due to slow dynamics and high peptide flexibility. Here, we present a prototype of the approach which combines global peptide docking using ClusPro PeptiDock and all-atom enhanced simulations using Gaussian accelerated molecular dynamics (GaMD). For three distinct model peptides, the lowest backbone root-mean-square deviations (RMSDs) of their bound conformations relative to X-ray structures obtained from PeptiDock were 3.3–4.8 Å, being medium quality predictions according to the Critical Assessment of PRediction of Interactions (CAPRI) criteria. GaMD simulations refined the peptide-protein complex structures with significantly reduced peptide backbone RMSDs of 0.6–2.7 Å, yielding two high quality (sub-angstrom) and one medium quality models. Furthermore, the GaMD simulations identified important low-energy conformational states and revealed the mechanism of peptide binding to the target proteins. Therefore, PeptiDock+GaMD is a promising approach for exploring peptide-protein interactions.
Collapse
Affiliation(s)
- Jinan Wang
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States
| | - Andrey Alekseenko
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, United States.,Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, NY, United States
| | - Dima Kozakov
- Laufer Center for Physical and Quantitative Biology, Stony Brook University, Stony Brook, NY, United States.,Department of Applied Mathematics and Statistics, Stony Brook University, Stony Brook, NY, United States
| | - Yinglong Miao
- Center for Computational Biology and Department of Molecular Biosciences, University of Kansas, Lawrence, KS, United States
| |
Collapse
|
6
|
Lapelosa M, Burrone O, Rocchia W. Specific Residue Interactions Regulate the Binding of Dengue Antigens to Broadly Neutralizing EDE Antibodies. ChemistryOpen 2018; 7:604-610. [PMID: 30151331 PMCID: PMC6099161 DOI: 10.1002/open.201800121] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2018] [Indexed: 11/06/2022] Open
Abstract
Antibodies binding to antigens present on the dengue virus (DENV) represent the main defense mechanism of the host organism against the pathogen. Among the antibodies elicited by DENV and that bind to DII of protein E, EDE1-C8 can bind all DENV serotypes. Our analysis reveals the key residues in this interaction as well as structurally conserved hydrogen bonds located at the binding interface. They stabilize the dengue antigen-antibody complex among the EDE1 group of antibodies (Abs). Combining structural alignments with molecular dynamics simulations in the EDE1 Abs, we identified the critical elements that provide a major energetic contribution to the association of antigens from protein E with Abs. We discuss possible molecular insights into the binding mechanism by using a surrogate molecular entity resembling the protein E that forms native salt bridges and hydrogen bonds, including inferences on the light of high-resolution crystal structures of dengue Fab complexes. Finally, the molecular determinants, the free energy profile, and the binding mechanism provide inspiration for potential strategies in protein engineering to design novel immunogens of protein E against DENV.
Collapse
Affiliation(s)
- Mauro Lapelosa
- CONCEPT LabDepartment of Drug Discovery and DevelopmentItalian Institute of TechnologyGenovaviaMorego 6016163Italy
| | - Oscar Burrone
- Department of Molecular ImmunologyICGEB-International Center for Genetic Engineering and BiotechnologyTriesteVia Padriciano 9934149Italy
| | - Walter Rocchia
- CONCEPT LabDepartment of Drug Discovery and DevelopmentItalian Institute of TechnologyGenovaviaMorego 6016163Italy
| |
Collapse
|