A protein alignment scoring system sensitive at all evolutionary distances

Summary Protein sequence alignments generally are constructed with the aid of a substitution matrix that specifies a score for aligning each pair of amino acids. Assuming a simple random protein model, it can be shown that any such matrix, when used for evaluating variable-length local alignments, is implicitly a log-odds matrix, with a specific probability distribution for amino acid pairs to which it is uniquely tailored. Given a model of protein evolution from which such distributions may be derived, a substitution matrix adapted to detecting relationships at any chosen evolutionary distance can be constructed. Because in a database search it generally is not known a priori what evolutionary distances will characterize the similarities found, it is necessary to employ an appropriate range of matrices in order not to overlook potential homologies. This paper formalizes this concept by defining a scoring system that is sensitive at all detectable evolutionary distances. The statistical behavior of this scoring system is analyzed, and it is shown that for a typical protein database search, estimating the originally unknown evolutionary distance appropriate to each alignment costs slightly over two bits of information, or somewhat less than a factor of five in statistical significance. A much greater cost may be incurred, however, if only a single substitution matrix, corresponding to the wrong evolutionary distance, is employed.     

Unbiased estimation of the rates of synonymous and nonsynonymous substitution

Summary The current convention in estimating the number of substitutions per synonymous site (K _S ) and per nonsynonymous site (K _A ) between two protein-coding genes is to count each twofold degenerate site as one-third synonymous and two-thirds nonsynonymous because one of the three possible changes at such a site is synonymous and the other two are nonsynonymous. This counting rule can considerably overestimate theK _S value because transitional mutations tend to occur more often than transversional mutations and because most transitional mutations at twofold degenerate sites are synonymous. A new method that gives unbiased estimates is proposed. An application of the new and the old method to 14 pairs of mouse and rat genes shows that the new method gives aK _S value very close to the number of substitutions per fourfold degenerate site whereas the old method gives a value 30% higher. Both methods give aK _A value close to the number of substitutions per nondegenerate site.     

Time and biosequences

In this paper we discuss and demonstrate the importance of several factors relative to the relationship between time and evolution of biosequences. In both quantitative and qualitative measurements of the genetic distances, the compositional constraints of the nucleotide sequences play a very important role. We demonstrate that when homologous sequences significantly differ in base composition we get erratic branching order and/or wrong evaluation of the evolutionary rates. We must consider that every gene may have a different evolutionary dynamic along its sequence, generally linked to its functional constraints; this too can seriously affect its clocklike behavior. We report some cases showing how these factors can affect the quantitative measurements of the genetic distances of biosequences. Presented at the NATO Advanced Research Workshop onGenome Organization and Evolution, Spetsai, Greece, 16–22 September 1992     

Sequential1H-NMR assignments of neurotoxin III from the sea anemoneHeteractis macrodactylus and structural comparison with related toxins

The complete sequence-specific assignment of resonances in the¹H-NMR spectrum of the polypeptide neurotoxin III (Hm III) from the sea anemoneHeteractis macrodactylus is described. Comparison of the chemical shifts and pattern of NOEs for Hm III with those for the related toxin Hp III fromHeteractis paumotensis, which differs only in the substitution of Asn for Tyr at position 11, shows that the overall secondary and tertiary structures are conserved. The largest differences in chemical shift caused by the substitution at position 11 are observed for the NH resonances of Arg-13, Thr-14, Ala-15, Leu-17, and Cys-26. The CH resonances influenced most are those of ASP-6, Gly-9, Leu-17, and Glu-42, while the most affected CH resonances are from Leu-17, Glu-28, and Lys-32. The absence of long-range NOEs to the aromatic ring of Tyr-11 as well as the lack of significant chemical shift effects on residues outside the loop comprising residues 7–16 confirm that this part of the loop makes no long-lived contacts with the rest of the molecule. The deviations from random coil shifts of Hm III are compared with those of the related anemone toxins Hp II, Hp III, and toxin I fromStichodactyla helianthus (Sh I). The similarity in deviations in chemical shift as a function of sequence position for these four toxins emphasizes the overall structural homology among these polypeptides.     

835型氨基酸分析仪分析中停电的处理方法

本文介绍了日立835-50型氨基酸分析仪分析过程中遇到停电时,仪器管路的冲洗方法。     

基于遗传算法的土地利用优化：NSGA-Ⅱ和NSGA-Ⅲ的对比研究

土地利用优化通常要兼顾不同群体的多种要求,理论上是复杂的超多目标(4个及以上)优化问题。但实际操作中却往往被简化为多目标(2—3个)优化问题,通过一种流行的多目标优化算法第Ⅱ代非支配排序遗传算法(NSGA-Ⅱ)求解。究其原因是对超多目标优化算法认知的缺失和与多目标优化算法理论对比的匮乏。对NSGA系列中应用最广泛的多目标优化算法NSGA-Ⅱ和最新提出、面向超多目标优化的算法NSGA-Ⅲ进行探究,从理论和实验两方面对Ⅲ和Ⅱ进行对比,从而探究二者进行土地利用优化时的优劣。在理论上,对比两种算法原理的异同。在实验中,分别设计多目标(3个目标)和超多目标(13个目标)土地利用优化问题,利用两种算法进行求解。对实验结果采用四层架构、六大指标进行全面评价,以对比两种算法的可用性。理论对比发现,两个算法只有种群多样性保护的方法不同,其中NSGA-Ⅲ是基于与固定的参考点的距离,而NSGA-Ⅱ则是基于相邻解间的距离。通过实验对比发现,NSGA-Ⅲ在超多目标优化时运算速度快,且产生的最优方案实用价值更高,NSGA-Ⅱ在算法的有效性方面更有优势。     

聚对苯二甲酸乙二醇酯水解酶检测方法的研究进展

聚对苯二甲酸乙二醇酯(polyethylene terephthalate,PET)是应用最广泛的合成聚酯之一。由于PET不易降解,在环境中积累,对陆地、水生生态系统以及人类健康构成严重威胁。基于生物酶催化的生物降解策略为PET回收利用提供了一种绿色途径,在过去20年间,已发现了多种PET水解酶,并通过蛋白质工程等手段来改善这些酶的降解性能,但是目前仍未找到适合大规模工业应用的PET水解酶。利用传统的检测方法筛选PET水解酶是一个缓慢而复杂的过程。为了促进PET酶法回收的工业化应用,需要研发高效的检测方法。近年来,研究人员开发了多种表征PET水解酶的分析方法。本文总结了可用于筛选PET水解酶的检测方法,如高效液相色谱法、紫外吸光度法和荧光激活液滴分选法等,并对其在筛选PET水解酶的应用方面进行了展望。     

Drosophila acidic ribosomal protein rpA2: Sequence and characterization

A cDNA encoding the Drosophila melanogaster acidic ribosomal protein rpA2 was cloned and sequenced. rpA2 is homologous to the Artemia salina acidic ribosomal protein eL12′. In situ hybridization to salivary gland polytene chromosomes localizes the rpA2 gene to band 21C. It is a single copy gene, with an mRNA of 0.8 kb. Two-dimensional gel electrophoresis of Drosophila ribosomal proteins followed by immuno-blotting showed that the rpA2 protein has an apparent relative mobility in SDS of 17 kD and an isoelectric point less than pH 5.0. Although the Drosophila gene rp21C may be the same as rpA2, the reported sequences differ. Comparisons of the aligned nucleotide sequences coding for the acidic ribosomal proteins rpA1 and rpA2 of Drosophila with those of other eukaryotes support the view of two separate, though closely related, groups of acidic proteins. Comparison with the Artemia homologues suggests that nucleotide identity may have been conserved by some constraint that acts in addition to the requirement for substantial similarity of amino acid sequences. © 1993 Wiley-Liss, Inc.