Human Immunodeficiency Virus Reverse Transcriptase and Protease Sequence Database: an expanded data model integrating natural language text and sequence analysis programs

The HIV Reverse Transcriptase and Protease Sequence Database is an on-line relational database that catalogs evolutionary and drug-related sequence variation in the human immunodeficiency virus (HIV) reverse transcriptase (RT) and protease enzymes, the molecular targets of anti-HIV therapy (http://hivdb.stanford.edu). The database contains a compilation of nearly all published HIV RT and protease sequences, including submissions from International Collaboration databases and sequences published in journal articles. Sequences are linked to data about the source of the sequence sample and the antiretroviral drug treatment history of the individual from whom the isolate was obtained. During the past year 3500 sequences have been added and the data model has been expanded to include drug susceptibility data on sequenced isolates. Database content has also been integrated with didactic text and the output of two sequence analysis programs.

The evolution of coloration and toxicity in the poison frog family (Dendrobatidae)

The poison frogs (family Dendrobatidae) are terrestrial anuran amphibians displaying a wide range of coloration and toxicity. These frogs generally have been considered to be aposematic, but relatively little research has been carried out to test the predictions of this hypothesis. Here we use a comparative approach to test one prediction of the hypothesis of aposematism: that coloration will evolve in tandem with toxicity. Recently, we developed a phylogenetic hypothesis of the evolutionary relationships among representative species of poison frogs, using sequences from three regions of mitochondrial DNA. In our analysis, we use that DNA-based phylogeny and comparative analysis of independent contrasts to investigate the correlation between coloration and toxicity in the poison frog family (Dendrobatidae). Information on the toxicity of different species was obtained from the literature. Two different measures of the brightness and extent of coloration were used. (i) Twenty-four human observers were asked to rank different photos of each different species in the analysis in terms of contrast to a leaf-littered background. (ii) Color photos of each species were scanned into a computer and a computer program was used to obtain a measure of the contrast of the colors of each species relative to a leaf-littered background. Comparative analyses of the results were carried out with two different models of character evolution: gradual change, with branch lengths proportional to the amount of genetic change, and punctuational change, with all change being associated with speciation events. Comparative analysis using either method or model indicated a significant correlation between the evolution of toxicity and coloration across this family. These results are consistent with the hypothesis that coloration in this group is aposematic.

Computation of Surface Electrical Potentials of Plant Cell Membranes : Correspondence to Published Zeta Potentials from Diverse Plant Sources

A Gouy-Chapman-Stern model has been developed for the computation of surface electrical potential (ψ0) of plant cell membranes in response to ionic solutes. The present model is a modification of an earlier version developed to compute the sorption of ions by wheat (Triticum aestivum L. cv Scout 66) root plasma membranes. A single set of model parameters generates values for ψ0 that correlate highly with published ζ potentials of protoplasts and plasma membrane vesicles from diverse plant sources. The model assumes ion binding to a negatively charged site (R− = 0.3074 μmol m−2) and to a neutral site (P0 = 2.4 μmol m−2) according to the reactions R− + IΖ ⇌ RIΖ−1 and P0 + IΖ ⇌ PIΖ, where IΖ represents an ion of charge Ζ. Binding constants for the negative site are 21,500 m−1 for H+, 20,000 m−1 for Al3+, 2,200 m−1 for La3+, 30 m−1 for Ca2+ and Mg2+, and 1 m−1 for Na+ and K+. Binding constants for the neutral site are 1/180 the value for binding to the negative site. Ion activities at the membrane surface, computed on the basis of ψ0, appear to determine many aspects of plant-mineral interactions, including mineral nutrition and the induction and alleviation of mineral toxicities, according to previous and ongoing studies. A computer program with instructions for the computation of ψ0, ion binding, ion concentrations, and ion activities at membrane surfaces may be requested from the authors.

Selection of Gβ Subunits with Point Mutations That Fail to Activate Specific Signaling Pathways In Vivo: Dissecting Cellular Responses Mediated by a Heterotrimeric G Protein in Dictyostelium discoideum

In Dictyostelium discoideum, a unique Gβ subunit is required for a G protein–coupled receptor system that mediates a variety of cellular responses. Binding of cAMP to cAR1, the receptor linked to the G protein G2, triggers a cascade of responses, including activation of adenylyl cyclase, gene induction, actin polymerization, and chemotaxis. Null mutations of the cAR1, Gα2, and Gβ genes completely impair all these responses. To dissect specificity in Gβγ signaling to downstream effectors in living cells, we screened a randomly mutagenized library of Gβ genes and isolated Gβ alleles that lacked the capacity to activate some effectors but retained the ability to regulate others. These mutant Gβ subunits were able to link cAR1 to G2, to support gene expression, and to mediate cAMP-induced actin polymerization, and some were able to mediate to chemotaxis toward cAMP. None was able to activate adenylyl cyclase, and some did not support chemotaxis. Thus, we separated in vivo functions of Gβγ by making point mutations on Gβ. Using the structure of the heterotrimeric G protein displayed in the computer program CHAIN, we examined the positions and the molecular interactions of the amino acids substituted in each of the mutant Gβs and analyzed the possible effects of each replacement. We identified several residues that are crucial for activation of the adenylyl cyclase. These residues formed an area that overlaps but is not identical to regions where bovine Gtβγ interacts with its regulators, Gα and phosducin.

Intron positions correlate with module boundaries in ancient proteins

We analyze the three-dimensional structure of proteins by a computer program that finds regions of sequence that contain module boundaries, defining a module as a segment of polypeptide chain bounded in space by a specific given distance. The program defines a set of “linker regions” that have the property that if an intron were to be placed into each linker region, the protein would be dissected into a set of modules all less than the specified diameter. We test a set of 32 proteins, all of ancient origin, and a corresponding set of 570 intron positions, to ask if there is a statistically significant excess of intron positions within the linker regions. For 28-Å modules, a standard size used historically, we find such an excess, with P < 0.003. This correlation is neither due to a compositional or sequence bias in the linker regions nor to a surface bias in intron positions. Furthermore, a subset of 20 introns, which can be putatively identified as old, lies even more explicitly within the linker regions, with P < 0.0003. Thus, there is a strong correlation between intron positions and three-dimensional structural elements of ancient proteins as expected by the introns-early approach. We then study a range of module diameters and show that, as the diameter varies, significant peaks of correlation appear for module diameters centered at 21.7, 27.6, and 32.9 Å. These preferred module diameters roughly correspond to predicted exon sizes of 15, 22, and 30 residues. Thus, there are significant correlations between introns, modules, and a quantized pattern of the lengths of polypeptide chains, which is the prediction of the “Exon Theory of Genes.”

GenEST, a powerful bidirectional link between cDNA sequence data and gene expression profiles generated by cDNA-AFLP

The release of vast quantities of DNA sequence data by large-scale genome and expressed sequence tag (EST) projects underlines the necessity for the development of efficient and inexpensive ways to link sequence databases with temporal and spatial expression profiles. Here we demonstrate the power of linking cDNA sequence data (including EST sequences) with transcript profiles revealed by cDNA-AFLP, a highly reproducible differential display method based on restriction enzyme digests and selective amplification under high stringency conditions. We have developed a computer program (GenEST) that predicts the sizes of virtual transcript-derived fragments (TDFs) of in silico-digested cDNA sequences retrieved from databases. The vast majority of the resulting virtual TDFs could be traced back among the thousands of TDFs displayed on cDNA-AFLP gels. Sequencing of the corresponding bands excised from cDNA-AFLP gels revealed no inconsistencies. As a consequence, cDNA sequence databases can be screened very efficiently to identify genes with relevant expression profiles. The other way round, it is possible to switch from cDNA-AFLP gels to sequences in the databases. Using the restriction enzyme recognition sites, the primer extensions and the estimated TDF size as identifiers, the DNA sequence(s) corresponding to a TDF with an interesting expression pattern can be identified. In this paper we show examples in both directions by analyzing the plant parasitic nematode Globodera rostochiensis. Various novel pathogenicity factors were identified by combining ESTs from the infective stage juveniles with expression profiles of ∼4000 genes in five developmental stages produced by cDNA-AFLP.

Highly specific protein sequence motifs for genome analysis

We present a method for discovering conserved sequence motifs from families of aligned protein sequences. The method has been implemented as a computer program called emotif (http://motif.stanford.edu/emotif). Given an aligned set of protein sequences, emotif generates a set of motifs with a wide range of specificities and sensitivities. emotif also can generate motifs that describe possible subfamilies of a protein superfamily. A disjunction of such motifs often can represent the entire superfamily with high specificity and sensitivity. We have used emotif to generate sets of motifs from all 7,000 protein alignments in the blocks and prints databases. The resulting database, called identify (http://motif.stanford.edu/identify), contains more than 50,000 motifs. For each alignment, the database contains several motifs having a probability of matching a false positive that range from 10−10 to 10−5. Highly specific motifs are well suited for searching entire proteomes, while generating very few false predictions. identify assigns biological functions to 25–30% of all proteins encoded by the Saccharomyces cerevisiae genome and by several bacterial genomes. In particular, identify assigned functions to 172 of proteins of unknown function in the yeast genome.