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.

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.

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.

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.

Analysis of the role of the Spitzenkörper in fungal morphogenesis by computer simulation of apical branching in Aspergillus niger

High-resolution video microscopy, image analysis, and computer simulation were used to study the role of the Spitzenkörper (Spk) in apical branching of ramosa-1, a temperature-sensitive mutant of Aspergillus niger. A shift to the restrictive temperature led to a cytoplasmic contraction that destabilized the Spk, causing its disappearance. After a short transition period, new Spk appeared where the two incipient apical branches emerged. Changes in cell shape, growth rate, and Spk position were recorded and transferred to the fungus simulator program to test the hypothesis that the Spk functions as a vesicle supply center (VSC). The simulation faithfully duplicated the elongation of the main hypha and the two apical branches. Elongating hyphae exhibited the growth pattern described by the hyphoid equation. During the transition phase, when no Spk was visible, the growth pattern was nonhyphoid, with consecutive periods of isometric and asymmetric expansion; the apex became enlarged and blunt before the apical branches emerged. Video microscopy images suggested that the branch Spk were formed anew by gradual condensation of vesicle clouds. Simulation exercises where the VSC was split into two new VSCs failed to produce realistic shapes, thus supporting the notion that the branch Spk did not originate by division of the original Spk. The best computer simulation of apical branching morphogenesis included simulations of the ontogeny of branch Spk via condensation of vesicle clouds. This study supports the hypothesis that the Spk plays a major role in hyphal morphogenesis by operating as a VSC—i.e., by regulating the traffic of wall-building vesicles in the manner predicted by the hyphoid model.

A two-year experience teaching computer literacy to first-year medical students using skill-based cohorts

Because it is widely accepted that providing information online will play a major role in both the teaching and practice of medicine in the near future, a short formal course of instruction in computer skills was proposed for the incoming class of students entering medical school at the State University of New York at Stony Brook. The syllabus was developed on the basis of a set of expected outcomes, which was accepted by the dean of medicine and the curriculum committee for classes beginning in the fall of 1997. Prior to their arrival, students were asked to complete a self-assessment survey designed to elucidate their initial skill base; the returned surveys showed students to have computer skills ranging from complete novice to that of a systems engineer. The classes were taught during the first three weeks of the semester to groups of students separated on the basis of their knowledge of and comfort with computers. Areas covered included computer basics, e-mail management, MEDLINE, and Internet search tools. Each student received seven hours of hands-on training followed by a test. The syllabus and emphasis of the classes were tailored to the initial skill base but the final test was given at the same level to all students. Student participation, test scores, and course evaluations indicated that this noncredit program was successful in achieving an acceptable level of comfort in using a computer for almost all of the student body.