Prediction of protein secondary structure content for the twilight zone sequences

Secondary protein structure carries information about local structural arrangements, which include three major conformations: alpha-helices, beta-strands, and coils. Significant majority of successful methods for prediction of the secondary structure is based on multiple sequence alignment. However, multiple alignment fails to provide accurate results when a sequence comes from the twilight zone, that is, it is characterized by low (<30%) homology. To this end, we propose a novel method for prediction of secondary structure content through comprehensive sequence representation, called PSSC-core. The method uses a multiple linear regression model and introduces a comprehensive feature-based sequence representation to predict amount of helices and strands for sequences from the twilight zone. The PSSC-core method was tested and compared with two other state-of-the-art prediction methods on a set of 2187 twilight zone sequences. The results indicate that our method provides better predictions for both helix and strand content. The PSSC-core is shown to provide statistically significantly better results when compared with the competing methods, reducing the prediction error by 5-7% for helix and 7-9% for strand content predictions. The proposed feature-based sequence representation uses a comprehensive set of physicochemical properties that are custom-designed for each of the helix and strand content predictions. It includes composition and composition moment vectors, frequency of tetra-peptides associated with helical and strand conformations, various property-based groups like exchange groups, chemical groups of the side chains and hydrophobic group, auto-correlations based on hydrophobicity, side-chain masses, hydropathy, and conformational patterns for beta-sheets. The PSSC-core method provides an alternative for predicting the secondary structure content that can be used to validate and constrain results of other structure prediction methods. At the same time, it also provides useful insight into design of successful protein sequence representations that can be used in developing new methods related to prediction of different aspects of the secondary protein structure.     

Characteristic amino acid combinations in olfactory G protein-coupled receptors

The human olfactory subgenome has recently been fully characterized with over 1000 genes. Although as many as two thirds of them are expected to be pseudogenes, it still leaves us with about half of all human G protein-coupled receptors being olfactory. It is therefore of great interest to characterize olfactory receptors with high precision. Usually it is done through sequence motifs that are not fully conserved, making an exact characterization difficult. In this paper, we propose a rule-based characterization of olfactory receptors derived from a multiple sequence alignment of human GPCRs. We show that just seven alignment sites are sufficient to characterize 99% of human olfactory GPCRs with one feature, a tyrosine at site 7.41, being of particular importance. We also show dependencies between sites near the extracellular and intracellular region of a membrane-embedded receptor, indicating that olfactory receptors are characterized by a combination of important residues in these two areas, whereas nonolfactory receptors tend to have residues of lower importance at the same sites.     

Linking stream habitats and spider distribution: spatial variations in trophic transfer across a forest–stream boundary

In headwater streams, many aquatic insects rely on terrestrial detritus, while their emergence from streams often subsidizes riparian generalist predators. However, spatial variations in such reciprocal trophic linkages remain poorly understood. The present study, conducted in a northern Japanese stream and the surrounding forest, showed that pool–riffle structure brought about heterogeneous distributions of detritus deposits and benthic aquatic insects. The resulting variations in aquatic insect emergence influenced the distributions of riparian web-building spiders. Pools with slow current stored greater amounts of detritus than riffles, allowing more benthic aquatic insects to develop in pools. The greater larval biomass in pools and greater tendency for riffle insects to drift into pools at metamorphosis resulted in an emergence rate of aquatic insects from pools that was some four to five times greater than from riffles. In the riparian forest, web-building spiders (Tetragnathidae and Linyphiidae) were distributed in accordance with the emergence rates of aquatic insects, upon which both spider groups heavily depended. Consequently, the riparian strips bordering pools had a density of tetragnathid spiders that was twice as high as that of the riparian strips adjacent to riffles. Moreover, although limitations of vegetation structure prevented the aggregation of linyphiid spiders around pools, linyphiid density normalized by shrub density was higher in habitats adjacent to pools than those adjacent to riffles. The results indicated that stream geomorphology, which affects the storage of terrestrial organic material and the export of such material to riparian forests via aquatic insect emergence, plays a role in determining the strength of terrestrial–aquatic linkages in headwater ecosystems.     

山东地区犬冠状病毒的分子流行病学调查及其基因型多样性分析

【目的】了解分析山东地区健康犬及腹泻犬中犬冠状病毒(CCoV)的分子流行病学及其基因型。【方法】采集2017年以来山东省宠物市场、流浪犬中心、犬舍、各地宠物医院等的健康和腹泻犬只的肛拭子样品,采用PCR方法检测CCoV及其CCoV-Ⅰ和CCoV-Ⅱ(CCoV-Ⅱa和CCoV-Ⅱb)基因分型情况,并对扩增出的M、S全基因进行测序分析。【结果】199份样品中,共检出CCoV阳性样品79份,健康犬中的检出率为40.2%(33/82),腹泻犬中的检出率为39.3%(46/117)。健康犬中CCoV-Ⅰ型与CCoV-Ⅱ型阳性检出率分别为78.8%(26/33)和51.5%(17/33);腹泻犬中CCoV-Ⅰ型与CCoV-Ⅱ型阳性检出率分别为39.13%(18/46)和80.4%(37/46)。基因型阳性比率分析表明,健康犬中单独感染CCoV-I型的比率最高,为48.5%(16/33);腹泻犬中单独感染CCoV-Ⅱa比率最高,为47.8%(22/46)。测序分析获得48株M基因序列,遗传进化分析表明,12株CCoV-Ⅰ型毒株中除1株外,其余均从健康犬中分离。36株CCoV-Ⅱ型毒株中除7株来自健康犬外,其余均为腹泻犬中获得。获得的3株S全基因均来自腹泻犬且均属于CCoV-Ⅱa亚型。【结论】山东地区犬群中CCoV检出率较高,说明CCoV广泛流行,健康犬和腹泻犬中均存在多重感染情况,健康犬中主要流行CCoV毒株为CCoV-Ⅰ型,腹泻犬中主要流行CCoV-Ⅱa型。     

临床分离志贺菌中CRISPR/Cas系统的分布及其与毒力基因的关系

【目的】了解临床分离志贺菌中CRISPR/Cas系统的分布特征并分析其与毒力基因的关系。【方法】以聚合酶链式反应(PCR)方法,采用10对引物分别对57株临床分离志贺菌中CRISPR1、cas2-cas1、cas6e-cas5、cas7、cse2、cse1-cas3基因和毒力基因ipaH、ial、ipaBCD、virA进行检测。对CRISPR1的PCR结果进行测序,并用CRISPR finder在线软件对CRISPR1基因座进行分析。通过卡方检验初步分析CRISPR/Cas系统与毒力基因的关系。【结果】测序结果显示,CRISPR1基因座中间隔序列数目较少且在不同菌株间一致性较高;57株志贺菌中,84.2% (48/57)的志贺菌中可检测到CRISPR/Cas系统,其中68.8% (33/48)的志贺菌中cas6e-cas5基因或(和) cse2基因中发现插入序列;毒力基因ipaH、ial、virA、ipaBCD的检出率依次为100%、100%、98.2%和87.7%;毒力基因ipaBCD的阳性率与活性CRISPR/Cas系统的分布无关(P>0.05)。【结论】CRISPR/Cas系统广泛存在于临床分离志贺菌中;部分cas基因中有插入序列;并未发现志贺菌中活性CRISPR/Cas系统与毒力基因的分布有关。     

Computational modeling of membrane proteins

The determination of membrane protein (MP) structures has always trailed that of soluble proteins due to difficulties in their overexpression, reconstitution into membrane mimetics, and subsequent structure determination. The percentage of MP structures in the protein databank (PDB) has been at a constant 1–2% for the last decade. In contrast, over half of all drugs target MPs, only highlighting how little we understand about drug‐specific effects in the human body. To reduce this gap, researchers have attempted to predict structural features of MPs even before the first structure was experimentally elucidated. In this review, we present current computational methods to predict MP structure, starting with secondary structure prediction, prediction of trans‐membrane spans, and topology. Even though these methods generate reliable predictions, challenges such as predicting kinks or precise beginnings and ends of secondary structure elements are still waiting to be addressed. We describe recent developments in the prediction of 3D structures of both α‐helical MPs as well as β‐barrels using comparative modeling techniques, de novo methods, and molecular dynamics (MD) simulations. The increase of MP structures has (1) facilitated comparative modeling due to availability of more and better templates, and (2) improved the statistics for knowledge‐based scoring functions. Moreover, de novo methods have benefited from the use of correlated mutations as restraints. Finally, we outline current advances that will likely shape the field in the forthcoming decade. Proteins 2015; 83:1–24. © 2014 Wiley Periodicals, Inc.     

An analysis approach to identify specific functional sites in orthologous proteins using sequence and structural information: Application to neuroserpin reveals regions that differentially regulate inhibitory activity

The analysis of sequence conservation is commonly used to predict functionally important sites in proteins. We have developed an approach that first identifies highly conserved sites in a set of orthologous sequences using a weighted substitution‐matrix‐based conservation score and then filters these conserved sites based on the pattern of conservation present in a wider alignment of sequences from the same family and structural information to identify surface‐exposed sites. This allows us to detect specific functional sites in the target protein and exclude regions that are likely to be generally important for the structure or function of the wider protein family. We applied our method to two members of the serpin family of serine protease inhibitors. We first confirmed that our method successfully detected the known heparin binding site in antithrombin while excluding residues known to be generally important in the serpin family. We next applied our sequence analysis approach to neuroserpin and used our results to guide site‐directed polyalanine mutagenesis experiments. The majority of the mutant neuroserpin proteins were found to fold correctly and could still form inhibitory complexes with tissue plasminogen activator (tPA). Kinetic analysis of tPA inhibition, however, revealed altered inhibitory kinetics in several of the mutant proteins, with some mutants showing decreased association with tPA and others showing more rapid dissociation of the covalent complex. Altogether, these results confirm that our sequence analysis approach is a useful tool that can be used to guide mutagenesis experiments for the detection of specific functional sites in proteins. Proteins 2015; 83:135–152. © 2014 Wiley Periodicals, Inc.     

Coiled‐coil length: Size does matter

Protein evolution is governed by processes that alter primary sequence but also the length of proteins. Protein length may change in different ways, but insertions, deletions and duplications are the most common. An optimal protein size is a trade‐off between sequence extension, which may change protein stability or lead to acquisition of a new function, and shrinkage that decreases metabolic cost of protein synthesis. Despite the general tendency for length conservation across orthologous proteins, the propensity to accept insertions and deletions is heterogeneous along the sequence. For example, protein regions rich in repetitive peptide motifs are well known to extensively vary their length across species. Here, we analyze length conservation of coiled‐coils, domains formed by an ubiquitous, repetitive peptide motif present in all domains of life, that frequently plays a structural role in the cell. We observed that, despite the repetitive nature, the length of coiled‐coil domains is generally highly conserved throughout the tree of life, even when the remaining parts of the protein change, including globular domains. Length conservation is independent of primary amino acid sequence variation, and represents a conservation of domain physical size. This suggests that the conservation of domain size is due to functional constraints. Proteins 2015; 83:2162–2169. © 2015 Wiley Periodicals, Inc.