Ab initio and template-based prediction of multi-class distance maps by two-dimensional recursive neural networks

Background

Prediction of protein structures from their sequences is still one of the open grand challenges of computational biology. Some approaches to protein structure prediction, especially ab initio ones, rely to some extent on the prediction of residue contact maps. Residue contact map predictions have been assessed at the CASP competition for several years now. Although it has been shown that exact contact maps generally yield correct three-dimensional structures, this is true only at a relatively low resolution (3–4 Å from the native structure). Another known weakness of contact maps is that they are generally predicted ab initio, that is not exploiting information about potential homologues of known structure.

Results

We introduce a new class of distance restraints for protein structures: multi-class distance maps. We show that C _α trace reconstructions based on 4-class native maps are significantly better than those from residue contact maps. We then build two predictors of 4-class maps based on recursive neural networks: one ab initio, or relying on the sequence and on evolutionary information; one template-based, or in which homology information to known structures is provided as a further input. We show that virtually any level of sequence similarity to structural templates (down to less than 10%) yields more accurate 4-class maps than the ab initio predictor. We show that template-based predictions by recursive neural networks are consistently better than the best template and than a number of combinations of the best available templates. We also extract binary residue contact maps at an 8 Å threshold (as per CASP assessment) from the 4-class predictors and show that the template-based version is also more accurate than the best template and consistently better than the ab initio one, down to very low levels of sequence identity to structural templates. Furthermore, we test both ab-initio and template-based 8 Å predictions on the CASP7 targets using a pre-CASP7 PDB, and find that both predictors are state-of-the-art, with the template-based one far outperforming the best CASP7 systems if templates with sequence identity to the query of 10% or better are available. Although this is not the main focus of this paper we also report on reconstructions of C _α traces based on both ab initio and template-based 4-class map predictions, showing that the latter are generally more accurate even when homology is dubious.

Conclusion

Accurate predictions of multi-class maps may provide valuable constraints for improved ab initio and template-based prediction of protein structures, naturally incorporate multiple templates, and yield state-of-the-art binary maps. Predictions of protein structures and 8 Å contact maps based on the multi-class distance map predictors described in this paper are freely available to academic users at the url http://distill.ucd.ie/.     

Riverscape genetics in brook lamprey: genetic diversity is less influenced by river fragmentation than by gene flow with the anadromous ecotype

Understanding the effect of human-induced landscape fragmentation on gene flow and evolutionary potential of wild populations has become a major concern. Here, we investigated the effect of riverscape fragmentation on patterns of genetic diversity in the freshwater resident European brook lamprey (Lampetra planeri) that has a low ability to pass obstacles to migration. We tested the hypotheses of (i) asymmetric gene flow following water current and (ii) an effect of gene flow with the closely related anadromous river lamprey (L. fluviatilis) ecotype on L. planeri genetic diversity. We genotyped 2472 individuals, including 225 L. fluviatilis, sampled from 81 sites upstream and downstream barriers to migration, in 29 western European rivers. Linear modelling revealed a strong positive relationship between genetic diversity and the distance from the river source, consistent with expected patterns of decreased gene flow into upstream populations. However, the presence of anthropogenic barriers had a moderate effect on spatial genetic structure. Accordingly, we found evidence for downstream-directed gene flow, supporting the hypothesis that barriers do not limit dispersal mediated by water flow. Downstream L. planeri populations in sympatry with L. fluviatilis displayed consistently higher genetic diversity. We conclude that genetic drift and slight downstream gene flow drive the genetic make-up of upstream L. planeri populations whereas gene flow between ecotypes maintains higher levels of genetic diversity in L. planeri populations sympatric with L. fluviatilis. We discuss the implications of these results for the design of conservation strategies of lamprey, and other freshwater organisms with several ecotypes, in fragmented dendritic river networks.Subject terms: Conservation biology, Ecological genetics, Evolutionary genetics, Genetic variation