Protein‐fold recognition using an improved single‐source K diverse shortest paths algorithm

Protein structure prediction, when construed as a fold recognition problem, is one of the most important applications of similarity search in bioinformatics. A new protein‐fold recognition method is reported which combines a single‐source K diverse shortest path (SSKDSP) algorithm with Enrichment of Network Topological Similarity (ENTS) algorithm to search a graphic feature space generated using sequence similarity and structural similarity metrics. A modified, more efficient SSKDSP algorithm is developed to improve the performance of graph searching. The new implementation of the SSKDSP algorithm empirically requires 82% less memory and 61% less time than the current implementation, allowing for the analysis of larger, denser graphs. Furthermore, the statistical significance of fold ranking generated from SSKDSP is assessed using ENTS. The reported ENTS‐SSKDSP algorithm outperforms original ENTS that uses random walk with restart for the graph search as well as other state‐of‐the‐art protein structure prediction algorithms HHSearch and Sparks‐X, as evaluated by a benchmark of 600 query proteins. The reported methods may easily be extended to other similarity search problems in bioinformatics and chemoinformatics. The SSKDSP software is available at http://compsci.hunter.cuny.edu/~leixie/sskdsp.html . Proteins 2016; 84:467–472. © 2016 Wiley Periodicals, Inc.     

基于蛋白互作网络识别非小细胞肺癌相关基因功能模块

本研究对非小细胞肺癌(non-small cell lung carcinoma,NSCLC)基因表达数据进行差异表达分析,并与蛋白质相互作用网络(PPIN)数据进行整合,进一步利用Heinz搜索算法识别NSCLC相关的基因功能模块,并对模块中的基因进行功能(GO term)和通路(KEGG)富集分析,旨在探究肺癌发病分子机制。蛋白互作网络分析得到一个包含96个基因和117个相互作用的功能模块,以及8个对NSCLC的发生和发展起到关键作用候选基因标志物。富集分析结果表明,这些基因主要富集于基因转录催化及染色质调控等生物学过程,并在基础转录因子、黏着连接、细胞周期、Wnt信号通路及HTLV-Ⅰ感染等生物学通路中发挥重要作用。本研究对非小细胞肺癌相关的基因和生物学通路进行预测,可用于肺癌的早期诊断和早期治疗,以降低肺癌死亡率。     

Dosimetric impact of statistical uncertainty on Monte Carlo dose calculation algorithm in volumetric modulated arc therapy using Monaco TPS for three different clinical cases

AimTo study the dosimetric impact of statistical uncertainty (SU) per plan on Monte Carlo (MC) calculation in Monaco? treatment planning system (TPS) during volumetric modulated arc therapy (VMAT) for three different clinical cases.BackgroundDuring MC calculation SU is an important factor to decide dose calculation accuracy and calculation time. It is necessary to evaluate optimal acceptance of SU for quality plan with reduced calculation time.Materials and methodsThree different clinical cases as the lung, larynx, and prostate treated using VMAT technique were chosen. Plans were generated with Monaco? V5.11 TPS with 2% statistical uncertainty. By keeping all other parameters constant, plans were recalculated by varying SU, 0.5%, 1%, 2%, 3%, 4%, and 5%. For plan evaluation, conformity index (CI), homogeneity index (HI), dose coverage to PTV, organ at risk (OAR) dose, normal tissue receiving dose ≥5 Gy and ≥10 Gy, integral dose (NTID), calculation time, gamma pass rate, calculation reproducibility and energy dependency were analyzed.ResultsCI and HI improve as SU increases from 0.5% to 5%. No significant dose difference was observed in dose coverage to PTV, OAR doses, normal tissue receiving dose ≥5 Gy and ≥10 Gy and NTID. Increase of SU showed decrease in calculation time, gamma pass rate and increase in PTV max dose. No dose difference was seen in calculation reproducibility and dependent on energy.ConclusionFor VMAT plans, SU can be accepted from 1% to 3% per plan with reduced calculation time without compromising plan quality and deliverability by accepting variations in point dose within the target.     

Assigning alleles to different loci in amplifications of duplicated loci

Duplicated loci, for example those associated with major histocompatibility complex (MHC) genes, often have similar DNA sequences that can be coamplified with a pair of primers. This results in genotyping difficulties and inaccurate analyses. Here, we present a method to assign alleles to different loci in amplifications of duplicated loci. This method simultaneously considers several factors that may each affect correct allele assignment. These are the sharing of identical alleles among loci, null alleles, copy number variation, negative amplification, heterozygote excess or heterozygote deficiency, and linkage disequilibrium. The possible multilocus genotypes are extracted from the alleles for each individual and weighted to estimate the allele frequencies. The likelihood of an allele configuration is calculated and is optimized with a heuristic algorithm. Monte‐Carlo simulations and three empirical MHC data sets are used as examples to evaluate the efficacy of our method under different conditions. Our new software, mhc‐typer V1.1, is freely available at https://github.com/huangkang1987/mhc-typer .     

Structural insights into the catalytic mechanism of 5-methylthioribose 1-phosphate isomerase

5-Methylthioribose 1-phosphate isomerase (M1Pi) is a crucial enzyme involved in the universally conserved methionine salvage pathway (MSP) where it is known to catalyze the conversion of 5-methylthioribose 1-phosphate (MTR-1-P) to 5-methylthioribulose 1-phosphate (MTRu-1-P) via a mechanism which remains unspecified till date. Furthermore, although M1Pi has a discrete function, it surprisingly shares high structural similarity with two functionally non-related proteins such as ribose-1,5-bisphosphate isomerase (R15Pi) and the regulatory subunits of eukaryotic translation initiation factor 2B (eIF2B). To identify the distinct structural features that lead to divergent functional obligations of M1Pi as well as to understand the mechanism of enzyme catalysis, the crystal structure of M1Pi from a hyperthermophilic archaeon Pyrococcus horikoshii OT3 was determined. A meticulous structural investigation of the dimeric M1Pi revealed the presence of an N-terminal extension and a hydrophobic patch absent in R15Pi and the regulatory α-subunit of eIF2B. Furthermore, unlike R15Pi in which a kink formation is observed in one of the helices, the domain movement of M1Pi is distinguished by a forward shift in a loop covering the active-site pocket. All these structural attributes contribute towards a hydrophobic microenvironment in the vicinity of the active site of the enzyme making it favorable for the reaction mechanism to commence. Thus, a hydrophobic active-site microenvironment in addition to the availability of optimal amino-acid residues surrounding the catalytic residues in M1Pi led us to propose its probable reaction mechanism via a cis-phosphoenolate intermediate formation.     

中国植物受威胁等级评估系统概述

濒危物种保护是生物多样性保护工作的重要组成部分, 而物种受威胁等级评估则是濒危物种保护的方向指引。经过多年的发展, 物种受威胁等级的评估由定性评估逐渐向定量评估为主、定性评估为辅的方向发展。本文综述了国内植物受威胁等级定量评估系统的研究进展, 同时介绍了国外较为成熟的IUCN红色名录评估系统、CITES评估系统、美国自然保育协会评估系统, 提出未来制定受威胁物种定量评估标准时要兼顾以下方面: (1)等级设置定义要明确、统一且合理; (2)评估标准应该定量化、客观且不冗余; (3)评估系统应该适应不同地理范围, 最好能同时表达出各范围的受威胁等级; (4)评估指标要包含物种动态信息, 能定量分析物种在过去或者未来的变化。此外, 本文认为国内的物种受威胁等级定量评估系统应该形成规范化的大纲, 加大宣传力度, 尽量将理论研究与具体的保护行动结合起来; 同时, 我国还应该采用全球广泛应用的受威胁等级评估系统获取物种受威胁等级, 将国内生物多样性保护工作纳入到全球范围中去。     

江苏省县域森林生态安全评价及空间计量分析

通过考察生态区位因素对森林生态安全的影响,建立评价指标体系,研究其空间相关性的内在效应机制,从而实现森林生态安全的评价与监测。以江苏省80个区县为研究对象,基于2000—2015年面板数据,运用熵权法、专家法及模糊物元法计算森林生态安全指数,然后结合气象类指标及地形类指标计算生态区位系数,再用该系数修正森林ESI,同时结合Arc GIS技术、空间相关性分析、SLM与SEM模型得出如下结论:(1)人口密度、单位面积能源消耗量、退耕还林面积占比等指标权重最大;(2)生态区位系数高值区主要分布在江苏南部少数地区,低值区主要分布在江苏东北部;(3)苏南地区森林生态安全状况整体好于苏北及中东部地区;(4) 2000—2015年,江苏省67.5%的区县森林ESI呈现出较明显下降趋势,反映出江苏省森林生态安全发展状况不容乐观;(5)江苏省县域森林ESI整体空间相关性显著(P≤0.01),但2000—2015年空间聚集程度有所下降,且Low-Low聚类显著性更强;(6)森林ESI在江苏省县域间为扩散效应与回流效应并存。     

Smart process development: Application of machine-learning and integrated process modeling for inclusion body purification processes

The development of a biopharmaceutical production process usually occurs sequentially, and tedious optimization of each individual unit operation is very time-consuming. Here, the conditions established as optimal for one-step serve as input for the following step. Yet, this strategy does not consider potential interactions between a priori distant process steps and therefore cannot guarantee for optimal overall process performance. To overcome these limitations, we established a smart approach to develop and utilize integrated process models using machine learning techniques and genetic algorithms. We evaluated the application of the data-driven models to explore potential efficiency increases and compared them to a conventional development approach for one of our development products. First, we developed a data-driven integrated process model using gradient boosting machines and Gaussian processes as machine learning techniques and a genetic algorithm as recommendation engine for two downstream unit operations, namely solubilization and refolding. Through projection of the results into our large-scale facility, we predicted a twofold increase in productivity. Second, we extended the model to a three-step model by including the capture chromatography. Here, depending on the selected baseline-process chosen for comparison, we obtained between 50% and 100% increase in productivity. These data show the successful application of machine learning techniques and optimization algorithms for downstream process development. Finally, our results highlight the importance of considering integrated process models for the whole process chain, including all unit operations.     

Selection of relevant features for EEG signal classification of schizophrenic patients

In this paper, EEG signals of 20 schizophrenic patients and 20 age-matched control participants are analyzed with the objective of determining the more informative channels and finally distinguishing the two groups. For each case, 22 channels of EEG were recorded. A two-stage feature selection algorithm is designed, such that, the more informative channels are first selected to enhance the discriminative information. Two methods, bidirectional search and plus-L minus-R (LRS) techniques are employed to select these informative channels. The interesting point is that most of selected channels are located in the temporal lobes (containing the limbic system) that confirm the neuro-phychological differences in these areas between the schizophrenic and normal participants. After channel selection, genetic algorithm (GA) is employed to select the best features from the selected channels. In this case, in addition to elimination of the less informative channels, the redundant and less discriminant features are also eliminated. A computationally fast algorithm with excellent classification results is obtained. Implementation of this efficient approach involves several features including autoregressive (AR) model parameters, band power, fractal dimension and wavelet energy. To test the performance of the final subset of features, classifiers including linear discriminant analysis (LDA) and support vector machine (SVM) are employed to classify the reduced feature set of the two groups. Using the bidirectional search for channel selection, a classification accuracy of 84.62% and 99.38% is obtained for LDA and SVM, respectively. Using the LRS technique for channel selection, a classification accuracy of 88.23% and 99.54% is also obtained for LDA and SVM, respectively. Finally, the results are compared and contrasted with two well-known methods namely, the single-stage feature selection (evolutionary feature selection) and principal component analysis (PCA)-based feature selection. The results show improved accuracy of classification in relatively low computational time with the two-stage feature selection.