PDB-REPRDB: a database of representative protein chains from the Protein Data Bank (PDB)

PDB-REPRDB is a database of representative protein chains from the Protein Data Bank (PDB). The previous version of PDB-REPRDB provided 48 representative sets, whose similarity criteria were predetermined, on the WWW. The current version is designed so that the user may obtain a quick selection of representative chains from PDB. The selection of representative chains can be dynamically configured according to the user’s requirement. The WWW interface provides a large degree of freedom in setting parameters, such as cut-off scores of sequence and structural similarity. One can obtain a representative list and classification data of protein chains from the system. The current database includes 20 457 protein chains from PDB entries (August 6, 2000). The system for PDB-REPRDB is available at the Parallel Protein Information Analysis system (PAPIA) WWW server (http://www.rwcp.or.jp/papia/).

The PDB data uniformity project

The Protein Data Bank (PDB; http://www.rcsb.org/pdb/) is the single worldwide archive of structural data of biological macromolecules. This paper describes the data uniformity project that is underway to address the inconsistency in PDB data.

Substrate-induced conformational change in a trimeric ornithine transcarbamoylase

The crystal structure of Escherichia coli ornithine transcarbamoylase (OTCase, EC 2.1.3.3) complexed with the bisubstrate analog N-(phosphonacetyl)-l-ornithine (PALO) has been determined at 2.8-Å resolution. This research on the structure of a transcarbamoylase catalytic trimer with a substrate analog bound provides new insights into the linkages between substrate binding, protein–protein interactions, and conformational change. The structure was solved by molecular replacement with the Pseudomonas aeruginosa catabolic OTCase catalytic trimer (Villeret, V., Tricot, C., Stalon, V. & Dideberg, O. (1995) Proc. Natl. Acad. Sci. USA 92, 10762–10766; Protein Data Bank reference pdb 1otc) as the model and refined to a crystallographic R value of 21.3%. Each polypeptide chain folds into two domains, a carbamoyl phosphate binding domain and an l-ornithine binding domain. The bound inhibitor interacts with the side chains and/or backbone atoms of Lys-53, Ser-55, Thr-56, Arg-57, Thr-58, Arg-106, His-133, Asn-167, Asp-231, Met-236, Leu-274, Arg-319 as well as Gln-82 and Lys-86 from an adjacent chain. Comparison with the unligated P. aeruginosa catabolic OTCase structure indicates that binding of the substrate analog results in closure of the two domains of each chain. As in E. coli aspartate transcarbamoylase, the 240s loop undergoes the largest conformational change upon substrate binding. The clinical implications for human OTCase deficiency are discussed.

Structural assignments to the Mycoplasma genitalium proteins show extensive gene duplications and domain rearrangements

The parasitic bacterium Mycoplasma genitalium has a small, reduced genome with close to a basic set of genes. As a first step toward determining the families of protein domains that form the products of these genes, we have used the multiple sequence programs psi-blast and geanfammer to match the sequences of the 467 gene products of M. genitalium to the sequences of the domains that form proteins of known structure [Protein Data Bank (PDB) sequences]. PDB sequences (274) match all of 106 M. genitalium sequences and some parts of another 85; thus, 41% of its total sequences are matched in all or part. The evolutionary relationships of the PDB domains that match M. genitalium are described in the structural classification of proteins (SCOP) database. Using this information, we show that the domains in the matched M. genitalium sequences come from 114 superfamilies and that 58% of them have arisen by gene duplication. This level of duplication is more than twice that found by using pairwise sequence comparisons. The PDB domain matches also describe the domain structure of the matched sequences: just over a quarter contain one domain and the rest have combinations of two or more domains.

The EMBL nucleotide sequence database

The EMBL Nucleotide Sequence Database (http://www.ebi.ac.uk/embl/) is maintained at the European Bioinformatics Institute (EBI) in an international collaboration with the DNA Data Bank of Japan (DDBJ) and GenBank at the NCBI (USA). Data is exchanged amongst the collaborating databases on a daily basis. The major contributors to the EMBL database are individual authors and genome project groups. Webin is the preferred web-based submission system for individual submitters, whilst automatic procedures allow incorporation of sequence data from large-scale genome sequencing centres and from the European Patent Office (EPO). Database releases are produced quarterly. Network services allow free access to the most up-to-date data collection via ftp, email and World Wide Web interfaces. EBI’s Sequence Retrieval System (SRS), a network browser for databanks in molecular biology, integrates and links the main nucleotide and protein databases plus many specialized databases. For sequence similarity searching a variety of tools (e.g. Blitz, Fasta, BLAST) are available which allow external users to compare their own sequences against the latest data in the EMBL Nucleotide Sequence Database and SWISS-PROT.

PDBsum: summaries and analyses of PDB structures

PDBsum is a web-based database providing a largely pictorial summary of the key information on each macromolecular structure deposited at the Protein Data Bank (PDB). It includes images of the structure, annotated plots of each protein chain’s secondary structure, detailed structural analyses generated by the PROMOTIF program, summary PROCHECK results and schematic diagrams of protein–ligand and protein–DNA interactions. RasMol scripts highlight key aspects of the structure, such as the protein’s domains, PROSITE patterns and protein–ligand interactions, for interactive viewing in 3D. Numerous links take the user to related sites. PDBsum is updated whenever any new structures are released by the PDB and is freely accessible via http://www.biochem.ucl.ac.uk/bsm/pdbsum.

A fully automatic evolutionary classification of protein folds: Dali Domain Dictionary version 3

The Dali Domain Dictionary (http://www.ebi.ac.uk/dali/domain) is a numerical taxonomy of all known structures in the Protein Data Bank (PDB). The taxonomy is derived fully automatically from measurements of structural, functional and sequence similarities. Here, we report the extension of the classification to match the traditional four hierarchical levels corresponding to: (i) supersecondary structural motifs (attractors in fold space), (ii) the topology of globular domains (fold types), (iii) remote homologues (functional families) and (iv) homologues with sequence identity above 25% (sequence families). The computational definitions of attractors and functional families are new. In September 2000, the Dali classification contained 10 531 PDB entries comprising 17 101 chains, which were partitioned into five attractor regions, 1375 fold types, 2582 functional families and 3724 domain sequence families. Sequence families were further associated with 99 582 unique homologous sequences in the HSSP database, which increases the number of effectively known structures several-fold. The resulting database contains the description of protein domain architecture, the definition of structural neighbours around each known structure, the definition of structurally conserved cores and a comprehensive library of explicit multiple alignments of distantly related protein families.

The RESID Database of protein structure modifications and the NRL-3D Sequence–Structure Database

The RESID Database is a comprehensive collection of annotations and structures for protein post-translational modifications including N-terminal, C-terminal and peptide chain cross-link modifications. The RESID Database includes systematic and frequently observed alternate names, Chemical Abstracts Service registry numbers, atomic formulas and weights, enzyme activities, taxonomic range, keywords, literature citations with database cross-references, structural diagrams and molecular models. The NRL-3D Sequence–Structure Database is derived from the three-dimensional structure of proteins deposited with the Research Collaboratory for Structural Bioinformatics Protein Data Bank. The NRL-3D Database includes standardized and frequently observed alternate names, sources, keywords, literature citations, experimental conditions and searchable sequences from model coordinates. These databases are freely accessible through the National Cancer Institute–Frederick Advanced Biomedical Computing Center at these web sites: http://www.ncifcrf.gov/RESID, http://www.ncifcrf.gov/ NRL-3D; or at these National Biomedical Research Foundation Protein Information Resource web sites: http://pir.georgetown.edu/pirwww/dbinfo/resid.html, http://pir.georgetown.edu/pirwww/dbinfo/nrl3d.html

Hydrophobic moments of protein structures: Spatially profiling the distribution

It is generally accepted that globular proteins fold with a hydrophobic core and a hydrophilic exterior. Might the spatial distribution of amino acid hydrophobicity exhibit common features? The hydrophobic profile detailing this distribution from the protein interior to exterior has been examined for 30 relatively diverse structures obtained from the Protein Data Bank, for 3 proteins of the 30S ribosomal subunit, and for a simple set of 14 decoys. A second-order hydrophobic moment has provided a simple measure of the spatial variation. Shapes of the calculated spatial profiles of all native structures have been found to be comparable. Consequently, profile shapes as well as particular profile features should assist in validating predicted protein structures and in discriminating between different protein-folding pathways. The spatial profiles of the 14 decoys are clearly distinguished from the profiles of their native structures.

Specific atomic groups and RNA helix geometry in acceptor stem recognition by a tRNA synthetase

Oligonucleotides that recapitulate the acceptor stems of tRNAs are substrates for aminoacylation by many tRNA synthetases in vitro, even though these substrates are missing the anticodon trinucleotides of the genetic code. In the case of tRNAAla a single acceptor stem G⋅U base pair at position 3·70 is essential, based on experiments where the wobble pair has been replaced by alternatives such as I⋅U, G⋅C, and A⋅U, among others. These experiments led to the conclusion that the minor-groove free 2-amino group (of guanosine) of the G⋅U wobble pair is essential for charging. Moreover, alanine-inserting tRNAs (amber suppressors) that replace G⋅U with mismatches such as G⋅A and C⋅A are partially active in vivo and can support growth of an Escherichia coli tRNAAla knockout strain, leading to the hypothesis that a helix irregularity and nucleotide functionalities are important for recognition. Herein we investigate the charging in vitro of oligonucleotide and full-length tRNA substrates that contain mismatches at the position of the G⋅U pair. Although most of these substrates have undetectable activity, G⋅A and C⋅A variants retain some activity, which is, nevertheless, reduced by at least 100-fold. Thus, the in vivo assays are much less sensitive to large changes in aminoacylation kinetic efficiency of 3·70 variants than is the in vitro assay system. Although these functional data do not clarify all of the details, it is now clear that specific atomic groups are substantially more important in determining kinetic efficiency than is a helical distortion. By implication, the activity of mutant tRNAs measured in the in vivo assays appears to be more dependent on factors other than aminoacylation kinetic efficiency.

Integration of multiple instructive cues by neural crest stem cells reveals cell-intrinsic biases in relative growth factor responsiveness

Growth factors can influence lineage determination of neural crest stem cells (NCSCs) in an instructive manner, in vitro. Because NCSCs are likely exposed to multiple signals in vivo, these findings raise the question of how stem cells would integrate such combined influences. Bone morphogenetic protein 2 (BMP2) promotes neuronal differentiation and glial growth factor 2 (GGF2) promotes glial differentiation; if NCSCs are exposed to saturating concentrations of both factors, BMP2 appears dominant. By contrast, if the cells are exposed to saturating concentrations of both BMP2 and transforming growth factor β1 (which promotes smooth muscle differentiation), the two factors appear codominant. Sequential addition experiments indicate that NCSCs require 48–96 hrs in GGF2 before they commit to a glial fate, whereas the cells commit to a smooth muscle fate within 24 hr in transforming growth factor β1. The delayed response to GGF2 does not reflect a lack of functional receptors; however, because the growth factor induces rapid mitogen-activated protein kinase phosphorylation in naive cells. Furthermore, GGF2 can attenuate induction of the neurogenic transcription factor mammalian achaete-scute homolog 1, by low doses of BMP2. This short-term antineurogenic influence of GGF2 is not sufficient for glial lineage commitment, however. These data imply that NCSCs exhibit cell-intrinsic biases in the timing and relative dosage sensitivity of their responses to instructive factors that influence the outcome of lineage decisions in the presence of multiple factors. The relative delay in glial lineage commitment, moreover, apparently reflects successive short-term and longer-term actions of GGF2. Such a delay may help to explain why glia normally differentiate after neurons, in vivo.

In vivo telomere dynamics of human hematopoietic stem cells

Aging in vivo and cell division in vitro are associated with telomere shortening. Several lines of evidence suggest that telomere length may be a good predictor of the long term replicative capacity of cells. To investigate the natural fate of chromosome telomeres of hematopoietic stem cells in vivo, we measured the telomere length of peripheral blood granulocytes from 11 fully engrafted bone marrow transplant recipients and from their respective donors. In 10 of 11 donor–recipient pairs, the telomere length was significantly reduced in the recipient and the extent of reduction correlated inversely with the number of nucleated cells infused. These data provide internally controlled in vivo evidence that, concomitantly with their proliferation, hematopoietic stem cells lose telomere length; it is possible that, as a result, their proliferative potential is reduced. These findings must be taken into account when developing new protocols in which few stem cells are used for bone marrow transplantation or for gene therapy.

Somatic mosaicism in Fanconi anemia: Evidence of genotypic reversion in lymphohematopoietic stem cells

Somatic mosaicism has been observed previously in the lymphocyte population of patients with Fanconi anemia (FA). To identify the cellular origin of the genotypic reversion, we examined each lymphohematopoietic and stromal cell lineage in an FA patient with a 2815–2816ins19 mutation in FANCA and known lymphocyte somatic mosaicism. DNA extracted from individually plucked peripheral blood T cell colonies and marrow colony-forming unit granulocyte–macrophage and burst-forming unit erythroid cells revealed absence of the maternal FANCA exon 29 mutation in 74.0%, 80.3%, and 86.2% of colonies, respectively. These data, together with the absence of the FANCA exon 29 mutation in Epstein–Barr virus-transformed B cells and its presence in fibroblasts, indicate that genotypic reversion, most likely because of back mutation, originated in a lymphohematopoietic stem cell and not solely in a lymphocyte population. Contrary to a predicted increase in marrow cellularity resulting from reversion in a hematopoietic stem cell, pancytopenia was progressive. Additional evaluations revealed a partial deletion of 11q in 3 of 20 bone marrow metaphase cells. By using interphase fluorescence in situ hybridization with an MLL gene probe mapped to band 11q23 to identify colony-forming unit granulocyte–macrophage and burst-forming unit erythroid cells with the 11q deletion, the abnormal clone was exclusive to colonies with the FANCA exon 29 mutation. Thus, we demonstrate the spontaneous genotypic reversion in a lymphohematopoietic stem cell. The subsequent development of a clonal cytogenetic abnormality in nonrevertant cells suggests that ex vivo correction of hematopoietic stem cells by gene transfer may not be sufficient for providing life-long stable hematopoiesis in patients with FA.

Regulation of the stem cell leukemia (SCL) gene: A tale of two fishes

The stem cell leukemia (SCL) gene encodes a tissue-specific basic helix–loop–helix (bHLH) protein with a pivotal role in hemopoiesis and vasculogenesis. Several enhancers have been identified within the murine SCL locus that direct reporter gene expression to subdomains of the normal SCL expression pattern, and long-range sequence comparisons of the human and murine SCL loci have identified additional candidate enhancers. To facilitate the characterization of regulatory elements, we have sequenced and analyzed 33 kb of the SCL genomic locus from the pufferfish Fugu rubripes, a species with a highly compact genome. Although the pattern of SCL expression is highly conserved from mammals to teleost fish, the genes flanking pufferfish SCL were unrelated to those known to flank both avian and mammalian SCL genes. These data suggest that SCL regulatory elements are confined to the region between the upstream and downstream flanking genes, a region of 65 kb in human and 8.5 kb in pufferfish. Consistent with this hypothesis, the entire 33-kb pufferfish SCL locus directed appropriate expression to hemopoietic and neural tissue in transgenic zebrafish embryos, as did a 10.4-kb fragment containing the SCL gene and extending to the 5′ and 3′ flanking genes. These results demonstrate the power of combining the compact genome of the pufferfish with the advantages that zebrafish provide for studies of gene regulation during development. Furthermore, the pufferfish SCL locus provides a powerful tool for the manipulation of hemopoiesis and vasculogenesis in vivo.

Thrombopoietic potential and serial repopulating ability of murine hematopoietic stem cells constitutively expressing interleukin 11.

Based on transplantation studies with bone marrow cultured under various conditions, a role of interleukin 11 (IL-11) in the self-renewal and/or the differentiation commitment of hematopoietic stem cells has been indicated. To better evaluate the in vivo effects of IL-11 on stem/progenitor cell biology, lethally irradiated mice were serially transplanted with bone marrow cells transduced with a defective retrovirus, termed MSCV-mIL-11, carrying the murine IL-11 (mIL-11) cDNA and the bacterial neomycin phosphotransferase (neo) gene. High serum levels (i.e., > 1 ng/ml) of mIL-11 in all (20/20) primary and 86% (12/14) of secondary long-term reconstituted mice, as well as 86% (12/14) of tertiary recipients examined at 6 weeks posttransplant, demonstrated persistence of vector expression subsequent to transduction of bone marrow precursors functionally definable as totipotent hematopoietic stem cells. In agreement with results obtained with human IL-11 in other myeloablation models, ectopic mIL-11 expression accelerated recovery of platelets, neutrophils, and, to some extent, total leukocytes while preferentially increasing peripheral platelet counts in fully reconstituted mice. When analyzed 5 months posttransplant, tertiary MSCV-mIL-11 recipients had a significantly greater percentage of G418-resistant colony-forming cells in their bone marrow compared with control MSCV animals. Collectively, these data show that persistent stimulation of platelet production by IL-11 is not detrimental to stem cell repopulating ability; rather, they suggest that IL-11 expression in vivo may have resulted in enhanced maintenance of the most primitive hematopoietic stem cell compartment. The prolonged expression achieved by the MSCV retroviral vector, despite the presence of a selectable marker, contrasts with the frequent transcriptional extinction observed with other retroviral vectors carrying two genes. These findings have potentially important implications for clinical bone marrow transplantation and gene therapy of the hematopoietic system.