Issues in incorporation semantic integrity in molecular biological object-oriented databases.

Issues critical to ensuring semantic integrity in molecular biological data collections have been identified and include complexity, exceptions, missing data, changing models, holism and integration, delocalized data, interoperability and nomenclature. This combination is peculiar to biology and presents some interesting problems as a result. Little is known about semantic checking in object-oriented databases in general, but because such technology appears highly suitable for modeling biological data, it is appropriate to examine the ways in which object-oriented technology can support this functionality. It is concluded that object-oriented technology will support semantic checking even in a complex domain like biology. We propose 10 guidelines for future work including ways of treating exceptional cases and 'positioning' of constraints in a schema.

Deducing protein structures using logic programming: exploiting minimum data of diverse types.

The extent to which a protein can be modeled from constraint data depends on the amount and quality of the data. This report quantifies a relationship between the amount of data and the achievable model resolution. In an information-theoretic framework the number of bits of information per residue needed to constrain a solution was calculated. The number of bits provided by different kinds of constraints was estimated from a tetrahedral lattice where all unique molecules of 6, 9, ..., 21 atoms were enumerated. Subsets of these molecules consistent with different constraint sets were then chosen, counted, and the root-mean-square distance between them calculated. This provided the desired relations. In a discrete system the number of possible models can be severely limited with relatively few constraints. An expert system that can model a protein from data of different types was built to illustrate the principle and was tested using known proteins as examples. C-alpha resolutions of 5 A are obtainable from 5 bits of information per amino acid and, in principle, from data that could be rapidly collected using standard biophysical techniques.

CDR3 length in antigen-specific immune receptors.

In both immunoglobulins (Ig) and T cell receptors (TCR), the rearrangement of V, D, and J region sequence elements during lymphocyte maturation creates an enormous degree of diversity in an area referred to as the complementarity determining region 3 (CDR3) loop. Variations in the particular V, D, and J elements used, precise points of recombination, and random nucleotide addition all lead to extensive length and sequence heterogeneity. CDR3 loops are often critical for antigen binding in Igs and appear to provide the principal peptide binding residues in TCRs. To better understand the physical and selective constraints on these sequences, we have compiled information on CDR3 size variation for Ig H, L (kappa and lambda) and TCR alpha, beta, gamma, and delta. Ig H and TCR delta CDR3s are the most variable in size and are significantly longer than L and gamma chains, respectively. In contrast, TCR alpha and beta chain distributions are highly constrained, with nearly identical average CDR3 lengths, and their length distributions are not altered by thymic selection. Perhaps most significantly, these CDR3 length profiles suggest that gamma/delta TCRs are more similar to Igs than to alpha/beta TCRs in their putative ligand binding region, and thus gamma/delta and alpha/beta T cells may have fundamentally different recognition properties.

Weighting in sequence space: a comparison of methods in terms of generalized sequences.

Four methods for weighting aligned biological sequences have recently appeared that differ mathematically, philosophically, and in their results. Thus, while there is consensus about the need to weight sequences, the method to use is contentious. A geometric analysis based on a continuous sequence space is presented that provides a common framework in which to compare the methods. It is concluded that there are two "best" methods. When the sequences are known to be phylogenetically related and a tree can be generated without introducing excessive stress into the data, the method of Altschul et al. [Altschul, S. F., Carroll, R. J. & Lipman, D. J. (1989) J. Mol. Biol. 207, 647-653] is appropriate. When the sequences are not known to be phylogenetically related or a tree cannot be produced without unduly distorting the distances between the sequences, a modification of the method of Sibbald and Argos [Sibbald, P. R. & Argos, P. (1990) J. Mol. Biol. 216, 813-818] is preferable.

PALM--a pattern language for molecular biology.

This paper presents a new pattern language, PALM, for describing patterns in molecular biology sequences. The language is intended for representing knowledge about such patterns in a declarative, clear and concise way. It is also shown that its expressive power enables the definition of any regular or context free language, and also higher languages in the Chomsky hierarchy by parameter attachment, variables and procedural attachment. It is also possible to define approximate patterns. The language is rigorously defined, and several examples of its use and expressive power are given.

Overseer: a nucleotide sequence searching tool.

Overseer is a computer program that searches databases of nucleic acid sequences for objects of interest to the user. Such objects may consist of any number of simpler building blocks such as repeats, palindromes or stem-loops, strings of particular bases with or without mismatches, etc. Written in standard Pascal, this program runs under Unix and VMS and should also run under other operating systems. A simple interface allows the user to generate interactively a file containing a description of the target to be found. The searching program runs non-interactively, processing the information from the file and searching the sequences. The results are output to a file. Search capabilities are quite flexible and the code is designed to be modified. Since the framework of the program is simple, adding new modules to search for new target types as the need arises is possible.

Identification of proteins in sequence databases from amino acid composition data.

Having obtained the amino acid composition of a protein, chemists and molecular biologists may wish to identify the protein from this data alone. In general such data will have errors associated with them and the length of the protein may be known only approximately or not at all. In this paper a method is described which enables searching of protein sequence databases for sequences or fragments of sequences which have a composition similar to the one being sought. Such searches are generally quite discriminating as shown by the examples provided. This method has been implemented as part of the computer program Scrutineer and is being freely distributed. It is simple to use.

Automated protein sequence pattern handling and PROSITE searching.

The protein sequence searching program Scrutineer has been modified to search for targets from a file. We are distributing a reformatted file of PROSITES which can be read by Scrutineer. In addition, Scrutineer still accepts targets typed in interactively but can now write them out in the format required as input. Since the input format is the same as the output format, target management and re-use is simple.

Weighting aligned protein or nucleic acid sequences to correct for unequal representation.

Aligned sequences from the same family (e.g. the haemoglobins) are seldom representative of the entire family. This is because (1) the sequence databases are heavily skewed toward a small number of organisms and (2) only a minute fraction of all the different family members have been sequenced. For many applications, such as using alignments or profiles to perform database searches for distantly related family members, such unequal representation requires correction. An algorithm to perform appropriate weighting of individual sequences is presented along with examples illustrating its efficacy.

The P-loop--a common motif in ATP- and GTP-binding proteins.

Many ATP- and GTP-binding proteins have a phosphate-binding loop (P-loop), the primary structure of which typically consists of a glycine-rich sequence followed by a conserved lysine and a serine or threonine. The three-dimensional structures of several ATP- and GTP-binding proteins containing P-loops have now been solved. In this review current knowledge of P-loops is discussed with the additional aim of illustrating the fascinating relationship between protein sequence, structure and function.

Trans-splicing of pre-mRNA is predicted to occur in a wide range of organisms including vertebrates.

Several known trans-splicing RNA structures were used to define a canonical trans-splicing structure which was then used to perform a computer search of the EMBL nucleotide database. In addition to most known trans-splicing structures, many putative new trans-splicing sites were detected. These were found in a broad range of organisms including the vertebrates. Control experiments indicate that the search predicts known false positives at a rate of only 20%. Trans-splicing may therefore be a very wide-spread phenomenon.

Scrutineer: a computer program that flexibly seeks and describes motifs and profiles in protein sequence databases.

Scrutineer is an interactive, user-friendly program designed to search for motifs, patterns and profiles in the Swissprot, Protein Identification Resource (PIR) or SeqDb protein sequence databases. Basic capabilities include (i) searches for strings of amino acids with multiple choices at a given position; (ii) searches for strings including variable-length segments and delocalized constraints; (iii) searches over subsets of a database or particular regions within each sequence (e.g. N-terminal one-third); (iv) searches involving secondary structure predictions, physicochemical characteristics, and the like; and (v) searches using aligned sequences as targets with various optional weighting schemes. The various search criteria and hits can be combined and complex targets located. Once the data are loaded into virtual memory, all occurrences in PIR release 22.0 (3.7 x 10(6) amino acids) of a given short string of amino acids (e.g. a hexamer) are found in approximately 36 s. Scrutineer can also describe the entire database, user-specified hits, user-defined regions of sequence and all hits. The source code and accompanying manual are being freely distributed.

Photosystem I-Mediated Regulation of Water Splitting in the Red Alga, Porphyra sanjuanensis.

The marine red alga, Porphyra sanjuanensis is found mainly in the high intertidal zone and at low tide subject to frequent and extreme water stress, often accompanied by high temperatures and light intensities. Such exposures can lead to severe desiccation which is accompanied by the progressive loss of photosynthetic activity. Even following the loss of more than 90% of the thallus water content the alga recovers rapidly when returned to seawater. This stress-induced, reversible inactivation of photosynthesis is believed to be a protective adaptation which prevents photodamage to the exposed alga. Effects of light, inhibitors of water splitting, and electron donors to PSI on variable fluorescence and water splitting suggest that activity of the oxygen evolving complex is regulated by the PSI-driven reduction of a component of intersystem electron transport.

How probable are antibody cross-reactions?

Antibodies to polypeptides are increasingly being used in research. Their specificity and tight but reversible binding make them ideal for applications such as identification of proteins, immunological quantification or purification, and peptide mapping. Antibodies are also used in medicine to deliver loads to specific sites in tissues, and in electron microscopy as heavy metal conjugates to locate antigens in thin sections. While these techniques depend on specificity of antibody binding, it is occasionally observed that cross-reactions occur. These cross-reactions are attributed to the existence of one or more antibody binding sites common to both polypeptides. It is important to know whether these cross-reactions are expected due to chance alone, or if they are improbable and likely due to some causative agent. Examples of such causative agents might include gene duplication events or convergence due to functional constraints. At the present time, good methods for predicting the probability and therefore the frequency of cross-reactions are unavailable. In this paper we apply some recently reported mathematical results to address the following questions: (1) What is the probability that polyclonal or monoclonal antibodies raised against a given polypeptide will cross-react with another polypeptide due to chance alone? (2) What is the probability that polyclonal or monoclonal antibodies raised against a given polypeptide will cross-react with one or more polypeptides in a pool of polypeptides? Approximate answers to these questions are presented for cases where amino acid compositions of linear polypeptides are known or unknown, but the amino acid sequence of one or more of the polypeptides is not known. Implications of the results for antibody use in protein research are discussed.

Copper in photosystem II: association with LHC II.

Photosystem II particles from spinach and barley contained 2.5 and 4.2 Cu per 300 chlorophylls respectively. This Cu was resistant to removal by EDTA. A large percentage of the PSII Cu in both plants is associated with the light-harvesting chlorophyll a/b protein, LHCII; 46% in barley and 76% in spinach. Several experiments have been performed to rule out the possibility that the Cu was introduced during the isolation procedures and to ensure that the Cu is associated with PSII. Since the PSII Cu is mainly associated with LHCII, it is unlikely that it is involved in O2 evolution.