Integrative modeling of diverse protein-peptide systems using CABS-dock.

The CABS model can be applied to a wide range of protein-protein and protein-peptide molecular modeling tasks, such as simulating folding pathways, predicting structures, docking, and analyzing the structural dynamics of molecular complexes. In this work, we use the CABS-dock tool in two diverse modeling tasks: 1) predicting the structures of amyloid protofilaments and 2) identifying cleavage sites in the peptide substrates of proteolytic enzymes. In the first case, simulations of the simultaneous docking of amyloidogenic peptides indicated that the CABS model can accurately predict the structures of amyloid protofilaments which have an in-register parallel architecture. Scoring based on a combination of symmetry criteria and estimated interaction energy values for bound monomers enables the identification of protofilament models that closely match their experimental structures for 5 out of 6 analyzed systems. For the second task, it has been shown that CABS-dock coarse-grained docking simulations can be used to identify the positions of cleavage sites in the peptide substrates of proteolytic enzymes. The cleavage site position was correctly identified for 12 out of 15 analyzed peptides. When combined with sequence-based methods, these docking simulations may lead to an efficient way of predicting cleavage sites in degraded proteins. The method also provides the atomic structures of enzyme-substrate complexes, which can give insights into enzyme-substrate interactions that are crucial for the design of new potent inhibitors.

Protein-Protein Docking with Large-Scale Backbone Flexibility Using Coarse-Grained Monte-Carlo Simulations.

Most of the protein-protein docking methods treat proteins as almost rigid objects. Only the side-chains flexibility is usually taken into account. The few approaches enabling docking with a flexible backbone typically work in two steps, in which the search for protein-protein orientations and structure flexibility are simulated separately. In this work, we propose a new straightforward approach for docking sampling. It consists of a single simulation step during which a protein undergoes large-scale backbone rearrangements, rotations, and translations. Simultaneously, the other protein exhibits small backbone fluctuations. Such extensive sampling was possible using the CABS coarse-grained protein model and Replica Exchange Monte Carlo dynamics at a reasonable computational cost. In our proof-of-concept simulations of 62 protein-protein complexes, we obtained acceptable quality models for a significant number of cases.

Protocols for All-Atom Reconstruction and High-Resolution Refinement of Protein-Peptide Complex Structures.

Structural characterizations of protein-peptide complexes may require further improvements. These may include reconstruction of missing atoms and/or structure optimization leading to higher accuracy models. In this work, we describe a workflow that generates accurate structural models of peptide-protein complexes starting from protein-peptide models in C-alpha representation generated using CABS-dock molecular docking. First, protein-peptide models are reconstructed from their C-alpha traces to all-atom representation using MODELLER. Next, they are refined using Rosetta FlexPepDock. The described workflow allows for reliable all-atom reconstruction of CABS-dock models and their further improvement to high-resolution models.

Protocols for Fast Simulations of Protein Structure Flexibility Using CABS-Flex and SURPASS.

Conformational flexibility of protein structures can play an important role in protein function. The flexibility is often studied using computational methods since experimental characterization can be difficult. Depending on protein system size, computational tools may require large computational resources or significant simplifications in the modeled systems to speed up calculations. In this work, we present the protocols for efficient simulations of flexibility of folded protein structures that use coarse-grained simulation tools of different resolutions: medium, represented by CABS-flex, and low, represented by SUPRASS. We test the protocols using a set of 140 globular proteins and compare the results with structure fluctuations observed in MD simulations, ENM modeling, and NMR ensembles. As demonstrated, CABS-flex predictions show high correlation to experimental and MD simulation data, while SURPASS is less accurate but promising in terms of future developments.

Computational reconstruction of atomistic protein structures from coarse-grained models.

Three-dimensional protein structures, whether determined experimentally or theoretically, are often too low resolution. In this mini-review, we outline the computational methods for protein structure reconstruction from incomplete coarse-grained to all atomistic models. Typical reconstruction schemes can be divided into four major steps. Usually, the first step is reconstruction of the protein backbone chain starting from the C-alpha trace. This is followed by side-chains rebuilding based on protein backbone geometry. Subsequently, hydrogen atoms can be reconstructed. Finally, the resulting all-atom models may require structure optimization. Many methods are available to perform each of these tasks. We discuss the available tools and their potential applications in integrative modeling pipelines that can transfer coarse-grained information from computational predictions, or experiment, to all atomistic structures.

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.

CABS-dock standalone: a toolbox for flexible protein-peptide docking.

SUMMARY: CABS-dock standalone is a multiplatform Python package for protein-peptide docking with backbone flexibility. The main feature of the CABS-dock method is its ability to simulate significant backbone flexibility of the entire protein-peptide system in a reasonable computational time. In the default mode, the package runs a simulation of fully flexible peptide searching for a binding site on the surface of a flexible protein receptor. The flexibility level of the molecules may be defined by the user. Furthermore, the CABS-dock standalone application provides users with full control over the docking simulation from the initial setup to the analysis of results. The standalone version is an upgrade of the original web server implementation-it introduces a number of customizable options, provides support for large-sized systems and offers a framework for deeper analysis of docking results. AVAILABILITY AND IMPLEMENTATION: CABS-dock standalone is distributed under the MIT licence, which is free for academic and non-profit users. It is implemented in Python and Fortran. The CABS-dock standalone source code, wiki with documentation and examples of use and installation instructions for Linux, macOS and Windows are available in the CABS-dock standalone repository at https://bitbucket.org/lcbio/cabsdock.

Modeling of Disordered Protein Structures Using Monte Carlo Simulations and Knowledge-Based Statistical Force Fields.

The description of protein disordered states is important for understanding protein folding mechanisms and their functions. In this short review, we briefly describe a simulation approach to modeling protein interactions, which involve disordered peptide partners or intrinsically disordered protein regions, and unfolded states of globular proteins. It is based on the CABS coarse-grained protein model that uses a Monte Carlo (MC) sampling scheme and a knowledge-based statistical force field. We review several case studies showing that description of protein disordered states resulting from CABS simulations is consistent with experimental data. The case studies comprise investigations of protein⁻peptide binding and protein folding processes. The CABS model has been recently made available as the simulation engine of multiscale modeling tools enabling studies of protein⁻peptide docking and protein flexibility. Those tools offer customization of the modeling process, driving the conformational search using distance restraints, reconstruction of selected models to all-atom resolution, and simulation of large protein systems in a reasonable computational time. Therefore, CABS can be combined in integrative modeling pipelines incorporating experimental data and other modeling tools of various resolution.

CABS-flex standalone: a simulation environment for fast modeling of protein flexibility.

SUMMARY: CABS-flex standalone is a Python package for fast simulations of protein structure flexibility. The package combines simulations of protein dynamics using CABS coarse-grained protein model with the reconstruction of selected models to all-atom representation and analysis of modeling results. CABS-flex standalone is designed to allow for command-line access to the CABS computations and complete control over simulation process. CABS-flex standalone is equipped with features such as: modeling of multimeric and large-size protein systems, contact map visualizations, analysis of similarities to the reference structure and configurable modeling protocol. For instance, the user may modify the simulation parameters, distance restraints, structural clustering scheme or all-atom reconstruction parameters. With these features CABS-flex standalone can be easily incorporated into other methodologies of structural biology. AVAILABILITY AND IMPLEMENTATION: CABS-flex standalone is distributed under the MIT license, which is free for academic and non-profit users. It is implemented in Python. CABS-flex source code, wiki with examples of use and installation instructions for Linux, macOS and Windows are available from the CABS-flex standalone repository at https://bitbucket.org/lcbio/cabsflex.

Protein-peptide docking using CABS-dock and contact information.

CABS-dock is a computational method for protein-peptide molecular docking that does not require predefinition of the binding site. The peptide is treated as fully flexible, while the protein backbone undergoes small fluctuations and, optionally, large-scale rearrangements. Here, we present a specific CABS-dock protocol that enhances the docking procedure using fragmentary information about protein-peptide contacts. The contact information is used to narrow down the search for the binding peptide pose to the proximity of the binding site. We used information on a single-chosen and randomly chosen native protein-peptide contact to validate the protocol on the peptiDB benchmark. The contact information significantly improved CABS-dock performance. The protocol has been made available as a new feature of the CABS-dock web server (at http://biocomp.chem.uw.edu.pl/CABSdock/). SHORT ABSTRACT: CABS-dock is a tool for flexible docking of peptides to proteins. In this article, we present a protocol for CABS-dock docking driven by information about protein-peptide contact(s). Using information on individual protein-peptide contacts allows to improve the accuracy of CABS-dock docking.

Modeling of Protein Structural Flexibility and Large-Scale Dynamics: Coarse-Grained Simulations and Elastic Network Models.

Fluctuations of protein three-dimensional structures and large-scale conformational transitions are crucial for the biological function of proteins and their complexes. Experimental studies of such phenomena remain very challenging and therefore molecular modeling can be a good alternative or a valuable supporting tool for the investigation of large molecular systems and long-time events. In this minireview, we present two alternative approaches to the coarse-grained (CG) modeling of dynamic properties of protein systems. We discuss two CG representations of polypeptide chains used for Monte Carlo dynamics simulations of protein local dynamics and conformational transitions, and highly simplified structure-based elastic network models of protein flexibility. In contrast to classical all-atom molecular dynamics, the modeling strategies discussed here allow the quite accurate modeling of much larger systems and longer-time dynamic phenomena. We briefly describe the main features of these models and outline some of their applications, including modeling of near-native structure fluctuations, sampling of large regions of the protein conformational space, or possible support for the structure prediction of large proteins and their complexes.

Role of Resultant Dipole Moment in Mechanical Dissociation of Biological Complexes.

Protein-peptide interactions play essential roles in many cellular processes and their structural characterization is the major focus of current experimental and theoretical research. Two decades ago, it was proposed to employ the steered molecular dynamics (SMD) to assess the strength of protein-peptide interactions. The idea behind using SMD simulations is that the mechanical stability can be used as a promising and an efficient alternative to computationally highly demanding estimation of binding affinity. However, mechanical stability defined as a peak in force-extension profile depends on the choice of the pulling direction. Here we propose an uncommon choice of the pulling direction along resultant dipole moment (RDM) vector, which has not been explored in SMD simulations so far. Using explicit solvent all-atom MD simulations, we apply SMD technique to probe mechanical resistance of ligand-receptor system pulled along two different vectors. A novel pulling direction-when ligand unbinds along the RDM vector-results in stronger forces compared to commonly used ligand unbinding along center of masses vector. Our observation that RDM is one of the factors influencing the mechanical stability of protein-peptide complex can be used to improve the ranking of binding affinities by using mechanical stability as an effective scoring function.

Kinetics and mechanical stability of the fibril state control fibril formation time of polypeptide chains: A computational study.

Fibril formation resulting from protein misfolding and aggregation is a hallmark of several neurodegenerative diseases such as Alzheimer's and Parkinson's diseases. Despite much progress in the understanding of the protein aggregation process, the factors governing fibril formation rates and fibril stability have not been fully understood. Using lattice models, we have shown that the fibril formation time is controlled by the kinetic stability of the fibril state but not by its energy. Having performed all-atom explicit solvent molecular dynamics simulations with the GROMOS43a1 force field for full-length amyloid beta peptides Aß40 and Aß42 and truncated peptides, we demonstrated that kinetic stability can be accessed via mechanical stability in such a way that the higher the mechanical stability or the kinetic stability, the faster the fibril formation. This result opens up a new way for predicting fibril formation rates based on mechanical stability that may be easily estimated by steered molecular dynamics.

CABS-flex 2.0: a web server for fast simulations of flexibility of protein structures.

Classical simulations of protein flexibility remain computationally expensive, especially for large proteins. A few years ago, we developed a fast method for predicting protein structure fluctuations that uses a single protein model as the input. The method has been made available as the CABS-flex web server and applied in numerous studies of protein structure-function relationships. Here, we present a major update of the CABS-flex web server to version 2.0. The new features include: extension of the method to significantly larger and multimeric proteins, customizable distance restraints and simulation parameters, contact maps and a new, enhanced web server interface. CABS-flex 2.0 is freely available at http://biocomp.chem.uw.edu.pl/CABSflex2.

Protein-peptide docking: opportunities and challenges.

Peptides have recently attracted much attention as promising drug candidates. Rational design of peptide-derived therapeutics usually requires structural characterization of the underlying protein-peptide interaction. Given that experimental characterization can be difficult, reliable computational tools are needed. In recent years, a variety of approaches have been developed for 'protein-peptide docking', that is, predicting the structure of the protein-peptide complex, starting from the protein structure and the peptide sequence, including variable degrees of information about the peptide binding site and/or conformation. In this review, we provide an overview of protein-peptide docking methods and outline their capabilities, limitations, and applications in structure-based drug design. Key challenges are also briefly discussed, such as modeling of large-scale conformational changes upon binding, scoring of predicted models, and optimal inclusion of varied types of experimental data and theoretical predictions into an integrative modeling process.

Coarse-Grained Modeling of the Interplay between Secondary Structure Propensities and Protein Fold Assembly.

We recently developed a new coarse-grained model of protein structure and dynamics [ Dawid et al. J. Chem. Theory Comput. 2017 , 13 ( 11 ), 5766 - 5779 ]. The model assumed a single bead representation of amino acid residues, where positions of such united residues were defined by centers of mass of four amino acid fragments. Replica exchange Monte Carlo sampling of the model chains provided good pictures of modeled structures and their dynamics. In its generic form the statistical knowledge-based force field of the model has been dedicated for single-domain globular proteins. Sequence-specific interactions are defined by three-letter secondary structure data. In the present work we demonstrate that different assignments and/or predictions of secondary structures are usually sufficient for enforcing cooperative formation of native-like folds of SURPASS chains for the majority of single-domain globular proteins. Simulations of native-like structure assembly for a representative set of globular proteins have shown that the accuracy of secondary structure data is usually not crucial for model performance, although some specific errors can strongly distort the obtained three-dimensional structures.

SURPASS Low-Resolution Coarse-Grained Protein Modeling.

Coarse-grained modeling of biomolecules has a very important role in molecular biology. In this work we present a novel SURPASS (Single United Residue per Pre-Averaged Secondary Structure fragment) model of proteins that can be an interesting alternative for existing coarse-grained models. The design of the model is unique and strongly supported by the statistical analysis of structural regularities characteristic for protein systems. Coarse-graining of protein chain structures assumes a single center of interactions per residue and accounts for preaveraged effects of four adjacent residue fragments. Knowledge-based statistical potentials encode complex interaction patterns of these fragments. Using the Replica Exchange Monte Carlo sampling scheme and a generic version of the SURPASS force field we performed test simulations of a representative set of single-domain globular proteins. The method samples a significant part of conformational space and reproduces protein structures, including native-like, with surprisingly good accuracy. Future extension of the SURPASS model on large biomacromolecular systems is briefly discussed.

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.