自动编码器方法的蛋白质二级结构预测

蛋白质二级结构预测是进行蛋白质三级结构研究的重要基础,氨基酸的编码方式对二级结构预测有一定的影响。本文应用了一种新的组合编码方式,即将基团编码与位置特异性打分矩阵(PSSM)进行组合的编码方式。本文中提出的基团编码是针对氨基酸的一种新的编码方式,基团编码是根据氨基酸内部组成来进行编码的,由42位属性组成。本文选取位置特异性打分矩阵(PSSM)中的Blosum62进化矩阵和新的基团编码进行组合,形成新的编码方式。然后对CB513和25pdb两组数据分别进行实验。本文中将采用贝叶斯分类器与自动编码器两种方法来对这种新的编码方式进行实验,然后比较这两种方法得到的两组数据的结果。可以很明显的发现采用自动编码器的实验结果要比使用贝叶斯分类器的结果要高出1.65%。在本文的实验中,可以提取特征的自动编码器的预测准确率更好。     

A genetic algorithm-based protocol for docking ensembles of small ligands using experimental restraints

Summary A genetic algorithm (GA) based method for docking ensembles of small, flexible ligands to receptor proteins using NMR-derived constraints is described. In this method, three translations and rotations of the ligand and the dihedral angles of the ligand are represented by binary strings and evolve under the genetic operators of cross-over, mutation, migration and selection. The fitness function for the selection process includes distance and dihedral restraints and a repulsive van der Waals term. The GA was applied to a three-atom model system as well as to the streptavidin-biotin complex using simulated intermolecular distance restraints. In both systems, the GA was able to obtain low-energy conformations when only a single binding site was simulated. Calculations were also performed using distance restraints from two distinct binding sites. In these simulations, the GA was able to obtain low-energy conformations corresponding to ligand molecules in each of the two sites. The inclusion of additional ligands in the ensemble did not result in an energetic benefit, confirming that only two ligand conformations were necessary to fulfill the distance restraints. This method allows for a direct investigation of the minimum number of ligand orientations necessary to fulfill experimental distance restraints, and simultaneously yields detailed structural information about each site.     

A dry season streamflow reconstruction of the critically endangered Formosan landlocked salmon habitat

This study presents a 259-year (1750–2008 CE) reconstruction for October to February (ONDJF) standardized streamflow index of the Chichiawan Stream, the sole habitat of the critically endangered Formosan landlocked salmon (Oncorhynchus formosanus) remaining natural population. The reconstruction was based on the earlywood ring-width variations of Taiwan Douglas-fir (Pseudotsuga wilsoniana), using ensemble empirical mode decomposition for chronology development. Results showed that a higher average ONDJF streamflow in the 1968–2008 period was positively associated with a wider Taiwan Douglas-fir earlywood ring. During that period, the average ONDJF streamflow was related to the Pacific Decadal Oscillation and El Niño-Southern Oscillation development. Spectral analysis indicated that the reconstructed index had significant peaks around 2∼7 years and less significant peaks around 20–35 years. Based on the derived and reconstructed indices, most of the average ONDJF streamflow between 1750 and 2008 was within the normal range. Furthermore, there was no statistical evidence to suggest that the flow regime of the last few decades had changed, which refuted the speculation that the dry season streamflow natural variability in the 1960s and 1970s might have contributed toward the critically endangered status of Formosan landlocked salmon. The results of this study will be useful for developing short- to mid-term strategies to protect the endangered species and manage water resources.     

Toward a molecular taxonomy

Summary The concepts of a super information source and ensemble averaging are used to estimate the amount of information stored in protein andt-RNA sequences. Specifically applied to cytochrome c and hemoglobins, information measures analogous to those found to be highly significant for DNA pair frequency data (D ₂ vs.R) by Gatlin (1968) prove to be extremely highly correlated with organism complexity. Super source stability and the possible taxonomic utility of the extraordinary clusterings obtained are discussed. A restriction on the construction of ancestral sequences and a possible handle on homology are also detailed.     

Homozygosity mapping identifies a GALK1 mutation as the cause of autosomal recessive congenital cataracts in 4 adult siblings

Objective

Monogenic congenital cataract is one of the most genetically heterogeneous ocular conditions with almost 30 different genes involved in its etiology. In adult patients, genotype–phenotype correlations are troubled by eye surgery during infancy and/or long-term ocular complications. Here, we describe the molecular diagnosis of GALK1 deficiency as the cause of autosomal recessive congenital cataract in a family from Costa Rica.

Methods

Four affected siblings were included in the study. All of them underwent eye surgery during the first decade but medical records were not available. Congenital cataract was diagnosed by report. Molecular analysis included genome wide homozygosity mapping using a 250 K SNP Affymetrix microarray followed by PCR amplification and direct nucleotide sequencing of candidate gene.

Results

Genome wide homozygosity mapping revealed a 6 Mb region of homozygosity shared by two affected siblings at 17q25. The GALK1 gene was included in this interval and direct sequencing of this gene revealed a homozygous c.1144C>T mutation (p.Q382*) in all four affected subjects.

Conclusions

This work demonstrates the utility of homozygosity mapping in the retrospective diagnosis of a family with congenital cataracts in which ocular surgery at early age, the lack of medical records, and the presence of long term eye complications, impeded a clear clinical diagnosis during the initial phases of evaluation.     

Ensemble modeling for strain development of l-lysine-producing Escherichia coli

One of the main strategies to improve the production of relevant metabolites has been the manipulation of single or multiple key genes in the metabolic pathways. This kind of strategy requires several rounds of experiments to identify enzymes that impact either yield or productivity. The use of mathematical tools to facilitate this process is desirable. In this work, we apply the Ensemble Modeling (EM) framework, which uses phenotypic data (effects of enzyme overexpression or genetic knockouts on the steady-state production rate) to screen for potential models capable of describe existing data and thus gaining insight to improve strains for l-lysine production. Described herein is a strategy to generate a set of kinetic models that describe a set of enzyme overexpression phenotypes previously determined in an Escherichia coli strain that produces increased levels of l-lysine in an industrial laboratory. This final ensemble of models captures the kinetic characteristics of the cell through screening of phenotypes after sequential overexpression of enzymes. Furthermore, these models demonstrate some predictive capability, as starting from the reference producing strain (overexpressing desensitized dihydrodipicolinate synthetase (dapA*)) this set of models is able to predict that the desensitization of aspartate kinase (lysC*) is the next rate-controlling step in the l-lysine pathway. Moreover, this set of models allows for the generation of further targets for testing, for example, phosphoenolpyruvate (Ppc), aspartate aminotransferase (AspC), and glutamate dehydrogenase (GdhA). This work demonstrates the usefulness, applicability, and scope that the Ensemble Modeling framework offers to build production strains.     

朴素贝叶斯判别的判别效用分析

通过模拟实验研究了朴素贝叶斯判别的独立性假定对误判率的影响.实验结果表明:对于两个二维正态总体,当两类的距离为1、2、3时,在真实密度下进行判别的误判率随类内特征间相关系数的增加而减小,而朴素贝叶斯判别的误判率与特征间相关系数变化无关,因此随着类内特征间相关系数的增加,朴素贝叶斯判别的误判率明显增加,当类内特征间相关系数均为0.8时,朴素贝叶斯判别的误判率分别提高43.8%、220.8%和785.7%;对于高维正态总体模拟实验显示出同样的规律.     

基于支持向量机融合网络的蛋白质折叠子识别研究

在不依赖于序列相似性的条件下,蛋白质折叠子识别是一种分析蛋白质结构的重要方法.提出了一种三层支持向量机融合网络,从蛋白质的氨基酸序列出发,对27类折叠子进行识别.融合网络使用支持向量机作为成员分类器,采用“多对多”的多类分类策略,将折叠子的6种特征分为主要特征和次要特征,构建了多个差异的融合方案,然后对这些融合方案进行动态选择得到最终决策.当分类之前难以确定哪些参与组合的特征种类能够使分类结果最好时,提供了一种可靠的解决方案来自动选择特征信息互补最大的组合,保证了最佳分类结果.最后,识别系统对独立测试样本的总分类精度达到61.04%.结果和对比表明,此方法是一种有效的折叠子识别方法.     

Proteins at Work: A COMBINED SMALL ANGLE X-RAY SCATTERING AND THEORETICAL DETERMINATION OF THE MULTIPLE STRUCTURES INVOLVED ON THE PROTEIN KINASE FUNCTIONAL LANDSCAPE*

C-terminal Src kinase (Csk) phosphorylates and down-regulates the Src family tyrosine kinases (SFKs). Crystallographic studies of Csk found an unusual arrangement of the SH2 and SH3 regulatory domains about the kinase core, forming a compact structure. However, recent structural studies of mutant Csk in the presence of an inhibitor indicate that the enzyme accesses an expanded structure. To investigate whether wt-Csk may also access open conformations we applied small angle x-ray scattering (SAXS). We find wt-Csk frequently occupies an extended conformation where the regulatory domains are removed from the kinase core. In addition, all-atom structure-based simulations indicate Csk occupies two free energy basins. These basins correspond to ensembles of distinct global conformations of Csk: a compact structure and an extended structure. The transitions between these structures are entropically driven and accessible via thermal fluctuations that break local interactions. We further characterized the ensemble by generating theoretical scattering curves for mixed populations of conformations from both basins and compared the predicted scattering curves to the experimental profile. This population-combination analysis is more consistent with the experimental data than any rigid model. It suggests that Csk adopts a broad ensemble of conformations in solution, populating extended conformations not observed in the crystal structure that may play an important role in the regulation of Csk. The methodology developed here is broadly applicable to biological macromolecules and will provide useful information about what ensembles of conformations are consistent with the experimental data as well as the ubiquitous dynamic reversible assembly processes inherent in biology.