The Mouse Gene Expression Database (GXD)

The Gene Expression Database (GXD) is a community resource of gene expression information for the laboratory mouse. By combining the different types of expression data, GXD aims to provide increasingly complete information about the expression profiles of genes in different mouse strains and mutants, thus enabling valuable insights into the molecular networks that underlie normal development and disease. GXD is integrated with the Mouse Genome Database (MGD). Extensive interconnections with sequence databases and with databases from other species, and the development and use of shared controlled vocabularies extend GXD’s utility for the analysis of gene expression information. GXD is accessible through the Mouse Genome Informatics web site at http://www.informatic s.jax.org/ or directly at http://www.informatics.jax.org/me nus/expression_menu.shtml.

Centromere mapping and orientation of the molecular linkage map of rice (Oryza sativa L.).

Rice has become a model cereal plant for molecular genetic research. Rice has the most comprehensive molecular linkage maps with more than 2000 DNA markers and shows synteny and colinearity with the maps of other cereal crops. Until now, however, no information was available about the positions of centromeres and arm locations of markers on the molecular linkage map. Secondary and telotrisomics were used to assign restriction fragment length polymorphism markers to specific chromosome arms and thereby to map the positions of centromeres. More than 170 restriction fragment length polymorphism markers were assigned to specific chromosome arms through gene dosage analysis using the secondary and telotrisomics and the centromere positions were mapped on all 12 linkage groups. The orientations of seven linkage groups were reversed to fit the "short arm on top" convention and the corrected map is presented.

Mapping the active site in vasoactive intestinal peptide to a core of four amino acids: Neuroprotective drug design

The understanding of the molecular mechanisms leading to peptide action entails the identification of a core active site. The major 28-aa neuropeptide, vasoactive intestinal peptide (VIP), provides neuroprotection. A lipophilic derivative with a stearyl moiety at the N-terminal and norleucine residue replacing the Met-17 was 100-fold more potent than VIP in promoting neuronal survival, acting at femtomolar–picomolar concentration. To identify the active site in VIP, over 50 related fragments containing an N-terminal stearic acid attachment and an amidated C terminus were designed, synthesized, and tested for neuroprotective properties. Stearyl-Lys-Lys-Tyr-Leu-NH2 (derived from the C terminus of VIP and the related peptide, pituitary adenylate cyclase activating peptide) captured the neurotrophic effects offered by the entire 28-aa parent lipophilic derivative and protected against β-amyloid toxicity in vitro. Furthermore, the 4-aa lipophilic peptide recognized VIP-binding sites and enhanced choline acetyltransferase activity as well as cognitive functions in Alzheimer’s disease-related in vivo models. Biodistribution studies following intranasal administration of radiolabeled peptide demonstrated intact peptide in the brain 30 min after administration. Thus, lipophilic peptide fragments offer bioavailability and stability, providing lead compounds for drug design against neurodegenerative diseases.

The underlying pathway structure of biochemical reaction networks

Bioinformatics is yielding extensive, and in some cases complete, genetic and biochemical information about individual cell types and cellular processes, providing the composition of living cells and the molecular structure of its components. These components together perform integrated cellular functions that now need to be analyzed. In particular, the functional definition of biochemical pathways and their role in the context of the whole cell is lacking. In this study, we show how the mass balance constraints that govern the function of biochemical reaction networks lead to the translation of this problem into the realm of linear algebra. The functional capabilities of biochemical reaction networks, and thus the choices that cells can make, are reflected in the null space of their stoichiometric matrix. The null space is spanned by a finite number of basis vectors. We present an algorithm for the synthesis of a set of basis vectors for spanning the null space of the stoichiometric matrix, in which these basis vectors represent the underlying biochemical pathways that are fundamental to the corresponding biochemical reaction network. In other words, all possible flux distributions achievable by a defined set of biochemical reactions are represented by a linear combination of these basis pathways. These basis pathways thus represent the underlying pathway structure of the defined biochemical reaction network. This development is significant from a fundamental and conceptual standpoint because it yields a holistic definition of biochemical pathways in contrast to definitions that have arisen from the historical development of our knowledge about biochemical processes. Additionally, this new conceptual framework will be important in defining, characterizing, and studying biochemical pathways from the rapidly growing information on cellular function.

Quantitative modeling of stochastic systems in molecular biology by using stochastic Petri nets

An integrated understanding of molecular and developmental biology must consider the large number of molecular species involved and the low concentrations of many species in vivo. Quantitative stochastic models of molecular interaction networks can be expressed as stochastic Petri nets (SPNs), a mathematical formalism developed in computer science. Existing software can be used to define molecular interaction networks as SPNs and solve such models for the probability distributions of molecular species. This approach allows biologists to focus on the content of models and their interpretation, rather than their implementation. The standardized format of SPNs also facilitates the replication, extension, and transfer of models between researchers. A simple chemical system is presented to demonstrate the link between stochastic models of molecular interactions and SPNs. The approach is illustrated with examples of models of genetic and biochemical phenomena where the UltraSAN package is used to present results from numerical analysis and the outcome of simulations.

Peptide design by artificial neural networks and computer-based evolutionary search

A technique for systematic peptide variation by a combination of rational and evolutionary approaches is presented. The design scheme consists of five consecutive steps: (i) identification of a “seed peptide” with a desired activity, (ii) generation of variants selected from a physicochemical space around the seed peptide, (iii) synthesis and testing of this biased library, (iv) modeling of a quantitative sequence-activity relationship by an artificial neural network, and (v) de novo design by a computer-based evolutionary search in sequence space using the trained neural network as the fitness function. This strategy was successfully applied to the identification of novel peptides that fully prevent the positive chronotropic effect of anti-β1-adrenoreceptor autoantibodies from the serum of patients with dilated cardiomyopathy. The seed peptide, comprising 10 residues, was derived by epitope mapping from an extracellular loop of human β1-adrenoreceptor. A set of 90 peptides was synthesized and tested to provide training data for neural network development. De novo design revealed peptides with desired activities that do not match the seed peptide sequence. These results demonstrate that computer-based evolutionary searches can generate novel peptides with substantial biological activity.

Self-assembly of fibronectin into fibrillar networks underneath dipalmitoyl phosphatidylcholine monolayers: Role of lipid matrix and tensile forces

The cell-mediated assembly of fibronectin (Fn) into fibrillar matrices is a complex multistep process that is incompletely understood because of the chemical complexity of the extracellular matrix and a lack of experimental control over molecular interactions and dynamic events. We have identified conditions under which Fn assembles into extended fibrillar networks after adsorption to a dipalmitoyl phosphatidylcholine (DPPC) monolayer in contact with physiological buffer. We propose a sequential model for the Fn assembly pathway, which involves the orientation of Fn underneath the lipid monolayer by insertion into the liquid expanded (LE) phase of DPPC. Attractive interactions between these surface-anchored proteins and the liquid condensed (LC) domains leads to Fn enrichment at domain edges. Spontaneous self-assembly into fibrillar networks, however, occurs only after expansion of the DPPC monolayer from the LC phase though the LC/LE phase coexistence. Upon monolayer expansion, the domain boundaries move apart while attractive interactions among Fn molecules and between Fn and domain edges produce a tensile force on the proteins that initiates fibril assembly. The resulting fibrils have been characterized in situ by using fluorescence and light-scattering microscopy. We have found striking similarities between fibrils produced under DPPC monolayers and those found on cellular surfaces, including their assembly pathways.

Molecular dynamics of MHC genesis unraveled by sequence analysis of the 1,796,938-bp HLA class I region

The intensely studied MHC has become the paradigm for understanding the architectural evolution of vertebrate multigene families. The 4-Mb human MHC (also known as the HLA complex) encodes genes critically involved in the immune response, graft rejection, and disease susceptibility. Here we report the continuous 1,796,938-bp genomic sequence of the HLA class I region, linking genes between MICB and HLA-F. A total of 127 genes or potentially coding sequences were recognized within the analyzed sequence, establishing a high gene density of one per every 14.1 kb. The identification of 758 microsatellite provides tools for high-resolution mapping of HLA class I-associated disease genes. Most importantly, we establish that the repeated duplication and subsequent diversification of a minimal building block, MIC-HCGIX-3.8–1-P5-HCGIV-HLA class I-HCGII, engendered the present-day MHC. That the currently nonessential HLA-F and MICE genes have acted as progenitors to today’s immune-competent HLA-ABC and MICA/B genes provides experimental evidence for evolution by “birth and death,” which has general relevance to our understanding of the evolutionary forces driving vertebrate multigene families.

Mapping regions containing binding residues within functional domains of Plasmodium vivax and Plasmodium knowlesi erythrocyte-binding proteins

Invasion of erythrocytes by malaria parasites is mediated by specific molecular interactions. Whereas Plasmodium vivax and Plasmodium knowlesi use the Duffy blood group antigen, Plasmodium falciparum uses sialic acid residues of glycophorin A as receptors to invade human erythrocytes. P. knowlesi uses the Duffy antigen as well as other receptors to invade rhesus erythrocytes by multiple pathways. Parasite ligands that bind these receptors belong to a family of erythrocyte-binding proteins (EBP). The EBP family includes the P. vivax and P. knowlesi Duffy-binding proteins, P. knowlesi β and γ proteins, which bind alternate receptors on rhesus erythrocytes, and P. falciparum erythrocyte-binding antigen (EBA-175), which binds sialic acid residues of human glycophorin A. Binding domains of each EBP lie in a conserved N-terminal cysteine-rich region, region II, which contains around 330 amino acids with 12 to 14 conserved cysteines. Regions containing binding residues have now been mapped within P. vivax and P. knowlesi β region II. Chimeric domains containing P. vivax region II sequences fused to P. knowlesi β region II sequences were expressed on the surface of COS cells and tested for binding to erythrocytes. Binding residues of P. vivax region II lie in a 170-aa stretch between cysteines 4 and 7, and binding residues of P. knowlesi β region II lie in a 53-aa stretch between cysteines 4 and 5. Mapping regions responsible for receptor recognition is an important step toward understanding the structural basis for the interaction of these parasite ligands with host receptors.

Topoisomerase II-mediated site-directed alkylation of DNA by psorospermin and its use in mapping other topoisomerase II poison binding sites

Psorospermin is a plant natural product that shows significant in vivo activity against P388 mouse leukemia. The molecular basis for this selectivity is unknown, although psorospermin has been demonstrated to intercalate into DNA and alkylate N7 of guanine. Significantly, the alkylation reactivity of psorospermin at specific sites on DNA increased 25-fold in the presence of topoisomerase II. In addition, psorospermin trapped the topoisomerase II-cleaved complex formation at the same site. These results imply that the efficacy of psorospermin is related to its interaction with the topoisomerase II–DNA complex. Because thermal treatment of (N7 guanine)–DNA adducts leads to DNA strand breakage, we were able to determine the site of alkylation of psorospermin within the topoisomerase II gate site and infer that intercalation takes place at the gate site between base pairs at the +1 and +2 positions. These results provide not only additional mechanistic information on the mode of action of the anticancer agent psorospermin but also structural insights into the design of an additional class of topoisomerase II poisons. Because the alkylation site for psorospermin in the presence of topoisomerase II can be assigned unambiguously and the intercalation site inferred, this drug is a useful probe for other topoisomerase poisons where the sites for interaction are less well defined.

Comparative mapping of Andropogoneae: Saccharum L. (sugarcane) and its relation to sorghum and maize

Comparative genetic maps of Papuan Saccharum officinarum L. (2n = 80) and S. robustum (2n = 80) were constructed by using single-dose DNA markers (SDMs). SDM-framework maps of S. officinarum and S. robustum were compared with genetic maps of sorghum and maize by way of anchor restriction fragment length polymorphism probes. The resulting comparisons showed striking colinearity between the sorghum and Saccharum genomes. There were no differences in marker order between S. officinarum and sorghum. Furthermore, there were no alterations in SDM order between S. officinarum and S. robustum. The S. officinarum and S. robustum maps also were compared with the map of the polysomic octoploid S. spontaneum ‘SES 208’ (2n = 64, x = 8), thus permitting relations to homology groups (“chromosomes”) of S. spontaneum to be studied. Investigation of transmission genetics in S. officinarum and S. robustum confirmed preliminary results that showed incomplete polysomy in these species. Because of incomplete polysomy, multiple-dose markers could not be mapped for lack of a genetic model for their segregation. To coalesce S. officinarum and S. robustum linkage groups into homology groups (composed of homologous pairing partners), they were compared with sorghum (2n = 20), which functioned as a synthetic diploid. Groupings suggested by comparative mapping were found to be highly concordant with groupings based on highly polymorphic restriction fragment length polymorphism probes detecting multiple SDMs. The resulting comparative maps serve as bridges to allow information from one Andropogoneae to be used by another, for breeding, ecology, evolution, and molecular biology.

Molecular Interaction Map of the Mammalian Cell Cycle Control and DNA Repair Systems

Eventually to understand the integrated function of the cell cycle regulatory network, we must organize the known interactions in the form of a diagram, map, and/or database. A diagram convention was designed capable of unambiguous representation of networks containing multiprotein complexes, protein modifications, and enzymes that are substrates of other enzymes. To facilitate linkage to a database, each molecular species is symbolically represented only once in each diagram. Molecular species can be located on the map by means of indexed grid coordinates. Each interaction is referenced to an annotation list where pertinent information and references can be found. Parts of the network are grouped into functional subsystems. The map shows how multiprotein complexes could assemble and function at gene promoter sites and at sites of DNA damage. It also portrays the richness of connections between the p53-Mdm2 subsystem and other parts of the network.

Molecular cloning and characterization of a murine pre-B-cell growth-stimulating factor/stromal cell-derived factor 1 receptor, a murine homolog of the human immunodeficiency virus 1 entry coreceptor fusin

Pre-B-cell growth-stimulating factor/stromal cell-derived factor 1 (PBSF/SDF-1) is a member of the CXC group of chemokines that is initially identified as a bone marrow stromal cell-derived factor and as a pre-B-cell stimulatory factor. Although most chemokines are thought to be inducible inflammatory mediators, PBSF/SDF-1 is essential for perinatal viability, B lymphopoiesis, bone marrow myelopoiesis, and cardiac ventricular septal formation, and it has chemotactic activities on resting lymphocytes and monocytes. In this paper, we have isolated a cDNA that encodes a seven transmembrane-spanning-domain receptor, designated pre-B-cell-derived chemokine receptor (PB-CKR) from a murine pre-B-cell clone, DW34. The deduced amino acid sequence has 90% identity with that of a HUMSTSR/fusin, a human immunodeficiency virus 1 (HIV-1) entry coreceptor. However, the second extracellular region has lower identity (67%) compared with HUMSTSR/fusin. PB-CKR is expressed during embryo genesis and in many organs and T cells of adult mice. Murine PBSF/SDF-1 induced an increase in intracellular free Ca2+ in DW34 cells and PB-CKR-transfected Chinese hamster ovary (CHO) cells, suggesting that PB-CKR is a functional receptor for murine PBSF/SDF-1. Murine PBSF/SDF-1 also induced Ca2+ influx in fusin-transfected CHO cells. On the other hand, considering previous results that HIV-1 does not enter murine T cells that expressed human CD4, PB-CKR may not support HIV-1 infection. Thus, PB-CKR will be an important tool for functional mapping of HIV-1 entry coreceptor fusin and for understanding the function of PBSF/SDF-1 further.

High-resolution mapping and isolation of a yeast artificial chromosome contig containing fw2.2: A major fruit weight quantitative trait locus in tomato

A high-resolution physical and genetic map of a major fruit weight quantitative trait locus (QTL), fw2.2, has been constructed for a region of tomato chromosome 2. Using an F2 nearly isogenic line mapping population (3472 individuals) derived from Lycopersicon esculentum (domesticated tomato) × Lycopersicon pennellii (wild tomato), fw2.2 has been placed near TG91 and TG167, which have an interval distance of 0.13 ± 0.03 centimorgan. The physical distance between TG91 and TG167 was estimated to be ≤ 150 kb by pulsed-field gel electrophoresis of tomato DNA. A physical contig composed of six yeast artificial chromosomes (YACs) and encompassing fw2.2 was isolated. No rearrangements or chimerisms were detected within the YAC contig based on restriction fragment length polymorphism analysis using YAC-end sequences and anchored molecular markers from the high-resolution map. Based on genetic recombination events, fw2.2 could be narrowed down to a region less than 150 kb between molecular markers TG91 and HSF24 and included within two YACs: YAC264 (210 kb) and YAC355 (300 kb). This marks the first time, to our knowledge, that a QTL has been mapped with such precision and delimited to a segment of cloned DNA. The fact that the phenotypic effect of the fw2.2 QTL can be mapped to a small interval suggests that the action of this QTL is likely due to a single gene. The development of the high-resolution genetic map, in combination with the physical YAC contig, suggests that the gene responsible for this QTL and other QTLs in plants can be isolated using a positional cloning strategy. The cloning of fw2.2 will likely lead to a better understanding of the molecular biology of fruit development and to the genetic engineering of fruit size characteristics.

Molecular cytogenetic dissection of human chromosomes 3 and 21 evolution

Chromosome painting in placental mammalians illustrates that genome evolution is marked by chromosomal synteny conservation and that the association of chromosomes 3 and 21 may be the largest widely conserved syntenic block known for mammals. We studied intrachromosomal rearrangements of the syntenic block 3/21 by using probes derived from chromosomal subregions with a resolution of up to 10–15 Mbp. We demonstrate that the rearrangements visualized by chromosome painting, mostly translocations, are only a fraction of the actual chromosomal changes that have occurred during evolution. The ancestral segment order for both primates and carnivores is still found in some species in both orders. From the ancestral primate/carnivore condition an inversion is needed to derive the pig homolog, and a fission of chromosome 21 and a pericentric inversion is needed to derive the Bornean orangutan condition. Two overlapping inversions in the chromosome 3 homolog then would lead to the chromosome form found in humans and African apes. This reconstruction of the origin of human chromosome 3 contrasts with the generally accepted scenario derived from chromosome banding in which it was proposed that only one pericentric inversion was needed. From the ancestral form for Old World primates (now found in the Bornean orangutan) a pericentric inversion and centromere shift leads to the chromosome ancestral for all Old World monkeys. Intrachromosomal rearrangements, as shown here, make up a set of potentially plentiful and informative markers that can be used for phylogenetic reconstruction and a more refined comparative mapping of the genome.