The estimation of statistical parameters for local alignment score distributions

The distribution of optimal local alignment scores of random sequences plays a vital role in evaluating the statistical significance of sequence alignments. These scores can be well described by an extreme-value distribution. The distribution’s parameters depend upon the scoring system employed and the random letter frequencies; in general they cannot be derived analytically, but must be estimated by curve fitting. For obtaining accurate parameter estimates, a form of the recently described ‘island’ method has several advantages. We describe this method in detail, and use it to investigate the functional dependence of these parameters on finite-length edge effects.

BAliBASE (Benchmark Alignment dataBASE): enhancements for repeats, transmembrane sequences and circular permutations

BAliBASE is specifically designed to serve as an evaluation resource to address all the problems encountered when aligning complete sequences. The database contains high quality, manually constructed multiple sequence alignments together with detailed annotations. The alignments are all based on three-dimensional structural superpositions, with the exception of the transmembrane sequences. The first release provided sets of reference alignments dealing with the problems of high variability, unequal repartition and large N/C-terminal extensions and internal insertions. Here we describe version 2.0 of the database, which incorporates three new reference sets of alignments containing structural repeats, trans­membrane sequences and circular permutations to evaluate the accuracy of detection/prediction and alignment of these complex sequences. BAliBASE can be viewed at the web site http://www-igbmc.u-strasbg.fr/BioInfo/BAliBASE2/index.html or can be downloaded from ftp://ftp-igbmc.u-strasbg.fr/pub/BAliBASE2/.

A database and tools for 3-D protein structure comparison and alignment using the Combinatorial Extension (CE) algorithm

The database reported here is derived using the Combinatorial Extension (CE) algorithm which compares pairs of protein polypeptide chains and provides a list of structurally similar proteins along with their structure alignments. Using CE, structure–structure alignments can provide insights into biological function. When a protein of known function is shown to be structurally similar to a protein of unknown function, a relationship might be inferred; a relationship not necessarily detectable from sequence comparison alone. Establishing structure–structure relationships in this way is of great importance as we enter an era of structural genomics where there is a likelihood of an increasing number of structures with unknown functions being determined. Thus the CE database is an example of a useful tool in the annotation of protein structures of unknown function. Comparisons can be performed on the complete PDB or on a structurally representative subset of proteins. The source protein(s) can be from the PDB (updated monthly) or uploaded by the user. CE provides sequence alignments resulting from structural alignments and Cartesian coordinates for the aligned structures, which may be analyzed using the supplied Compare3D Java applet, or downloaded for further local analysis. Searches can be run from the CE web site, http://cl.sdsc.edu/ce.html, or the database and software downloaded from the site for local use.

PALI—a database of Phylogeny and ALIgnment of homologous protein structures

PALI (release 1.2) contains three-dimensional (3-D) structure-dependent sequence alignments as well as structure-based phylogenetic trees of homologous protein domains in various families. The data set of homologous protein structures has been derived by consulting the SCOP database (release 1.50) and the data set comprises 604 families of homologous proteins involving 2739 protein domain structures with each family made up of at least two members. Each member in a family has been structurally aligned with every other member in the same family (pairwise alignment) and all the members in the family are also aligned using simultaneous super­position (multiple alignment). The structural alignments are performed largely automatically, with manual interventions especially in the cases of distantly related proteins, using the program STAMP (version 4.2). Every family is also associated with two dendrograms, calculated using PHYLIP (version 3.5), one based on a structural dissimilarity metric defined for every pairwise alignment and the other based on similarity of topologically equivalent residues. These dendrograms enable easy comparison of sequence and structure-based relationships among the members in a family. Structure-based alignments with the details of structural and sequence similarities, superposed coordinate sets and dendrograms can be accessed conveniently using a web interface. The database can be queried for protein pairs with sequence or structural similarities falling within a specified range. Thus PALI forms a useful resource to help in analysing the relationship between sequence and structure variation at a given level of sequence similarity. PALI also contains over 653 ‘orphans’ (single member families). Using the web interface involving PSI_BLAST and PHYLIP it is possible to associate the sequence of a new protein with one of the families in PALI and generate a phylogenetic tree combining the query sequence and proteins of known 3-D structure. The database with the web interfaced search and dendrogram generation tools can be accessed at http://pa uling.mbu.iisc.ernet.in/~pali.

STACK: Sequence Tag Alignment and Consensus Knowledgebase

STACK is a tool for detection and visualisation of expressed transcript variation in the context of developmental and pathological states. The datasystem organises and reconstructs human transcripts from available public data in the context of expression state. The expression state of a transcript can include developmental state, pathological association, site of expression and isoform of expressed transcript. STACK consensus transcripts are reconstructed from clusters that capture and reflect the growing evidence of transcript diversity. The comprehensive capture of transcript variants is achieved by the use of a novel clustering approach that is tolerant of sub-sequence diversity and does not rely on pairwise alignment. This is in contrast with other gene indexing projects. STACK is generated at least four times a year and represents the exhaustive processing of all publicly available human EST data extracted from GenBank. This processed information can be explored through 15 tissue-specific categories, a disease-related category and a whole-body index and is accessible via WWW at http://www.sanbi.ac.za/Dbases.html. STACK represents a broadly applicable resource, as it is the only reconstructed transcript database for which the tools for its generation are also broadly available (http://www.sanbi.ac.za/CODES).

Active site mutant transgene confers tolerance to human β-glucuronidase without affecting the phenotype of MPS VII mice

Mucopolysaccharidosis type VII (MPS VII; Sly syndrome) is an autosomal recessive lysosomal storage disorder due to an inherited deficiency of β-glucuronidase. A naturally occurring mouse model for this disease was discovered at The Jackson Laboratory and shown to be due to homozygosity for a 1-bp deletion in exon 10 of the gus gene. The murine model MPS VII (gusmps/mps) has been very well characterized and used extensively to evaluate experimental strategies for lysosomal storage diseases, including bone marrow transplantation, enzyme replacement therapy, and gene therapy. To enhance the value of this model for enzyme and gene therapy, we produced a transgenic mouse expressing the human β-glucuronidase cDNA with an amino acid substitution at the active site nucleophile (E540A) and bred it onto the MPS VII (gusmps/mps) background. We demonstrate here that the mutant mice bearing the active site mutant human transgene retain the clinical, morphological, biochemical, and histopathological characteristics of the original MPS VII (gusmps/mps) mouse. However, they are now tolerant to immune challenge with human β-glucuronidase. This “tolerant MPS VII mouse model” should be useful for preclinical trials evaluating the effectiveness of enzyme and/or gene therapy with the human gene products likely to be administered to human patients with MPS VII.

Requirement for natural killer T (NKT) cells in the induction of allograft tolerance

In this study, we investigated the role of Vα14 natural killer T (NKT) cells in transplant immunity. The ability to reject allografts was not significantly different between wild-type (WT) and Vα14 NKT cell-deficient mice. However, in models in which tolerance was induced against cardiac allografts by blockade of lymphocyte function-associated antigen-1/intercellular adhesion molecule-1 or CD28/B7 interactions, long-term acceptance of the grafts was observed only in WT but not Vα14 NKT cell-deficient mice. Adoptive transfer with Vα14 NKT cells restored long-term acceptance of allografts in Vα14 NKT cell-deficient mice. The critical role of Vα14 NKT cells to mediate immunosuppression was also observed in vitro in mixed lymphocyte cultures in which lymphocyte function-associated antigen-1/intercellular adhesion molecule-1 or CD28/B7 interactions were blocked. Experiments using IL-4- or IFN-γ-deficient mice suggested a critical contribution of IFN-γ to the Vα14 NKT cell-mediated allograft acceptance in vivo. These results indicate a critical contribution of Vα14 NKT cells to the induction of allograft tolerance and provide a useful model to investigate the regulatory role of Vα14 NKT cells in various immune responses.

Peyer's patches are required for oral tolerance to proteins

To clarify the role of Peyer's patches in oral tolerance induction, BALB/c mice were treated in utero with lymphotoxin β-receptor Ig fusion protein to generate mice lacking Peyer's patches. When these Peyer's patch-null mice were fed 25 mg of ovalbumin (OVA) before systemic immunization, OVA-specific IgG Ab responses in serum and spleen were seen, in marked contrast to low responses in OVA-fed normal mice. Further, high T-cell-proliferative- and delayed-type hypersensitivity responses were seen in Peyer's patch-null mice given oral OVA before systemic challenge. Higher levels of CD4+ T-cell-derived IFN-γ, IL-4, IL-5, and IL-10 syntheses were noted in Peyer's patch-null mice fed OVA, whereas OVA-fed normal mice had suppressed cytokine levels. In contrast, oral administration of trinitrobenzene sulfonic acid (TNBS) to Peyer's patch-null mice resulted in reduced TNBS-specific serum Abs and splenic B cell antitrinitrophenyl Ab-forming cell responses after skin painting with picryl chloride. Further, when delayed-type hypersensitivity and splenic T cell proliferative responses were examined, Peyer's patch-null mice fed TNBS were unresponsive to hapten. Peyer's patch-null mice fed trinitrophenyl-OVA failed to induce systemic unresponsiveness to hapten or protein. These findings show that organized Peyer's patches are required for oral tolerance to proteins, whereas haptens elicit systemic unresponsiveness via the intestinal epithelial cell barrier.

Glutamic acid decarboxylase and glutamate receptor changes during tolerance and dependence to benzodiazepines

Protracted administration of diazepam elicits tolerance, whereas discontinuation of treatment results in signs of dependence. Tolerance to the anticonvulsant action of diazepam is present in an early phase (6, 24, and 36 h) but disappears in a late phase (72–96 h) of withdrawal. In contrast, signs of dependence such as decrease in open-arm entries on an elevated plus-maze and increased susceptibility to pentylenetetrazol-induced seizures were apparent 96 h (but not 12, 24, or 48 h) after diazepam withdrawal. During the first 72 h of withdrawal, tolerance is associated with changes in the expression of GABAA (γ-aminobutyric acid type A) receptor subunits (decrease in γ2 and α1; increase in α5) and with an increase of mRNA expression of the most abundant form of glutamic acid decarboxylase (GAD), GAD67. In contrast, dl-α-amino-3-hydroxy-5-methylisoxazole-4-propionic acid (AMPA) receptor GluR1 subunit mRNA and cognate protein, which are normal during the early phase of diazepam withdrawal, increase by approximately 30% in cortex and hippocampus in association with the appearance of signs of dependence 96 h after diazepam withdrawal. Immunohistochemical studies of GluR1 subunit expression with gold-immunolabeling technique reveal that the increase of GluR1 subunit protein is localized to layer V pyramidal neurons and their apical dendrites in the cortex, and to pyramidal neurons and in their dendritic fields in hippocampus. The results suggest an involvement of GABA-mediated processes in the development and maintenance of tolerance to diazepam, whereas excitatory amino acid-related processes (presumably via AMPA receptors) may be involved in the expression of signs of dependence after withdrawal.

ParAlign: a parallel sequence alignment algorithm for rapid and sensitive database searches

There is a need for faster and more sensitive algorithms for sequence similarity searching in view of the rapidly increasing amounts of genomic sequence data available. Parallel processing capabilities in the form of the single instruction, multiple data (SIMD) technology are now available in common microprocessors and enable a single microprocessor to perform many operations in parallel. The ParAlign algorithm has been specifically designed to take advantage of this technology. The new algorithm initially exploits parallelism to perform a very rapid computation of the exact optimal ungapped alignment score for all diagonals in the alignment matrix. Then, a novel heuristic is employed to compute an approximate score of a gapped alignment by combining the scores of several diagonals. This approximate score is used to select the most interesting database sequences for a subsequent Smith–Waterman alignment, which is also parallelised. The resulting method represents a substantial improvement compared to existing heuristics. The sensitivity and specificity of ParAlign was found to be as good as Smith–Waterman implementations when the same method for computing the statistical significance of the matches was used. In terms of speed, only the significantly less sensitive NCBI BLAST 2 program was found to outperform the new approach. Online searches are available at http://dna.uio.no/search/

A Selaginella lepidophylla Trehalose-6-Phosphate Synthase Complements Growth and Stress-Tolerance Defects in a Yeast tps1 Mutant1

The accumulation of the disaccharide trehalose in anhydrobiotic organisms allows them to survive severe environmental stress. A plant cDNA, SlTPS1, encoding a 109-kD protein, was isolated from the resurrection plant Selaginella lepidophylla, which accumulates high levels of trehalose. Protein-sequence comparison showed that SlTPS1 shares high similarity to trehalose-6-phosphate synthase genes from prokaryotes and eukaryotes. SlTPS1 mRNA was constitutively expressed in S. lepidophylla. DNA gel-blot analysis indicated that SlTPS1 is present as a single-copy gene. Transformation of a Saccharomyces cerevisiae tps1Δ mutant disrupted in the ScTPS1 gene with S. lepidophylla SlTPS1 restored growth on fermentable sugars and the synthesis of trehalose at high levels. Moreover, the SlTPS1 gene introduced into the tps1Δ mutant was able to complement both deficiencies: sensitivity to sublethal heat treatment at 39°C and induced thermotolerance at 50°C. The osmosensitive phenotype of the yeast tps1Δ mutant grown in NaCl and sorbitol was also restored by the SlTPS1 gene. Thus, SlTPS1 protein is a functional plant homolog capable of sustaining trehalose biosynthesis and could play a major role in stress tolerance in S. lepidophylla.

Transgenically Enhanced Sorbitol Synthesis Facilitates Phloem Boron Transport and Increases Tolerance of Tobacco to Boron Deficiency1

The mobility of elements within plants contributes to a plant species' tolerance of nutrient deficiencies in the soil. The genetic manipulation of within-plant nutrient movement may therefore provide a means to enhance plant growth under conditions of variable soil nutrient availability. In these experiments tobacco (Nicotiana tabacum) was engineered to synthesize sorbitol, and the resultant effect on phloem mobility of boron (B) was determined. In contrast to wild-type tobacco, transgenic tobacco plants containing sorbitol exhibit a marked increase in within-plant B mobility and a resultant increase in plant growth and yield when grown with limited or interrupted soil B supply. Growth of transgenic tobacco could be maintained by reutilization of B present in mature tissues or from B supplied as a foliar application to mature leaves. In contrast, B present in mature leaves of control tobacco lines could not be used to provide the B requirements for new plant growth. 10B-labeling experiments verified that B is phloem mobile in transgenic tobacco but is immobile in control lines. These results demonstrate that the transgenic enhancement of within-plant nutrient mobility is a viable approach to improve plant tolerance of nutrient stress.

Arabidopsis Roots and Shoots Have Different Mechanisms for Hypoxic Stress Tolerance

Arabidopsis has inducible responses for tolerance of O2 deficiency. Plants previously exposed to 5% O2 were more tolerant than the controls to hypoxic stress (0.1% O2 for 48 h) in both roots and shoots, but hypoxic acclimation did not improve tolerance to anoxia (0% O2). The acclimation of shoots was not dependent on the roots: increased shoot tolerance was observed when the roots of the plants were removed. An adh (alcohol dehydrogenase) null mutant did not show acclimation of the roots but retained the shoot survival response. Abscisic acid treatment also differentiated the root and shoot responses; pretreatment induced root survival in hypoxic stress conditions (0.1% O2) but did not induce any increase in the survival of shoots. Cycloheximide blocked both root and shoot acclimation, indicating that both acclimation mechanisms are dependent on protein synthesis.

Cellular mechanisms of neuropathic pain, morphine tolerance, and their interactions

Compelling evidence has accumulated over the last several years from our laboratory, as well as others, indicating that central hyperactive states resulting from neuronal plastic changes within the spinal cord play a critical role in hyperalgesia associated with nerve injury and inflammation. In our laboratory, chronic constriction injury of the common sciatic nerve, a rat model of neuropathic pain, has been shown to result in activation of central nervous system excitatory amino acid receptors and subsequent intracellular cascades including protein kinase C translocation and activation, nitric oxide production, and nitric oxide-activated poly(ADP ribose) synthetase activation. Similar cellular mechanisms also have been implicated in the development of tolerance to the analgesic effects of morphine. A recently observed phenomenon, the development of “dark neurons,” is associated with both chronic constriction injury and morphine tolerance. A site of action involved in both hyperalgesia and morphine tolerance is in the superficial laminae of the spinal cord dorsal horn. These observations suggest that hyperalgesia and morphine tolerance may be interrelated at the level of the superficial laminae of the dorsal horn by common neural substrates that interact at the level of excitatory amino acid receptor activation and subsequent intracellular events. The demonstration of interrelationships between neural mechanisms underlying hyperalgesia and morphine tolerance may lead to a better understanding of the neurobiology of these two phenomena in particular and pain in general. This knowledge may also provide a scientific basis for improved pain management with opiate analgesics.