Camps 2.0: exploring the sequence and structure space of prokaryotic, eukaryotic, and viral membrane proteins

Structural bioinformatics of membrane proteins is still in its infancy, and the picture of their fold space is only beginning to emerge. Because only a handful of three-dimensional structures are available, sequence comparison and structure prediction remain the main tools for investigating sequence-structure relationships in membrane protein families. Here we present a comprehensive analysis of the structural families corresponding to α-helical membrane proteins with at least three transmembrane helices. The new version of our CAMPS database (CAMPS 2.0) covers nearly 1300 eukaryotic, prokaryotic, and viral genomes. Using an advanced classification procedure, which is based on high-order hidden Markov models and considers both sequence similarity as well as the number of transmembrane helices and loop lengths, we identified 1353 structurally homogeneous clusters roughly corresponding to membrane protein folds. Only 53 clusters are associated with experimentally determined three-dimensional structures, and for these clusters CAMPS is in reasonable agreement with structure-based classification approaches such as SCOP and CATH. We therefore estimate that ～1300 structures would need to be determined to provide a sufficient structural coverage of polytopic membrane proteins. CAMPS 2.0 is available at http://webclu.bio.wzw.tum.de/CAMPS2.0/.     

PROSO II--a new method for protein solubility prediction

Many fields of science and industry depend on efficient production of active protein using heterologous expression in Escherichia coli. The solubility of proteins upon expression is dependent on their amino acid sequence. Prediction of solubility from sequence is therefore highly valuable. We present a novel machine-learning-based model called PROSO II which makes use of new classification methods and growth in experimental data to improve coverage and accuracy of solubility predictions. The classification algorithm is organized as a two-layered structure in which the output of a primary Parzen window model for sequence similarity and a logistic regression classifier of amino acid k-mer composition serve as input for a second-level logistic regression classifier. Compared with previously published research our model is trained on five times more data than used by any other method before (82 000 proteins). When tested on a separate holdout set not used at any point of method development our server attained the best results in comparison with other currently available methods: accuracy 75.4%, Matthew's correlation coefficient 0.39, sensitivity 0.731, specificity 0.759, gain (soluble) 2.263. In summary, due to utilization of cutting edge machine learning technologies combined with the largest currently available experimental data set the PROSO II server constitutes a substantial improvement in protein solubility predictions. PROSO II is available at http://mips.helmholtz-muenchen.de/prosoII.     

Phylogenetic web profiler

SUMMARY: Phylogenetic Web Profiler (PWP) is a web-based service designed to perform phylogenetic profiling of proteins against genomes. The current version offers a selection of 63 completed genomes and available plasmids as annotated in the PEDANT genome database. Unlike currently available applications, this tool offers several choices of ortholog prediction parameters including E-value cutoff, percent length difference tolerance, and annotation similarity. Additional features include tight integration with the PEDANT database and tools to analyze properties of predicted proteins. PWP should prove very useful for the analysis of functional-linkage between proteins.     

Replacement or exclusion of the B-branch bacteriopheophytin in the purple bacterial reaction centre: The HB cofactor is not required for assembly or core function of the Rhodobacter sphaeroides complex

All of the membrane-embedded cofactors of the purple bacterial reaction centre have well-defined functional or structural roles, with the exception of the bacteriopheophytin (H_B) located approximately half-way across the membrane on the so-called inactive- or B-branch of cofactors. Sequence alignments indicate that this bacteriochlorin cofactor is a conserved feature of purple bacterial reaction centres, and a pheophytin is also found at this position in the Photosystem-II reaction centre. Possible structural or functional consequences of replacing the H_B bacteriopheophytin by bacteriochlorophyll were investigated in the Rhodobacter sphaeroides reaction centre through mutagenesis of residue Leu L185 to His (LL185H). Results from absorbance spectroscopy indicated that the LL185H mutant assembled with a bacteriochlorophyll at the H_B position, but this did not affect the capacity of the reaction centre to support photosynthetic growth, or change the kinetics of charge separation along the A-branch of cofactors. It was also found that mutation of residue Ala M149 to Trp (AM149W) caused the reaction centre to assemble without an H_B bacteriochlorin, demonstrating that this cofactor is not required for correct assembly of the reaction centre. The absence of a cofactor at this position did not affect the capacity of the reaction centre to support photosynthetic growth, or the kinetics of A-branch electron transfer. A combination of X-ray crystallography and FTIR difference spectroscopy confirmed that the H_B cofactor was absent in the AM149W mutant, and that this had not produced any significant disturbance of the adjacent ubiquinol reductase (Q_B) site. The data are discussed with respect to possible functional roles of the H_B bacteriopheophytin, and we conclude that the reason(s) for conservation of a bacteriopheophytin cofactor at this position in purple bacterial reaction centres are likely to be different from those underlying conservation of a pheophytin at the analogous position in Photosystem-II.     

Mechanisms of the changes in the rate and nature of reparative regeneration of the aortic endothelium after repeated injuries

An experimental model of repeated cryogenic lesions in the rat abdominal aorta endothelium with a concentrically decreasing area of the defect has been worked out. In reendothelization of every successive defect participate the cells of the endothelial layer that is formed anew after the preceding lesions. As a result of repeated lesions the regeneration rate increases by 1.2 times. By means of scanning radioautography after 3H-thymidine++ administration the index of the labelled nuclei of the endothelial cells (EC) has been demonstrated to increase by 2.2 times. In the experiments with gamma-irradiation of the vessel before the last cryodestruction (this results in blockade of proliferation) an increased rate of EC migration has been revealed. Repeated lesions also produce an increase in the heteromorphism degree of the endothelial layer. This is mainly manifested as appearance of multinuclear EC clusters. These changes in rate and in character of the endothelium regeneration can be determined by the mechanisms similar to clonal proliferative senescence of EC in vitro.     

Changes in three-dimensional structure of the internal elastic membrane of rat aorta after its mechanical injury and regeneration

By means of scanning electron microscopy of chemically extracted preparations the dynamics of changes in the three-dimensional structure of the internal elastic membrane (IEM) of the rat aorta has been investigated after its lesion by a vascular clip. The mechanical lesion results in rupture of the IEM along external borders of the instrument lips and in crushing of its central part. A niche is formed, along its periphery it is surrounded with a practically intact IEM. The aorta regeneration is accompanied with neoelastogenesis. In the center of the niche newly formed elastic structures appear later, but the IEM reparation develops more actively. During the neoelastogenesis some stages are distinguished: at first separate elastic fibers appear, which then anastomosed, uniting into fasciculi and laminae. It is supposed that the IEM-restoration at regeneration depends on synthetic activity of smooth myocytes and on the other hand, the changes in the IEM structure can regulate their migration and metabolism.     

Prediction of β-turns and β-turn types by a novel bidirectional Elman-type recurrent neural network with multiple output layers (MOLEBRNN)

Prediction of β-turns from amino acid sequences has long been recognized as an important problem in structural bioinformatics due to their frequent occurrence as well as their structural and functional significance. Because various structural features of proteins are intercorrelated, secondary structure information has been often employed as an additional input for machine learning algorithms while predicting β-turns. Here we present a novel bidirectional Elman-type recurrent neural network with multiple output layers (MOLEBRNN) capable of predicting multiple mutually dependent structural motifs and demonstrate its efficiency in recognizing three aspects of protein structure: β-turns, β-turn types, and secondary structure. The advantage of our method compared to other predictors is that it does not require any external input except for sequence profiles because interdependencies between different structural features are taken into account implicitly during the learning process. In a sevenfold cross-validation experiment on a standard test dataset our method exhibits the total prediction accuracy of 77.9% and the Mathew's Correlation Coefficient of 0.45, the highest performance reported so far. It also outperforms other known methods in delineating individual turn types. We demonstrate how simultaneous prediction of multiple targets influences prediction performance on single targets. The MOLEBRNN presented here is a generic method applicable in a variety of research fields where multiple mutually depending target classes need to be predicted. Availability: http://webclu.bio.wzw.tum.de/predator-web/.     

Characterization of a thaumarchaeal symbiont that drives incomplete nitrification in the tropical sponge Ianthella basta

Marine sponges represent one of the few eukaryotic groups that frequently harbour symbiotic members of the Thaumarchaeota, which are important chemoautotrophic ammonia-oxidizers in many environments. However, in most studies, direct demonstration of ammonia-oxidation by these archaea within sponges is lacking, and little is known about sponge-specific adaptations of ammonia-oxidizing archaea (AOA). Here, we characterized the thaumarchaeal symbiont of the marine sponge Ianthella basta using metaproteogenomics, fluorescence in situ hybridization, qPCR and isotope-based functional assays. ‘Candidatus Nitrosospongia ianthellae’ is only distantly related to cultured AOA. It is an abundant symbiont that is solely responsible for nitrite formation from ammonia in I. basta that surprisingly does not harbour nitrite-oxidizing microbes. Furthermore, this AOA is equipped with an expanded set of extracellular subtilisin-like proteases, a metalloprotease unique among archaea, as well as a putative branched-chain amino acid ABC transporter. This repertoire is strongly indicative of a mixotrophic lifestyle and is (with slight variations) also found in other sponge-associated, but not in free-living AOA. We predict that this feature as well as an expanded and unique set of secreted serpins (protease inhibitors), a unique array of eukaryotic-like proteins, and a DNA-phosporothioation system, represent important adaptations of AOA to life within these ancient filter-feeding animals.     

Co-evolving residues in membrane proteins

MOTIVATION: The analysis of co-evolving residues has been exhaustively evaluated for the prediction of intramolecular amino acid contacts in soluble proteins. Although a variety of different methods for the detection of these co-evolving residues have been developed, the fraction of correctly predicted contacts remained insufficient for their reliable application in the construction of structural models. Membrane proteins, which constitute between one-fourth and one-third of all proteins in an organism, were only considered in few individual case studies. RESULTS: We present the first general study of correlated mutations in alpha-helical membrane proteins. Using seven different prediction algorithms, we extracted co-evolving residues for 14 membrane proteins having a solved 3D structure. On average, distances between correlated pairs of residues lying on different transmembrane segments were found to be significantly smaller compared to a random prediction. Covariation of residues was frequently found in direct sequence neighborhood to helix-helix contacts. Based on the results obtained from individual prediction methods, we constructed a consensus prediction for every protein in the dataset that combines obtained correlations from different prediction algorithms and simultaneously removes likely false positives. Using this consensus prediction, 53% of all predicted residue pairs were found within one helix turn of an observed helix-helix contact. Based on the combination of co-evolving residues detected with the four best prediction algorithms, interacting helices could be predicted with a specificity of 83% and sensitivity of 42%. AVAILABILITY: http://webclu.bio.wzw.tum.de/helixcorr/