A database and tools for 3-D protein structure comparison and alignment using the Combinatorial Extension (CE) algorithm

The database reported here is derived using the Combinatorial Extension (CE) algorithm which compares pairs of protein polypeptide chains and provides a list of structurally similar proteins along with their structure alignments. Using CE, structure–structure alignments can provide insights into biological function. When a protein of known function is shown to be structurally similar to a protein of unknown function, a relationship might be inferred; a relationship not necessarily detectable from sequence comparison alone. Establishing structure–structure relationships in this way is of great importance as we enter an era of structural genomics where there is a likelihood of an increasing number of structures with unknown functions being determined. Thus the CE database is an example of a useful tool in the annotation of protein structures of unknown function. Comparisons can be performed on the complete PDB or on a structurally representative subset of proteins. The source protein(s) can be from the PDB (updated monthly) or uploaded by the user. CE provides sequence alignments resulting from structural alignments and Cartesian coordinates for the aligned structures, which may be analyzed using the supplied Compare3D Java applet, or downloaded for further local analysis. Searches can be run from the CE web site, http://cl.sdsc.edu/ce.html, or the database and software downloaded from the site for local use.

Delineating developmental and metabolic pathways in vivo by expression profiling using the RIKEN set of 18,816 full-length enriched mouse cDNA arrays

We have systematically characterized gene expression patterns in 49 adult and embryonic mouse tissues by using cDNA microarrays with 18,816 mouse cDNAs. Cluster analysis defined sets of genes that were expressed ubiquitously or in similar groups of tissues such as digestive organs and muscle. Clustering of expression profiles was observed in embryonic brain, postnatal cerebellum, and adult olfactory bulb, reflecting similarities in neurogenesis and remodeling. Finally, clustering genes coding for known enzymes into 78 metabolic pathways revealed a surprising coordination of expression within each pathway among different tissues. On the other hand, a more detailed examination of glycolysis revealed tissue-specific differences in profiles of key regulatory enzymes. Thus, by surveying global gene expression by using microarrays with a large number of elements, we provide insights into the commonality and diversity of pathways responsible for the development and maintenance of the mammalian body plan.

Roles for genomic imprinting and the zygotic genome in placental development

The placenta contains several types of feto-maternal interfaces where zygote-derived cells interact with maternal cells or maternal blood for the promotion of fetal growth and viability. The genetic factors regulating the interactions between different cell types within feto-maternal interfaces and the relative contributions of the maternal and zygotic genomes are poorly understood. Genomic imprinting, the epigenetic process responsible for parental origin-dependent functional differences between homologous chromosomes, has been proposed to contribute to these events. Previous studies showed that mouse conceptuses with an absence of imprinted differences between the two copies of chromosome 12 (upon paternal inheritance of both copies) die late in gestation and have a variety of defects, including placentomegaly. Here we examined the role of chromosome 12 imprinting in these placentae in more detail. We show that the spatial interactions between different cell types within feto-maternal interfaces are defective and identify abnormal behaviors in both zygote-derived and maternal cells that are attributed to the genome of the zygote but not the mother. These include compromised invasion of the maternal decidualized endometrium and the central maternal artery situated within it by zygote-derived trophoblast, abnormalities in the wall of the central maternal artery, and defects within the zygote-derived cellular layer of the labyrinth, which is in direct contact with maternal blood. These findings demonstrate multiple roles for chromosome 12 imprinting in the placenta that have not previously been associated with imprinting effects. They provide insights into the function of imprinting in placental development and have evolutionary and clinical implications.

Characterization of a high-voltage-activated IA current with a role in spike timing and locomotor pattern generation

Transient A-type K+ channels (IA) in neurons have been implicated in the delay of the spike onset and the decrease in the firing frequency. Here we have characterized biophysically and pharmacologically an IA current in lamprey locomotor network neurons that is activated by suprathreshold depolarization and is specifically blocked by catechol at 100 μM. The biophysical properties of this current are similar to the mammalian Kv3.4 channel. The role of the IA current both in single neuron firing and in locomotor pattern generation was analyzed. The IA current facilitates Na+ channel recovery from inactivation and thus sustains repetitive firing. The role of the IA current in motor pattern generation was examined by applying catechol during fictive locomotion induced by N-methyl-d-aspartate. Blockade of this current increased the locomotor burst frequency and decreased the firing of motoneurons. Although an alternating motor pattern could still be generated, the cycle duration was less regular, with ventral roots bursts failing on some cycles. Our results thus provide insights into the contribution of a high-voltage-activated IA current to the regulation of firing properties and motor coordination in the lamprey spinal cord.

Molecular mechanisms of cocaine reward: Combined dopamine and serotonin transporter knockouts eliminate cocaine place preference

Cocaine blocks uptake by neuronal plasma membrane transporters for dopamine (DAT), serotonin (SERT), and norepinephrine (NET). Cocaine reward/reinforcement has been linked to actions at DAT or to blockade of SERT. However, knockouts of neither DAT, SERT, or NET reduce cocaine reward/reinforcement, leaving substantial uncertainty about cocaine's molecular mechanisms for reward. Conceivably, the molecular bases of cocaine reward might display sufficient redundancy that either DAT or SERT might be able to mediate cocaine reward in the other's absence. To test this hypothesis, we examined double knockout mice with deletions of one or both copies of both the DAT and SERT genes. These mice display viability, weight gain, histologic features, neurochemical parameters, and baseline behavioral features that allow tests of cocaine influences. Mice with even a single wild-type DAT gene copy and no SERT copies retain cocaine reward/reinforcement, as measured by conditioned place-preference testing. However, mice with no DAT and either no or one SERT gene copy display no preference for places where they have previously received cocaine. The serotonin dependence of cocaine reward in DAT knockout mice is thus confirmed by the elimination of cocaine place preference in DAT/SERT double knockout mice. These results provide insights into the brain molecular targets necessary for cocaine reward in knockout mice that develop in their absence and suggest novel strategies for anticocaine medication development.

A demographic approach to selection

The concepts of demography provide a means of combining the ecological approach to population growth with the genetical approach to natural selection. We have utilized the demographic theory of natural selection developed by Norton and Charlesworth to analyze life history schedules of births and deaths for populations of genotypes in Drosophila pseudoobscura. Our populations illustrate a stable genetic equilibrium, an unstable genetic equilibrium, and a case of no equilibrium. We have estimated population growth rates and Darwinian fitnesses of the genotypes and have explored the role of population growth in determining natural selection. The age-specific rates of births and deaths provide insights into components of selection. Both viability and fertility are important components in our populations.

Gene genealogies and population variation in plants

Early in the development of plant evolutionary biology, genetic drift, fluctuations in population size, and isolation were identified as critical processes that affect the course of evolution in plant species. Attempts to assess these processes in natural populations became possible only with the development of neutral genetic markers in the 1960s. More recently, the application of historically ordered neutral molecular variation (within the conceptual framework of coalescent theory) has allowed a reevaluation of these microevolutionary processes. Gene genealogies trace the evolutionary relationships among haplotypes (alleles) with populations. Processes such as selection, fluctuation in population size, and population substructuring affect the geographical and genealogical relationships among these alleles. Therefore, examination of these genealogical data can provide insights into the evolutionary history of a species. For example, studies of Arabidopsis thaliana have suggested that this species underwent rapid expansion, with populations showing little genetic differentiation. The new discipline of phylogeography examines the distribution of allele genealogies in an explicit geographical context. Phylogeographic studies of plants have documented the recolonization of European tree species from refugia subsequent to Pleistocene glaciation, and such studies have been instructive in understanding the origin and domestication of the crop cassava. Currently, several technical limitations hinder the widespread application of a genealogical approach to plant evolutionary studies. However, as these technical issues are solved, a genealogical approach holds great promise for understanding these previously elusive processes in plant evolution.

Direct Evidence for Rapid Degradation of Bacillus thuringiensis Toxin mRNA as a Cause of Poor Expression in Plants1

It is well established that the expression of Bacillus thuringiensis (B.t.) toxin genes in higher plants is severely limited at the mRNA level, but the cause remains controversial. Elucidating whether mRNA accumulation is limited transcriptionally or posttranscriptionally could contribute to effective gene design as well as provide insights about endogenous plant gene-expression mechanisms. To resolve this controversy, we compared the expression of an A/U-rich wild-type cryIA(c) gene and a G/C-rich synthetic cryIA(c) B.t.-toxin gene under the control of identical 5′ and 3′ flanking sequences. Transcriptional activities of the genes were equal as determined by nuclear run-on transcription assays. In contrast, mRNA half-life measurements demonstrated directly that the wild-type transcript was markedly less stable than that encoded by the synthetic gene. Sequences that limit mRNA accumulation were located at more than one site within the coding region, and some appeared to be recognized in Arabidopsis but not in tobacco (Nicotiana tabacum). These results support previous observations that some A/U-rich sequences can contribute to mRNA instability in plants. Our studies further indicate that some of these sequences may be differentially recognized in tobacco cells and Arabidopsis.

BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} acts as a phosphate analog in proteins phosphorylated on aspartate: Structure of a BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} complex with phosphoserine phosphatase

Protein phosphoaspartate bonds play a variety of roles. In response regulator proteins of two-component signal transduction systems, phosphorylation of an aspartate residue is coupled to a change from an inactive to an active conformation. In phosphatases and mutases of the haloacid dehalogenase (HAD) superfamily, phosphoaspartate serves as an intermediate in phosphotransfer reactions, and in P-type ATPases, also members of the HAD family, it serves in the conversion of chemical energy to ion gradients. In each case, lability of the phosphoaspartate linkage has hampered a detailed study of the phosphorylated form. For response regulators, this difficulty was recently overcome with a phosphate analog, BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document}, which yields persistent complexes with the active site aspartate of their receiver domains. We now extend the application of this analog to a HAD superfamily member by solving at 1.5-Å resolution the x-ray crystal structure of the complex of BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} with phosphoserine phosphatase (PSP) from Methanococcus jannaschii. The structure is comparable to that of a phosphoenzyme intermediate: BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} is bound to Asp-11 with the tetrahedral geometry of a phosphoryl group, is coordinated to Mg2+, and is bound to residues surrounding the active site that are conserved in the HAD superfamily. Comparison of the active sites of BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document}⋅PSP and BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document}⋅CeY, a receiver domain/response regulator, reveals striking similarities that provide insights into the function not only of PSP but also of P-type ATPases. Our results indicate that use of BeF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} \begin{equation*}{\mathrm{_{3}^{-}}}\end{equation*}\end{document} for structural studies of proteins that form phosphoaspartate linkages will extend well beyond response regulators.

Conserved promoter elements in the CYP6B gene family suggest common ancestry for cytochrome P450 monooxygenases mediating furanocoumarin detoxification.

Despite the fact that Papilio glaucus and Papilio polyxenes share no single hostplant species, both species feed to varying extents on hostplants that contain furanocoumarins. P. glaucus contains two nearly identical genes, CYP6B4v2 and CYP6B5v1, and P. polyxenes contains two related genes, CYP6B1v3 and CYP6B3v2. Except for CYP6B3v2, the substrate specificity of which has not yet been defined, each of the encoded cytochrome P450 monooxygenases (P450s) metabolizes an array of linear furanocoumarins. All four genes are transcriptionally induced in larvae by exposure to the furanocoumarin xanthotoxin; several are also induced by other furanocoumarins. Comparisons of the organizational structures of these genes indicate that all have the same intron/exon arrangement. Sequences in the promoter regions of the P. glaucus CYP6B4v2/CYP6B5v1 genes and the P. polyxenes CYP6B3v2 gene are similar but not identical to the -146 to -97 region of CYP6B1v3 gene, which contains a xanthotoxin-responsive element (XRE-xan) important for basal and xanthotoxin-inducible transcription of CYP6B1v3. Complements of the xenobiotic-responsive element (XRE-AhR) in the dioxin-inducible human and rat CYP1A1 genes also exist in all four promoters, suggesting that these genes may be regulated by dioxin. Antioxidant-responsive elements (AREs) in mouse and rat glutathione S-transferase genes and the Barbie box element (Bar) in the bacterial CYP102 gene exist in the CYP6B1v3, CYP6B4v2, and CYP6B5v1 promoters. Similarities in the protein sequences, intron positions, and xanthotoxin- and xenobiotic-responsive promoter elements indicate that these insect CYP6B genes are derived from a common ancestral gene. Evolutionary comparisons between these P450 genes are the first available for a group of insect genes transcriptionally regulated by hostplant allelochemicals and provide insights into the process by which insects evolve specialized feeding habits.

The ocular albinism type 1 gene product is a membrane glycoprotein localized to melanosomes.

Ocular albinism type 1 (OA1) is an inherited disorder characterized by severe reduction of visual acuity, photophobia, and retinal hypopigmentation. Ultrastructural examination of skin melanocytes and of the retinal pigment epithelium reveals the presence of macromelanosomes, suggesting a defect in melanosome biogenesis. The gene responsible for OA1 is exclusively expressed in pigment cells and encodes a predicted protein of 404 aa displaying several putative transmembrane domains and sharing no similarities with previously identified molecules. Using polyclonal antibodies we have identified the endogenous OA1 protein in retinal pigment epithelial cells, in normal human melanocytes and in various melanoma cell lines. Two forms of the OA1 protein were identified by Western analysis, a 60-kDa glycoprotein and a doublet of 48 and 45 kDa probably corresponding to unglycosylated precursor polypeptides. Upon subcellular fractionation and phase separation with the nonionic detergent Triton X-114, the OA1 protein segregated into the melanosome-rich fraction and behaved as an authentic integral membrane protein. Immunofluorescence and immunogold analyses on normal human melanocytes confirmed the melanosomal membrane localization of the endogenous OA1 protein, consistent with its possible involvement in melanosome biogenesis. The identification of a novel melanosomal membrane protein involved in a human disease will provide insights into the mechanisms that control the cell-specific pathways of subcellular morphogenesis.

A molecular mechanism for the effect of lithium on development.

Lithium, one of the most effective drugs for the treatment of bipolar (manic-depressive) disorder, also has dramatic effects on morphogenesis in the early development of numerous organisms. How lithium exerts these diverse effects is unclear, but the favored hypothesis is that lithium acts through inhibition of inositol monophosphatase (IMPase). We show here that complete inhibition of IMPase has no effect on the morphogenesis of Xenopus embryos and present a different hypothesis to explain the broad action of lithium. Our results suggest that lithium acts through inhibition of glycogen synthase kinase-3 beta (GSK-3 beta), which regulates cell fate determination in diverse organisms including Dictyostelium, Drosophila, and Xenopus. Lithium potently inhibits GSK-3 beta activity (Ki = 2 mM), but is not a general inhibitor of other protein kinases. In support of this hypothesis, lithium treatment phenocopies loss of GSK-3 beta function in Xenopus and Dictyostelium. These observations help explain the effect of lithium on cell-fate determination and could provide insights into the pathogenesis and treatment of bipolar disorder.

Long-circulating bacteriophage as antibacterial agents.

The increased prevalence of multidrug-resistant bacterial pathogens motivated us to attempt to enhance the therapeutic efficacy of bacteriophages. The therapeutic application of phages as antibacterial agents was impeded by several factors: (i) the failure to recognize the relatively narrow host range of phages; (ii) the presence of toxins in crude phage lysates; and (iii) a lack of appreciation for the capacity of mammalian host defense systems, particularly the organs of the reticuloendothelial system, to remove phage particles from the circulatory system. In our studies involving bacteremic mice, the problem of the narrow host range of phage was dealt with by using selected bacterial strains and virulent phage specific for them. Toxin levels were diminished by purifying phage preparations. To reduce phage elimination by the host defense system, we developed a serial-passage technique in mice to select for phage mutants able to remain in the circulatory system for longer periods of time. By this approach we isolated long-circulating mutants of Escherichia coli phage lambda and of Salmonella typhimurium phage P22. We demonstrated that the long-circulating lambda mutants also have greater capability as antibacterial agents than the corresponding parental strain in animals infected with lethal doses of bacteria. Comparison of the parental and mutant lambda capsid proteins revealed that the relevant mutation altered the major phage head protein E. The use of toxin-free, bacteria-specific phage strains, combined with the serial-passage technique, may provide insights for developing phage into therapeutically effective antibacterial agents.

Comparative interleukin (IL-2)/interferon IFN-gamma and IL-4/IL-10 responses during acute infection of macaques inoculated with attenuated nef-truncated or pathogenic SICmac251 virus.

Comparison of immune responses to infection by a pathogenic or a nonpathogenic immunodeficiency virus in macaques may provide insights into pathogenetic events leading to simian AIDS. This work is aimed at exploring cytokine expression during infection by simian immunodeficiency virus (SIV). We used semiquantitative reverse transcription-PCR to monitor interleukin (IL)-2/interferon (IFN)-gamma (Th1-like), and IL-4/IL-10 (Th2-like) expression in unmanipulated peripheral blood mononuclear cells (PBMCs), during the acute phase of infection of eight cynomolgus macaques (Macaca fascicularis) with a pathogenic primary isolate of SIVmac251 (full-length nef), and of four other cynomolgus macaques by an attenuated molecular clone of SIVmac251 (nef-truncated). All the monkeys became infected, as clearly shown by the presence of infected PBMCs and by seroconversion. Nevertheless, PBMC-associated virus loads and p27 antigenemia in monkeys infected by the attenuated virus clone remained lower than those observed in animals infected with the pathogenic SIVmac251 isolate. A rise of IL-10 mRNA expression occurred in both groups of monkeys coincident with the peak of viral replication. In monkeys infected with the pathogenic SIVmac251, IL-2, IL-4, and IFN-gamma mRNAs were either weakly detectable or undetectable. On the contrary, animals infected by the attenuated virus exhibited an overexpression of these cytokine mRNAs during the first weeks after inoculation. The lack of expression of these cytokines in monkeys infected with the pathogenic primary isolate may reflect early immunodeficiency.

Three-dimensional structure of recombinant human osteogenic protein 1: structural paradigm for the transforming growth factor beta superfamily.

We report the three-dimensional structure of osteogenic protein 1 (OP-1, also known as bone morphogenetic protein 7) to 2.8-A resolution. OP-1 is a member of the transforming growth factor beta (TGF-beta) superfamily of proteins and is able to induce new bone formation in vivo. Members of this superfamily share sequence similarity in their C-terminal regions and are implicated in embryonic development and adult tissue repair. Our crystal structure makes possible the structural comparison between two members of the TGF-beta superfamily. We find that although there is limited sequence identity between OP-1 and TGF-beta 2, they share a common polypeptide fold. These results establish a basis for proposing the OP-1/TGF-beta 2 fold as the primary structural motif for the TGF-beta superfamily as a whole. Detailed comparison of the OP-1 and TGF-beta 2 structures has revealed striking differences that provide insights into how these growth factors interact with their receptors.