Application of docking methodologies to modeled proteins

Protein docking is essential for structural characterization of protein interactions. Besides providing the structure of protein complexes, modeling of proteins and their complexes is important for understanding the fundamental principles and specific aspects of protein interactions. The accuracy of protein modeling, in general, is still less than that of the experimental approaches. Thus, it is important to investigate the applicability of docking techniques to modeled proteins. We present new comprehensive benchmark sets of protein models for the development and validation of protein docking, as well as a systematic assessment of free and template-based docking techniques on these sets. As opposed to previous studies, the benchmark sets reflect the real case modeling/docking scenario where the accuracy of the models is assessed by the modeling procedure, without reference to the native structure (which would be unknown in practical applications). We also expanded the analysis to include docking of protein pairs where proteins have different structural accuracy. The results show that, in general, the template-based docking is less sensitive to the structural inaccuracies of the models than the free docking. The near-native docking poses generated by the template-based approach, typically, also have higher ranks than those produces by the free docking (although the free docking is indispensable in modeling the multiplicity of protein interactions in a crowded cellular environment). The results show that docking techniques are applicable to protein models in a broad range of modeling accuracy. The study provides clear guidelines for practical applications of docking to protein models.     

Accounting for global protein deformability during protein-protein and protein-ligand docking

Computational docking methods are valuable tools aimed to simplify the costly process of drug development and improvement. Most current approaches assume a rigid receptor structure to allow virtual screening of large numbers of possible ligands and putative binding sites on a receptor molecule. However, inclusion of receptor flexibility can be of critical importance since binding of a ligand can lead to changes in the receptor protein conformation that are sterically necessary to accommodate a ligand. Recent approaches to efficiently account for receptor flexibility during docking simulations are reviewed. In particular, accounting efficiently for global conformational changes of the protein backbone during docking is a still challenging unsolved problem. An approximate method has recently been suggested that is based on relaxing the receptor conformation during docking in pre-calculated soft collective degrees of freedom (M. Zacharias, Rapid protein-ligand docking using soft modes from molecular dynamics simulations to account for protein deformability: binding of FK506 to FKBP, Proteins: Struct., Funct., Genet. 54 (2004) 759-767). Test applications on protein-protein docking and on docking the inhibitor staurosporine to the apo-form of cAMP-dependent protein kinase A catalytic domain indicate significant improvement of docking results compared to rigid docking at a very modest computational demand. Accounting for receptor conformational changes in pre-calculated global degrees of freedom might offer a promising route to improve systematic docking screening simulations.     

Application of molecular modeling to analysis of inhibition of kinesin motor proteins of the BimC subfamily by monastrol and related compounds

Application of molecular modeling approaches has potential to contribute to rational drug design. These approaches may be especially useful when attempting to elucidate the structural features associated with novel drug targets. In this study, molecular docking and molecular dynamics were applied to studies of inhibition of the human motor protein denoted HsEg5 and other homologues in the BimC subfamily. These proteins are essential for mitosis, so compounds that inhibit their activity may have potential as anticancer therapeutics. The discovery of a small-molecule cell-permeable inhibitor, monastrol, has stimulated research in this area. Interestingly, monastrol is reported to inhibit the human and Xenopus forms of Eg5, but not those from Drosophila and Aspergillus. In this study, homology modeling was used to generate models of the Xenopus, Drosophila, and Aspergillus homologues, using the crystal structure of the human protein in complex with monastrol as a template. A series of known inhibitors was docked into each of the homologues, and the differences in binding energies were consistent with reported experimental data. Molecular dynamics revealed significant changes in the structure of the Aspergillus homologue that may contribute to its relative insensitivity to monastrol and related compounds.     

Flexible ligand-flexible protein docking in protein kinase systems

Quite a few reviews on molecular docking have already appeared. This mini-review focuses on methods that incorporate protein flexibility in docking rather than those that treat protein targets as rigid molecules. This is still a challenging problem but there are encouraging recent advances. These methods will be reviewed particularly in light of their applications to protein kinases and phosphatases. In addition to obtaining correct docking pose, recent developments on exploring docking pathways are also highlighted.     

Molecular modeling of the interaction between heparan sulfate and cellular growth factors: bringing pieces together

Heparan sulfate is a polysaccharide belonging to the glycaminoglycan family. It interacts with numerous proteins of the extracellular matrix, in particular cellular growth factors. The number of experimental protein-heparin sulfate complexes obtained by crystallography or nuclear magnetic resonance is limited. Alternatively, computational approaches can be employed. Generally, they restrain the conformation of the glycosidic rings and linkages in order to reduce the complexity of the problem. Modeling the interaction between protein and heparan sulfate is indeed challenging because of the large size of the fragment needed for a strong binding, the flexibility brought by the glycosidic rings and linkages and the high density of negative charges. We propose a two-step method based on molecular docking and molecular dynamics simulation. Molecular docking allows exploring the positioning of a rigid heparin sulfate fragment on the protein surface. Molecular dynamics refine selected docking models by explicitly representing solvent molecules and not restraining the polysaccharide backbone. The interaction of a hexamer of heparin sulfate was studied in interaction with fibroblast growth factor 2 and stromal cell-derived factor 1α. This approach shed light on the plasticity of the growth factors interacting with heparan sulfate. This approach can be extended to the study of other protein/glycosaminoglycan complexes.     

Molecular docking of balanol to dynamics snapshots of protein kinase A

Even if the structure of a receptor has been determined experimentally, it may not be a conformation to which a ligand would bind when induced fit effects are significant. Molecular docking using such a receptor structure may thus fail to recognize a ligand to which the receptor can bind with reasonable affinity. Here, we examine one way to alleviate this problem by using an ensemble of receptor conformations generated from a molecular dynamics simulation for molecular docking. Two molecular dynamics simulations were conducted to generate snapshots for protein kinase A: one with the ligand bound, the other without. The ligand, balanol, was then docked to conformations of the receptors presented by these trajectories. The Lamarckian genetic algorithm in Autodock [Goodsell et al. J Mol Recognit 1996;9(1):1-5; Morris et al. J Comput Chem 1998;19(14):1639-1662] was used in the docking. Three ligand models were used: rigid, flexible, and flexible with torsional potentials. When the snapshots were taken from the molecular dynamics simulation of the protein-ligand complex, the correct docking structure could be recovered easily by the docking algorithm in all cases. This was an easier case for challenging the docking algorithm because, by using the structure of the protein in a protein-ligand complex, one essentially assumed that the protein already had a pocket to which the ligand can fit well. However, when the snapshots were taken from the ligand-free protein simulation, which is more useful for a practical application when the structure of the protein-ligand complex is not known, several clusters of structures were found. Of the 10 docking runs for each snapshot, at least one structure was close to the correctly docked structure when the flexible-ligand models were used. We found that a useful way to identify the correctly docked structure was to locate the structure that appeared most frequently as the lowest energy structure in the docking experiments to different snapshots.     

Optimized grid-based protein-protein docking as a global search tool followed by incorporating experimentally derivable restraints

Unbound protein docking, or the computational prediction of the structure of a protein complex from the structures of its separated components, is of importance but still challenging. A practical approach toward reliable results for unbound docking is to incorporate experimentally derived information with computation. To this end, truly systematic search of the global docking space is desirable. The fast Fourier transform (FFT) docking is a systematic search method with high computational efficiency. However, by using FFT to perform unbound docking, possible conformational changes upon binding must be treated implicitly. To better accommodate the implicit treatment of conformational flexibility, we develop a rational approach to optimize "softened" parameters for FFT docking. In connection with the increased "softness" of the parameters in this global search step, we use a revised rule to select candidate models from the search results. For complexes designated as of low and medium difficulty for unbound docking, these adaptations of the original FTDOCK program lead to substantial improvements of the global search results. Finally, we show that models resulted from FFT-based global search can be further filtered with restraints derivable from nuclear magnetic resonance (NMR) chemical shift perturbation or mutagenesis experiments, leading to a small set of models that can be feasibly refined and evaluated using computationally more expensive methods and that still include high-ranking near-native conformations.     

DynaDock: A new molecular dynamics‐based algorithm for protein–peptide docking including receptor flexibility

Molecular docking programs play an important role in drug development and many well‐established methods exist. However, there are two situations for which the performance of most approaches is still not satisfactory, namely inclusion of receptor flexibility and docking of large, flexible ligands like peptides. In this publication a new approach is presented for docking peptides into flexible receptors. For this purpose a two step procedure was developed: first, the protein–peptide conformational space is scanned and approximate ligand poses are identified and second, the identified ligand poses are refined by a new molecular dynamics‐based method, optimized potential molecular dynamics (OPMD). The OPMD approach uses soft‐core potentials for the protein–peptide interactions and applies a new optimization scheme to the soft‐core potential. Comparison with refinement results obtained by conventional molecular dynamics and a soft‐core scaling approach shows significant improvements in the sampling capability for the OPMD method. Thus, the number of starting poses needed for successful refinement is much lower than for the other methods. The algorithm was evaluated on 15 protein–peptide complexes with 2–16mer peptides. Docking poses with peptide RMSD values <2.10 Å from the equilibrated experimental structures were obtained in all cases. For four systems docking into the unbound receptor structures was performed, leading to peptide RMSD values <2.12 Å. Using a specifically fitted scoring function in 11 of 15 cases the best scoring poses featured a peptide RMSD ≤2.10 Å. Proteins 2010. © 2009 Wiley‐Liss, Inc.     

Incorporating replacement free energy of binding‐site waters in molecular docking

Binding‐site water molecules play a crucial role in protein‐ligand recognition, either being displaced upon ligand binding or forming water bridges to stabilize the complex. However, rigorously treating explicit binding‐site waters is challenging in molecular docking, which requires to fully sample ensembles of waters and to consider the free energy cost of replacing waters. Here, we describe a method to incorporate structural and energetic properties of binding‐site waters into molecular docking. We first developed a solvent property analysis (SPA) program to compute the replacement free energies of binding‐site water molecules by post‐processing molecular dynamics trajectories obtained from ligand‐free protein structure simulation in explicit water. Next, we implemented a distance‐dependent scoring term into DOCK scoring function to take account of the water replacement free energy cost upon ligand binding. We assessed this approach in protein targets containing important binding‐site waters, and we demonstrated that our approach is reliable in reproducing the crystal binding geometries of protein‐ligand‐water complexes, as well as moderately improving the ligand docking enrichment performance. In addition, SPA program (free available to academic users upon request) may be applied in identifying hot‐spot binding‐site residues and structure‐based lead optimization. Proteins 2014; 82:1765–1776. © 2014 Wiley Periodicals, Inc.     

Ensemble docking of multiple protein structures: considering protein structural variations in molecular docking

One approach to incorporate protein flexibility in molecular docking is the use of an ensemble consisting of multiple protein structures. Sequentially docking each ligand into a large number of protein structures is computationally too expensive to allow large-scale database screening. It is challenging to achieve a good balance between docking accuracy and computational efficiency. In this work, we have developed a fast, novel docking algorithm utilizing multiple protein structures, referred to as ensemble docking, to account for protein structural variations. The algorithm can simultaneously dock a ligand into an ensemble of protein structures and automatically select an optimal protein structure that best fits the ligand by optimizing both ligand coordinates and the conformational variable m, where m represents the m-th structure in the protein ensemble. The docking algorithm was validated on 10 protein ensembles containing 105 crystal structures and 87 ligands in terms of binding mode and energy score predictions. A success rate of 93% was obtained with the criterion of root-mean-square deviation <2.5 A if the top five orientations for each ligand were considered, comparable to that of sequential docking in which scores for individual docking are merged into one list by re-ranking, and significantly better than that of single rigid-receptor docking (75% on average). Similar trends were also observed in binding score predictions and enrichment tests of virtual database screening. The ensemble docking algorithm is computationally efficient, with a computational time comparable to that for docking a ligand into a single protein structure. In contrast, the computational time for the sequential docking method increases linearly with the number of protein structures in the ensemble. The algorithm was further evaluated using a more realistic ensemble in which the corresponding bound protein structures of inhibitors were excluded. The results show that ensemble docking successfully predicts the binding modes of the inhibitors, and discriminates the inhibitors from a set of noninhibitors with similar chemical properties. Although multiple experimental structures were used in the present work, our algorithm can be easily applied to multiple protein conformations generated by computational methods, and helps improve the efficiency of other existing multiple protein structure(MPS)-based methods to accommodate protein flexibility.     

FitEM2EM--tools for low resolution study of macromolecular assembly and dynamics

Studies of the structure and dynamics of macromolecular assemblies often involve comparison of low resolution models obtained using different techniques such as electron microscopy or atomic force microscopy. We present new computational tools for comparing (matching) and docking of low resolution structures, based on shape complementarity. The matched or docked objects are represented by three dimensional grids where the value of each grid point depends on its position with regard to the interior, surface or exterior of the object. The grids are correlated using fast Fourier transformations producing either matches of related objects or docking models depending on the details of the grid representations. The procedures incorporate thickening and smoothing of the surfaces of the objects which effectively compensates for differences in the resolution of the matched/docked objects, circumventing the need for resolution modification. The presented matching tool FitEM2EMin successfully fitted electron microscopy structures obtained at different resolutions, different conformers of the same structure and partial structures, ranking correct matches at the top in every case. The differences between the grid representations of the matched objects can be used to study conformation differences or to characterize the size and shape of substructures. The presented low-to-low docking tool FitEM2EMout ranked the expected models at the top.     

Unravelling the dynamics of organisms in a changing world using ecological modelling

Understanding and predicting the dynamics of organisms is a central objective in ecology and conservation biology, and modelling provides a solution to tackling this problem. However, the complex nature of ecological systems means that for a thorough understanding of ecological dynamics at hierarchical scales, a set of modeling approaches need to be adopted. This review illustrates how modelling approaches can be used to understand the dynamics of organisms in applied ecological problems, focussing on mechanistic models at a local scale and statistical models at a broad scale. Mechanistic models incorporate ecological processes explicitly and thus are likely to be robust under novel conditions. Models based on behavioural decisions by individuals represent a typical example of the successful application of mechanistic models to applied problems. Considering the data-hungry nature of such mechanistic models, model complexity and parameterisation need to be explored further for a quick and widespread implementation of this model type. For broad-scale phenomena, statistical models play an important role in dealing with problems that are often inherent in data. Examples include models for quantifying population trends from long-term, large-scale data and those for comparative methods of extinction risk. Novel statistical approaches also allow mechanistic models to be parameterised using readily obtained data at a macro scale. In conclusion, the complementary use and improvement of multiple model types, the increased use of novel model parameterisation, the examination of model transferability and the achievement of wider biodiversity information availability are key challenges for the effective use of modelling in applied ecological problems.     

Biased coarse-grained molecular dynamics simulation approach for flexible fitting of X-ray structure into cryo electron microscopy maps

Several approaches have been introduced to interpret, in terms of high-resolution structure, low-resolution structural data as obtained from cryo-EM. As conformational changes are often observed in biological molecules, these techniques need to take into account the flexibility of proteins. Flexibility has been described in terms of movement between rigid domains and between rigid secondary structure elements, which present some limitations for studying dynamical properties. Normal mode analysis has also been used, but is limited to medium resolution data. All-atom molecular dynamics fitting techniques are more appropriate to fit structures into higher-resolution data as full protein flexibility is considered, but are cumbersome in terms of computational time. Here, we introduce a coarse-grained approach; a Go-model was used to represent biological molecules, combined with biased molecular dynamics to reproduce accurately conformational transitions. Illustrative examples on simulated data are shown. Accurate fittings can be obtained for resolution ranging from 5 to 20 Å. The approach was also tested on experimental data of Elongation Factor G and Escherichia coli RNA polymerase, where its validity is compared to previous models obtained from different techniques. This comparison demonstrates that quantitative flexible techniques, as opposed to manual docking, need to be considered to interpret low-resolution data.     

Automatic prediction of flexible regions improves the accuracy of protein-protein docking models

Computational models of protein-protein docking that incorporate backbone flexibility can predict perturbations of the backbone and side chains during docking and produce protein interaction models with atomic accuracy. Most previous models usually predefine flexible regions by visually comparing the bound and unbound structures. In this paper, we propose a general method to automatically identify the flexible hinges for domain assembly and the flexible loops for loop refinement, in addition to predicting the corresponding movements of the identified active residues. We conduct experiments to evaluate performance of our approach on two test sets. Comparison of results on test set I between algorithms with and without prediction of flexible regions demonstrate the superior recovery of energy funnels in many target interactions using the new loop refinement model. In addition, our decoys are superior for each target. Indeed, the total number of satisfactory models is almost double that of other programs. The results on test set II docking tests produced by our domain assembly method also show encouraging results. Of the three targets examined, one exhibits energy funnel and the best models of the other two targets all meet the conditions of acceptable accuracy. Results demonstrate that the automatic prediction of flexible backbone regions can greatly improve the performance of protein-protein docking models.     

Implicit solvent models for flexible protein-protein docking by molecular dynamics simulation

The suitability of three implicit solvent models for flexible protein-protein docking by procedures using molecular dynamics simulation is investigated. The three models are (i) the generalized Born (GB) model implemented in the program AMBER6.0; (ii) a distance-dependent dielectric (DDD) model; and (iii) a surface area-dependent model that we have parameterized and call the NPSA model. This is a distance-dependent dielectric model modified by neutralizing the ionizable side-chains and adding a surface area-dependent solvation term. These solvent models were first tested in molecular dynamics simulations at 300 K of the native structures of barnase, barstar, segment B1 of protein G, and three WW domains. These protein structures display a range of secondary structure contents and stabilities. Then, to investigate the performance of the implicit solvent models in protein docking, molecular dynamics simulations of barnase/barstar complexation, as well as PIN1 WW domain/peptide complexation, were conducted, starting from separated unbound structures. The simulations show that the NPSA model has significant advantages over the DDD and GB models in maintaining the native structures of the proteins and providing more accurate docked complexes.     

Docking proteins and peptides under evolutionary constraints in Critical Assessment of PRediction of Interactions rounds 38 to 45

Computational structural prediction of macromolecular interactions is a fundamental tool toward the global understanding of cellular processes. The Critical Assessment of PRediction of Interactions (CAPRI) community-wide experiment provides excellent opportunities for blind testing computational docking methods and includes original targets, thus widening the range of docking applications. Our participation in CAPRI rounds 38 to 45 enabled us to expand the way we include evolutionary information in structural predictions beyond our standard free docking InterEvDock pipeline. InterEvDock integrates a coarse-grained potential that accounts for interface coevolution based on joint multiple sequence alignments of two protein partners (co-alignments). However, even though such co-alignments could be built for none of the CAPRI targets in rounds 38 to 45, including host-pathogen and protein-oligosaccharide complexes and a redesigned interface, we identified multiple strategies that can be used to incorporate evolutionary constraints, which helped us to identify the most likely macromolecular binding modes. These strategies include template-based modeling where only local adjustments should be applied when query-template sequence identity is above 30% and larger perturbations are needed below this threshold; covariation-based structure prediction for individual protein partners; and the identification of evolutionarily conserved and structurally recurrent anchoring interface motifs. Overall, we submitted correct predictions among the top 5 models for 12 out of 19 interface challenges, including four High- and five Medium-quality predictions. Our top 20 models included correct predictions for three out of the five targets we missed in the top 5, including two targets for which misleading biological data led us to downgrade correct free docking models.     

Towards the prediction of protein interaction partners using physical docking

Deciphering the whole network of protein interactions for a given proteome (‘interactome’) is the goal of many experimental and computational efforts in Systems Biology. Separately the prediction of the structure of protein complexes by docking methods is a well‐established scientific area. To date, docking programs have not been used to predict interaction partners. We provide a proof of principle for such an approach. Using a set of protein complexes representing known interactors in their unbound form, we show that a standard docking program can distinguish the true interactors from a background of 922 non‐redundant potential interactors. We additionally show that true interactions can be distinguished from non‐likely interacting proteins within the same structural family. Our approach may be put in the context of the proposed ‘funnel‐energy model’; the docking algorithm may not find the native complex, but it distinguishes binding partners because of the higher probability of favourable models compared with a collection of non‐binders. The potential exists to develop this proof of principle into new approaches for predicting interaction partners and reconstructing biological networks.     

Computational approaches for protein function prediction: A combined strategy from multiple sequence alignment to molecular docking-based virtual screening

The functional characterization of proteins represents a daily challenge for biochemical, medical and computational sciences. Although finally proved on the bench, the function of a protein can be successfully predicted by computational approaches that drive the further experimental assays. Current methods for comparative modeling allow the construction of accurate 3D models for proteins of unknown structure, provided that a crystal structure of a homologous protein is available. Binding regions can be proposed by using binding site predictors, data inferred from homologous crystal structures, and data provided from a careful interpretation of the multiple sequence alignment of the investigated protein and its homologs. Once the location of a binding site has been proposed, chemical ligands that have a high likelihood of binding can be identified by using ligand docking and structure-based virtual screening of chemical libraries. Most docking algorithms allow building a list sorted by energy of the lowest energy docking configuration for each ligand of the library. In this review the state-of-the-art of computational approaches in 3D protein comparative modeling and in the study of protein–ligand interactions is provided. Furthermore a possible combined/concerted multistep strategy for protein function prediction, based on multiple sequence alignment, comparative modeling, binding region prediction, and structure-based virtual screening of chemical libraries, is described by using suitable examples. As practical examples, Abl-kinase molecular modeling studies, HPV-E6 protein multiple sequence alignment analysis, and some other model docking-based characterization reports are briefly described to highlight the importance of computational approaches in protein function prediction.     

Evolutionary models of extended phenotypes

A variety of theoretical models incorporate phenotypes expressed in the external environment, but a core question is whether such traits generate dynamics that alter evolution. This has proven to be a challenging and controversial proposition. However, several recent modelling frameworks provide insight: indirect genetic effect (IGE) models, niche construction models, and evolutionary feedback models. These distinct approaches converge upon the observation that gene action at a distance generates feedback that expands the range of trait values and evolutionary rates that we should expect to observe in empirical studies. Such conceptual replication provides solid evidence that traits with extended effects have important evolutionary consequences, but more empirical work is needed to evaluate the predictive power of different modelling approaches.