A flexible method to align large numbers of biological sequences

Summary A method for the alignment of two or more biological sequences is described. The method is a direct extension of the method of Taylor (1987) incorporating a consensus sequence approach and allows considerable freedom in the control of the clustering of the sequences. At one extreme this is equivalent to the earlier method (Taylor 1987), whereas at the other, the clustering approaches the binary method of Feng and Doolittle (1987). Such freedom allows the program to be adapted to particular problems, which has the important advantage of resulting in considerable savings in computer time, allowing very large problems to be tackled. Besides a detailed analysis of the alignment of the cytochrome c superfamily, the clustering and alignment of the PIR sequence data bank (3500 sequences approx.) is described.     

The SKS of the KMSKS signature of class I aminoacyl-tRNA synthetases corresponds to the GKT/S sequence characteristic of the ATP-binding site of many proteins

On the basis of protein modification studies and primary structure comparison, we propose that the SKS sequence within the KMSKS signature of the class 1 aminoacyl-tRNA synthetases corresponds to the GKT(or S) sequence considered as a signature of the nucleotide triphosphate-binding site of many proteins.     

Designing repeat proteins: modular leucine-rich repeat protein libraries based on the mammalian ribonuclease inhibitor family

We present a novel approach to design repeat proteins of the leucine-rich repeat (LRR) family for the generation of libraries of intracellular binding molecules. From an analysis of naturally occurring LRR proteins, we derived the concept to assemble repeat proteins with randomized surface positions from libraries of consensus repeat modules. As a guiding principle, we used the mammalian ribonuclease inhibitor (RI) family, which comprises cytosolic LRR proteins known for their extraordinary affinities to many RNases. By aligning the amino acid sequences of the internal repeats of human, pig, rat, and mouse RI, we derived a first consensus sequence for the characteristic alternating 28 and 29 amino acid residue A-type and B-type repeats. Structural considerations were used to replace all conserved cysteine residues, to define less conserved positions, and to decide where to introduce randomized amino acid residues. The so devised consensus RI repeat library was generated at the DNA level and assembled by stepwise ligation to give libraries of 2-12 repeats. Terminal capping repeats, known to shield the continuous hydrophobic core of the LRR domain from the surrounding solvent, were adapted from human RI. In this way, designed LRR protein libraries of 4-14 LRRs (equivalent to 130-415 amino acid residues) were obtained. The biophysical analysis of randomly chosen library members showed high levels of soluble expression in the Escherichia coli cytosol, monomeric behavior as characterized by gel-filtration, and alpha-helical CD spectra, confirming the success of our design approach.     

Retroviral Oligonucleotide Distributions Correlate with Biased Nucleotide Compositions of Retrovirus Sequences,Suggesting a Duplicative Stepwise Molecular Evolution

A computer-assisted analysis was made of 24 complete nucleotide sequences selected from the vertebrate retroviruses to represent the ten viral groups. The conclusions of this analysis extend and strengthen the previously made hypothesis on the Moloney murine leukemia virus: The evolution of the nucleotide sequence appears to have occurred mainly through at least three overlapping levels of duplication: (1) The distributions of overrepresented (3–6)-mers are consistent with the universal rule of a trend toward TG/CT excess and with the persistence of a certain degree of symmetry between the two strands of DNA. This suggests one or several original tandemly repeated sequences and some inverted duplications. (2) The existence of two general core consensuses at the level of these (3–6)-mers supports the hypothesis of a common evolutionary origin of vertebrate retroviruses. Consensuses more specific to certain sequences are compatible with phylogenetic trees established independently. The consensuses could correspond to intermediary evolutionary stages. (3) Most of the (3–6)-mers with a significantly higher than average frequency appear to be internally repeated (with monomeric or oligomeric internal iterations) and seem to be at least partly the cause of the bias observed by other researchers at the level of retroviral nucleotide composition. They suggest a third evolutionary stage by slippage-like stepwise local duplications. Received: 3 January 1996 / Accepted: 27 March 1996     

Carrying capacity in a heterogeneous environment with habitat connectivity

A large body of theory predicts that populations diffusing in heterogeneous environments reach higher total size than if non‐diffusing, and, paradoxically, higher size than in a corresponding homogeneous environment. However, this theory and its assumptions have not been rigorously tested. Here, we extended previous theory to include exploitable resources, proving qualitatively novel results, which we tested experimentally using spatially diffusing laboratory populations of yeast. Consistent with previous theory, we predicted and experimentally observed that spatial diffusion increased total equilibrium population abundance in heterogeneous environments, with the effect size depending on the relationship between r and K. Refuting previous theory, however, we discovered that homogeneously distributed resources support higher total carrying capacity than heterogeneously distributed resources, even with species diffusion. Our results provide rigorous experimental tests of new and old theory, demonstrating how the traditional notion of carrying capacity is ambiguous for populations diffusing in spatially heterogeneous environments.     

Consensus and comprehensive linkage maps of bovine chromosome 24

This study describes development of a consensus genetic linkage map of bovine chromosome 24 (BTA24). Eight participating laboratories contributed data for 58 unique markers including a total of 25 409 meioses. Eighteen markers, which were typed in more than one reference population, were used as potential anchors to generate a consensus framework map. The framework map contained 16 loci ordered with odds greater than 1000:1 and spanned 79.3 cM. Remaining markers were included in a comprehensive map relative to these anchors. The resulting BTA24 comprehensive map was 98.3 cM in length. Average marker intervals were 6.1 and 2.5 cM for framework and comprehensive maps, respectively. Marker order was generally consistent with previously reported BTA24 linkage maps. Only one discrepancy was found when comparing the comprehensive map with the published USDA-MARC linkage map. Integration of genetic information from different maps provides a high-resolution BTA24 linkage map.     

A broad range of Fab stabilities within a host of therapeutic IgGs

Although the functional properties of IgGs are well known, little has been published concerning the stability of whole IgG molecules. Stability is, however, a requirement for the development of antibodies for therapeutic or diagnostic applications. The hypervariable antigen-binding region (Fv) is responsible for stability variations between IgGs of identical subclass. To determine the range of stabilities that may be expected for human(ized) antibodies, differential scanning calorimetry was performed on 17 human(ized) antibodies from various in-house programs. The antigen-binding fragments (Fabs) of these antibodies exhibited thermal unfolding transitions with midpoints (T(M)s) varying from 57 to 82 degrees C. Antibodies with very low Fab stabilities were found to aggregate and express poorly. Fab instability was often associated with high levels of uncommonly observed amino acids or CDR loop lengths particularly within the variable heavy chain domain. Overall, the study provides a thermostability range for IgGs and suggests possible stability guidelines for developing antibody diagnostics or therapeutics.     

Mapping hidden potential identity elements by computing the average discriminating power of individual tRNA positions

The recently published discrete mathematical method, extended consensus partition (ECP), identifies nucleotide types at each position that are strictly absent from a given sequence set, while occur in other sets. These are defined as discriminating elements (DEs). In this study using the ECP approach, we mapped potential hidden identity elements that discriminate the 20 different tRNA identities. We filtered the tDNA data set for the obligatory presence of well-established tRNA features, and then separately for each identity set, the presence of already experimentally identified strictly present identity elements. The analysis was performed on the three kingdoms of life. We determined the number of DE, e.g. the number of sets discriminated by the given position, for each tRNA position of each tRNA identity set. Then, from the positional DE numbers obtained from the 380 pairwise comparisons of the 20 identity sets, we calculated the average excluding value (AEV) for each tRNA position. The AEV provides a measure on the overall discriminating power of each position. Using a statistical analysis, we show that positional AEVs correlate with the number of already identified identity elements. Positions having high AEV but lacking published identity elements predict hitherto undiscovered tRNA identity elements.     

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.     

pH对毕赤酵母表达重组人复合a干扰素的降解影响

使用Pichia pastoris表达重组人复合a干扰素(cIFN)会发生降解、聚合等不均一表达的现象。在5 L发酵罐中考察了不同诱导pH对cIFN表达产生降解的影响, 结果发现在适合酵母生长的pH 3.0~7.0范围内, 当诱导pH为4.0~5.0时, cIFN不均一表达现象最少, 生物活性达到2.5×108 IU/mL。通过测定发酵液中总蛋白酶活和细胞活性寻找了cIFN降解出现的原因：发现低诱导pH下细胞死亡率升高释放更多酶系, 高诱导pH下蛋白酶活性明显增大, 两者都使蛋白酶作用加强, 加剧cIFN的降解; 特别是诱导pH为7.0时, 适宜的pH使蛋白酶酶活陡升, 将cIFN完全降解。