Grouping of amino acids and recognition of protein structurally conserved regions by reduced alphabets of amino acids

Sequence alignment is a common method for finding protein structurally conserved/similar regions. However, sequence alignment is often not accurate if sequence identities between to-be-aligned sequences are less than 30%. This is because that for these sequences, different residues may play similar structural roles and they are incorrectly aligned during the sequence alignment using substitution matrix consisting of 20 types of residues. Based on the similarity of physicochemical features, residues can be clustered into a few groups. Using such simplified alphabets, the complexity of protein sequences is reduced and at the same time the key information encoded in the sequences remains. As a result, the accuracy of sequence alignment might be improved if the residues are properly clustered. Here, by using a database of aligned protein structures (DAPS), a new clustering method based on the substitution scores is proposed for the grouping of residues, and substitution matrices of residues at different levels of simplification are constructed. The validity of the reduced alphabets is confirmed by relative entropy analysis. The reduced alphabets are applied to recognition of protein structurally conserved/similar regions by sequence alignment. The results indicate that the accuracy or efficiency of sequence alignment can be improved with the optimal reduced alphabet with N around 9. Supported by the National Natural Science Foundation of China (Grant Nos. 90403120, 10474041 and 10021001) and the Nonlinear Project (973) of the NSM     

蛋白质序列复杂性简化与非比对序列分析

非比对序列分析是最新发展的一种序列分析方法,具有计算效率高并适用于分析低相似性的序列,已成功用于DNA的序列分析中.但是由于蛋白质序列的复杂性,非比对序列分析对于蛋白质序列分析的准确度却不高.用将20种天然氨基酸残基归类的方法,简化了蛋白质序列的复杂性,并运用到对蛋白质的非比对序列分析中,有效地提高了序列分析的准确性.     

A workbench for multiple alignment construction and analysis

Multiple sequence alignment can be a useful technique for studying molecular evolution, as well as for analyzing relationships between structure or function and primary sequence. We have developed for this purpose an interactive program, MACAW (Multiple Alignment Construction and Analysis Workbench), that allows the user to construct multiple alignments by locating, analyzing, editing, and combining "blocks" of aligned sequence segments. MACAW incorporates several novel features. (1) Regions of local similarity are located by a new search algorithm that avoids many of the limitations of previous techniques. (2) The statistical significance of blocks of similarity is evaluated using a recently developed mathematical theory. (3) Candidate blocks may be evaluated for potential inclusion in a multiple alignment using a variety of visualization tools. (4) A user interface permits each block to be edited by moving its boundaries or by eliminating particular segments, and blocks may be linked to form a composite multiple alignment. No completely automatic program is likely to deal effectively with all the complexities of the multiple alignment problem; by combining a powerful similarity search algorithm with flexible editing, analysis and display tools, MACAW allows the alignment strategy to be tailored to the problem at hand.     

Conformational Study Based on Simulated Annealing and 1H NMR of Peptides Comprising the Consensus Sequence Arg5-Asp-Val-Arg-Gly9: Effects of the Substitution Ser 12 by Ala 12 or of 3 Residues Deletion at the N-Terminus

Summary

A conformational search by simulated annealing has been performed on two peptides derivated from the tetradecapeptide used to isolate the Xenopus laevis skin maturation RXVRG-endoprotease. The Ala 12 derivative, obtained by substitution in the hydrophobic C terminal fragment and the undecapeptide 4–14, obtained by deletion of an acidic rich tripeptide, were studied. No unique structure has been found for the tetradecapeptide Ala 12. This structural disorganization could explain the loss of activity of the endoprotease towards the subsituted peptide. For the undecapeptide, two different models in accordance with the NMR data were found. The conformational differences between these two models are locat ed in the consensus sequence and in each case an hairpin-like conformation is observed. These results could be related to the enhanced cleavage activity of the maturation enzyme. The obtained structures are also compared with those of the original tetradecapeptide.     

Structure‐based classification of FAD binding sites: A comparative study of structural alignment tools

A total of six different structural alignment tools (TM‐Align, TriangleMatch, CLICK, ProBis, SiteEngine and GA‐SI) were assessed for their ability to perform two particular tasks: (i) discriminating FAD (flavin adenine dinucleotide) from non‐FAD binding sites, and (ii) performing an all‐to‐all comparison on a set of 883 FAD binding sites for the purpose of classifying them. For the first task, the consistency of each alignment method was evaluated, showing that every method is able to distinguish FAD and non‐FAD binding sites with a high Matthews correlation coefficient. Additionally, GA‐SI was found to provide alignments different from those of the other approaches. The results obtained for the second task revealed more significant differences among alignment methods, as reflected in the poor correlation of their results and highlighted clearly by the independent evaluation of the structural superimpositions generated by each method. The classification itself was performed using the combined results of all methods, using the best result found for each comparison of binding sites. A number of different clustering methods (Single‐linkage, UPGMA, Complete‐linkage, SPICKER and k‐Means clustering) were also used. The groups of similar binding sites (proteins) or clusters generated by the best performing method were further analyzed in terms of local sequence identity, local structural similarity and conservation of analogous contacts with the FAD ligands. Each of the clusters was characterized by a unique set of structural features or patterns, demonstrating that the groups generated truly reflect the structural diversity of FAD binding sites. Proteins 2016; 84:1728–1747. © 2016 Wiley Periodicals, Inc.     

Paired natural cysteine mutation mapping: aid to constraining models of protein tertiary structure.

This paper discusses the benefit of mapping paired cysteine mutation patterns as a guide to identifying the positions of protein disulfide bonds. This information can facilitate the computer modeling of protein tertiary structure. First, a simple, paired natural-cysteine-mutation map is presented that identifies the positions of putative disulfide bonds in protein families. The method is based on the observation that if, during the process of evolution, a disulfide-bonded cysteine residue is not conserved, then it is likely that its counterpart will also be mutated. For each target protein, protein databases were searched for the primary amino acid sequences of all known members of distinct protein families. Primary sequence alignment was carried out using PileUp algorithms in the GCG package. To search for correlated mutations, we listed only the positions where cysteine residues were highly conserved and emphasized the mutated residues. In proteins of known three-dimensional structure, a striking pattern of paired cysteine mutations correlated with the positions of known disulfide bridges. For proteins of unknown architecture, the mutation maps showed several positions where disulfide bridging might occur.     

BCL::Align—Sequence alignment and fold recognition with a custom scoring function online

BCL::Align is a multiple sequence alignment tool that utilizes the dynamic programming method in combination with a customizable scoring function for sequence alignment and fold recognition. The scoring function is a weighted sum of the traditional PAM and BLOSUM scoring matrices, position-specific scoring matrices output by PSI-BLAST, secondary structure predicted by a variety of methods, chemical properties, and gap penalties. By adjusting the weights, the method can be tailored for fold recognition or sequence alignment tasks at different levels of sequence identity. A Monte Carlo algorithm was used to determine optimized weight sets for sequence alignment and fold recognition that most accurately reproduced the SABmark reference alignment test set. In an evaluation of sequence alignment performance, BCL::Align ranked best in alignment accuracy (Cline score of 22.90 for sequences in the Twilight Zone) when compared with Align-m, ClustalW, T-Coffee, and MUSCLE. ROC curve analysis indicates BCL::Align's ability to correctly recognize protein folds with over 80% accuracy. The flexibility of the program allows it to be optimized for specific classes of proteins (e.g. membrane proteins) or fold families (e.g. TIM-barrel proteins). BCL::Align is free for academic use and available online at http://www.meilerlab.org/.     

氨基酸高产菌株创制的关键技术

氨基酸是构成蛋白质的基本单元,被广泛应用于食品、医药、饲料和化工等多个领域。利用微生物细胞工厂生产氨基酸,具备原料可再生、生产过程条件温和、产品纯度高、环境污染小等优点,能够助力实现碳中和。借助代谢工程和合成生物学技术,对大肠杆菌(Escherichia coli)和谷氨酸棒杆菌(Corynebacterium glutamicum)进行定向设计、改造与优化,创制了高产氨基酸的微生物细胞工厂,实现了支链氨基酸、天冬氨酸族氨基酸、谷氨酸族氨基酸和芳香族氨基酸的生物炼制。本文对高产氨基酸的大肠杆菌细胞工厂和谷氨酸棒杆菌细胞工厂创制过程进行分析,以期对高性能微生物细胞工厂的创制提供参考。     

The role of pattern databases in sequence analysis

In the wake of the numerous now-fruitful genome projects, we are entering an era rich in biological data. The field of bioinformatics is poised to exploit this information in increasingly powerful ways, but the abundance and growing complexity both of the data and of the tools and resources required to analyse them are threatening to overwhelm us. Databases and their search tools are now an essential part of the research environment. However, the rate of sequence generation and the haphazard proliferation of databases have made it difficult to keep pace with developments. In an age of information overload, researchers want rapid, easy-to-use, reliable tools for functional characterisation of newly determined sequences. But what are those tools? How do we access them? Which should we use? This review focuses on a particular type of database that is increasingly used in the task of routine sequence analysis--the so-called pattern database. The paper aims to provide an overview of the current status of pattern databases in common use, outlining the methods behind them and giving pointers on their diagnostic strengths and weaknesses.     

A comparison of fluorescamine and o-phthaldialdehyde as effective blocking reagents in protein sequence analyses by the Beckman sequencer

Use of o-phthaldialdehyde to chemically reduce the newly generated amino termini responsible for the progressively increasing background during an extended amino acid sequence analysis in a liquid phase sequencer has been described. The results have been compared with Fluram blocking using apomyoglobin and rabbit C-reactive protein as standard and unknown samples, respectively.