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/.     

Five colour morphs and three new species of Gyptis (Hesionidae, Annelida) under a jetty in Edithburgh, South Australia

We report five different colour morphs of the hesionid polychaete genus Gyptis co-occurring in a small area in shallow water under Edithburgh jetty, South Australia. The five morphs cannot be separated using standard morphological features, but phylogenetic analyses of sequence data (COI and ITS1) unequivocally show that three species are present, introduced as Gyptis simpsonorum , new species, G. paucilineata , new species and G. polymorpha , new species. Gyptis simpsonorum has a speckled pigmentation pattern and G. paucilineata a few transverse lines on specific segments. Gyptis polymorpha is polymorphic with three different, distinct pigmentation patterns, either as dense transverse lines, as a thin longitudinal, mid-dorsal line, or as an uniformly dark dorsum. The speckled pigmentation likely represents the plesiomorphic condition, and G. propinqua is the closest known relative to these new species.     

Entomotoxic action of Sambucus nigra agglutinin i in Acyrthosiphon pisum aphids and Spodoptera exigua caterpillars through caspase‐3‐like‐dependent apoptosis

In this project, the toxicity and mechanism of action of the ricin‐B‐related lectin SNA‐I from elderberry (Sambucus nigra) in the pea aphid (Acyrthosiphon pisum) and the beet armyworm (Spodoptera exigua), two important pest insects in agriculture, were studied. SNA‐I is a chimeric lectin belonging to the class of ribosome‐inactivating proteins and consists of an A‐chain with N‐glycosidase activity and a carbohydrate‐binding B‐chain. Incorporation of 2 mg/ml of SNA‐I in the diet of neonates and adults of A. pisum caused 40–46% mortality within 2 days, while in third instars of S. exigua, the larval biomass was significantly reduced by 12% after feeding for 3 days on a diet containing 5 mg/g of SNA‐I. Interestingly, extracts of the (mid)gut of treated A. pisum and S. exigua demonstrated DNA fragmentation and this was accompanied with an increase in caspase‐3‐like activity. The involvement of cell death or apoptosis in the entomotoxicity of SNA‐I through induction of caspase‐3‐like activity was also confirmed by addition of the permeable caspase‐3 inhibitor III in the diet, leading to a rescue of the treated aphid neonates. Finally, similar to the chimeric lectin SNA‐I, the hololectin SNA‐II, consisting of two carbohydrate‐binding B‐chains caused high mortality to neonate A. pisum aphids with an LC₅₀ of 1.59 mg/ml, suggesting that the entomotoxic action of the lectins under study mainly relies on their carbohydrate‐binding activity. © 2010 Wiley Periodicals, Inc.     

Recombinant human factor VIII-specific affinity ligands selected from phage-displayed combinatorial libraries of protein A.

Factor VIII-specific affibodies were selected from phage displayed libraries constructed by combinatorial mutagenesis of an alpha helical bacterial receptor domain derived from staphylococcal protein A. Bead-immobilized recombinant human factor VIII (rVIII) (80 and 90 kDa chains) protein was used during competitive biopannings in the presence of free 80-kDa chain protein, resulting in the selection of several binders that showed dissociation constants (Kd) in the range 100-200 nM as determined by biosensor analyses. One variant (Z[rVIII:3], 90-kDa chain specific) was further characterized in small-scale affinity chromatography experiments, and showed efficient and selective recovery of biologically active rVIII from Chinese hamster ovary cell supernatant-derived feed stocks. The purity of the enriched rVIII was comparable with rVIII material purified by immunoaffinity chromatography using a 90-kDa chain-specific monoclonal antibody. Interestingly, epitope mapping showed that the monoclonal antibody and the affibody ligand competed for the same or at least overlapping epitopes on rVIII. In addition, the Z[rVIII:3] variant was produced by peptide synthesis with a C-terminal cysteine to enable directed coupling to solid supports. This 59-residue protein was analyzed by circular dichroism and showed a secondary structure content similar to that of the parental Z domain used as scaffold. In biosensor studies, the synthetic affibody was immobilized recruiting the C-terminal cysteine residue, and demonstrated to bind both recombinantly produced and plasma-derived factor VIII. From a secondary library, constructed by re-randomization of relevant positions identified after alignment of the first-generation variants, a panel of affinity-improved second-generation affibodies were selected of which one clone showed a dissociation constant (Kd) for rVIII of 5 nM. Several of these variants also showed higher apparent binding efficiencies towards rVIII when analyzed as immobilized ligands in biosensor experiments. Taken together, the results suggest that affibody ligands produced by bacterial or synthetic routes could be of interest as an alternative to monoclonal antibodies in purification processes or as diagnostic or monitoring tools.     

An Improved Negative Staining Technique Using A Thin Quartz Membrane as Sample Support

A negative staining technique is presented baaed on the use of 40-60 nm quartz membrane supported by a silicon grid. The quartz membrane is fabricated by thermal growth of silicon dioxide on a silicon substrate followed by an anisotropic silicon etching step giving rectangular holes in the silicon substrate. The hydrophilic membrane is shown to be ideally suited for negative staining due to its spreading characteristics, homogeneity, heat resistance and mechanical stability. Micrographs of phage λ are presented showing the detailed structure of the tail. A simple method of calculating the number of adsorbed particles based on diffusion limited association is also presented.     

Molecular coevolution among cryptically simple expansion segments of eukaryotic 26S/28S rRNAs

The set of "expansion segments" of any eukaryotic 26S/28S ribosomal RNA (rRNA) gene is responsible for the bulk of the difference in length between the prokaryotic 23S rRNA gene and the eukaryotic 26S/28S rRNA gene. The expansion segments are also responsible for interspecific fluctuations in length during eukaryotic evolution. They show a consistent bias in base composition in any species; for example, they are AT rich in Drosophila melanogaster and GC rich in vertebrate species. Dot-matrix comparisons of sets of expansion segments reveal high similarities between members of a set within any 28S rRNA gene of a species, in contrast to the little or spurious similarity that exists between sets of expansion segments from distantly related species. Similarities among members of a set of expansion segments within any 28S rRNA gene cannot be accounted for by their base-compositional bias alone. In contrast, no significant similarity exists within a set of "core" segments (regions between expansion segments) of any 28S rRNA gene, although core segments are conserved between species. The set of expansion segments of a 26S/28S gene is coevolving as a unit in each species, at the same time as the family of 28S rRNA genes, as a whole, is undergoing continual homogenization, making all sets of expansion segments from all ribosomal DNA (rDNA) arrays in a species similar in sequence. Analysis of DNA simplicity of 26S/28S rRNA genes shows a direct correlation between significantly high relative simplicity factors (RSFs) and sequence similarity among a set of expansion segments. A similar correlation exists between RSF values, overall rDNA lengths, and the lengths of individual expansion segments. Such correlations suggest that most length fluctuations reflect the gain and loss of simple sequence motifs by slippage-like mechanisms. We discuss the molecular coevolution of expansion segments, which takes place against a background of slippage-like and unequal crossing-over mechanisms of turnover that are responsible for the accumulation of interspecific differences in rDNA sequences.      

Functional dedifferentiation of parathyroid cells results both from a defective regulation and action of cytoplasmic Ca2+

Monolayer culture of bovine parathyroid cells for 24 hours resulted in a right-shift of the dose-effect relationships for Ca²⁺-inhibition of parathyroid hormone (PTH) release and the dependence of the cytoplasmic Ca²⁺ concentration (Ca²⁺) on extracellular Ca²⁺ as well as in a less suppressible hormone release. After 4 days of culture, hormone secretion was almost non-suppressible and Ca _i ²⁺ increased poorly in response to a rise in extracelluiar Ca²⁺. Ionomycin, a Ca²⁺ ionophore, raised Ca _i ²⁺ , but there was only a small inhibition of PTH release and the correlation between Ca _i ²⁺ and secretion was weak. A deteriorated Ca _i ²⁺ regulation and a decreased inhibitory action of cytoplasmic Ca²⁺ on PTH release were also found in ceils from human parathyroid adenomas. Functional dedifferentiation of the parathyroid cell thus results from both defective regulation and action of cytoplasmic Ca²⁺.