A new topological method to measure protein structure similarity

A method for the quantitative evaluation of structural similarity between protein pairs is developed that makes use of a Delaunay-based topological mapping. The result of the mapping is a three-dimensional array which is representative of the global structural topology and whose elements can be used to construe an integral scoring scheme. This scoring scheme was tested for its dependence on the protein length difference in a pairwise comparison, its ability to provide a reasonable means for structural similarity comparison within a family of structural neighbors of similar length, and its sensitivity to the differences in protein conformation. It is shown that such a topological evaluation of similarity is capable of providing insight into these points of interest. Protein structure comparison using the method is computationally efficient and the topological scores, although providing different information about protein similarity, correlate well with the distance root-mean-square deviation values calculated by rigid-body structural alignment.     

Structure and patterns of sequence variation in the mitochondrial DNA control region of the great cats

Mitochondrial DNA control region structure and variation were determined in the five species of the genus Panthera. Comparative analyses revealed two hypervariable segments, a central conserved region, and the occurrence of size and sequence heteroplasmy. As observed in the domestic cat, but not commonly seen in other animals, two repetitive sequence arrays (RS-2 with an 80-bp motif and RS-3 with a 6-10-bp motif) were identified. The 3' ends of RS-2 and RS-3 were highly conserved among species, suggesting that these motifs have different functional constraints. Control region sequences provided improved phylogenetic resolution grouping the sister taxa lion (Panthera leo) and leopard (Panthera pardus), with the jaguar (Panthera onca).     

Structural elucidation of the EPS of slime producing Brevundimonas vesicularis sp. isolated from a paper machine

The slime forming bacteria Brevundimonas vesicularis sp. was isolated from a paper mill and its EPS was produced on laboratory scale. After production, the exopolysaccharide (EPS) was purified and analysed for its purity and homogeneity, HPSEC revealed one distinct population with a molecular mass of more than 2,000 kDa. The protein content was around 9 w/w%. The sample was analysed to determine its chemical structure. The EPS was found to consist of rhamnose, glucose, galacturonic acid and glucuronic acid. Due to the presence of uronic acids the molar ratio between the four sugars found varies from 3:5:2:4 by sugar composition analyses after methanolysis to 1:1:1:1 found by NMR. A repeating unit with a molecular mass of 678 Da was confirmed by MALDI-TOF mass spectrometry after mild acid treatment. 13C and 1H hetero- and homonuclear 2D NMR spectroscopy of the native and partial hydrolysed EPS revealed a repeating unit, no non-sugar substituents were present.     

Faecal genetic analysis to determine the presence and distribution of elusive carnivores: design and feasibility for the Iberian lynx

Noninvasive methods using genetic markers have been suggested as ways to overcome difficulties associated with documenting the presence of elusive species. We present and assess a novel, reliable and effective molecular genetic technique for the unequivocal genetic identification of faeces from the endangered Iberian lynx (Lynx pardinus). From mitochondrial DNA (mtDNA) cytochrome b and D-loop region sequences, we designed four species-specific primers (for products 130-161 bp long) that were considered to be likely to amplify degraded DNA. We compared two DNA extraction methods, various DNA amplification conditions and the robustness and specificity of the primer pairs with 87 lynx samples from 5 potentially different lynx populations and with 328 samples of other carnivore species. The utility of the identification technique was tested with faeces of different ages, with faeces from controlled field experiments, and with faeces collected from locales with possible lynx populations from throughout the state of Andalusia, Spain (8052 km2). Faecal mtDNA extraction was more efficient using PBS wash of the faeces instead of a faeces homogenate. Our assay increased from 92.6 to 99% efficiency with a second amplification and a reduction in template concentration to overcome polymerase chain reaction (PCR) inhibition. Our assay never produced false positives, and correctly identified all lynx faeces. Of 252 faeces samples of unknown species collected throughout Andalusia, 26.6% (from three different areas) were classified as Iberian lynx, 1.4% showed evidence of PCR inhibition and 1.2% were of uncertain origin. This method has proven to be a reliable technique that can be incorporated into large-scale surveys of Iberian lynx populations and exemplifies an approach that can easily be extended to other species.     

Use of SAMPL for a study of DNA polymorphism,genetic diversity and possible gene tagging in bread wheat

Selective Amplification of Microsatellite Polymorphic Loci (SAMPL) technology was used in bread wheat for the first time for a study of genetic diversity, genotype identification and gene tagging. The diversity studies involved 55 wheat genotypes and two SAMPL primer pairs (SAMPL-6 and SAMPL-7, each with a M-CAG primer), which together gave 43 polymorphic bands out of a total of 87 SAMPL bands. The average polymorphic information content (PIC) of SAMPL primers was 0.221 and that of SAMPL markers was 0.264. The marker index of SAMPL markers was 9.61. The genetic similarity (GS) coefficients for 1,485 pairs of genotypes ranged from 0.35 to 0.96 with an average of 0.65. A dendrogram was prepared on the basis of a similarity matrix using the UPGMA algorithm, which corresponded well with the results of principal component analysis (PCA). From a total of 55 genotypes, 54 could be distinguished using the SAMPL banding patterns of both primers. For gene tagging, 568 bands from a total of 1,185 SAMPL bands detected polymorphism between each of the three pairs of parents differing for grain protein content (GPC), pre-harvest sprouting tolerance (PHST) and grain weight (GW). An association of six bands with GPC, of seven bands with PHST and four bands with GW was observed using bulked segregant analysis (BSA). Received: 5 April 2001 / Accepted: 17 May 2001     

Crystal structure and nucleotide sequence of an anionic trypsin from chum salmon (Oncorhynchus keta) in comparison with Atlantic salmon (Salmo salar) and bovine trypsin

The nucleotide sequence and crystal structure of chum salmon trypsin (CST) are now reported. The cDNA isolated from the pyloric caeca of chum salmon encodes 222 amino acid residues, the same number of residues as the anionic Atlantic salmon trypsin (AST), but one residue less than bovine beta-trypsin (BT). The net charge on CST determined from the sum of all charged amino acid side-chains is -3. There are 79 sequence differences between CST and BT, but only seven sequence differences between CST and AST. Anionic CST isolated from pyloric caeca has also been purified and crystallized; the structure of the CST-benzamidine complex has been determined to 1.8A resolution. The overall tertiary structure of CST is similar to that of AST and BT, but some differences are observed among the three trypsins. The most striking difference is at the C terminus of CST, where the expected last two residues are absent. The absence of these residues likely increases the flexibility of CST by the loss of important interactions between the N and C-terminal domains. Similarly, the lack of Tyr151 in CST (when compared with BT) allows more space for Gln192 in the active site thereby increasing substrate accessibility to the binding pocket. Lys152 in CST also adopts the important role of stabilizing the loop from residue 142 to 153. These observations on CST provide a complementary view of a second cold-adapted trypsin, which in comparison with the structures of AST and BT, suggest a structural basis for differences in enzymatic activity between enzymes from cold-adapted species and mammals.     

Quantifying the similarities within fold space

We have used GRATH, a graph-based structure comparison algorithm, to map the similarities between the different folds observed in the CATH domain structure database. Statistical analysis of the distributions of the fold similarities has allowed us to assess the significance for any similarity. Therefore we have examined whether it is best to represent folds as discrete entities or whether, in fact, a more accurate model would be a continuum wherein folds overlap via common motifs. To do this we have introduced a new statistical measure of fold similarity, termed gregariousness. For a particular fold, gregariousness measures how many other folds have a significant structural overlap with that fold, typically comprising 40% or more of the larger structure. Gregarious folds often contain commonly occurring super-secondary structural motifs, such as beta-meanders, greek keys, alpha-beta plait motifs or alpha-hairpins, which are matching similar motifs in other folds. Apart from one example, all the most gregarious folds matching 20% or more of the other folds in the database, are alpha-beta proteins. They also occur in highly populated architectural regions of fold space, adopting sandwich-like arrangements containing two or more layers of alpha-helices and beta-strands.Domains that exhibit a low gregariousness, are those that have very distinctive folds, with few common motifs or motifs that are packed in unusual arrangements. Most of the superhelices exhibit low gregariousness despite containing some commonly occurring super-secondary structural motifs. In these folds, these common motifs are combined in an unusual way and represent a small proportion of the fold (<10%). Our results suggest that fold space may be considered as continuous for some architectural arrangements (e.g. alpha-beta sandwiches), in that super-secondary motifs can be used to link neighbouring fold groups. However, in other regions of fold space much more discrete topologies are observed with little similarity between folds.     

Comparative Molecular Evolution of Primary (<Emphasis Type="BoldItalic">Buchnera</Emphasis>) and Secondary Symbionts of Aphids Based on Two Protein-Coding Genes

A+T content, phylogenetic relationships, codon usage, evolutionary rates, and ratio of synonymous versus non-synonymous substitutions have been studied in partial sequences of the atpD and aroQ/pheA genes of primary (Buchnera) and secondary symbionts of aphids and a set of selected non-symbiotic bacteria, belonging to the five subdivisions of the Proteobacteria. Compared to the homologous genes of the last group, both genes belonging to Buchnera behave in a similar way, showing a higher A+T content, forming a monophyletic group, a loss in codon bias, especially in third base position, an evolutionary acceleration and an increase in the number of non-synonymous substitutions, confirming previous results reported elsewhere for other genes. When available, these properties have been partly observed with the secondary symbionts, but with values that are intermediate between Buchnera and free living Proteobacteria. They show high A+T content, but not as high as Buchnera, a non-solved phylogenetic position between Buchnera, and the other γ-Proteobacteria, a loss in codon bias, again not as high as in Buchnera and a significant evolutionary acceleration in the case of the three atpD genes, but not when considering aroQ/pheA genes. These results give support to the hypothesis that they are symbionts at different stages of the symbiotic accommodation to the host.     

The PRINTS database: a resource for identification of protein families

The PRINTS database houses a collection of protein fingerprints, which may be used to assign family and functional attributes to uncharacterised sequences, such as those currently emanating from the various genome-sequencing projects. The April 2002 release includes 1,700 family fingerprints, encoding approximately 10,500 motifs, covering a range of globular and membrane proteins, modular polypeptides and so on. Fingerprints are groups of conserved motifs that, taken together, provide diagnostic protein family signatures. They derive much of their potency from the biological context afforded by matching motif neighbours; this makes them at once more flexible and powerful than single-motif approaches. The technique further departs from other pattern-matching methods by readily allowing the creation of fingerprints at superfamily-, family- and subfamily-specific levels, thereby allowing more fine-grained diagnoses. Here, we provide an overview of the method of protein fingerprinting and how the results of fingerprint analyses are used to build PRINTS and its relational cousin, PRINTS-S.     

Analysis of bilateral inverse symmetry in whole bacterial chromosomes

The positions of the 64 DNA tri-nucleotides (triplets) along the Borrelia burgdorferi chromosome were determined and cumulative position plots (CPP) were obtained. Analysis of CPP for complementary triplets revealed close correlations in complementary triplet frequencies (CTF) between opposing leading and lagging strands. Such bilateral inverse symmetry (BIS) applied also to complementary mono- and di-nucleotides and to some >3 n-tuples. At the level of individual bases BIS explains Chargaff's second parity rule for whole bacterial chromosomes. Using shuffled control sequences we show that single-base BIS was not the source of higher-order BIS. Analysis of CTF in 45 other chromosomes suggests that BIS is a general property of eubacteria. BIS at the various levels may be due to the very similar numbers of codons used in chromosomal halves. Evolutionarily, BIS could have resulted from asymmetric substitution of bases combined with genetic rearrangements. However, the provocative theoretical alternative of whole-genome inverse duplication is here considered.