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

MicroProteinDB: A database to provide knowledge on sequences, structures and function of ncRNA-derived microproteins

Omics-based technologies have revolutionized our comprehension of microproteins encoded by ncRNAs, revealing their abundant presence and pivotal roles within complex functional landscapes. Here, we developed MicroProteinDB (http://bio-bigdata.hrbmu.edu.cn/MicroProteinDB), which offers and visualizes the extensive knowledge to aid retrieval and analysis of computationally predicted and experimentally validated microproteins originating from various ncRNA types. Employing prediction algorithms grounded in diverse deep learning approaches, MicroProteinDB comprehensively documents the fundamental physicochemical properties, secondary and tertiary structures, interactions with functional proteins, family domains, and inter-species conservation of microproteins. With five major analytical modules, it will serve as a valuable knowledge for investigating ncRNA-derived microproteins.

Bioinformatics tools for the study of bioactive peptides from vegetal sources: evolution and future perspectives

Bioactive peptides from vegetal sources have been shown to have functional properties as anti-inflammatory, antioxidant, antihypertensive or antidiabetic capacity. For this reason, they have been proposed as an interesting and promising alternative to improve human health. In recent years, the numerous advances in the bioinformatics field for in silico prediction have speeded up the discovery of bioactive peptides, also reducing the associated costs when using an integrated approach between the classical and bioinformatics discovery. This review aims to provide an overview of the evolution, limitations and latest advances in the field of bioinformatics and computational tools, and specifically make a critical and comprehensive insight into computational techniques used to study the mechanism of interaction that allows the explanation of plant bioactive peptide functionality. In particular, molecular docking is considered key to explain the different functionalities that have been previously identified. The assumptions to simplify such a high complex environment implies a degree of uncertainty that can only be guaranteed and validated by in vitro or in vivo studies, however, the combination of databases, software and bioinformatics applications with the classical approach has become a promising procedure for the study of bioactive peptides.

Computational studies and synthesis of 131iodine-labeled nocardiotide A analogs as a peptide-based theragnostic radiopharmaceutical ligand for cancer targeting SSTR2

Radiolabeled peptides belong to a highly specific group of radiotracers used in oncology, particularly for diagnostics and cancer therapy. With the notable advantages of high binding affinity and selectivity to cancer cells, they have proven to be very useful in nuclear medicine. As a result, efforts have been focused on discovering new peptide sequences for radiopeptide preparation. Nocardiotide A, a cyclic hexapeptide comprising the amino acids cyclo-Trp-Ile-Trp-Leu-Val-Ala (cWIWLVA) isolated from Nocardiopsis sp., has shown significant cytotoxicity against cancer cells, rendering it a suitable candidate for the process. Therefore, the present study aimed to design a stable and effective radiopeptide by labeling nocardiotide A with iodine-131 (131I), ensuring that its affinity to SSTR2 is not compromised. In silico study showed that structural modification of nocardiotide A labeled with 131iodine exhibited good affinity value, forming hydrogen bonds with key residues, such as Q.102 and T.194, which are essential in SSTR2. Based on the results, cyclic hexapeptides of cWIWLYA were selected for further synthesis, and its peptide product was confirmed by the presence of an ionic molecule peak m/z [M + Na]+ 855.4332 (yield, 25.60%). In vitro tests conducted on cWIWLYA showed that cWIWLYA can bind to HeLa cancer cells. Radiopeptide synthesis was initiated with radiolabeling of cWIWLYA by 131I using the chloramine-T method that showed a radiochemical yield of 93.37%. Non-radioactive iodine labeling reaction showed that iodination was successful, which detected the presence of di-iodinated peptide (I2-cWIWLYA) with m/z [M + Na]+ 1107.1138. In summary, a radiopeptide derived from nocardiotide A showed great potential for further development as a diagnostic and therapeutic agent in cancer treatment.

Machine Learning Advances in Predicting Peptide/Protein-Protein Interactions Based on Sequence Information for Lead Peptides Discovery

Peptides have shown increasing advantages and significant clinical value in drug discovery and development. With the development of high-throughput technologies and artificial intelligence (AI), machine learning (ML) methods for discovering new lead peptides have been expanded and incorporated into rational drug design. Predictions of peptide-protein interactions (PepPIs) and protein-protein interactions (PPIs) are both opportunities and challenges in computational biology, which will help to better understand the mechanisms of disease and provide the impetus for the discovery of lead peptides. This paper comprehensively reviews computational models for PepPI and PPI predictions. It begins with an introduction of various databases of peptide ligands and target proteins. Then it discusses data formats and feature representations for proteins and peptides. Furthermore, classical ML methods and emerging deep learning (DL) methods that can be used to train prediction models of PepPI and PPI are classified into four categories, and their advantages and disadvantages are analyzed. To assess the relative performance of different models, different validation protocols and evaluation indexes are discussed. The goal of this review is to help researchers quickly get started to develop computational frameworks using these integrated resources and eventually promote the discovery of lead peptides.

Peptide-derived coordination frameworks for biomimetic and selective separation

Peptide-derived metal-organic frameworks (PMOFs) have emerged as a class of biomimetic materials with attractive performances in analytical and bioanalytical chemistry. The incorporation of biomolecule peptides gives the frameworks conformational flexibility, guest adaptability, built-in chirality, and molecular recognition ability, which greatly accelerate the applications of PMOFs in enantiomeric separation, affinity separation, and the enrichment of bioactive species from complicated samples. This review focuses on the recent advances in the engineering and applications of PMOFs in selective separation. The unique biomimetic size-, enantio-, and affinity-selective performances for separation are discussed along with the chemical structures and functions of MOFs and peptides. Updates of the applications of PMOFs in adaptive separation of small molecules, chiral separation of drug molecules, and affinity isolation of bioactive species are summarized. Finally, the promising future and remaining challenges of PMOFs for selective separation of complex biosamples are discussed.

Computational design of dynamic receptor-peptide signaling complexes applied to chemotaxis

Engineering protein biosensors that sensitively respond to specific biomolecules by triggering precise cellular responses is a major goal of diagnostics and synthetic cell biology. Previous biosensor designs have largely relied on binding structurally well-defined molecules. In contrast, approaches that couple the sensing of flexible compounds to intended cellular responses would greatly expand potential biosensor applications. Here, to address these challenges, we develop a computational strategy for designing signaling complexes between conformationally dynamic proteins and peptides. To demonstrate the power of the approach, we create ultrasensitive chemotactic receptor-peptide pairs capable of eliciting potent signaling responses and strong chemotaxis in primary human T cells. Unlike traditional approaches that engineer static binding complexes, our dynamic structure design strategy optimizes contacts with multiple binding and allosteric sites accessible through dynamic conformational ensembles to achieve strongly enhanced signaling efficacy and potency. Our study suggests that a conformationally adaptable binding interface coupled to a robust allosteric transmission region is a key evolutionary determinant of peptidergic GPCR signaling systems. The approach lays a foundation for designing peptide-sensing receptors and signaling peptide ligands for basic and therapeutic applications.

Glycyl- l-histidyl- l-lysine-Cu2+ rescues cigarette smoking-induced skeletal muscle dysfunction via a sirtuin 1-dependent pathway

BACKGROUND

Skeletal muscle dysfunction is an important co-morbidity in patients with chronic obstructive pulmonary disease (COPD) and is significantly associated with increased mortality. Oxidative stress has been demonstrated an important trigger for COPD-related skeletal muscle dysfunction. Glycine-histidine-lysine (GHK) is an active tripeptide, which is a normal component of human plasma, saliva, and urine; promotes tissue regeneration; and acts as an anti-inflammatory and antioxidant properties. The purpose of this study was to determine whether GHK is involved in COPD-related skeletal muscle dysfunction.

METHODS

The plasma GHK level in patients with COPD (n = 9) and age-paired healthy subjects (n = 11) were detected using reversed-phase high-performance liquid chromatography. The complex GHK with Cu (GHK-Cu) was used in in vitro (C2C12 myotubes) and in vivo experiments (cigarette smoking [CS]-exposure mouse model) to explore the involvement of GHK in CS-induced skeletal muscle dysfunction.

RESULTS

Compared with healthy control, plasma GHK levels were decreased in patients with COPD (70.27 ± 38.87 ng/mL vs. 133.0 ± 54.54 ng/mL, P = 0.009). And plasma GHK levels in patients with COPD were associated with pectoralis muscle area (R = 0.684, P = 0.042), inflammatory factor TNF-α (R = -0.696, P = 0.037), and antioxidative stress factor SOD2 (R = 0.721, P = 0.029). GHK-Cu was found to rescue CSE-induced skeletal muscle dysfunction in C2C12 myotubes, as evidenced by increased expression of myosin heavy chain, reduced expression of MuRF1 and atrogin-1, elevated mitochondrial content, and enhanced resistance to oxidative stress. In CS-induced muscle dysfunction C57BL/6 mice, GHK-Cu treatment (0.2 and 2 mg/kg) reduces CS-induced muscle mass loss (skeletal muscle weight (1.19 ± 0.09% vs. 1.29 ± 0.06%, 1.40 ± 0.05%; P < 0.05) and muscle cross-sectional area elevated (1055 ± 552.4 μm² vs. 1797 ± 620.9 μm² , 2252 ± 534.0 μm² ; P < 0.001), and also rescues CS-induced muscle weakness, indicated by improved grip strength (175.5 ± 36.15 g vs. 257.6 ± 37.98 g, 339.1 ± 72.22 g; P < 0.01). Mechanistically, GHK-Cu directly binds and activates SIRT1(the binding energy was -6.1 kcal/mol). Through activating SIRT1 deacetylation, GHK-Cu inhibits FoxO3a transcriptional activity to reduce protein degradation, deacetylates Nrf2 and contribute to its action of reducing oxidative stress by generation of anti-oxidant enzymes, increases PGC-1α expression to promote mitochondrial function. Finally, GHK-Cu could protect mice against CS-induced skeletal muscle dysfunction via SIRT1.

CONCLUSIONS

Plasma glycyl- l-histidyl- l-lysine level in patients with chronic obstructive pulmonary disease was significantly decreased and was significantly associated with skeletal muscle mass. Exogenous administration of glycyl- l-histidyl- l-lysine-Cu²⁺ could protect against cigarette smoking-induced skeletal muscle dysfunction via sirtuin 1.

Yeast Protein Asf1 Possesses Modulating Activity towards Protein Kinase CK2

Protein kinase CK2 plays an important role in cell survival and protects regulatory proteins from caspase-mediated degradation during apoptosis. The consensus sequence of proteins phosphorylated by CK2 contains a cluster of acidic amino acids around the phosphorylation site. The poly-acidic sequence in yeast protein Asf1 is similar to the acidic loop in CK2β, which possesses a regulatory function. We observed that the overexpression of Asf1 in yeast cells influences cell growth. Experiments performed in vitro and in vivo indicate that yeast protein Asf1 inhibits protein kinase CK2. Our data suggest that each CK2 isoform might be regulated in a different way. Deletion of the amino or carboxyl end of Asf1 reveals that the acidic cluster close to the C-terminus is responsible for the activation or inhibition of CK2 activity.

Protocol for iterative optimization of modified peptides bound to protein targets

Peptides are commonly used as therapeutic agents. However, they suffer from easy degradation and instability. Replacing natural by non-natural amino acids can avoid these problems, and potentially improve the affinity towards the target protein. Here, we present a computational pipeline to optimize peptides based on adding non-natural amino acids while improving their binding affinity. The workflow is an iterative computational evolution algorithm, inspired by the PARCE protocol, that performs single-point mutations on the peptide sequence using modules from the Rosetta framework. The modifications can be guided based on the structural properties or previous knowledge of the biological system. At each mutation step, the affinity to the protein is estimated by sampling the complex conformations and applying a consensus metric using various open protein-ligand scoring functions. The mutations are accepted based on the score differences, allowing for an iterative optimization of the initial peptide. The sampling/scoring scheme was benchmarked with a set of protein-peptide complexes where experimental affinity values have been reported. In addition, a basic application using a known protein-peptide complex is also provided. The structure- and dynamic-based approach allows users to optimize bound peptides, with the option to personalize the code for further applications. The protocol, called mPARCE, is available at: https://github.com/rochoa85/mPARCE/.

LRH-1 (liver receptor homolog-1) derived affinity peptide ligand to inhibit interactions between β-catenin and LRH-1 in pancreatic cancer cells: from computational design to experimental validation

Poor prognosis, rapid progression and the lack of an effective treatment make pancreatic cancer one of the most lethal malignancies. Recent studies point to a role for liver receptor homolog-1 (LRH-1) in pathogenesis of pancreatic cancer and suggest prevention of the β-catenin/LRH-1 complex formation as a potential strategy for inhibition of the pancreas cancer cells progression. In the current investigation, we have followed a biomimetic strategy and designed an affinity peptide with sequence DEMEEPQQTE to inhibit formation of the β-catenin/LRH-1 complex. Quantitative real-time PCR experiments on the AsPC-1 pancreatic metastatic cells showed that the peptide has an inhibitory effect on the Wnt signaling proliferation line by reducing the expression levels of the CCND1, CCNE1, and MYC genes. Furthermore, the increased expression level of BAX gene showed that AsPC-1 cells were directed to the apoptosis pathway. At last, POU5F1, KLF4, and CD44 gene expression levels suggested that the peptide has an inhibitory effect on the stemness feature of the AsPC-1 cells. Here, we introduced a novel peptide inhibitor targeting an important protein-protein interaction, the β-catenin/LRH-1 complex, which may provide highly promising starting points for subsequent drug design. Communicated by Ramaswamy H. Sarma.

Identification of a Potential mRNA-based Vaccine Candidate against the SARS-CoV-2 Spike Glycoprotein: A Reverse Vaccinology Approach

The emergence of the novel coronavirus (SARS-CoV-2) in December 2019 has generated a devastating global consequence which makes the development of a rapidly deployable, effective and safe vaccine candidate an imminent global health priority. The design of most vaccine candidates has been directed at the induction of antibody responses against the trimeric spike glycoprotein of SARS-CoV-2, a class I fusion protein that aids ACE2 (angiotensin-converting enzyme 2) receptor binding. A variety of formulations and vaccinology approaches are being pursued for targeting the spike glycoprotein, including simian and human replication-defective adenoviral vaccines, subunit protein vaccines, nucleic acid vaccines and whole-inactivated SARS-CoV-2. Here, we directed a reverse vaccinology approach towards the design of a nucleic acid (mRNA-based) vaccine candidate. The "YLQPRTFLL" peptide sequence (position 269-277) which was predicted to be a B cell epitope and likewise a strong binder of the HLA*A-0201 was selected for the design of the vaccine candidate, having satisfied series of antigenicity assessments. Through the codon optimization protocol, the nucleotide sequence for the vaccine candidate design was generated and targeted at the human toll-like receptor 7 (TLR7). Bioinformatics analyses showed that the sequence "UACCUGCAGCCGCGUACCUUCCUGCUG" exhibited a strong affinity and likewise was bound to a stable cavity in the TLR7 pocket. This study is therefore expected to contribute to the research efforts directed at securing definitive preventive measures against the SARS-CoV-2 infection.

Predicting Protein-Peptide Complex Structures by Accounting for Peptide Flexibility and the Physicochemical Environment

Predicting protein-peptide complex structures is crucial to the understanding of a vast variety of peptide-mediated cellular processes and to peptide-based drug development. Peptide flexibility and binding mode ranking are the two major challenges for protein-peptide complex structure prediction. Peptides are highly flexible molecules, and therefore, brute-force modeling of peptide conformations of interest in protein-peptide docking is beyond current computing power. Inspired by the fact that the protein-peptide binding process is like protein folding, we developed a novel strategy, named MDockPeP2, which tries to address these challenges using physicochemical information embedded in abundant monomeric proteins with an exhaustive search strategy, in combination with an integrated global search and a local flexible minimization method. Only the peptide sequence and the protein crystal structure are required. The method was systemically assessed using a newly constructed structural database of 89 nonredundant protein-peptide complexes with the peptide sequence length ranging from 5 to 29 in which about half of the peptides are longer than 15 residues. MDockPeP2 yielded a total success rate of 58.4% (70.8, 79.8%) for the bound docking (i.e., with the bound receptor and fully flexible peptides) and 19.0% (44.8, 70.7%) for the challenging unbound docking when top 10 (100, 1000) models were considered for each prediction. MDockPeP2 achieved significantly higher success rates on two other datasets, peptiDB and LEADS-PEP, which contain only short- and medium-size peptides (≤ 15 residues). For peptiDB, our method obtained a success rate of 62.0% for the bound docking and 35.9% for the unbound docking when the top 10 models were considered. For LEADS-PEP, MDockPeP2 achieved a success rate of 69.8% when the top 10 models were considered. The program is available at https://zougrouptoolkit.missouri.edu/mdockpep2/download.html.

Divide-and-link peptide docking: a fragment-based peptide docking protocol

Protein-peptide interactions are crucial for various important cellular regulations, and are also a basis for understanding protein-protein interactions, protein folding and peptide drug design. Due to the limited structural data obtained using experimental methods, it is necessary to predict protein-peptide interaction modes using computational methods. In the present work, we designed a fragment-based docking protocol, Divide-and-Link Peptide Docking (DLPepDock), to predict protein-peptide binding modes. This protocol contains the following steps: dividing the peptide into fragments and separately docking the fragments using a third-party small molecular docking tool, linking the docked fragmental poses to form the whole peptide conformations via fragmental coordinate transformation using our in-house program, removing unreasonable poses according to several geometrical filters, extracting representative conformations after clustering for further minimization using the steepest descent and conjugation gradient methods based on a full-atom molecular force field and finally scoring using the MM/PBSA binding energy calculation implemented in Amber. When tested on the LEADS-PEP benchmark data set of 26 diverse complexes with peptides of 6-12 residues, FlexPepDock ab initio and AutoDock CrankPep achieved superior results. DLPepDock performed better than the other 15 docking protocols implemented in nine docking programs (HPepDock, DockThor, rDock, Glide, LeDock, AutoDock, AutoDock Vina, Surflex, and GOLD). The Linux scripts to call the third-party tools and run all the calculations.

Design of Monovalent and Chimeric Tetravalent Dengue Vaccine Using an Immunoinformatics Approach

An immunoinformatics technique was used to predict a monovalent amide immunogen candidate capable of producing therapeutic antibodies as well as a potent immunogen candidate capable of acting as a universal vaccination against all dengue fever virus serotypes. The capsid protein is an attractive goal for anti-DENV due to its position in the dengue existence cycle. The widely accessible immunological data, advances in antigenic peptide prediction using reverse vaccinology, and the introduction of molecular docking in immunoinformatics have directed vaccine manufacturing. The C-proteins of DENV-1-4 serotypes were known as antigens to assist with logical design. Binding epitopes for TC cells, TH cells, and B cells is predicted from structural dengue virus capsid proteins. Each T cell epitope of C-protein integrated with a B cell as a templet was used as a vaccine and produce antibodies in contrast to serotype of the dengue virus. A chimeric tetravalent vaccine was created by combining four vaccines, each representing four dengue serotypes, to serve as a standard vaccine candidate for all four Sero groups. The LKRARNRVS, RGFRKEIGR, KNGAIKVLR, and KAINVLRGF from dengue 1, dengue 2, dengue 3, and dengue 4 epitopes may be essential immunotherapeutic representatives for controlling outbreaks.

PepFun: Open Source Protocols for Peptide-Related Computational Analysis

Peptide research has increased during the last years due to their applications as biomarkers, therapeutic alternatives or as antigenic sub-units in vaccines. The implementation of computational resources have facilitated the identification of novel sequences, the prediction of properties, and the modelling of structures. However, there is still a lack of open source protocols that enable their straightforward analysis. Here, we present PepFun, a compilation of bioinformatics and cheminformatics functionalities that are easy to implement and customize for studying peptides at different levels: sequence, structure and their interactions with proteins. PepFun enables calculating multiple characteristics for massive sets of peptide sequences, and obtaining different structural observables derived from protein-peptide complexes. In addition, random or guided library design of peptide sequences can be customized for screening campaigns. The package has been created under the python language based on built-in functions and methods available in the open source projects BioPython and RDKit. We present two tutorials where we tested peptide binders of the MHC class II and the Granzyme B protease.

Plant Antimicrobial Peptides: State of the Art, In Silico Prediction and Perspectives in the Omics Era

Even before the perception or interaction with pathogens, plants rely on constitutively guardian molecules, often specific to tissue or stage, with further expression after contact with the pathogen. These guardians include small molecules as antimicrobial peptides (AMPs), generally cysteine-rich, functioning to prevent pathogen establishment. Some of these AMPs are shared among eukaryotes (eg, defensins and cyclotides), others are plant specific (eg, snakins), while some are specific to certain plant families (such as heveins). When compared with other organisms, plants tend to present a higher amount of AMP isoforms due to gene duplications or polyploidy, an occurrence possibly also associated with the sessile habit of plants, which prevents them from evading biotic and environmental stresses. Therefore, plants arise as a rich resource for new AMPs. As these molecules are difficult to retrieve from databases using simple sequence alignments, a description of their characteristics and in silico (bioinformatics) approaches used to retrieve them is provided, considering resources and databases available. The possibilities and applications based on tools versus database approaches are considerable and have been so far underestimated.

Protein-Protein Interactions Mediated by Intrinsically Disordered Protein Regions Are Enriched in Missense Mutations

Because proteins are fundamental to most biological processes, many genetic diseases can be traced back to single nucleotide variants (SNVs) that cause changes in protein sequences. However, not all SNVs that result in amino acid substitutions cause disease as each residue is under different structural and functional constraints. Influential studies have shown that protein–protein interaction interfaces are enriched in disease-associated SNVs and depleted in SNVs that are common in the general population. These studies focus primarily on folded (globular) protein domains and overlook the prevalent class of protein interactions mediated by intrinsically disordered regions (IDRs). Therefore, we investigated the enrichment patterns of missense mutation-causing SNVs that are associated with disease and cancer, as well as those present in the healthy population, in structures of IDR-mediated interactions with comparisons to classical globular interactions. When comparing the different categories of interaction interfaces, division of the interface regions into solvent-exposed rim residues and buried core residues reveal distinctive enrichment patterns for the various types of missense mutations. Most notably, we demonstrate a strong enrichment at the interface core of interacting IDRs in disease mutations and its depletion in neutral ones, which supports the view that the disruption of IDR interactions is a mechanism underlying many diseases. Intriguingly, we also found an asymmetry across the IDR interaction interface in the enrichment of certain missense mutation types, which may hint at an increased variant tolerance and urges further investigations of IDR interactions.

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

MOTIVATION

Protein-peptide interactions mediate a wide variety of cellular and biological functions. Methods for predicting these interactions have garnered a lot of interest over the past few years, as witnessed by the rapidly growing number of peptide-based therapeutic molecules currently in clinical trials. The size and flexibility of peptides has shown to be challenging for existing automated docking software programs.

RESULTS

Here we present AutoDock CrankPep or ADCP in short, a novel approach to dock flexible peptides into rigid receptors. ADCP folds a peptide in the potential field created by the protein to predict the protein-peptide complex. We show that it outperforms leading peptide docking methods on two protein-peptide datasets commonly used for benchmarking docking methods: LEADS-PEP and peptiDB, comprised of peptides with up to 15 amino acids in length. Beyond these datasets, ADCP reliably docked a set of protein-peptide complexes containing peptides ranging in lengths from 16 to 20 amino acids. The robust performance of ADCP on these longer peptides enables accurate modeling of peptide-mediated protein-protein interactions and interactions with disordered proteins.

AVAILABILITY AND IMPLEMENTATION

ADCP is distributed under the LGPL 2.0 open source license and is available at http://adcp.scripps.edu. The source code is available at https://github.com/ccsb-scripps/ADCP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Coarse-grained (hybrid) integrative modeling of biomolecular interactions

The computational modeling field has vastly evolved over the past decades. The early developments of simplified protein systems represented a stepping stone towards establishing more efficient approaches to sample intricated conformational landscapes. Downscaling the level of resolution of biomolecules to coarser representations allows for studying protein structure, dynamics and interactions that are not accessible by classical atomistic approaches. The combination of different resolutions, namely hybrid modeling, has also been proved as an alternative when mixed levels of details are required. In this review, we provide an overview of coarse-grained/hybrid models focusing on their applicability in the modeling of biomolecular interactions. We give a detailed list of ready-to-use modeling software for studying biomolecular interactions allowing various levels of coarse-graining and provide examples of complexes determined by integrative coarse-grained/hybrid approaches in combination with experimental information.

The emerging role of computational design in peptide macrocycle drug discovery

Drug discovery is a laborious process with rising cost per new drug. Peptide macrocycles are promising therapeutics, though conformational flexibility can reduce target affinity and specificity. Recent computational advancements address this problem by enabling rational design of rigidly folded peptide macrocycles.

AREAS COVERED

This review summarizes currently approved peptide macrocycle therapeutics and discusses advantages of mesoscale drugs over small molecules or protein therapeutics. It describes the history, rationale, and state of the art of computational tools, such as Rosetta, that allow the design of rigidly structured peptide macrocycles. The emerging pipeline for designing peptide macrocycle drugs is described, including current challenges in designing permeable molecules that can emulate the chameleonic behavior of natural macrocycles. Prospects for reducing computational cost and improving accuracy with emerging computational technologies are also discussed.

EXPERT OPINION

To embrace computational design of peptide macrocycle drugs, we must shift current attitudes regarding the role of computation in drug discovery, and move beyond Lipinski's rules. This technology has the potential to shift failures to earlier in silico stages of the drug discovery process, improving success rates in costly clinical trials. Given the available tools, now is the time for drug developers to incorporate peptide macrocycle design into drug discovery pipelines.

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

Protein-peptide docking, which predicts the complex structure between a protein and a peptide, is a valuable computational tool in peptide therapeutics development and the mechanistic investigation of peptides involved in cellular processes. Although current peptide docking approaches are often able to sample near-native peptide binding modes, correctly identifying those near-native modes from decoys is still challenging because of the extremely high complexity of the peptide binding energy landscape. In this study, we have developed an efficient postdocking rescoring protocol using a combined scoring function of knowledge-based ITScorePP potentials and physics-based MM-GBSA energies. Tested on five benchmark/docking test sets, our postdocking strategy showed an overall significantly better performance in binding mode prediction and score-rmsd correlation than original docking approaches. Specifically, our postdocking protocol outperformed original docking approaches with success rates of 15.8 versus 10.5% for pepATTRACT on the Global_57 benchmark, 5.3 versus 5.3% for CABS-dock on the Global_57 benchmark, 17.0 versus 11.3% for FlexPepDock on the LEADS-PEP data set, 40.3 versus 33.9% for HPEPDOCK on the Local_62 benchmark, and 64.2 versus 52.8% for HPEPDOCK on the LEADS-PEP data set when the top prediction was considered. These results demonstrated the efficacy and robustness of our postdocking protocol.

Proteasome Inhibitor Drugs

Proteasomes are large, multicatalytic protein complexes that cleave cellular proteins into peptides. There are many distinct forms of proteasomes that differ in catalytically active subunits, regulatory subunits, and associated proteins. Proteasome inhibitors are an important class of drugs for the treatment of multiple myeloma and mantle cell lymphoma, and they are being investigated for other diseases. Bortezomib (Velcade) was the first proteasome inhibitor to be approved by the US Food and Drug Administration. Carfilzomib (Kyprolis) and ixazomib (Ninlaro) have recently been approved, and more drugs are in development. While the primary mechanism of action is inhibition of the proteasome, the downstream events that lead to selective cell death are not entirely clear. Proteasome inhibitors have been found to affect protein turnover but at concentrations that are much higher than those achieved clinically, raising the possibility that some of the effects of proteasome inhibitors are mediated by other mechanisms.

Flexible docking of peptides to proteins using CABS-dock

Molecular docking of peptides to proteins can be a useful tool in the exploration of the possible peptide binding sites and poses. CABS-dock is a method for protein-peptide docking that features significant conformational flexibility of both the peptide and the protein molecules during the peptide search for a binding site. The CABS-dock has been made available as a web server and a standalone package. The web server is an easy to use tool with a simple web interface. The standalone package is a command-line program dedicated to professional users. It offers a number of advanced features, analysis tools and support for large-sized systems. In this article, we outline the current status of the CABS-dock method, its recent developments, applications, and challenges ahead.

PepPro: A Nonredundant Structure Data Set for Benchmarking Peptide-Protein Computational Docking

We present a nonredundant benchmark, coined PepPro, for testing peptide-protein docking algorithms. Currently, PepPro contains 89 nonredundant experimentally determined peptide-protein complex structures, with peptide sequence lengths ranging from 5 to 30 amino acids. The benchmark covers peptides with distinct secondary structures, including helix, partial helix, a mixture of helix and β-sheet, β-sheet formed through binding, β-sheet formed through self-folding, and coil. In addition, unbound proteins' structures are provided for 58 complexes and can be used for testing the ability of a docking algorithm handling the conformational changes of proteins during the binding process. PepPro should benefit the docking community for the development and improvement of peptide docking algorithms. The benchmark is available at http://zoulab.dalton.missouri.edu/PepPro_benchmark. © 2019 Wiley Periodicals, Inc.

HPEPDOCK: a web server for blind peptide-protein docking based on a hierarchical algorithm

Protein-peptide interactions are crucial in many cellular functions. Therefore, determining the structure of protein-peptide complexes is important for understanding the molecular mechanism of related biological processes and developing peptide drugs. HPEPDOCK is a novel web server for blind protein-peptide docking through a hierarchical algorithm. Instead of running lengthy simulations to refine peptide conformations, HPEPDOCK considers the peptide flexibility through an ensemble of peptide conformations generated by our MODPEP program. For blind global peptide docking, HPEPDOCK obtained a success rate of 33.3% in binding mode prediction on a benchmark of 57 unbound cases when the top 10 models were considered, compared to 21.1% for pepATTRACT server. HPEPDOCK also performed well in docking against homology models and obtained a success rate of 29.8% within top 10 predictions. For local peptide docking, HPEPDOCK achieved a high success rate of 72.6% on a benchmark of 62 unbound cases within top 10 predictions, compared to 45.2% for HADDOCK peptide protocol. Our HPEPDOCK server is computationally efficient and consumed an average of 29.8 mins for a global peptide docking job and 14.2 mins for a local peptide docking job. The HPEPDOCK web server is available at http://huanglab.phys.hust.edu.cn/hpepdock/.

Modified Potential Functions Result in Enhanced Predictions of a Protein Complex by All-Atom Molecular Dynamics Simulations, Confirming a Stepwise Association Process for Native Protein-Protein Interactions

The relative prevalence of native protein-protein interactions (PPIs) are the cornerstone for understanding the structure, dynamics and mechanisms of function of protein complexes. In this study, we develop a scheme for scaling the protein-water interaction in the CHARMM36 force field, in order to better fit the solvation free energy of amino acids side-chain analogues. We find that the molecular dynamics simulation with the scaled force field, CHARMM36s, as well as a recently released version, CHARMM36m, effectively improve on the overly sticky association of proteins, such as ubiquitin. We investigate the formation of a heterodimer protein complex between the SAM domains of the EphA2 receptor and the SHIP2 enzyme by performing a combined total of 48 μs simulations with the different potential functions. While the native SAM heterodimer is only predicted at a low rate of 6.7% with the original CHARMM36 force field, the yield is increased to 16.7% with CHARMM36s, and to 18.3% with CHARMM36m. By analyzing the 25 native SAM complexes formed in the simulations, we find that their formation involves a preorientation guided by Coulomb interactions, consistent with an electrostatic steering mechanism. In 12 cases, the complex could directly transform to the native protein interaction surfaces with only small adjustments in domain orientation. In the other 13 cases, orientational and/or translational adjustments are needed to reach the native complex. Although the tendency for non-native complexes to dissociate has nearly doubled with the modified potential functions, a dissociation followed by a reassociation to the correct complex structure is still rare. Instead, the remaining non-native complexes undergo configurational changes/surface searching, which, however, rarely leads to native structures on a time scale of 250 ns. These observations provide a rich picture of the mechanisms of protein-protein complex formation and suggest that computational predictions of native complex PPIs could be improved further.

Affordable Microfluidic Bead-Sorting Platform for Automated Selection of Porous Particles Functionalized with Bioactive Compounds

The ability to rapidly and accurately evaluate bioactive compounds immobilized on porous particles is crucial in the discovery of drugs, diagnostic reagents, ligands, and catalysts. Existing options for solid phase screening of bioactive compounds, while highly effective and well established, can be cost-prohibitive for proof-of-concept and early stage work, limiting its applicability and flexibility in new research areas. Here, we present a low-cost microfluidics-based platform enabling automated screening of small porous beads from solid-phase peptide libraries with high sensitivity and specificity, to identify leads with high binding affinity for a biological target. The integration of unbiased computer assisted image processing and analysis tools, provided the platform with the flexibility of sorting through beads with distinct fluorescence patterns. The customized design of the microfluidic device helped with handling beads with different diameters (~100-300 µm). As a microfluidic device, this portable novel platform can be integrated with a variety of analytical instruments to perform screening. In this study, the system utilizes fluorescence microscopy and unsupervised image analysis, and can operate at a sorting speed of up to 125 beads/hr (~3.5 times faster than a trained operator) providing >90% yield and >90% bead sorting accuracy. Notably, the device has proven successful in screening a model solid-phase peptide library by showing the ability to select beads carrying peptides binding a target protein (human IgG).

In Silico Peptide Ligation: Iterative Residue Docking and Linking as a New Approach to Predict Protein-Peptide Interactions

Peptide–protein interactions are corner-stones of living functions involved in essential mechanisms, such as cell signaling. Given the difficulty of obtaining direct experimental structural biology data, prediction of those interactions is of crucial interest for the rational development of new drugs, notably to fight diseases, such as cancer or Alzheimer’s disease. Because of the high flexibility of natural unconstrained linear peptides, prediction of their binding mode in a protein cavity remains challenging. Several theoretical approaches have been developed in the last decade to address this issue. Nevertheless, improvements are needed, such as the conformation prediction of peptide side-chains, which are dependent on peptide length and flexibility. Here, we present a novel in silico method, Iterative Residue Docking and Linking (IRDL), to efficiently predict peptide–protein interactions. In order to reduce the conformational space, this innovative method splits peptides into several short segments. Then, it uses the performance of intramolecular covalent docking to rebuild, sequentially, the complete peptide in the active site of its protein target. Once the peptide is constructed, a rescoring step is applied in order to correctly rank all IRDL solutions. Applied on a set of 11 crystallized peptide–protein complexes, the IRDL method shows promising results, since it is able to retrieve experimental binding conformations with a Root Mean Square Deviation (RMSD) below 2 Å in the top five ranked solutions. For some complexes, IRDL method outperforms two other docking protocols evaluated in this study. Hence, IRDL is a new tool that could be used in drug design projects to predict peptide–protein interactions.

Interaction of Insulin-Like Growth Factor-Binding Protein 3 With Hyaluronan and Its Regulation by Humanin and CD44

Insulin-like growth factor-binding protein-3 (IGFBP-3) belongs to a family of IGF-binding proteins. Humanin is a peptide known to bind residues 215-232 of mature IGFBP-3 in the C-terminal region of the protein. This region of IGFBP-3 was shown earlier to bind certain glycosaminoglycans including hyaluronan (HA). Here, we characterized the binding affinities of the IGFBP-3 protein and peptide (²¹⁵-KKGFYKKKQCRPSKGRKR-²³²) to HA and to humanin and found that HA binds with a weaker affinity to this region than does humanin. Either HA or humanin could bind to this IGFBP-3 segment, but not simultaneously. The HA receptor, CD44, blocked HA binding to IGFBP-3 but had no effect on binding of humanin to either IGFBP-3 or its peptide. Upon incubation of HA with CD44 and either IGFBP-3 protein or peptide, humanin was effective at binding and sequestering IGFBP-3 or peptide, thereby enabling access of CD44 to HA. We show that IGFBP-3 and humanin in the medium of A549 lung cancer cells can immunoprecipitate in a complex. However, the fraction of IGFBP-3 in the medium that is able to bind HA was not complexed with humanin suggesting that HA binding to the 215-232 segment renders it inaccessible for binding to humanin. Moreover, while the cytotoxic effects of IGFBP-3 on cell viability were reversed by humanin, blocking HA-CD44 interaction with an anti-CD44 antibody in combination with IGFBP-3 did not have an additive negative effect on cell viability suggesting that IGFBP-3 exerts its cytotoxic effects on cell survival through a mechanism that depends on HA-CD44 interactions.

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

Given the importance of peptide-mediated protein interactions in cellular processes, protein-peptide docking has received increasing attention. Here, we have developed a Hierarchical flexible Peptide Docking approach through fast generation and ensemble docking of peptide conformations, which is referred to as HPepDock. Tested on the LEADS-PEP benchmark data set of 53 diverse complexes with peptides of 3-12 residues, HPepDock performed significantly better than the 11 docking protocols of five small-molecule docking programs (DOCK, AutoDock, AutoDock Vina, Surflex, and GOLD) in predicting near-native binding conformations. HPepDock was also evaluated on the 19 bound/unbound and 10 unbound/unbound protein-peptide complexes of the Glide SP-PEP benchmark and showed an overall better performance than Glide SP-PEP+MM-GBSA and FlexPepDock in both bound and unbound docking. HPepDock is computationally efficient, and the average running time for docking a peptide is ∼15 min with the range from about 1 min for short peptides to around 40 min for long peptides.

Heterodimer Binding Scaffolds Recognition via the Analysis of Kinetically Hot Residues

Physical interactions between proteins are often difficult to decipher. The aim of this paper is to present an algorithm that is designed to recognize binding patches and supporting structural scaffolds of interacting heterodimer proteins using the Gaussian Network Model (GNM). The recognition is based on the (self) adjustable identification of kinetically hot residues and their connection to possible binding scaffolds. The kinetically hot residues are residues with the lowest entropy, i.e., the highest contribution to the weighted sum of the fastest modes per chain extracted via GNM. The algorithm adjusts the number of fast modes in the GNM's weighted sum calculation using the ratio of predicted and expected numbers of target residues (contact and the neighboring first-layer residues). This approach produces very good results when applied to dimers with high protein sequence length ratios. The protocol's ability to recognize near native decoys was compared to the ability of the residue-level statistical potential of Lu and Skolnick using the Sternberg and Vakser decoy dimers sets. The statistical potential produced better overall results, but in a number of cases its predicting ability was comparable, or even inferior, to the prediction ability of the adjustable GNM approach. The results presented in this paper suggest that in heterodimers at least one protein has interacting scaffold determined by the immovable, kinetically hot residues. In many cases, interacting proteins (especially if being of noticeably different sizes) either behave as a rigid lock and key or, presumably, exhibit the opposite dynamic behavior. While the binding surface of one protein is rigid and stable, its partner's interacting scaffold is more flexible and adaptable.

General Prediction of Peptide-MHC Binding Modes Using Incremental Docking: A Proof of Concept

The class I major histocompatibility complex (MHC) is capable of binding peptides derived from intracellular proteins and displaying them at the cell surface. The recognition of these peptide-MHC (pMHC) complexes by T-cells is the cornerstone of cellular immunity, enabling the elimination of infected or tumoral cells. T-cell-based immunotherapies against cancer, which leverage this mechanism, can greatly benefit from structural analyses of pMHC complexes. Several attempts have been made to use molecular docking for such analyses, but pMHC structure remains too challenging for even state-of-the-art docking tools. To overcome these limitations, we describe the use of an incremental meta-docking approach for structural prediction of pMHC complexes. Previous methods applied in this context used specific constraints to reduce the complexity of this prediction problem, at the expense of generality. Our strategy makes no assumption and can potentially be used to predict binding modes for any pMHC complex. Our method has been tested in a re-docking experiment, reproducing the binding modes of 25 pMHC complexes whose crystal structures are available. This study is a proof of concept that incremental docking strategies can lead to general geometry prediction of pMHC complexes, with potential applications for immunotherapy against cancer or infectious diseases.

Enriching Peptide Libraries for Binding Affinity and Specificity Through Computationally Directed Library Design

Peptide reagents with high affinity or specificity for their target protein interaction partner are of utility for many important applications. Optimization of peptide binding by screening large libraries is a proven and powerful approach. Libraries designed to be enriched in peptide sequences that are predicted to have desired affinity or specificity characteristics are more likely to yield success than random mutagenesis. We present a library optimization method in which the choice of amino acids to encode at each peptide position can be guided by available experimental data or structure-based predictions. We discuss how to use analysis of predicted library performance to inform rounds of library design. Finally, we include protocols for more complex library design procedures that consider the chemical diversity of the amino acids at each peptide position and optimize a library score based on a user-specified input model.

Structure-based Methods for Binding Mode and Binding Affinity Prediction for Peptide-MHC Complexes

Understanding the mechanisms involved in the activation of an immune response is essential to many fields in human health, including vaccine development and personalized cancer immunotherapy. A central step in the activation of the adaptive immune response is the recognition, by T-cell lymphocytes, of peptides displayed by a special type of receptor known as Major Histocompatibility Complex (MHC). Considering the key role of MHC receptors in T-cell activation, the computational prediction of peptide binding to MHC has been an important goal for many immunological applications. Sequence- based methods have become the gold standard for peptide-MHC binding affinity prediction, but structure-based methods are expected to provide more general predictions (i.e., predictions applicable to all types of MHC receptors). In addition, structural modeling of peptide-MHC complexes has the potential to uncover yet unknown drivers of T-cell activation, thus allowing for the development of better and safer therapies. In this review, we discuss the use of computational methods for the structural modeling of peptide-MHC complexes (i.e., binding mode prediction) and for the structure-based prediction of binding affinity.

Efficient conformational ensemble generation of protein-bound peptides

Conformation generation of protein-bound peptides is critical for the determination of protein–peptide complex structures. Despite significant progress in conformer generation of small molecules, few methods have been developed for modeling protein-bound peptide conformations. Here, we have developed a fast de novo peptide modeling algorithm, referred to as MODPEP, for conformational sampling of protein-bound peptides. Given a sequence, MODPEP builds the peptide 3D structure from scratch by assembling amino acids or helix fragments based on constructed rotamer and helix libraries. The MODPEP algorithm was tested on a diverse set of 910 experimentally determined protein-bound peptides with 3–30 amino acids from the PDB and obtained an average accuracy of 1.90 Å when 200 conformations were sampled for each peptide. On average, MODPEP obtained a success rate of 74.3% for all the 910 peptides and ≥ 90% for short peptides with 3–10 amino acids in reproducing experimental protein-bound structures. Comparative evaluations of MODPEP with three other conformer generation methods, PEP-FOLD3, RDKit, and Balloon, have also been performed in both accuracy and success rate. MODPEP is fast and can generate 100 conformations for less than one second. The fast MODPEP will be beneficial for large-scale de novo modeling and docking of peptides. The MODPEP program and libraries are available for download at http://huanglab.phys.hust.edu.cn/.

Modeling EphB4-EphrinB2 protein-protein interaction using flexible docking of a short linear motif

Background

Many protein–protein interactions are mediated by a short linear motif. Usually, amino acid sequences of those motifs are known or can be predicted. It is much harder to experimentally characterize or predict their structure in the bound form. In this work, we test a possibility of using flexible docking of a short linear motif to predict the interaction interface of the EphB4-EphrinB2 complex (a system extensively studied for its significance in tumor progression).

Methods

In the modeling, we only use knowledge about the motif sequence and experimental structures of EphB4-EphrinB2 complex partners. The proposed protocol enables efficient modeling of significant conformational changes in the short linear motif fragment during molecular docking simulation. For the docking simulations, we use the CABS-dock method for docking fully flexible peptides to flexible protein receptors (available as a server at http://biocomp.chem.uw.edu.pl/CABSdock/). Based on the docking result, the protein–protein complex is reconstructed and refined.

Results

Using this novel protocol, we obtained an accurate EphB4-EphrinB2 interaction model.

Conclusions

The results show that the CABS-dock method may be useful as the primary docking tool in specific protein–protein docking cases similar to EphB4-EphrinB2 complex—that is, where a short linear motif fragment can be identified.

A New Mixed All-Atom/Coarse-Grained Model: Application to Melittin Aggregation in Aqueous Solution

We introduce a new mixed resolution, all-atom/coarse-grained approach (AACG), for modeling peptides in aqueous solution and apply it to characterizing the aggregation of melittin. All of the atoms in peptidic components are represented, while a single site is used for each water molecule. With the full flexibility of the peptide retained, our AACG method achieves speedups by a factor of 3-4 for CPU time reduction and another factor of roughly 7 for diffusion. An Ewald treatment permits the inclusion of long-range electrostatic interactions. These characteristics fit well with the requirements for studying peptide association and aggregation, where the system sizes and time scales require considerable computational resources with all-atom models. In particular, AACG is well suited for biologics since changes in peptide shape and long-range electrostatics may play an important role. The application of AACG to melittin, a 26-residue peptide with a well-known propensity to aggregate in solution, serves as an initial demonstration of this technology for studying peptide aggregation. We observed the formation of melittin aggregates during our simulations and characterized the time-evolution of aggregate size distribution, buried surface areas, and residue contacts. Key interactions including π-cation and π-stacking involving TRP19 were also examined. Our AACG simulations demonstrated a clear salt effect and a moderate temperature effect on aggregation and support the molten globule model of melittin aggregates. As a showcase, this work illustrates the useful role for AACG in investigations of peptide aggregation and its potential to guide formulation and design of biologics.

FlexPepDock lessons from CAPRI peptide-protein rounds and suggested new criteria for assessment of model quality and utility

CAPRI rounds 28 and 29 included, for the first time, peptide-receptor targets of three different systems, reflecting increased appreciation of the importance of peptide-protein interactions. The CAPRI rounds allowed us to objectively assess the performance of Rosetta FlexPepDock, one of the first protocols to explicitly include peptide flexibility in docking, accounting for peptide conformational changes upon binding. We discuss here successes and challenges in modeling these targets: we obtain top-performing, high-resolution models of the peptide motif for cases with known binding sites but there is a need for better modeling of flanking regions, as well as better selection criteria, in particular for unknown binding sites. These rounds have also provided us the opportunity to reassess the success criteria, to better reflect the quality of a peptide-protein complex model. Using all models submitted to CAPRI, we analyze the correlation between current classification criteria and the ability to retrieve critical interface features, such as hydrogen bonds and hotspots. We find that loosening the backbone (and ligand) RMSD threshold, together with a restriction on the side chain RMSD measure, allows us to improve the selection of high-accuracy models. We also suggest a new measure to assess interface hydrogen bond recovery, which is not assessed by the current CAPRI criteria. Finally, we find that surprisingly much can be learned from rather inaccurate models about binding hotspots, suggesting that the current status of peptide-protein docking methods, as reflected by the submitted CAPRI models, can already have a significant impact on our understanding of protein interactions. Proteins 2017; 85:445-462. © 2016 Wiley Periodicals, Inc.

Modeling disordered protein interactions from biophysical principles

Disordered protein-protein interactions (PPIs), those involving a folded protein and an intrinsically disordered protein (IDP), are prevalent in the cell, including important signaling and regulatory pathways. IDPs do not adopt a single dominant structure in isolation but often become ordered upon binding. To aid understanding of the molecular mechanisms of disordered PPIs, it is crucial to obtain the tertiary structure of the PPIs. However, experimental methods have difficulty in solving disordered PPIs and existing protein-protein and protein-peptide docking methods are not able to model them. Here we present a novel computational method, IDP-LZerD, which models the conformation of a disordered PPI by considering the biophysical binding mechanism of an IDP to a structured protein, whereby a local segment of the IDP initiates the interaction and subsequently the remaining IDP regions explore and coalesce around the initial binding site. On a dataset of 22 disordered PPIs with IDPs up to 69 amino acids, successful predictions were made for 21 bound and 18 unbound receptors. The successful modeling provides additional support for biophysical principles. Moreover, the new technique significantly expands the capability of protein structure modeling and provides crucial insights into the molecular mechanisms of disordered PPIs.

A substantial fraction of the proteins encoded in genomes are intrinsically disordered proteins (IDPs), which lack a single stable structure in the native state. IDPs serve many functions including mediating protein-protein interactions (PPIs). Such disordered PPIs are prevalent in important regulatory pathways, including many interactions of the tumor suppressor protein p53. To elucidate the molecular mechanisms of disordered PPIs, obtaining tertiary structure information is essential; however, they are difficult to study with experimental techniques and existing computational protein-protein and protein-peptide modeling methods are unable to model disordered PPIs. Here we present a novel computational method for modeling the structure of disordered PPIs, which is the first of this sort. The method, IDP-LZerD, is designed to follow a known biophysical picture of the mechanism of how IDPs interact with structured proteins. IDP-LZerD successfully modeled the majority of disordered PPIs tested. This technique opens up new possibilities for structural studies of IDPs and their interactions.

Application of the ATTRACT Coarse-Grained Docking and Atomistic Refinement for Predicting Peptide-Protein Interactions

Peptide-protein interactions are abundant in the cell and form an important part of the interactome. Large-scale modeling of peptide-protein complexes requires a fully blind approach; i.e., simultaneously predicting the peptide-binding site and the peptide conformation to high accuracy. Here, we present one of the first fully blind peptide-protein docking protocols, pepATTRACT. It combines a coarse-grained ensemble docking search of the entire protein surface with two stages of atomistic flexible refinement. pepATTRACT yields high-quality predictions for 70 % of the cases when tested on a large benchmark of peptide-protein complexes. This performance in fully blind mode is similar to state-of-the-art local docking approaches that use information on the location of the binding site. Limiting the search to the peptide-binding region, the resulting pepATTRACT-local approach further improves the performance. Docking scripts for pepATTRACT and pepATTRACT-local can be generated via a web interface at www.attract.ph.tum.de/peptide.html . Here, we explain how to set up a docking run with the pepATTRACT web interface and demonstrate its usage by an application on binding of disordered regions from tumor suppressor p53 to a partner protein.

AnchorDock for Blind Flexible Docking of Peptides to Proteins

Due to increasing interest in peptides as signaling modulators and drug candidates, several methods for peptide docking to their target proteins are under active development. The "blind" docking problem, where the peptide-binding site on the protein surface is unknown, presents one of the current challenges in the field. AnchorDock protocol was developed by Ben-Shimon and Niv to address this challenge.This protocol narrows the docking search to the most relevant parts of the conformational space. This is achieved by pre-folding the free peptide and by computationally detecting anchoring spots on the surface of the unbound protein. Multiple flexible simulated annealing molecular dynamics (SAMD) simulations are subsequently carried out, starting from pre-folded peptide conformations, constrained to the various precomputed anchoring spots.Here, AnchorDock is demonstrated using two known protein-peptide complexes. A PDZ-peptide complex provides a relatively easy case due to the relatively small size of the protein, and a typical peptide conformation and binding region; a more challenging example is a complex between USP7_N-term and a p53-derived peptide, where the protein is larger, and the peptide conformation and a binding site are generally assumed to be unknown. AnchorDock returned native-like solutions ranked first and third for the PDZ and USP7 complexes, respectively. We describe the procedure step by step and discuss possible modifications where applicable.

Detection of Peptide-Binding Sites on Protein Surfaces Using the Peptimap Server

Peptide-mediated interactions are of primordial importance to the cell, and the structure of such interaction provides an important starting point for their further characterization. In many cases, the structure of the peptide-protein complex has not been solved by experiment, and modeling tools need to be applied to generate structural models of the interaction. PeptiMap is a protocol that identifies the peptide-binding site when only the structure of the receptor is known, but no information about where the peptide binds is available. This is achieved by mapping the surface for solvents to identify ligand-binding sites, similar in approach to ANCHORMAP in which amino acids are mapped. Peptimap is a free open access web-based server. It can be accessed at http://peptimap.cluspro.org .

Template-Based Prediction of Protein-Peptide Interactions by Using GalaxyPepDock

We introduce a web server called GalaxyPepDock that predicts protein-peptide interactions based on templates. With the continuously increasing size of the protein structure database, the probability of finding related proteins for templates is increasing. GalaxyPepDock takes a protein structure and a peptide sequence as input and returns protein-peptide complex structures as output. Templates for protein-peptide complex structures are selected from the structure database considering similarity to the target protein structure and to putative protein-peptide interactions as estimated by protein structure alignment and peptide sequence alignment. Complex structures are then built from the template structures by template-based modeling. By further structure refinement that performs energy-based optimization, structural aspects that are missing in the template structures or that are not compatible with the given protein and peptide are refined. During the refinement, flexibilities of both protein and peptide induced by binding are considered. The atomistic protein-peptide interactions predicted by GalaxyPepDock can offer important clues for designing new peptides with desired binding properties.

Modeling Peptide-Protein Structure and Binding Using Monte Carlo Sampling Approaches: Rosetta FlexPepDock and FlexPepBind

Many signaling and regulatory processes involve peptide-mediated protein interactions, i.e., the binding of a short stretch in one protein to a domain in its partner. Computational tools that generate accurate models of peptide-receptor structures and binding improve characterization and manipulation of known interactions, help to discover yet unknown peptide-protein interactions and networks, and bring into reach the design of peptide-based drugs for targeting specific systems of medical interest.Here, we present a concise overview of the Rosetta FlexPepDock protocol and its derivatives that we have developed for the structure-based characterization of peptide-protein binding. Rosetta FlexPepDock was built to generate precise models of protein-peptide complex structures, by effectively addressing the challenge of the considerable conformational flexibility of the peptide. Rosetta FlexPepBind is an extension of this protocol that allows characterizing peptide-binding affinities and specificities of various biological systems, based on the structural models generated by Rosetta FlexPepDock. We provide detailed descriptions and guidelines for the usage of these protocols, and on a specific example, we highlight the variety of different challenges that can be met and the questions that can be answered with Rosetta FlexPepDock.