The first genome from the third domain of life

Protocells and LUCA: Transport of substances from first physicochemical principles

Models of the transport of substances in protocells are considered from first physicochemical principles. Functional similarities and differences in the transport systems of archaea, cyanobacteria, E. coli, and diatoms have been analyzed. Based on the selection of the most important transport systems, a model of transport of substances through the membrane of the last universal common ancestor, LUCA, was constructed. Models of isotope separation in protocells were considered. Based on the proposed models, the difference in isotope concentrations in rocks can be predicted, which can serve as an indicator of the presence of life in the early stages of evolution. Mechanisms of energy conversion for the simplest forms of directed motion in protocells are considered. A special stage in the evolution of protocells is proposed - the minimal mobile cell.

Biocuration at the Saccharomyces genome database

Saccharomyces Genome Database is an online resource dedicated to managing information about the biology and genetics of the model organism, yeast (Saccharomyces cerevisiae). This information is derived primarily from scientific publications through a process of human curation that involves manual extraction of data and their organization into a comprehensive system of knowledge. This system provides a foundation for further analysis of experimental data coming from research on yeast as well as other organisms. In this review we will demonstrate how biocuration and biocurators add a key component, the biological context, to our understanding of how genes, proteins, genomes and cells function and interact. We will explain the role biocurators play in sifting through the wealth of biological data to incorporate and connect key information. We will also discuss the many ways we assist researchers with their various research needs. We hope to convince the reader that manual curation is vital in converting the flood of data into organized and interconnected knowledge, and that biocurators play an essential role in the integration of scientific information into a coherent model of the cell.

From microbial gene essentiality to novel antimicrobial drug targets

Background

Bacterial respiratory tract infections, mainly caused by Streptococcus pneumoniae, Haemophilus influenzae and Moraxella catarrhalis are among the leading causes of global mortality and morbidity. Increased resistance of these pathogens to existing antibiotics necessitates the search for novel targets to develop potent antimicrobials.

Result

Here, we report a proof of concept study for the reliable identification of potential drug targets in these human respiratory pathogens by combining high-density transposon mutagenesis, high-throughput sequencing, and integrative genomics. Approximately 20% of all genes in these three species were essential for growth and viability, including 128 essential and conserved genes, part of 47 metabolic pathways. By comparing these essential genes to the human genome, and a database of genes from commensal human gut microbiota, we identified and excluded potential drug targets in respiratory tract pathogens that will have off-target effects in the host, or disrupt the natural host microbiota. We propose 249 potential drug targets, 67 of which are targets for 75 FDA-approved antimicrobials and 35 other researched small molecule inhibitors. Two out of four selected novel targets were experimentally validated, proofing the concept.

Conclusion

Here we have pioneered an attempt in systematically combining the power of high-density transposon mutagenesis, high-throughput sequencing, and integrative genomics to discover potential drug targets at genome-scale. By circumventing the time-consuming and expensive laboratory screens traditionally used to select potential drug targets, our approach provides an attractive alternative that could accelerate the much needed discovery of novel antimicrobials.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-958) contains supplementary material, which is available to authorized users.

The Little Fly that Could: Wizardry and Artistry of Drosophila Genomics

For more than 100 years now, the fruit fly Drosophila melanogaster has been at the forefront of our endeavors to unlock the secrets of the genome. From the pioneering studies of chromosomes and heredity by Morgan and his colleagues, to the generation of fly models for human disease, Drosophila research has been at the forefront of genetics and genomics. We present a broad overview of some of the most powerful genomics tools that keep Drosophila research at the cutting edge of modern biomedical research.

The reference genome sequence of Saccharomyces cerevisiae: then and now

The genome of the budding yeast Saccharomyces cerevisiae was the first completely sequenced from a eukaryote. It was released in 1996 as the work of a worldwide effort of hundreds of researchers. In the time since, the yeast genome has been intensively studied by geneticists, molecular biologists, and computational scientists all over the world. Maintenance and annotation of the genome sequence have long been provided by the Saccharomyces Genome Database, one of the original model organism databases. To deepen our understanding of the eukaryotic genome, the S. cerevisiae strain S288C reference genome sequence was updated recently in its first major update since 1996. The new version, called "S288C 2010," was determined from a single yeast colony using modern sequencing technologies and serves as the anchor for further innovations in yeast genomic science.

A yeast-based recombination assay for homing endonuclease activity

Homing endonucleases (HEs) are natural enzymes that cleave long DNA target with a high specificity and trigger homologous recombination at the exact site of the break. Such mechanisms can thus be used for all the applications covered today by the generic name of "genome engineering": targeted sequence insertion, removal, or editing. However, before being able to address those applications, the engineering of HEs must be mastered so that any potential target would be efficiently and specifically recognized and cleaved. Working on the I-CreI model, we have developed a very powerful platform to generate HEs with new tailored specificity. We have put in place the first in vivo, functional, high throughput assay to generate I-CreI variants and measure their activity. We use semi-rational design combined with proprietary in silico predictions to design and synthesize I-CreI mutants that are tested for their capacity to induce homologous recombination in a yeast cell. The process has been standardized and robotized so that we can generate thousands of I-CreI derivatives, characterize their cleavage profile, and deliver them for further applications in the research, therapeutic, or agrobusiness fields.

Potassium and sodium transport in non-animal cells: the Trk/Ktr/HKT transporter family

Bacterial Trk and Ktr, fungal Trk and plant HKT form a family of membrane transporters permeable to K(+) and/or Na(+) and characterized by a common structure probably derived from an ancestral K(+) channel subunit. This transporter family, specific of non-animal cells, displays a large diversity in terms of ionic permeability, affinity and energetic coupling (H(+)-K(+) or Na(+)-K(+) symport, K(+) or Na(+) uniport), which might reflect a high need for adaptation in organisms living in fluctuating or dilute environments. Trk/Ktr/HKT transporters are involved in diverse functions, from K(+) or Na(+) uptake to membrane potential control, adaptation to osmotic or salt stress, or Na(+) recirculation from shoots to roots in plants. Structural analyses of bacterial Ktr point to multimeric structures physically interacting with regulatory subunits. Elucidation of Trk/Ktr/HKT protein structures along with characterization of mutated transporters could highlight functional and evolutionary relationships between ion channels and transporters displaying channel-like features.

Membrane transporters and drug development: relevance to pharmacogenomics, nutrigenomics, epigenetics, and systems biology

The study of membrane transporters may result in breakthroughs in the discovery of new drugs and the development of safer drugs. Membrane transporters are essential for fundamental cellular functions and normal physiological processes. These molecules influence drug absorption and distribution and play key roles in drug therapeutic effects. A primary goal of current research in drug discovery and development is to fully understand the interactions between transporters and drugs at both the system levels in the human body and the individual level for personalized therapy. Systematic studies of membrane transporters will help in not only better understanding of diseases from the systems biology point of view but also better drug design and development. The exploration of both pharmacogenomics and systems biology in transporters is necessary to connect individuals' genetic profiles with systematic drug responses in the human body. Understanding of gene-diet interactions and the effects of epigenetic changes on transporter gene expression may help improve clinical drug efficacy. The integration of pharmacogenomics, nutrigenomics, epigenetics, and systems biology may enable us to move from disease treatment to disease prevention and optimal health. The key issues in such integrative understanding include the correlations between structure and function, genotype and phenotype, and systematic interactions among transporters, other proteins, nutrients, drugs, and the environment. The exploration in these key issues may ultimately contribute to personalized medicine with high efficacy but less toxicity.

Overview of subcellular fractionation procedures for the yeast Saccharomyces cerevisiae

This unit provides an overview of centrifugation-based fractionation procedures adapted for the yeast Saccharomyces cerevisiae. The goals, merits, limitations, and critical parameters of are discussed in order to facilitate the development of subcellular fractionation strategies. Topics include yeast cell lysate preparation, differential velocity centrifugation, density gradient centrifugation, and the analysis of subcellular fractions.

The Very Early Stages of Biological Evolution and the Nature of the Last Common Ancestor of the Three Major Cell Domains

Quantitative estimates of the gene complement of the last common ancestor of all extant organisms, that is, the cenancestor, may be hindered by ancient horizontal gene transfer events and polyphyletic gene losses, as well as by biases in genome databases and methodological artifacts. Nevertheless, most reports agree that the last common ancestor resembled extant prokaryotes. A significant number of the highly conserved genes are sequences involved in the synthesis, degradation, and binding of RNA, including transcription and translation. Although the gene complement of the cenancestor includes sequences that may have originated in different epochs, the extraordinary conservation of RNA-related sequences supports the hypothesis that the last common ancestor was an evolutionary outcome of the so-called RNA/protein world. The available evidence suggests that the cenancestor was not a hyperthermophile, but it is currently not possible to assess its ecological niche or its mode of energy acquisition and carbon sources.

Question 7: comparative genomics and early cell evolution: a cautionary methodological note

Inventories of the gene content of the last common ancestor (LCA), i.e., the cenancestor, include sequences that may have undergone horizontal transfer events, as well as sequences that have originated in different pre-cenancestral epochs. However, the universal distribution of highly conserved genes involved in RNA metabolism provide insights into early stages of cell evolution during which RNA played a much more conspicuous biological role, and is consistent with the hypothesis that extant living systems were preceded by an RNA/protein world. Insights into the traits of primitive entities from which the LCA evolved may be derived from the analysis of paralogous gene families, including those formed by sequences that resulted from internal elongation events. Three major types of paralogous gene families can be recognized. The importance of this grouping for understanding the traits of early cells is discussed.

The calpain system

The calpain system originally comprised three molecules: two Ca2+-dependent proteases, mu-calpain and m-calpain, and a third polypeptide, calpastatin, whose only known function is to inhibit the two calpains. Both mu- and m-calpain are heterodimers containing an identical 28-kDa subunit and an 80-kDa subunit that shares 55-65% sequence homology between the two proteases. The crystallographic structure of m-calpain reveals six "domains" in the 80-kDa subunit: 1). a 19-amino acid NH2-terminal sequence; 2). and 3). two domains that constitute the active site, IIa and IIb; 4). domain III; 5). an 18-amino acid extended sequence linking domain III to domain IV; and 6). domain IV, which resembles the penta EF-hand family of polypeptides. The single calpastatin gene can produce eight or more calpastatin polypeptides ranging from 17 to 85 kDa by use of different promoters and alternative splicing events. The physiological significance of these different calpastatins is unclear, although all bind to three different places on the calpain molecule; binding to at least two of the sites is Ca2+ dependent. Since 1989, cDNA cloning has identified 12 additional mRNAs in mammals that encode polypeptides homologous to domains IIa and IIb of the 80-kDa subunit of mu- and m-calpain, and calpain-like mRNAs have been identified in other organisms. The molecules encoded by these mRNAs have not been isolated, so little is known about their properties. How calpain activity is regulated in cells is still unclear, but the calpains ostensibly participate in a variety of cellular processes including remodeling of cytoskeletal/membrane attachments, different signal transduction pathways, and apoptosis. Deregulated calpain activity following loss of Ca2+ homeostasis results in tissue damage in response to events such as myocardial infarcts, stroke, and brain trauma.

Rhesus monkey simian immunodeficiency virus infection as a model for assessing the role of selenium in AIDS

The objective of this study was to determine whether simian immunodeficiency virus (SIV) infection of macaques could be used as a model system to assess the role of selenium in AIDS. Plasma and serum selenium levels were determined by standard assays in monkeys before and after inoculation of SIV. SIV-infected cells or cells expressing the HIV Tat protein were labeled with 75Se, and protein extracts were prepared and electrophoresed to analyze selenoprotein expression. Total tRNA was isolated from CEMx174 cells infected with SIV or from KK1 cells infected with HIV, and selenocysteine tRNA isoforms were characterized by reverse phase chromatography. SIV-infected monkeys show a decrease in blood selenium levels similar to that observed in AIDS with development of SAIDS. Cells infected with SIV in vitro exhibit reduced selenoprotein levels and an accumulation of small molecular weight selenium compounds relative to uninfected cells. Examination of the selenocysteine tRNA isoforms in HIV-infected KK1 cells or SIV-infected CEMx174 cells reveals an isoform distribution characteristic of selenium-deficient cells. Furthermore, transfection of Jurkat E6 cells with the Tat gene selectively altered selenoprotein synthesis, with GPX4 and Sep15 being the most inhibited and TR1 the most enhanced. Taken together, the data show that monkeys infected with SIV in vivo and cells infected with SIV in vitro will provide appropriate models for investigating the mechanism(s) responsible for reduced selenium levels that accompany the progression of AIDS in HIV disease.

Global perspectives on proteins: comparing genomes in terms of folds, pathways and beyond

The sequencing of complete genomes provides us with a global view of all the proteins in an organism. Proteomic analysis can be done on a purely sequence-based level, with a focus on finding homologues and grouping them into families and clusters of orthologs. However, incorporating protein structure into this analysis provides valuable simplification; it allows one to collect together very distantly related sequences, thus condensing the proteome into a minimal number of 'parts.' We describe issues related to surveying proteomes in terms of structural parts, including methods for fold assignment and formats for comparisons (eg top-10 lists and whole-genome trees), and show how biases in the databases and in sampling can affect these surveys. We illustrate our main points through a case study on the unique protein properties evident in many thermophile genomes (eg more salt bridges). Finally, we discuss metabolic pathways as an even greater simplification of genomes. In comparison to folds these allow the organization of many more genes into coherent systems, yet can nevertheless be understood in many of the same terms.

Molecular evolution of the nontandemly repeated genes of the histone 3 multigene family

In some species, histone gene clusters consist of tandem arrays of each type of histone gene, whereas in other species the genes may be clustered but not arranged in tandem. In certain species, however, histone genes are found scattered across several different chromosomes. This study examines the evolution of histone 3 (H3) genes that are not arranged in large clusters of tandem repeats. Although H3 amino acid sequences are highly conserved both within and between species, we found that the nucleotide sequence divergence at synonymous sites is high, indicating that purifying selection is the major force for maintaining H3 amino acid sequence homogeneity over long-term evolution. In cases where synonymous-site divergence was low, recent gene duplication appeared to be a better explanation than gene conversion. These results, and other observations on gene inactivation, organization, and phylogeny, indicated that these H3 genes evolve according to a birth-and-death process under strong purifying selection. Thus, we found little evidence to support previous claims that all H3 proteins, regardless of their genome organization, undergo concerted evolution. Further analyses of the structure of H3 proteins revealed that the histones of higher eukaryotes might have evolved from a replication-independent-like H3 gene.

A new moss genetics: targeted mutagenesis in Physcomitrella patens

The potential of moss as a model system to study plant biology is associated with their relatively simple developmental pattern that nevertheless resembles the basic organization of the body plan of land plants, the direct access to cell-lineage analysis, their similar responses to plant growth factors and environmental stimuli as those observed in other land plants, and the dominance of the gametophyte in the life cycle that facilitates genetic approaches. Transformation studies in the moss Physcomitrella patens have revealed a totally unique feature for plants, i.e., that foreign DNA sequences integrate in the genome preferentially at targeted locations by homologous recombination, enabling for the first time in plants the application of the powerful molecular genetic approaches used routinely in bacteria, yeast, and since 1989, the mouse embryonic stem cells. This article reviews our current knowledge of Physcomitrella patens transformation and its unique suitability for functional genomic studies.

Open systems: panoramic views of gene expression

Since their development in the early 1990s, differential gene expression (DGE) technologies have been applied to a multitude of biological challenges, both for the purpose of basic biological research and as a valuable tool for the discovery and development of pharmaceuticals. In this review we survey a class of DGE technologies collectively referred to as 'open' architecture systems. These technologies are distinct from the 'closed' DGE technologies (quantitative PCR, chip technologies), in that no pre-existing biological or sequence information is necessary and they are applicable to any species. Examples of open systems include GeneCalling, SAGE, TOGA, READS, and their progenitor DGE technologies, differential display and cDNA representational difference analysis. We review these technologies and summarize a specific application using GeneCalling for novel gene discovery. Additionally, the significance of data management and experimental design in this new age of expression analysis is discussed.

Evolution of the Na-P(i) cotransport systems

Membrane transport systems for P(i) transport are key elements in maintaining homeostasis of P(i) in organisms as diverse as bacteria and human. Two Na-P(i) cotransporter families with well-described functional properties in vertebrates, namely NaPi-II and NaPi-III, show conserved structural features with prokaryotic origin. A clear vertical relationship can be established among the mammalian protein family NaPi-III, a homologous system in C. elegans, the yeast system Pho89, and the bacterial P(i) transporter Pit. An alternative lineage connects the mammalian NaPi-II-related transporters with homologous proteins from Caenorhabditis elegans and Vibrio cholerae. The present review focuses on the molecular evolution of the NaPi-II protein family. Preliminary results indicate that the NaPi-II homologue cloned from V. cholerae is indeed a functional P(i) transporter when expressed in Xenopus oocytes. The closely related NaPi-II isoforms NaPi-IIa and NaPi-IIb are responsible for regulated epithelial Na-dependent P(i) transport in all vertebrates. Most species express two different NaPi-II proteins with the exception of the flounder and Xenopus laevis, which rely on only a single isoform. Using an RT-PCR-based approach with degenerate primers, we were able to identify NaPi-II-related mRNAs in a variety of vertebrates from different families. We hypothesize that the original NaPi-IIb-related gene was duplicated early in vertebrate development. The appearance of NaPi-IIa correlates with the development of the mammalian nephron.

Genetic research and nutritional individuality

Recent genetic research builds on a base established over the last century by physicians and nutritional scientists, who introduced the concept of biochemical individuality and documented its significance for understanding a wide variety of problems in human health. Current comparative genomic investigations on a variety of organisms (Haemophilus influenzae, Saccharomyces cerevisiae, Caenorhabditis elegans, Drosophila melanogaster, Homo sapiens) have established the existence of numerous orthologs (proteins in different organisms that show significant sequence similarities over 80% of their lengths), suggesting significant conservation of structure and probably some of function as well. At the same time, molecular comparisons among individuals within our own species show the existence of abundant molecular variants, many of which have been shown to have functional significance in nutritional and related metabolic contexts. The combination of biochemical individuality and known functional utilities of allelic variants should converge to create a situation in which nutritional optima can be specified as part of comprehensive lifestyle prescriptions tailored to the needs of each person.

Characterization of 16 novel human genes showing high similarity to yeast sequences

The entire set of open reading frames (ORFs) of Saccharomyces cerevisiae has been used to perform systematic similarity searches against nucleic acid and protein databases: with the aim of identifying interesting homologies between yeast and mammalian genes. Many similarities were detected: mostly with known genes. However: several yeast ORFs were only found to match human partial sequence tags: indicating the presence of human transcripts still uncharacterized that have a homologous counterpart in yeast. About 30 such transcripts were further studied and named HUSSY (human sequence similar to yeast). The 16 most interesting are presented in this paper along with their sequencing and mapping data. As expected: most of these genes seem to be involved in basic metabolic and cellular functions (lipoic acid biosynthesis: ribulose-5-phosphate-3-epimerase: glycosyl transferase: beta-transducin: serine-threonine-kinase: ABC proteins: cation transporters). Genes related to RNA maturation were also found (homologues to DIM1: ROK1-RNA-elicase and NFS1). Furthermore: five novel human genes were detected (HUSSY-03: HUSSY-22: HUSSY-23: HUSSY-27: HUSSY-29) that appear to be homologous to yeast genes whose function is still undetermined. More information on this work can be obtained at the website http://grup.bio.unipd.it/hussy

Differences between predicted and observed sequences in Saccharomyces cerevisiae

We recently studied the protein composition of a Saccharomyces cerevisiae wine yeast strain (K310) of enological interest. About 2,500 spots of 8-250 kDa observed molecular mass were resolved by two-dimensional gel electrophoresis. Experimental molecular masses and isoelectric points were calculated for most of them. Twenty-seven proteins were subjected to Edman microsequencing. N-terminal sequences of 12/27 proteins were determined, whereas internal sequences of 6/27 proteins were obtained following in situ proteolysis. Comparison between the experimental data and those reported in the SWISS-PROT database revealed some differences between genotypic and phenotypic sequences. These are indicative of the changes a protein can undergo with respect to the primary structure coded by the genomic DNA. Our results highlight the need to complement genomic analysis with detailed proteomics in order to refine the vast amount of information provided by DNA sequencing and to find an exact correlation between genome and proteome.

Partial deletion of a loop region in the high affinity K+ transporter HKT1 changes ionic permeability leading to increased salt tolerance

HKT1 is a high affinity K(+) transporter protein that is a member of a large superfamily of transporters found in plants, bacteria, and fungi. These transporters are primarily involved in K(+) uptake and are energized by Na(+) or H(+). HKT1 is energized by Na(+) but also mediates low affinity Na(+) uptake and may therefore be a pathway for Na(+) uptake, which is toxic to plants. The aim of this study was to identify regions of HKT1 that are involved in K(+)/Na(+) selectivity and alter the amino acid composition in those regions to increase the ionic selectivity of the transporter. A highly charged loop was identified, and two deletions were created that resulted in the removal of charged and uncharged amino acids. The functional changes caused by the deletions were studied in yeast and Xenopus oocytes. The deletions improved the K(+)/Na(+) selectivity of the transporter and increased the salt tolerance of the yeast cells in which they were expressed. In light of recent structural models of members of this symporter superfamily, it was necessary to determine the orientation of this highly charged loop. Introduction of an epitope tag allowed us to demonstrate that this loop faces the outside of the membrane where it is likely to facilitate the interaction with cations such as K(+) and Na(+). This study has identified an important structural feature in HKT1 that in part determines its K(+)/Na(+) selectivity. Understanding the structural basis of the functional characteristics in transporters such as HKT1 may have important implications for increasing the salt tolerance of higher plants.

One of two genes encoding glycyl-tRNA synthetase in Saccharomyces cerevisiae provides mitochondrial and cytoplasmic functions

In the yeast Saccharomyces cerevisiae, two genes (GRS1 and GRS2) encode glycyl-tRNA synthetase (GlyRS1 and GlyRS2, respectively). 59% of the sequence of GlyRS2 is identical to that of GlyRS1. Others have proposed that GRS1 and GRS2 encode the cytoplasmic and mitochondrial enzymes, respectively. In this work, we show that GRS1 encodes both functions, whereas GRS2 is dispensable. In addition, both cytoplasmic and mitochondrial phenotypes of the knockout allele of GRS1 in S. cerevisiae are complemented by the expression of the only known gene for glycyl-tRNA synthetase in Schizosaccharomyces pombe. Thus, a single gene for glycyl-tRNA synthetase likely encodes both cytoplasmic and mitochondrial activities in most or all yeast. Phylogenetic analysis shows that GlyRS2 is a predecessor of all yeast GlyRS homologues. Thus, GRS1 appears to be the result of a duplication of GRS2, which itself is pseudogene-like.

Distribution and abundance of microsatellites in the yeast genome can Be explained by a balance between slippage events and point mutations

We fit a Markov chain model of microsatellite evolution introduced by Kruglyak et al. to data on all di-, tri-, and tetranucleotide repeats in the yeast genome. Our results suggest that many features of the distribution of abundance and length of microsatellites can be explained by this simple model, which incorporates a competition between slippage events and base pair substitutions, with no need to invoke selection or constraints on the lengths. Our results provide some new information on slippage rates for individual repeat motifs, which suggest that AT-rich trinucleotide repeats have higher slippage rates. As our model predicts, we found that many repeats were adjacent to shorter repeats of the same motif. However, we also found a significant tendency of microsatellites of different motifs to cluster.

Searching for drug targets in microbial genomes

Comparative analysis of the complete genome sequences of 10 bacterial pathogens available in the public databases offers the first insights into the drug discovery approaches of the near future. Genes that are conserved in different genomes often turn out to be essential, which makes them attractive targets for new broad-spectrum antibiotics. Subtractive genome analysis reveals the genes that are conserved in all or most of the pathogenic bacteria but not in eukaryotes; these are the most obvious candidates for drug targets. Species-specific genes, on the other hand, may offer the possibility to design drugs against a particular, narrow group of pathogens.

Genome evolution and the evolution of exon-shuffling--a review

Recent studies on the genomes of protists, plants, fungi and animals confirm that the increase in genome size and gene number in different eukaryotic lineages is paralleled by a general decrease in genome compactness and an increase in the number and size of introns. It may thus be predicted that exon-shuffling has become increasingly significant with the evolution of larger, less compact genomes. To test the validity of this prediction, we have analyzed the evolutionary distribution of modular proteins that have clearly evolved by intronic recombination. The results of this analysis indicate that modular multidomain proteins produced by exon-shuffling are restricted in their evolutionary distribution. Although such proteins are present in all major groups of metazoa from sponges to chordates, there is practically no evidence for the presence of related modular proteins in other groups of eukaryotes. The biological significance of this difference in the composition of the proteomes of animals, fungi, plants and protists is best appreciated when these modular proteins are classified with respect to their biological function. The majority of these proteins can be assigned to functional categories that are inextricably linked to multicellularity of animals, and are of absolute importance in permitting animals to function in an integrated fashion: constituents of the extracellular matrix, proteases involved in tissue remodelling processes, various proteins of body fluids, membrane-associated proteins mediating cell-cell and cell-matrix interactions, membrane associated receptor proteins regulating cell cell communications, etc. Although some basic types of modular proteins seem to be shared by all major groups of metazoa, there are also groups of modular proteins that appear to be restricted to certain evolutionary lineages. In summary, the results suggest that exon-shuffling acquired major significance at the time of metazoan radiation. It is interesting to note that the rise of exon-shuffling coincides with a spectacular burst of evolutionary creativity: the Big Bang of metazoan radiation. It seems probable that modular protein evolution by exon-shuffling has contributed significantly to this accelerated evolution of metazoa, since it facilitated the rapid construction of multidomain extracellular and cell surface proteins that are indispensable for multicellularity.

Structural models of the KtrB, TrkH, and Trk1,2 symporters based on the structure of the KcsA K(+) channel

Three-dimensional computer modeling is used to further investigate the hypothesis forwarded in the accompanying paper of an evolutionary relationship between four related families of K(+) sympoter proteins and the superfamily of K(+) channel proteins. Atomic-scale models are developed for the transmembrane regions of one member from each of the three more distinct symporter families, i.e., a TrkH protein from Escherichia coli, a KtrB protein from Aquifex aeolicus, and a Trk1,2 protein from Schizosaccharomyces pombe. The portions of the four consecutive M1-P-M2 motifs in the symporters that can be aligned with K(+) channel sequences are modeled directly from the recently determined crystal structure of the KcsA K(+) channel from Streptomyces lividans. The remaining portions are developed using our previously accumulated theoretical modeling criteria and principles. Concurrently, the use of these criteria and principles is further supported by the now verified predictions of our previous K(+) channel modeling efforts and the degree to which they are satisfied by the known structure of the KcsA protein. Thus the observed ability of the portions of the symporter models derived from the KcsA crystal structure to also satisfy the theoretical modeling criteria provides additional support for an evolutionary link with K(+) channel proteins. Efforts to further satisfy the criteria and principles suggest that the symporter proteins from fungi and plants (i.e., Trk1,2 and HKT1) form dimeric and/or tetrameric complexes in the membrane. Furthermore, analysis of the atomic-scale models in relation to the sequence conservation within and between the protein families suggests structural details for previously proposed mechanisms for the linked symport of K(+) with Na(+) and H(+). Suggestions are also given for experiments to test these structures and hypotheses.

Evolutionary relationship between K(+) channels and symporters

The hypothesis is presented that at least four families of putative K(+) symporter proteins, Trk and KtrAB from prokaryotes, Trk1,2 from fungi, and HKT1 from wheat, evolved from bacterial K(+) channel proteins. Details of this hypothesis are organized around the recently determined crystal structure of a bacterial K(+) channel: i. e., KcsA from Streptomyces lividans. Each of the four identical subunits of this channel has two fully transmembrane helices (designated M1 and M2), plus an intervening hairpin segment that determines the ion selectivity (designated P). The symporter sequences appear to contain four sequential M1-P-M2 motifs (MPM), which are likely to have arisen from gene duplication and fusion of the single MPM motif of a bacterial K(+) channel subunit. The homology of MPM motifs is supported by a statistical comparison of the numerical profiles derived from multiple sequence alignments formed for each protein family. Furthermore, these quantitative results indicate that the KtrAB family of symporters has remained closest to the single-MPM ancestor protein. Strong sequence evidence is also found for homology between the cytoplasmic C-terminus of numerous bacterial K(+) channels and the cytoplasm-resident TrkA and KtrA subunits of the Trk and KtrAB symporters, which in turn are homologous to known dinucleotide-binding domains of other proteins. The case for homology between bacterial K(+) channels and the four families of K(+) symporters is further supported by the accompanying manuscript, in which the patterns of residue conservation are demonstrated to be similar to each other and consistent with the known 3D structure of the KcsA K(+) channel.

The physiological contribution of Acinetobacter PcaK, a transport system that acts upon protocatechuate, can be masked by the overlapping specificity of VanK

VanK is the fourth member of the ubiquitous major facilitator superfamily of transport proteins to be identified that, together with PcaK, BenK, and MucK, contributes to aromatic catabolism in Acinetobacter sp. strain ADP1. VanK and PcaK have overlapping specificity for p-hydroxybenzoate and, most clearly, for protocatechuate: inactivation of both proteins severely impairs growth with protocatechuate, and the activity of either protein alone can mask the phenotype associated with inactivation of its homolog. Furthermore, vanK pcaK double-knockout mutants appear completely unable to grow in liquid culture with the hydroaromatic compound quinate, although such cells on plates convert quinate to protocatechuate, which then accumulates extracellularly and is readily visible as purple staining. This provides genetic evidence that quinate is converted to protocatechuate in the periplasm and is in line with the early argument that quinate catabolism should be physically separated from aromatic amino acid biosynthesis in the cytoplasm so as to avoid potential competition for intermediates common to both pathways. Previous studies of aromatic catabolism in Acinetobacter have taken advantage of the ability to select directly strains that contain a spontaneous mutation blocking the beta-ketoadipate pathway and preventing the toxic accumulation of carboxymuconate. By using this procedure, strains with a mutation in structural or regulatory genes blocking degradation of vanillate, p-hydroxybenzoate, or protocatechuate were selected. In this study, the overlapping specificity of the VanK and PcaK permeases was exploited to directly select strains with a mutation in either vanK or pcaK. Spontaneous mutations identified in vanK include a hot spot for frameshift mutation due to contraction of a G6 mononucleotide repeat as well as point mutations producing amino acid substitutions useful for analysis of VanK structure and function. Preliminary second-site suppression analysis using transformation-facilitated PCR mutagenesis in one VanK mutant gave results similar to those using LacY, the prototypic member of the major facilitator superfamily, consistent with the two proteins having a similar mechanism of action. The selection for transport mutants described here for Acinetobacter may also be applicable to Pseudomonas putida, where the PcaK permease has an additional role in chemotaxis.

Change to alanine of one out of four selectivity filter glycines in KtrB causes a two orders of magnitude decrease in the affinities for both K+ and Na+ of the Na+ dependent K+ uptake system KtrAB from Vibrio alginolyticus

KtrAB from Vibrio alginolyticus is a recently described new type of high affinity bacterial K+ uptake system. Its activity assayed in an Escherichia coli K+ uptake negative mutant depended on Na+ ions (Km of 40 microM). Subunit KtrB contains four putative P-loops. The selectivity filter from each P-loop contains a conserved glycine residue. Residue Gly-290 from the third P-loop selectivity filter in KtrB was exchanged for Ala, Ser or Asp. KtrB variants Ser-290 and Asp-290 were without activity. In contrast, KtrB variant Ala-290 was still active. This variant transported K+ with a two orders of magnitude decrease in apparent affinity for both K+ and Na+ with little effect on Vmax.

Proteins connecting the nuclear pore complex with the nuclear interior

While much has been learned in recent years about the movement of soluble transport factors across the nuclear pore complex (NPC), comparatively little is known about intranuclear trafficking. We isolated the previously identified Saccharomyces protein Mlp1p (myosin-like protein) by an assay designed to find nuclear envelope (NE) associated proteins that are not nucleoporins. We localized both Mlp1p and a closely related protein that we termed Mlp2p to filamentous structures stretching from the nucleoplasmic face of the NE into the nucleoplasm, similar to the homologous vertebrate and Drosophila Tpr proteins. Mlp1p can be imported into the nucleus by virtue of a nuclear localization sequence (NLS) within its COOH-terminal domain. Overexpression experiments indicate that Mlp1p can form large structures within the nucleus which exclude chromatin but appear highly permeable to proteins. Remarkably, cells harboring a double deletion of MLP1 and MLP2 were viable, although they showed a slower net rate of active nuclear import and faster passive efflux of a reporter protein. Our data indicate that the Tpr homologues are not merely NPC-associated proteins but that they can be part of NPC-independent, peripheral intranuclear structures. In addition, we suggest that the Tpr filaments could provide chromatin-free conduits or tracks to guide the efficient translocation of macromolecules between the nucleoplasm and the NPC.

Patterns of protein-fold usage in eight microbial genomes: a comprehensive structural census

Eight microbial genomes are compared in terms of protein structure. Specifically, yeast, H. influenzae, M. genitalium, M. jannaschii, Synechocystis, M. pneumoniae, H. pylori, and E. coli are compared in terms of patterns of fold usage-whether a given fold occurs in a particular organism. Of the approximately 340 soluble protein folds currently in the structure databank (PDB), 240 occur in at least one of the eight genomes, and 30 are shared amongst all eight. The shared folds are depleted in allhelical structure and enriched in mixed helix-sheet structure compared to the folds in the PDB. The top-10 most common of the shared 30 are enriched in superfolds, uniting many non-homologous sequence families, and are especially similar in overall architecture-eight having helices packed onto a central sheet. They are also very different from the common folds in the PBD, highlighting databank biases. Folds can be ranked in terms of expression as well as genome duplication. In yeast the top-10 most highly expressed folds are considerably different from the most highly duplicated folds. A tree can be constructed grouping genomes in terms of their shared folds. This has a remarkably similar topology to more conventional classifications, based on very different measures of relatedness. Finally, folds of membrane proteins can be analyzed through transmembrane-helix (TM) prediction. All the genomes appear to have similar usage patterns for these folds, with the occurrence of a particular fold falling off rapidly with increasing numbers of TM-elements, according to a "Zipf-like" law. This implies there are no marked preferences for proteins with particular numbers of TM-helices (e.g. 7-TM) in microbial genomes.

Constructing multigenome views of whole microbial genomes

We have designed and implemented a system to carry out cross-genome comparisons of open reading frames (ORFs) from multiple genomes. This implementation includes a genome profiling system that allows us to explore pairwise comparisons at different levels of match similarity and ask biologically motivated queries involving number and identity of ORFs, their function, functional category, distribution in genomes or in biological domains, and statistics on their matches and match families. This analysis required precise definition of new classification terms and concepts. We define the terms genomic signature, summary signature, biologic domain signature, domain class, match level, match family, and extended match family, then use these terms to define concepts, including genomically universal proteins and proteins characteristics of sets of genomes. We initiate an analysis based on automated FASTA (Pearson, 1996) comparison of 22,419 conceptually translated protein sequences from nine microbial genomes.

Sulfur reduction by human erythrocytes

Washed human erythrocytes incubated with glucose and S8 and purged with N2 produced H2S at a nearly constant rate of 170 mumol (L cells)-1 min-1, which continued for several hours. In sealed vials up to 25 mM HS- accumulated. Glucose caused the fastest H2S production, although either lactate or glycerol could support slower rates. When glucose was added without S8, anoxic H2S production nonetheless occurred at about 1.5% of the maximum rate, after 24 hr totaling 0.5 mmol H2S (L cells)-1, suggesting the presence of endogenous reducible sulfur. Anaerobic conditions were not required, since oxygenated cells produced H2S from S8 at 80% of the anoxic rate. Using cell lysates, production of H2S occurred after addition of either glutathione, NADH, or NADPH. The observations suggest possible physiological roles for H2S as an electron carrier, and are consistent with an evolutionary relationship between eukaryotic cytoplasm and sulfur-reducing Archaea.

Nutritional management of enteropathy

The prognosis for nutritional management of enteropathy in children is good when the enteropathy is reversible with the use of a food elimination diet, such as cow's-milk-sensitive enteropathy, but is poor when enteropathy is irreversible, such as microvillous atrophy. However, nutritional management is central to the care of all children with small intestinal enteropathy. Enteral nutrition (provision of liquid formula diets by mouth or by tube) is possible in most cases, but in some children with intractable diarrhea, parenteral nutrition needs to supplement enteral feeding. The choice of enteral feeding ranges from elemental to partial hydrolysate.

Findings emerging from complete microbial genome sequences

Sixteen microorganisms, including one eukaryote, four archaeons, and 11 eubacteria, have been completely sequenced and published. More than 50 genomes are scheduled to be completed by the year 2000. This explosive growth of information is forcing change in many scientific disciplines (e.g. bioinformatics and molecular genetics), spawning new fields, and even changing the way scientific information is used and shared. Novel, global genome sequence comparisons seem slow to appear but the infrastructure for these projects is being built, and we expect exciting developments in the near future.

Comparing genomes in terms of protein structure: surveys of a finite parts list

We give an overview of the emerging field of structural genomics, describing how genomes can be compared in terms of protein structure. As the number of genes in a genome and the total number of protein folds are both quite limited, these comparisons take the form of surveys of a finite parts list, similar in respects to demographic censuses. Fold surveys have many similarities with other whole-genome characterizations, e.g., analyses of motifs or pathways. However, structure has a number of aspects that make it particularly suitable for comparing genomes, namely the way it allows for the precise definition of a basic protein module and the fact that it has a better defined relationship to sequence similarity than does protein function. An essential requirement for a structure survey is a library of folds, which groups the known structures into 'fold families.' This library can be built up automatically using a structure comparison program, and we described how important objective statistical measures are for assessing similarities within the library and between the library and genome sequences. After building the library, one can use it to count the number of folds in genomes, expressing the results in the form of Venn diagrams and 'top-10' statistics for shared and common folds. Depending on the counting methodology employed, these statistics can reflect different aspects of the genome, such as the amount of internal duplication or gene expression. Previous analyses have shown that the common folds shared between very different microorganisms, i.e., in different kingdoms, have a remarkably similar structure, being comprised of repeated strand-helix-strand super-secondary structure units. A major difficulty with this sort of 'fold-counting' is that only a small subset of the structures in a complete genome are currently known and this subset is prone to sampling bias. One way of overcoming biases is through structure prediction, which can be applied uniformly and comprehensively to a whole genome. Various investigators have, in fact, already applied many of the existing techniques for predicting secondary structure and transmembrane (TM) helices to the recently sequenced genomes. The results have been consistent: microbial genomes have similar fractions of strands and helices even though they have significantly different amino acid composition. The fraction of membrane proteins with a given number of TM helices falls off rapidly with more TM elements, approximately according to a Zipf law. This latter finding indicates that there is no preference for the highly studied 7-TM proteins in microbial genomes. Continuously updated tables and further information pertinent to this review are available over the web at http://bioinfo.mbb.yale.edu/genome.

Structure-function relationships in OxlT, the oxalate/formate transporter of Oxalobacter formigenes. Topological features of transmembrane helix 11 as visualized by site-directed fluorescent labeling

Analysis of hydropathy suggests that in OxlT, the oxalate/formate antiporter of Oxalobacter formigenes, lysine 355 is within transmembrane helix no. 11. To test this idea, we used single-cysteine, histidine-tagged OxlT variants to study the organization of a 30-residue segment (residues 344-373) containing this region. Topology was examined by probing the A345C and A370C proteins with Oregon Green maleimide carboxylic acid, an impermeant and fluorescent thiol-reactive agent. Examination of purified protein showed that only A370C was fluorescent after treating intact cells with the probe, while both proteins were modified in tests with isolated membrane ghosts. In addition, labeling of A370C, but not A345C, was blocked when external cysteines were protected with the impermeant and nonfluorescent agent, methanethiosulfonate ethyltrimethylammonium. These findings confirm that A345 faces the cytoplasm, while A370C faces the periplasm. A similar study focused on 13 single-cysteine variants positioned throughout the target segment. That work revealed a striking discontinuity in reactivity toward Oregon Green maleimide; cysteines within a 10-residue central core (residues 351-360) were not labeled when membranes were probed, but were readily modified after protein denaturation. We suggest this core resides within the lipid bilayer, unavailable to an impermeant reporter. Since this region includes position 355, we also suggest that lysine 355 lies within the OxlT hydrophobic sector, where it may facilitate the binding and translocation of the anionic substrates, oxalate and formate.