AutoDock CrankPep: combining folding and docking to predict protein-peptide complexes

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

MDM2 and MDMX are potential targets for p53-dependent cancer therapy. Peptides are key in cellular immunology and oncology, and cyclic peptides generally have higher half-life than their linear counterparts. However, prediction of cyclic peptide-protein binding is challenging with normal molecular simulation approaches because of high peptide flexibility. Here, we used global peptide docking, normal molecular dynamics, Gaussian accelerated molecular dynamics (GaMD), two-dimensional (2D) potential of mean force (PMF) profiles, and solvated interaction energy (SIE) techniques to investigate the interactions of MDM2/MDMX with three N-to-C-terminal cyclic peptide-based inhibitors. We determined the possible cyclic peptide-MDM2/MDMX complex structures via 2D PMF profiles and SIE calculations. Our findings increase the accuracy of peptide-protein structural prediction, which may facilitate cyclic peptide drug design. Advancements in the computational methods and computing power may further aid in addressing the challenges in cyclic peptide drug design. Communicated by Ramaswamy H. Sarma.

GM-DockZn: a geometry matching-based docking algorithm for zinc proteins

Motivation

Molecular docking is a widely used technique for large-scale virtual screening of the interactions between small-molecule ligands and their target proteins. However, docking methods often perform poorly for metalloproteins due to additional complexity from the three-way interactions among amino-acid residues, metal ions and ligands. This is a significant problem because zinc proteins alone comprise about 10% of all available protein structures in the protein databank. Here, we developed GM-DockZn that is dedicated for ligand docking to zinc proteins. Unlike the existing docking methods developed specifically for zinc proteins, GM-DockZn samples ligand conformations directly using a geometric grid around the ideal zinc-coordination positions of seven discovered coordination motifs, which were found from the survey of known zinc proteins complexed with a single ligand.

Results

GM-DockZn has the best performance in sampling near-native poses with correct coordination atoms and numbers within the top 50 and top 10 predictions when compared to several state-of-the-art techniques. This is true not only for a non-redundant dataset of zinc proteins but also for a homolog set of different ligand and zinc-coordination systems for the same zinc proteins. Similar superior performance of GM-DockZn for near-native-pose sampling was also observed for docking to apo-structures and cross-docking between different ligand complex structures of the same protein. The highest success rate for sampling nearest near-native poses within top 5 and top 1 was achieved by combining GM-DockZn for conformational sampling with GOLD for ranking. The proposed geometry-based sampling technique will be useful for ligand docking to other metalloproteins.

Availability and implementation

GM-DockZn is freely available at www.qmclab.com/ for academic users.

Supplementary information

Supplementary data are available at Bioinformatics online.

Targeting Tumors Using Peptides

To penetrate solid tumors, low molecular weight (Mw < 10 KDa) compounds have an edge over antibodies: their higher penetration because of their small size. Because of the dense stroma and high interstitial fluid pressure of solid tumors, the penetration of higher Mw compounds is unfavored and being small thus becomes an advantage. This review covers a wide range of peptidic ligands—linear, cyclic, macrocyclic and cyclotidic peptides—to target tumors: We describe the main tools to identify peptides experimentally, such as phage display, and the possible chemical modifications to enhance the properties of the identified peptides. We also review in silico identification of peptides and the most salient non-peptidic ligands in clinical stages. We later focus the attention on the current validated ligands available to target different tumor compartments: blood vessels, extracelullar matrix, and tumor associated macrophages. The clinical advances and failures of these ligands and their therapeutic conjugates will be discussed. We aim to present the reader with the state-of-the-art in targeting tumors, by using low Mw molecules, and the tools to identify new ligands.

Highly Flexible Ligand Docking: Benchmarking of the DockThor Program on the LEADS-PEP Protein-Peptide Data Set

Protein-peptide interactions play a crucial role in many cellular and biological functions, which justify the increasing interest in the development of peptide-based drugs. However, predicting experimental binding modes and affinities in protein-peptide docking remains a great challenge for most docking programs due to some particularities of this class of ligands, such as the high degree of flexibility. In this paper, we present the performance of the DockThor program on the LEADS-PEP data set, a benchmarking set composed of 53 diverse protein-peptide complexes with peptides ranging from 3 to 12 residues and with up to 51 rotatable bonds. The DockThor performance for pose prediction on redocking studies was compared with some state-of-the-art docking programs that were also evaluated on the LEADS-PEP data set, AutoDock, AutoDock Vina, Surflex, GOLD, Glide, rDock, and DINC, as well as with the task-specific docking protocol HPepDock. Our results indicate that DockThor could dock 40% of the cases with an overall backbone RMSD below 2.5 Å when the top-scored docking pose was considered, exhibiting similar results to Glide and outperforming other protein-ligand docking programs, whereas rDock and HPepDock achieved superior results. Assessing the docking poses closest to the crystal structure (i.e., best-RMSD pose), DockThor achieved a success rate of 60% in pose prediction. Due to the great overall performance of handling peptidic compounds, the DockThor program can be considered as suitable for docking highly flexible and challenging ligands, with up to 40 rotatable bonds. DockThor is freely available as a virtual screening Web server at https://www.dockthor.lncc.br/ .

AutoGridFR: Improvements on AutoDock Affinity Maps and Associated Software Tools

Precomputed affinity maps are used by AutoDock to efficiently describe rigid biomolecules called receptors in automated docking. These maps greatly speed up the docking process and allow users to experiment with the forcefield. Here, we present AutoGridFR (AGFR): a software tool facilitating the calculation of these maps. We describe a new version of the AutoSite algorithm that improves the description of binding pockets automatically detected on receptors, and an algorithm for adding affinity gradients which help search methods optimize solution using fewer evaluations of the scoring functions. AGFR supports the calculation of maps for various advanced docking techniques such as covalent docking, hydrated docking, and docking with flexible receptor sidechains. Maps are stored in a single file along with metadata supporting data provenance, reproducibility, and facilitating their management. Finally, maps can be calculated from the command line or through a modern graphical user interface which also supports their visualization. © 2019 Wiley Periodicals, Inc.

Modeling Electrostatic Force in Protein-Protein Recognition

Electrostatic interactions are important for understanding molecular interactions, since they are long-range interactions and can guide binding partners to their correct binding positions. To investigate the role of electrostatic forces in molecular recognition, we calculated electrostatic forces between binding partners separated at various distances. The investigation was done on a large set of 275 protein complexes using recently developed DelPhiForce tool and in parallel, evaluating the total electrostatic force via electrostatic association energy. To accomplish the goal, we developed a method to find an appropriate direction to move one chain of protein complex away from its bound position and then calculate the corresponding electrostatic force as a function of separation distance. It is demonstrated that at large distances between the partners, the electrostatic force (magnitude and direction) is consistent among the protocols used and the main factors contributing to it are the net charge of the partners and their interfaces. However, at short distances, where partners form specific pair-wise interactions or de-solvation penalty becomes significant, the outcome depends on the precise balance of these factors. Based on the electrostatic force profile (force as a function of distance), we group the cases into four distinctive categories, among which the most intriguing is the case termed "soft landing." In this case, the electrostatic force at large distances is favorable assisting the partners to come together, while at short distance it opposes binding, and thus slows down the approach of the partners toward their physical binding.

Docking Flexible Cyclic Peptides with AutoDock CrankPep

While a new therapeutic cyclic peptide is approved nearly every year, docking large macrocycles has remained challenging. Here, we present a new version of our peptide docking software AutoDock CrankPep (ADCP), extended to dock peptides cyclized through their backbone and/or side chain disulfide bonds. We show that within the top 10 solutions, ADCP identifies the proper interactions for 71% of a data set of 38 complexes, thus making it a useful tool for rational peptide-based drug design.