Recent Advances in Computational Models for the Study of Protein-Peptide Interactions

ProtTrans and multi-window scanning convolutional neural networks for the prediction of protein-peptide interaction sites

This study delves into the prediction of protein-peptide interactions using advanced machine learning techniques, comparing models such as sequence-based, standard CNNs, and traditional classifiers. Leveraging pre-trained language models and multi-view window scanning CNNs, our approach yields significant improvements, with ProtTrans standing out based on 2.1 billion protein sequences and 393 billion amino acids. The integrated model demonstrates remarkable performance, achieving an AUC of 0.856 and 0.823 on the PepBCL Set_1 and Set_2 datasets, respectively. Additionally, it attains a Precision of 0.564 in PepBCL Set 1 and 0.527 in PepBCL Set 2, surpassing the performance of previous methods. Beyond this, we explore the application of this model in cancer therapy, particularly in identifying peptide interactions for selective targeting of cancer cells, and other fields. The findings of this study contribute to bioinformatics, providing valuable insights for drug discovery and therapeutic development.

Correlation of Experimental and Calculated Inhibition Constants of Protease Inhibitor Complexes

Predicting the potency of inhibitors is key to in silico screening of promising synthetic or natural compounds. Here we describe a predictive workflow that provides calculated inhibitory values, which concord well with empirical data. Calculations of the free interaction energy ΔG with the YASARA plugin FoldX were used to derive inhibition constants Ki from PDB coordinates of protease-inhibitor complexes. At the same time, corresponding KD values were obtained from the PRODIGY server. These results correlated well with the experimental values, particularly for serine proteases. In addition, analyses were performed for inhibitory complexes of cysteine and aspartic proteases, as well as of metalloproteases, whereby the PRODIGY data appeared to be more consistent. Based on our analyses, we calculated theoretical Ki values for trypsin with sunflower trypsin inhibitor (SFTI-1) variants, which yielded the more rigid Pro14 variant, with probably higher potency than the wild-type inhibitor. Moreover, a hirudin variant with an Arg1 and Trp3 is a promising basis for novel thrombin inhibitors with high potency. Further examples from antibody interaction and a cancer-related effector-receptor system demonstrate that our approach is applicable to protein interaction studies beyond the protease field.

Free Energy-Based Computational Methods for the Study of Protein-Peptide Binding Equilibria

This chapter discusses the theory and application of physics-based free energy methods to estimate protein-peptide binding free energies. It presents a statistical mechanics formulation of molecular binding, which is then specialized in three methodologies: (1) alchemical absolute binding free energy estimation with implicit solvation, (2) alchemical relative binding free energy estimation with explicit solvation, and (3) potential of mean force binding free energy estimation. Case studies of protein-peptide binding application taken from the recent literature are discussed for each method.

Large-Scale Molecular Dynamics Simulations of Cellular Compartments

Molecular dynamics or MD simulation is gradually maturing into a tool for constructing in vivo models of living cells in atomistic details. The feasibility of such models is bolstered by integrating the simulations with data from microscopic, tomographic and spectroscopic experiments on exascale supercomputers, facilitated by the use of deep learning technologies. Over time, MD simulation has evolved from tens of thousands of atoms to over 100 million atoms comprising an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium. In this chapter, we present a step-by-step outline for preparing, executing and analyzing such large-scale MD simulations of biological systems that are essential to life processes. All scripts are provided via GitHub.

An affinity prediction approach for the ligand of E3 ligase Cbl-b and an insight into substrate binding pattern

Protein-protein interactions (PPIs) are essentially fundamental to all cellular processes, so that developing small molecule inhibitors of PPIs have great significance despite representing a huge challenge. Studying PPIs with the help of peptide motifs could obtain the structural information and reference significance to reduce the difficulty in the development of small molecules. Computational methods are powerful tools to characterize peptide-protein interactions, especially molecular dynamics simulation and binding free energy calculation. Here, we established an affinity prediction model suitable for Casitas B lymphoma-b (Cbl-b) and phosphorylated motif system. According to the affinity data set of multiple truncated peptides, the force field, solvent model, and internal dielectric constant of molecular mechanics/generalized Born surface area (MM/GBSA) method were optimized. Further, we predicted the affinity of the rationally designed new sequences through this model and obtained a new 6-mer motif with a 7-fold increase in affinity and the comprehensive structure-activity relationship. Moreover, we proposed an insight of unexpected activity of the truncated 5-mer peptide and revealed the possible binding mode of the new highly active 6-mer motif by extended simulation. Our results showed that the activity enhancement of the truncated peptide was caused by the acetyl-mediated conformation change. The side chain of Arg and pTyr in the 6-mer motif co-occupied the site p1 to form numerous hydrogen bond interactions and increased hydrophobic interaction formed with Tyr266, leading to the higher affinity. The present work provided a reference to investigate the PPI of Cbl-b and phosphorylated substrates and guided the development of Cbl-b inhibitors.

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

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.

Implementing and Assessing an Alchemical Method for Calculating Protein-Protein Binding Free Energy

Protein-protein binding is fundamental to most biological processes. It is important to be able to use computation to accurately estimate the change in protein-protein binding free energy due to mutations in order to answer biological questions that would be experimentally challenging, laborious, or time-consuming. Although nonrigorous free-energy methods are faster, rigorous alchemical molecular dynamics-based methods are considerably more accurate and are becoming more feasible with the advancement of computer hardware and molecular simulation software. Even with sufficient computational resources, there are still major challenges to using alchemical free-energy methods for protein-protein complexes, such as generating hybrid structures and topologies, maintaining a neutral net charge of the system when there is a charge-changing mutation, and setting up the simulation. In the current study, we have used the pmx package to generate hybrid structures and topologies, and a double-system/single-box approach to maintain the net charge of the system. To test the approach, we predicted relative binding affinities for two protein-protein complexes using a nonequilibrium alchemical method based on the Crooks fluctuation theorem and compared the results with experimental values. The method correctly identified stabilizing from destabilizing mutations for a small protein-protein complex, and a larger, more challenging antibody complex. Strong correlations were obtained between predicted and experimental relative binding affinities for both protein-protein systems.

Bioinformatics and Biosimulations as Toolbox for Peptides and Peptidomimetics Design: Where Are We?

Peptides and peptidomimetics are strongly re-emerging as amenable candidates in the development of therapeutic strategies against a plethora of pathologies. In particular, these molecules are extremely suitable to treat diseases in which a major role is played by protein-protein interactions (PPIs). Unlike small organic compounds, peptides display both a high degree of specificity avoiding secondary off-targets effects and a relatively low degree of toxicity. Further advantages are provided by the possibility to easily conjugate peptides to functionalized nanoparticles, so improving their delivery and cellular uptake. In many cases, such molecules need to assume a specific three-dimensional conformation that resembles the bioactive one of the endogenous ligand. To this end, chemical modifications are introduced in the polypeptide chain to constrain it in a well-defined conformation, and to improve the drug-like properties. In this context, a successful strategy for peptide/peptidomimetics design and optimization is to combine different computational approaches ranging from structural bioinformatics to atomistic simulations. Here, we review the computational tools for peptide design, highlighting their main features and differences, and discuss selected protocols, among the large number of methods available, used to assess and improve the stability of the functional folding of the peptides. Finally, we introduce the simulation techniques employed to predict the binding affinity of the designed peptides for their targets.

Quantitative mapping of protein-peptide affinity landscapes using spectrally encoded beads

Transient, regulated binding of globular protein domains to Short Linear Motifs (SLiMs) in disordered regions of other proteins drives cellular signaling. Mapping the energy landscapes of these interactions is essential for deciphering and perturbing signaling networks but is challenging due to their weak affinities. We present a powerful technology (MRBLE-pep) that simultaneously quantifies protein binding to a library of peptides directly synthesized on beads containing unique spectral codes. Using MRBLE-pep, we systematically probe binding of calcineurin (CN), a conserved protein phosphatase essential for the immune response and target of immunosuppressants, to the PxIxIT SLiM. We discover that flanking residues and post-translational modifications critically contribute to PxIxIT-CN affinity and identify CN-binding peptides based on multiple scaffolds with a wide range of affinities. The quantitative biophysical data provided by this approach will improve computational modeling efforts, elucidate a broad range of weak protein-SLiM interactions, and revolutionize our understanding of signaling networks.

Large-scale docking predicts that sORF-encoded peptides may function through protein-peptide interactions in Arabidopsis thaliana

Several recent studies indicate that small Open Reading Frames (sORFs) embedded within multiple eukaryotic non-coding RNAs can be translated into bioactive peptides of up to 100 amino acids in size. However, the functional roles of the 607 Stress Induced Peptides (SIPs) previously identified from 189 Transcriptionally Active Regions (TARs) in Arabidopsis thaliana remain unclear. To provide a starting point for functional annotation of these plant-derived peptides, we performed a large-scale prediction of peptide binding sites on protein surfaces using coarse-grained peptide docking. The docked models were subjected to further atomistic refinement and binding energy calculations. A total of 530 peptide-protein pairs were successfully docked. In cases where a peptide encoded by a TAR is predicted to bind at a known ligand or cofactor-binding site within the protein, it can be assumed that the peptide modulates the ligand or cofactor-binding. Moreover, we predict that several peptides bind at protein-protein interfaces, which could therefore regulate the formation of the respective complexes. Protein-peptide binding analysis further revealed that peptides employ both their backbone and side chain atoms when binding to the protein, forming predominantly hydrophobic interactions and hydrogen bonds. In this study, we have generated novel predictions on the potential protein-peptide interactions in A. thaliana, which will help in further experimental validation.

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.

Assessment of a Single Decoupling Alchemical Approach for the Calculation of the Absolute Binding Free Energies of Protein-Peptide Complexes

The computational modeling of peptide inhibitors to target protein-protein binding interfaces is growing in interest as these are often too large, too shallow, and too feature-less for conventional small molecule compounds. Here, we present a rare successful application of an alchemical binding free energy method for the calculation of converged absolute binding free energies of a series of protein-peptide complexes. Specifically, we report the binding free energies of a series of cyclic peptides derived from the LEDGF/p75 protein to the integrase receptor of the HIV1 virus. The simulations recapitulate the effect of mutations relative to the wild-type binding motif of LEDGF/p75, providing structural, energetic and dynamical interpretations of the observed trends. The equilibration and convergence of the calculations are carefully analyzed. Convergence is aided by the adoption of a single-decoupling alchemical approach with implicit solvation, which circumvents the convergence difficulties of conventional double-decoupling protocols. We hereby present the single-decoupling methodology and critically evaluate its advantages and limitations. We also discuss some of the challenges and potential pitfalls of binding free energy calculations for complex molecular systems which have generally limited their applicability to the quantitative study of protein-peptide binding equilibria.

Amyloid-Like β-Aggregates as Force-Sensitive Switches in Fungal Biofilms and Infections

Cellular aggregation is an essential step in the formation of biofilms, which promote fungal survival and persistence in hosts. In many of the known yeast cell adhesion proteins, there are amino acid sequences predicted to form amyloid-like β-aggregates. These sequences mediate amyloid formation in vitro. In vivo, these sequences mediate a phase transition from a disordered state to a partially ordered state to create patches of adhesins on the cell surface. These β-aggregated protein patches are called adhesin nanodomains, and their presence greatly increases and strengthens cell-cell interactions in fungal cell aggregation. Nanodomain formation is slow (with molecular response in minutes and the consequences being evident for hours), and strong interactions lead to enhanced biofilm formation. Unique among functional amyloids, fungal adhesin β-aggregation can be triggered by the application of physical shear force, leading to cellular responses to flow-induced stress and the formation of robust biofilms that persist under flow. Bioinformatics analysis suggests that this phenomenon may be widespread. Analysis of fungal abscesses shows the presence of surface amyloids in situ, a finding which supports the idea that phase changes to an amyloid-like state occur in vivo. The amyloid-coated fungi bind the damage-associated molecular pattern receptor serum amyloid P component, and there may be a consequential modulation of innate immune responses to the fungi. Structural data now suggest mechanisms for the force-mediated induction of the phase change. We summarize and discuss evidence that the sequences function as triggers for protein aggregation and subsequent cellular aggregation, both in vitro and in vivo.

DINC 2.0: A New Protein-Peptide Docking Webserver Using an Incremental Approach

Molecular docking is a standard computational approach to predict binding modes of protein-ligand complexes by exploring alternative orientations and conformations of the ligand (i.e., by exploring ligand flexibility). Docking tools are largely used for virtual screening of small drug-like molecules, but their accuracy and efficiency greatly decays for ligands with more than 10 flexible bonds. This prevents a broader use of these tools to dock larger ligands, such as peptides, which are molecules of growing interest in cancer research. To overcome this limitation, our group has previously proposed a meta-docking strategy, called DINC, to predict binding modes of large ligands. By incrementally docking overlapping fragments of a ligand, DINC allowed predicting binding modes of peptide-based inhibitors of transcription factors involved in cancer. Here, we describe DINC 2.0, a revamped version of the DINC webserver with enhanced capabilities and a more user-friendly interface. DINC 2.0 allows docking ligands that were previously too challenging for DINC, such as peptides with more than 25 flexible bonds. The webserver is freely accessible at http://dinc.kavrakilab.org, together with additional documentation and video tutorials. Our team will provide continuous support for this tool and is working on extending its applicability to other challenging fields, such as personalized immunotherapy against cancer. Cancer Res; 77(21); e55-57. ©2017 AACR.

Molecular Mechanisms of Drug Action: X-ray Crystallography at the Basis of Structure-based and Ligand-based Drug Design

Biological systems are recognized for their complexity and diversity and yet we sometimes manage to cure disease via the administration of small chemical drug molecules. At first, active ingredients were found accidentally and at that time there did not seem a need to understand the molecular mechanism of drug functioning. However, the urge to develop new drugs, the discovery of multipurpose characteristics of some drugs, and the necessity to remove unwanted secondary drug effects, incited the pharmaceutical sector to rationalize drug design. This did not deliver success in the years directly following its conception, but it drove the evolution of biochemical and biophysical techniques to enable the characterization of molecular mechanisms of drug action. Functional and structural data generated by biochemists and structural biologists became a valuable input for computational biologists, chemists and bioinformaticians who could extrapolate in silico, based on variations in the structural aspects of the drug molecules and their target. This opened up new avenues with much improved predictive power because of a clearer perception of the role and impact of structural elements in the intrinsic affinity and specificity of the drug for its target. In this chapter, we review how crystal structures can initiate structure-based drug design in general.

Free Energy of Binding and Mechanism of Interaction for the MEEVD-TPR2A Peptide-Protein Complex

The association between the MEEVD C-terminal peptide from the heat shock protein 90 (Hsp90) and tetratricopeptide repeat A (TPR2A) domain of the heat shock organizing protein (Hop) is a useful prototype to study the fundamental molecular details about the Hop-Hsp90 interaction. We study here the mechanism of binding/unbinding and compute the standard binding free energy and potential of mean force for the association of the MEEVD peptide to the TPR2A domain using the Adaptive Biasing Force (ABF) methodology. We observe conformational changes of the peptide and the protein receptor induced by binding. We measure the binding free energy of -8.4 kcal/mol, which is consistent with experimental estimates. The simulations achieve multiple unbinding and rebinding events along a consistent pathway connecting the binding site to solvent. The MEEVD peptide slowly dissociates disrupting the hydrogen bonds first, then tilting on the side while preserving the interaction with the side chain of residue Asp 5 of the peptide. After this initial displacement, the peptide completely dissociates and moves into the solvent. Rebinding of the MEEVD peptide from the solvent to the receptor binding site occurs slowly through the portal of entry. Unbinding and rebinding go through intermediate states characterized by the peptide interacting with a lateral helix, helix A1, of the receptor with mainly Asp 5, Val 4, and Glu 3 of the peptide. This newly discovered intermediate structure is characterized by numerous contacts with the receptor which lead to complete formation of the bound complex. The structure of the bound complex obtained after rebinding is structurally very similar to the crystal structure of the complex (0.48 Å root-mean square deviation). The residues Asp 5, Val 4, and Glu 3 adopt conformations and intermolecular contacts with excellent structural similarity with the native ones. Finally, the dissociation and reassociation of MEEVD induce hydration/dehydration transitions, which provide insights on the role of desolvation and solvation processes in protein-peptide binding.

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.