Guiding protein docking with geometric and evolutionary information

Structural modeling of molecular assemblies promises to improve our understanding of molecular interactions and biological function. Even when focusing on modeling structures of protein dimers from knowledge of monomeric native structure, docking two rigid structures onto one another entails exploring a large configurational space. This paper presents a novel approach for docking protein molecules and elucidating native-like configurations of protein dimers. The approach makes use of geometric hashing to focus the docking of monomeric units on geometrically complementary regions through rigid-body transformations. This geometry-based approach improves the feasibility of searching the combined configurational space. The search space is narrowed even further by focusing the sought rigid-body transformations around molecular surface regions composed of amino acids with high evolutionary conservation. This condition is based on recent findings, where analysis of protein assemblies reveals that many functional interfaces are significantly conserved throughout evolution. Different search procedures are employed in this work to search the resulting narrowed configurational space. A proof-of-concept energy-guided probabilistic search procedure is also presented. Results are shown on a broad list of 18 protein dimers and additionally compared with data reported by other labs. Our analysis shows that focusing the search around evolutionary-conserved interfaces results in lower lRMSDs.     

Scoring docked conformations generated by rigid-body protein-protein docking

Rigid-body methods, particularly Fourier correlation techniques, are very efficient for docking bound (co-crystallized) protein conformations using measures of surface complementarity as the target function. However, when docking unbound (separately crystallized) conformations, the method generally yields hundreds of false positive structures with good scores but high root mean square deviations (RMSDs). This paper describes a two-step scoring algorithm that can discriminate near-native conformations (with less than 5 A RMSD) from other structures. The first step includes two rigid-body filters that use the desolvation free energy and the electrostatic energy to select a manageable number of conformations for further processing, but are unable to eliminate all false positives. Complete discrimination is achieved in the second step that minimizes the molecular mechanics energy of the retained structures, and re-ranks them with a combined free-energy function which includes electrostatic, solvation, and van der Waals energy terms. After minimization, the improved fit in near-native complex conformations provides the free-energy gap required for discrimination. The algorithm has been developed and tested using docking decoys, i.e., docked conformations generated by Fourier correlation techniques. The decoy sets are available on the web for testing other discrimination procedures. Proteins 2000;40:525-537.     

Fragment-Based flexible ligand docking by evolutionary optimization.

A new computational approach for the efficient docking of flexible ligands in a rigid protein is presented. It exploits the binding modes of functional groups determined by an exhaustive search with solvation. The search in ligand conformational space is performed by a genetic algorithm whose scoring function approximates steric effects and intermolecular hydrogen bonds. Ligand conformations generated by the genetic algorithm are docked in the protein binding site by optimizing the fit of their fragments to optimal positions of chemically related functional groups. We show that the use of optimal binding modes of molecular fragments allows to dock known inhibitors with about ten rotatable bonds in the active site of the uncomplexed and complexed conformations of thrombin and HIV-1 protease.     

Protein models docking benchmark 2

Structural characterization of protein–protein interactions is essential for our ability to understand life processes. However, only a fraction of known proteins have experimentally determined structures. Such structures provide templates for modeling of a large part of the proteome, where individual proteins can be docked by template‐free or template‐based techniques. Still, the sensitivity of the docking methods to the inherent inaccuracies of protein models, as opposed to the experimentally determined high‐resolution structures, remains largely untested, primarily due to the absence of appropriate benchmark set(s). Structures in such a set should have predefined inaccuracy levels and, at the same time, resemble actual protein models in terms of structural motifs/packing. The set should also be large enough to ensure statistical reliability of the benchmarking results. We present a major update of the previously developed benchmark set of protein models. For each interactor, six models were generated with the model‐to‐native C^α RMSD in the 1 to 6 Å range. The models in the set were generated by a new approach, which corresponds to the actual modeling of new protein structures in the “real case scenario,” as opposed to the previous set, where a significant number of structures were model‐like only. In addition, the larger number of complexes (165 vs. 63 in the previous set) increases the statistical reliability of the benchmarking. We estimated the highest accuracy of the predicted complexes (according to CAPRI criteria), which can be attained using the benchmark structures. The set is available at http://dockground.bioinformatics.ku.edu . Proteins 2015; 83:891–897. © 2015 Wiley Periodicals, Inc.     

Scoring optimisation of unbound protein-protein docking including protein binding site predictions

The prediction of the structure of the protein-protein complex is of great importance to better understand molecular recognition processes. During systematic protein-protein docking, the surface of a protein molecule is scanned for putative binding sites of a partner protein. The possibility to include external data based on either experiments or bioinformatic predictions on putative binding sites during docking has been systematically explored. The external data were included during docking with a coarse-grained protein model and on the basis of force field weights to bias the docking search towards a predicted or known binding region. The approach was tested on a large set of protein partners in unbound conformations. The significant improvement of the docking performance was found if reliable data on the native binding sites were available. This was possible even if data for single key amino acids at a binding interface are included. In case of binding site predictions with limited accuracy, only modest improvement compared with unbiased docking was found. The optimisation of the protocol to bias the search towards predicted binding sites was found to further improve the docking performance resulting in approximately 40% acceptable solutions within the top 10 docking predictions compared with 22% in case of unbiased docking of unbound protein structures.     

GEMDOCK: a generic evolutionary method for molecular docking

We have developed an evolutionary approach for flexible ligand docking. This approval, GEMDOCK, uses a Generic Evolutionary Method for molecular DOCKing and an empirical scoring function. The former combines both discrete and continuous global search strategies with local search strategies to speed up convergence, whereas the latter results in rapid recognition of potential ligands. GEMDOCK was tested on a diverse data set of 100 protein-ligand complexes from the Protein Data Bank. In 79% of these complexes, the docked lowest energy ligand structures had root-mean-square derivations (RMSDs) below 2.0 A with respect to the corresponding crystal structures. The success rate increased to 85% if the structure water molecules were retained. We evaluated GEMDOCK on two cross-docking experiments in which each ligand of a protein ensemble was docked into each protein of the ensemble. Seventy-six percent of the docked structures had RMSDs below 2.0 A when the ligands were docked into foreign structures. We analyzed and validated GEMDOCK with respect to various search spaces and scoring functions, and found that if the scoring function was perfect, then the predicted accuracy was also essentially perfect. This study suggests that GEMDOCK is a useful tool for molecular recognition and may be used to systematically evaluate and thus improve scoring functions.     

Defining fitness in evolutionary models

The analysis of evolutionary models requires an appropriate definition for fitness. In this paper, I review such definitions in relation to the five major dimensions by which models may be described, namely (i) finite versus infinite (or very large) population size, (ii) type of environment (constant, fixed length, temporally stochastic, temporally predictable, spatially stochastic, spatially predictable and social environment), (iii) density-independent or density-dependent, (iv) inherent population dynamics (equilibrium, cyclical and chaotic), and (v) frequency-dependent or independent. In simple models, the Malthusian parameter ‘r’ or the net reproductive rate R ₀ may be satisfactory, but once density-dependence or complex population dynamics is introduced the invasion exponent should be used. Defining fitness in a social environment or when there is frequency-dependence requires special consideration.     

Optimizing long intrinsic disorder predictors with protein evolutionary information

Protein existing as an ensemble of structures, called intrinsically disordered, has been shown to be responsible for a wide variety of biological functions and to be common in nature. Here we focus on improving sequence-based predictions of long (>30 amino acid residues) regions lacking specific 3-D structure by means of four new neural-network-based Predictors Of Natural Disordered Regions (PONDRs): VL3, VL3H, VL3P, and VL3E. PONDR VL3 used several features from a previously introduced PONDR VL2, but benefitted from optimized predictor models and a slightly larger (152 vs. 145) set of disordered proteins that were cleaned of mislabeling errors found in the smaller set. PONDR VL3H utilized homologues of the disordered proteins in the training stage, while PONDR VL3P used attributes derived from sequence profiles obtained by PSI-BLAST searches. The measure of accuracy was the average between accuracies on disordered and ordered protein regions. By this measure, the 30-fold cross-validation accuracies of VL3, VL3H, and VL3P were, respectively, 83.6 +/- 1.4%, 85.3 +/- 1.4%, and 85.2 +/- 1.5%. By combining VL3H and VL3P, the resulting PONDR VL3E achieved an accuracy of 86.7 +/- 1.4%. This is a significant improvement over our previous PONDRs VLXT (71.6 +/- 1.3%) and VL2 (80.9 +/- 1.4%). The new disorder predictors with the corresponding datasets are freely accessible through the web server at http://www.ist.temple.edu/disprot.     

Observability in dynamic evolutionary models

In the paper observability problems are considered in basic dynamic evolutionary models for sexual and asexual populations. Observability means that from the (partial) knowledge of certain phenotypic characteristics the whole evolutionary process can be uniquely recovered. Sufficient conditions are given to guarantee observability for both sexual and asexual populations near an evolutionarily stable state.     

Best alpha-helical transmembrane protein topology predictions are achieved using hidden Markov models and evolutionary information

Methods that predict the topology of helical membrane proteins are standard tools when analyzing any proteome. Therefore, it is important to improve the performance of such methods. Here we introduce a novel method, PRODIV-TMHMM, which is a profile-based hidden Markov model (HMM) that also incorporates the best features of earlier HMM methods. In our tests, PRODIV-TMHMM outperforms earlier methods both when evaluated on "low-resolution" topology data and on high-resolution 3D structures. The results presented here indicate that the topology could be correctly predicted for approximately two-thirds of all membrane proteins using PRODIV-TMHMM. The importance of evolutionary information for topology prediction is emphasized by the fact that compared with using single sequences, the performance of PRODIV-TMHMM (as well as two other methods) is increased by approximately 10 percentage units by the use of homologous sequences. On a more general level, we also show that HMM-based (or similar) methods perform superiorly to methods that focus mainly on identification of the membrane regions.     

Advancing metabolic models with kinetic information

1. Download : Download high-res image (145KB)
2. Download : Download full-size image

     

EADock: docking of small molecules into protein active sites with a multiobjective evolutionary optimization

In recent years, protein-ligand docking has become a powerful tool for drug development. Although several approaches suitable for high throughput screening are available, there is a need for methods able to identify binding modes with high accuracy. This accuracy is essential to reliably compute the binding free energy of the ligand. Such methods are needed when the binding mode of lead compounds is not determined experimentally but is needed for structure-based lead optimization. We present here a new docking software, called EADock, that aims at this goal. It uses an hybrid evolutionary algorithm with two fitness functions, in combination with a sophisticated management of the diversity. EADock is interfaced with the CHARMM package for energy calculations and coordinate handling. A validation was carried out on 37 crystallized protein-ligand complexes featuring 11 different proteins. The search space was defined as a sphere of 15 A around the center of mass of the ligand position in the crystal structure, and on the contrary to other benchmarks, our algorithm was fed with optimized ligand positions up to 10 A root mean square deviation (RMSD) from the crystal structure, excluding the latter. This validation illustrates the efficiency of our sampling strategy, as correct binding modes, defined by a RMSD to the crystal structure lower than 2 A, were identified and ranked first for 68% of the complexes. The success rate increases to 78% when considering the five best ranked clusters, and 92% when all clusters present in the last generation are taken into account. Most failures could be explained by the presence of crystal contacts in the experimental structure. Finally, the ability of EADock to accurately predict binding modes on a real application was illustrated by the successful docking of the RGD cyclic pentapeptide on the alphaVbeta3 integrin, starting far away from the binding pocket.     

Yeasts as models in evolutionary biology

A report of the workshop 'Evolutionary and Environmental Genomics of Yeasts', Heidelberg, Germany, 1-5 October 2008.     

Prediction of distant residue contacts with the use of evolutionary information

In this work we present a novel correlated mutations analysis (CMA) method that is significantly more accurate than previously reported CMA methods. Calculation of correlation coefficients is based on physicochemical properties of residues (predictors) and not on substitution matrices. This results in reliable prediction of pairs of residues that are distant in protein sequence but proximal in its three dimensional tertiary structure. Multiple sequence alignments (MSA) containing a sequence of known structure for 127 families from PFAM database have been selected so that all major protein architectures described in CATH classification database are represented. Protein sequences in the selected families were filtered so that only those evolutionarily close to the target protein remain in the MSA. The average accuracy obtained for the alpha beta class of proteins was 26.8% of predicted proximal pairs with average improvement over random accuracy (IOR) of 6.41. Average accuracy is 20.6% for the mainly beta class and 14.4% for the mainly alpha class. The optimum correlation coefficient cutoff (cc cutoff) was found to be around 0.65. The first predictor, which correlates to hydrophobicity, provides the most reliable results. The other two predictors give good predictions which can be used in conjunction to those of the first one. When stricter cc cutoff is chosen, the average accuracy increases significantly (38.76% for alpha beta class), but the trade off is a smaller number of predictions. The use of solvent accessible area estimations for filtering false positives out of the predictions is promising.     

Protected polymorphism and evolutionary stability in pleiotropic models with trait-specific dominance

When alleles have pleiotropic effects on a number of quantitative traits, the degree of dominance between a pair of alleles can be different for each trait. Such trait-specific dominance has been studied previously in models for the maintenance of genetic variation by antagonistic effects of an allele on two fitness components. By generalizing these models to an arbitrary number of fitness components or other phenotypic traits with different degrees of dominance, I show that genetic polymorphism is generally impossible without antagonistic fitness effects of different traits and without trait-specific dominance. I also investigate dominance and pleiotropy from a more long-term evolutionary perspective, allowing for the study of general ecological scenarios, and I discuss the effects of trait-specific dominance on evolutionary stability criteria. When selection is mainly directional and only trait-specific dominance and antagonism cause the emergence of polymorphism, then these polymorphisms can be overtaken by single mutants again, such that they are probably short-lived on an evolutionary time scale. Near evolutionarily singular points where directional selection is absent, trait-specific dominance and overdominance facilitate the emergence of polymorphism and cause evolutionary divergence in some cases. An important outcome of these models is that trait-specific dominance allows for the emergence of genetic polymorphisms without a selective disadvantage for heterozygotes. This removes the scope for the evolution of assortative mate choice and affects dominance modification. Sympatric speciation by disruptive ecological selection requires this heterozygote disadvantage in order to evolve, and therefore it becomes less plausible if the emergence of genetic polymorphism usually occurs via trait-specific dominance and antagonistic effects.     

Recombination estimation under complex evolutionary models with the coalescent composite-likelihood method

The composite-likelihood estimator (CLE) of the population recombination rate considers only sites with exactly two alleles under a finite-sites mutation model (McVean, G. A. T., P. Awadalla, and P. Fearnhead. 2002. A coalescent-based method for detecting and estimating recombination from gene sequences. Genetics 160:1231-1241). While in such a model the identity of alleles is not considered, the CLE has been shown to be robust to minor misspecification of the underlying mutational model. However, there are many situations where the putative mutation and demographic history can be quite complex. One good example is rapidly evolving pathogens, like HIV-1. First we evaluated the performance of the CLE and the likelihood permutation test (LPT) under more complex, realistic models, including a general time reversible (GTR) substitution model, rate heterogeneity among sites (Gamma), positive selection, population growth, population structure, and noncontemporaneous sampling. Second, we relaxed some of the assumptions of the CLE allowing for a four-allele, GTR + Gamma model in an attempt to use the data more efficiently. Through simulations and the analysis of real data, we concluded that the CLE is robust to severe misspecifications of the substitution model, but underestimates the recombination rate in the presence of exponential growth, population mixture, selection, or noncontemporaneous sampling. In such cases, the use of more complex models slightly increases performance in some occasions, especially in the case of the LPT. Thus, our results provide for a more robust application of the estimation of recombination rates.     

Developmental quantitative genetic models of evolutionary change

Discussions about evolutionary change in developmental processes or morphological structures are predicated on specific quantitative genetic models whose parameters predict whether evolutionary change can occur, its relative rate and direction, and if correlated change will occur in other related and unrelated structures. The appropriate genetic model should reflect the relevant genetical and developmental biology of the organisms, yet be simple enough in its parameters so that deductions can be made and hypotheses tested. As a consequence, the choice of the most appropriate genetic model for polygenically controlled traits is a complex tissue and the eventual choice of model is often a compromise between completeness of the model and computational expediency. Herein, we discuss several developmental quantitative genetic models for the evolution of development and morphology. The models range from the classical direct effects model to complex epigenetic models. Further, we demonstrate the algebraic equivalency of the Cowley and Atchley epigenetic model and Wagner's developmental mapping model. Finally, we propose a new multivariate model for continuous growth trajectories. The relative efficacy of these various models for understanding evolutionary change in developmental and morphological traits is discussed. © 1994 Wiley-Liss, Inc.     

Applying population-genetic models in theoretical evolutionary epidemiology

Much of the existing theory for the evolutionary biology of infectious diseases uses an invasion analysis approach. In this Ideas and Perspectives article, we suggest that techniques from theoretical population genetics can also be profitably used to study the evolutionary epidemiology of infectious diseases. We highlight four ways in which population-genetic models provide benefits beyond those provided by most invasion analyses: (i) they can make predictions about the rate of pathogen evolution; (ii) they explicitly draw out the mechanistic way in which the epidemiological dynamics feed into evolutionary change, and thereby provide new insights into pathogen evolution; (iii) they can make predictions about the evolutionary consequences of non-equilibrium epidemiological dynamics; (iv) they can readily incorporate the effects of multiple host dynamics, and thereby account for phenomena such as immunological history and/or host co-evolution.