On signal sequence polymorphisms and diseases of distribution.

We report a previously unappreciated property of the signals that target organelle-specific proteins to their subcellular sites of action. Such targeting sequences are shown to be polymorphic. We discovered this polymorphism when we cloned the mitochondrial manganese-containing superoxide dismutase from cell lines of normal individuals and patients with genetic diseases of premature aging and compared their sequences to each other and to those previously reported. The polymorphism consists of a single nucleotide change in the region of the DNA that encodes the signal sequence such that either an alanine or valine is present. Subsequently, eight cell lines were analyzed and all three possible combinations of the two signal sequences were observed. Such signal sequence polymorphisms could result in diseases of distribution, where essential proteins are not properly targeted, thereby leading to absolute or relative deficiencies of critical enzymes within specific cellular compartments. Progeria and related syndromes may be diseases of distribution.

Sequence motifs in adenoviral DNA block immune activation by stimulatory CpG motifs

Unmethylated CpG dinucleotides in particular base contexts (CpG-S motifs) are relatively common in bacterial DNA but are rare in vertebrate DNA. B cells and monocytes have the ability to detect such CpG-S motifs that trigger innate immune defenses with production of Th1-like cytokines. Despite comparable levels of unmethylated CpG dinucleotides, DNA from serotype 12 adenovirus is immune-stimulatory, but serotype 2 is nonstimulatory and can even inhibit activation by bacterial DNA. In type 12 genomes, the distribution of CpG-flanking bases is similar to that predicted by chance. However, in type 2 adenoviral DNA the immune stimulatory CpG-S motifs are outnumbered by a 15- to 30-fold excess of CpG dinucleotides in clusters of direct repeats or with a C on the 5′ side or a G on the 3′ side. Synthetic oligodeoxynucleotides containing these putative neutralizing (CpG-N) motifs block immune activation by CpG-S motifs in vitro and in vivo. Eliminating 52 of the 134 CpG-N motifs present in a DNA vaccine markedly enhanced its Th1-like function in vivo, which was increased further by the addition of CpG-S motifs. Thus, depending on the CpG motif, prokaryotic DNA can be either immune-stimulatory or neutralizing. These results have important implications for understanding microbial pathogenesis and molecular evolution and for the clinical development of DNA vaccines and gene therapy vectors.

Crossover isomer bias is the primary sequence-dependent property of immobilized Holliday junctions

Recombination of genes is essential to the evolution of genetic diversity, the segregation of chromosomes during cell division, and certain DNA repair processes. The Holliday junction, a four-arm, four-strand branched DNA crossover structure, is formed as a transient intermediate during genetic recombination and repair processes in the cell. The recognition and subsequent resolution of Holliday junctions into parental or recombined products appear to be critically dependent on their three-dimensional structure. Complementary NMR and time-resolved fluorescence resonance energy transfer experiments on immobilized four-arm DNA junctions reported here indicate that the Holliday junction cannot be viewed as a static structure but rather as an equilibrium mixture of two conformational isomers. Furthermore, the distribution between the two possible crossover isomers was found to depend on the sequence in a manner that was not anticipated on the basis of previous low-resolution experiments.

Imaging of H217O distribution in the brain of a live rat by using proton-detected 17O MRI

Imaging of H217O has a number of important applications. Mapping the distribution of H217O produced by oxidative metabolism of 17O-enriched oxygen gas may lead to a new method of metabolic functional imaging; regional cerebral blood flow also can be measured by measuring the H217O distribution after the injection of 17O-enriched physiological saline solution. Previous studies have proposed a method for indirect detection of 17O. The method is based on the shortening of the proton T2 in H217O solutions, caused by the residual 17O-1H scalar coupling and transferred to the bulk water via fast chemical exchange. It has been shown that the proton T2 of H217O solutions can be restored to that of H216O by irradiating the resonance frequency of the 17O nucleus. The indirect 17O image thus is obtained by taking the difference between two T2-weighted spin-echo images: one acquired after irradiation of the 17O resonance and one acquired without irradiation. It also has been established that, at relatively low concentrations of H217O, the indirect method yields an image that quantitatively reflects the H217O distribution in the sample. The method is referred to as PRIMO (proton imaging of oxygen). In this work, we show in vivo proton images of the H217O distribution in a rat brain after an i.v. injection of H217O-enriched physiological saline solution. Implementing the indirect detection method in an echo-planar imaging sequence enabled obtaining H217O images with good spatial and temporal resolution of few seconds.

Mitochondrial carbonic anhydrase CA VB: Differences in tissue distribution and pattern of evolution from those of CA VA suggest distinct physiological roles

A cDNA for a second mouse mitochondrial carbonic anhydrase (CA) called CA VB was identified by homology to the previously characterized murine CA V, now called CA VA. The full-length cDNA encodes a 317-aa precursor that contains a 33-aa classical mitochondrial leader sequence. Comparison of products expressed from cDNAs for murine CA VB and CA VA in COS cells revealed that both expressed active CAs that localized in mitochondria, and showed comparable activities in crude extracts and in mitochondria isolated from transfected COS cells. Northern blot analyses of total RNAs from mouse tissues and Western blot analyses of mouse tissue homogenates showed differences in tissue-specific expression between CA VB and CA VA. CA VB was readily detected in most tissues, while CA VA expression was limited to liver, skeletal muscle, and kidney. The human orthologue of murine CA VB was recently reported also. Comparison of the CA domain sequence of human CA VB with that reported here shows that the CA domains of CA VB are much more highly conserved between mouse and human (95% identity) than the CA domains of mouse and human CA VAs (78% identity). Analysis of phylogenetic relationships between these and other available human and mouse CA isozyme sequences revealed that mammalian CA VB evolved much more slowly than CA VA, accepting amino acid substitutions at least 4.5 times more slowly since each evolved from its respective human–mouse ancestral gene around 90 million years ago. Both the differences in tissue distribution and the much greater evolutionary constraints on CA VB sequences suggest that CA VB and CA VA have evolved to assume different physiological roles.

Point mutations and sequence variability in proteins: Redistributions of preexisting populations

Here we study the effect of point mutations in proteins on the redistributions of the conformational substates. We show that regardless of the location of a mutation in the protein structure and of its type, the observed movements of the backbone recur largely at the same positions in the structures. Despite the different interactions that are disrupted and formed by the residue substitution, not only are the conformations very similar, but the regions that move are also the same, regardless of their sequential or spatial distance from the mutation. This observation leads us to conclude that, apart from some extreme cases, the details of the interactions are not critically important in determining the protein conformation or in specifying which parts of the protein would be more prone to take on different local conformations in response to changes in the sequence. This finding further illustrates why proteins manifest a robustness toward many mutational events. This nonuniform distribution of the conformer population is consistently observed in a variety of protein structural types. Topology is critically important in determining folding pathways, kinetics, building block cutting, and anatomy trees. Here we show that topology is also very important in determining which regions of the protein structure will respond to sequence changes, regardless of the sequential or spatial location of the mutation.

A Di-Leucine Sequence and a Cluster of Acidic Amino Acids Are Required for Dynamic Retention in the Endosomal Recycling Compartment of Fibroblasts

Insulin-regulated aminopeptidase (IRAP), a transmembrane aminopeptidase, is dynamically retained within the endosomal compartment of fibroblasts. The characteristics of this dynamic retention are rapid internalization from the plasma membrane and slow recycling back to the cell surface. These specialized trafficking kinetics result in <15% of IRAP on the cell surface at steady state, compared with 35% of the transferrin receptor, another transmembrane protein that traffics between endosomes and the cell surface. Here we demonstrate that a 29-amino acid region of IRAP's cytoplasmic domain (residues 56–84) is necessary and sufficient to promote trafficking characteristic of IRAP. A di-leucine sequence and a cluster of acidic amino acids within this region are essential elements of the motif that slows IRAP recycling. Rapid internalization requires any two of three distinct motifs: M15,16, DED64–66, and LL76,77. The DED and LL sequences are part of the motif that regulates recycling, demonstrating that this motif is bifunctional. In this study we used horseradish peroxidase quenching of fluorescence to demonstrate that IRAP is dynamically retained within the transferrin receptor-containing general endosomal recycling compartment. Therefore, our data demonstrate that motifs similar to those that determine targeting among distinct membrane compartments can also regulate the rate of transport of proteins from endosomal compartments. We propose a model for dynamic retention in which IRAP is transported from the general endosomal recycling compartment in specialized, slowly budding recycling vesicles that are distinct from those that mediate rapid recycling back to the surface (e.g., transferrin receptor-containing transport vesicles). It is likely that the dynamic retention of IRAP is an example of a general mechanism for regulating the distribution of proteins between the surface and interior of cells.

A unified statistical framework for sequence comparison and structure comparison

We present an approach for assessing the significance of sequence and structure comparisons by using nearly identical statistical formalisms for both sequence and structure. Doing so involves an all-vs.-all comparison of protein domains [taken here from the Structural Classification of Proteins (scop) database] and then fitting a simple distribution function to the observed scores. By using this distribution, we can attach a statistical significance to each comparison score in the form of a P value, the probability that a better score would occur by chance. As expected, we find that the scores for sequence matching follow an extreme-value distribution. The agreement, moreover, between the P values that we derive from this distribution and those reported by standard programs (e.g., blast and fasta validates our approach. Structure comparison scores also follow an extreme-value distribution when the statistics are expressed in terms of a structural alignment score (essentially the sum of reciprocated distances between aligned atoms minus gap penalties). We find that the traditional metric of structural similarity, the rms deviation in atom positions after fitting aligned atoms, follows a different distribution of scores and does not perform as well as the structural alignment score. Comparison of the sequence and structure statistics for pairs of proteins known to be related distantly shows that structural comparison is able to detect approximately twice as many distant relationships as sequence comparison at the same error rate. The comparison also indicates that there are very few pairs with significant similarity in terms of sequence but not structure whereas many pairs have significant similarity in terms of structure but not sequence.

Occurrence of two somatostatin variants in the frog brain: characterization of the cDNAs, distribution of the mRNAs, and receptor-binding affinities of the peptides.

In tetrapods, only one gene encoding a somatostatin precursor has been identified so far. The present study reports the characterization of the cDNA clones that encode two distinct somatostatin precursors in the brain of the frog Rana ridibunda. The cDNAs were isolated by using degenerate oligonucleotides based on the sequence of the central region of somatostatin to screen a frog brain cDNA library. One of the cDNAs encodes a 115-amino acid protein (prepro-somatostatin-14; PSS1) that exhibits a high degree of structural similarity with the mammalian somatostatin precursor. The other cDNA encodes a 103-amino acid protein (prepro-[Pro2, Met13]somatostatin-14; PSS2) that contains the sequence of the somatostatin analog (peptide SS2) at its C terminus, but does not exhibit appreciable sequence similarity with PSS1 in the remaining region. In situ hybridization studies indicate differential expression of the PSS1 and PSS2 genes in the septum, the lateral part of the pallium, the amygdaloid complex, the posterior nuclei of the thalamus, the ventral hypothalamic nucleus, the torus semicircularis and the optic tectum. The somatostatin variant SS2 was significantly more potent (4-6 fold) than somatostatin itself in displacing [125I-Tyr0, D-Trp8] somatostatin-14 from its specific binding sites. The present study indicates that the two somatostatin variants could exert different functions in the frog brain and pituitary. These data also suggest that distinct genes encoding somatostatin variants may be expressed in the brain of other tetrapods.

Molecular identity and cellular distribution of advanced glycation endproduct receptors: relationship of p60 to OST-48 and p90 to 80K-H membrane proteins.

Advanced glycation endproducts (AGEs) are derivatives of nonenzymatic reactions between sugars and protein or lipids, and together with AGE-specific receptors are involved in numerous pathogenic processes associated with aging and hyperglycemia. Two of the known AGE-binding proteins isolated from rat liver membranes, p60 and p90, have been partially sequenced. We now report that the N-terminal sequence of p60 exhibits 95% identity to OST-48, a 48-kDa member of the oligosaccharyltransferase complex found in microsomal membranes, while sequence analysis of p90 revealed 73% and 85% identity to the N-terminal and internal sequences, respectively, of human 80K-H, a 80- to 87-kDa protein substrate for protein kinase C. AGE-ligand and Western analyses of purified oligosaccharyltransferase complex, enriched rough endoplasmic reticulum, smooth endoplasmic reticulum, and plasma membranes from rat liver or RAW 264.7 macrophages yielded a single protein of approximately 50 kDa recognized by both anti-p60 and anti-OST-48 antibodies, and also exhibited AGE-specific binding. Immunoprecipitated OST-48 from rat rough endoplasmic reticulum fractions exhibited both AGE binding and immunoreactivity to an anti-p60 antibody. Immune IgG raised to recombinant OST-48 and 80K-H inhibited binding of AGE-bovine serum albumin to cell membranes in a dose-dependent manner. Immunostaining and flow cytometry demonstrated the surface expression of OST-48 and 80K-H on numerous cell types and tissues, including mononuclear, endothelial, renal, and brain neuronal and glial cells. We conclude that the AGE receptor components p60 and p90 are identical to OST-48, and 80K-H, respectively, and that they together contribute to the processing of AGEs from extra- and intracellular compartments and in the cellular responses associated with these pathogenic substances.

Structure, tissue distribution, and chromosomal localization of the prepronociceptin gene.

Nociceptin (orphanin FQ), the newly discovered natural agonist of opioid receptor-like (ORL1) receptor, is a neuropeptide that is endowed with pronociceptive activity in vivo. Nociceptin is derived from a larger precursor, prepronociceptin (PPNOC), whose human, mouse, and rat genes we have now isolated. The PPNOC gene is highly conserved in the three species and displays organizational features that are strikingly similar to those of the genes of preproenkephalin, preprodynorphin, and preproopiomelanocortin, the precursors to endogenous opioid peptides, suggesting the four genes belong to the same family-i.e., have a common evolutionary origin. The PPNOC gene encodes a single copy of nociceptin as well as of other peptides whose sequence is strictly conserved across murine and human species; hence it is likely to be neurophysiologically significant. Northern blot analysis shows that the PPNOC gene is predominantly transcribed in the central nervous system (brain and spinal cord) and, albeit weakly, in the ovary, the sole peripheral organ expressing the gene. By using a radiation hybrid cell line panel, the PPNOC gene was mapped to the short arm of human chromosome 8 (8p21), between sequence-tagged site markers WI-5833 and WI-1172, in close proximity of the locus encoding the neurofilament light chain NEFL. Analysis of yeast artificial chromosome clones belonging to the WC8.4 contig covering the 8p21 region did not allow to detect the presence of the gene on these yeast artificial chromosomes, suggesting a gap in the coverage within this contig.

Primary structure and tissue distribution of the orphanin FQ precursor.

The heptadecapeptide orphanin FQ (OFQ) is a recently discovered neuropeptide that exhibits structural features reminiscent of the opioid peptides and that is an endogenous ligand to a G protein-coupled receptor sequentially related to the opioid receptors. We have cloned both the human and rat cDNAs encoding the OFQ precursor proteins, to investigate whether the sequence relationships existing between the opioid and OFQ systems are also found at the polypeptide precursor level, in particular whether the OFQ precursor would encode several bioactive peptides as do the opioid precursors, and to study the regional distribution of OFQ sites of synthesis. The entire precursor protein displays structural homology to the opioid peptide precursors, especially preprodynorphin and preproenkephalin. The predicted amino acid sequence of the OFQ precursor contains a putative signal peptide and one copy of the OFQ sequence flanked by pairs of basic amino acid residues. Carboxyl-terminal to the OFQ sequence, the human and rat precursors contain a stretch of 28 amino acids that is 100% conserved and thus may encode novel bioactive peptides. Two peptides derived from this stretch were synthesized but were found to be unable to activate the OFQ receptor, suggesting that if they are produced in vivo, these peptides would likely recognize receptors different from the OFQ receptor. To begin analyzing the sites of OFQ mRNA synthesis, Northern analysis of human and rat tissues were carried out and showed that the OFQ precursor mRNA is mainly expressed in the brain. In situ hybridization of rat brain slices demonstrated a regional distribution pattern of the OFQ precursor mRNA, which is distinct from that of the opioid peptide precursors. These data confirm that the OFQ system differs from the opioid system at the molecular level, although the OFQ and opioid precursors may have arisen from a common ancestral gene.

cDNA structure, tissue distribution, and chromosomal localization of rat PC7, a novel mammalian proprotein convertase closest to yeast kexin-like proteinases.

By using reverse transcription-coupled PCR on rat anterior pituitary RNA, we isolated a 285-bp cDNA coding for a novel subtilisin/kexin-like protein convertase (PC), called rat (r) PC7. By screening rat spleen and PC12 cell lambda gt11 cDNA libraries, we obtained a composite 3.5-kb full-length cDNA sequence of rPC7. The open reading frame codes for a prepro-PC with a 36-amino acid signal peptide, a 104-amino acid prosegment ending with a cleavable RAKR sequence, and a 747-amino acid type I membrane-bound glycoprotein, representing the mature form of this serine proteinase. Phylogenetic analysis suggests that PC7 represents the most divergent enzyme of the mammalian convertase family and that it is the closest member to the yeast convertases krp and kexin. Northern blot analyses demonstrated a widespread expression with the richest source of rPC7 mRNA being the colon and lymphoid-associated tissues. In situ hybridization revealed a distinctive tissue distribution that sometimes overlaps with that of furin, suggesting that PC7 has widespread proteolytic functions. The gene for PC7 (Pcsk7) was mapped to mouse chromosome 9 by linkage analysis of an interspecific backcross DNA panel.

Evolution of chlorophyll and bacteriochlorophyll: the problem of invariant sites in sequence analysis.

Competing hypotheses seek to explain the evolution of oxygenic and anoxygenic processes of photosynthesis. Since chlorophyll is less reduced and precedes bacteriochlorophyll on the modern biosynthetic pathway, it has been proposed that chlorophyll preceded bacteriochlorophyll in its evolution. However, recent analyses of nucleotide sequences that encode chlorophyll and bacteriochlorophyll biosynthetic enzymes appear to provide support for an alternative hypothesis. This is that the evolution of bacteriochlorophyll occurred earlier than the evolution of chlorophyll. Here we demonstrate that the presence of invariant sites in sequence datasets leads to inconsistency in tree building (including maximum-likelihood methods). Homologous sequences with different biological functions often share invariant sites at the same nucleotide positions. However, different constraints can also result in additional invariant sites unique to the genes, which have specific and different biological functions. Consequently, the distribution of these sites can be uneven between the different types of homologous genes. The presence of invariant sites, shared by related biosynthetic genes as well as those unique to only some of these genes, has misled the recent evolutionary analysis of oxygenic and anoxygenic photosynthetic pigments. We evaluate an alternative scheme for the evolution of chlorophyll and bacteriochlorophyll.

Prediction of protein folding class using global description of amino acid sequence.

We present a method for predicting protein folding class based on global protein chain description and a voting process. Selection of the best descriptors was achieved by a computer-simulated neural network trained on a data base consisting of 83 folding classes. Protein-chain descriptors include overall composition, transition, and distribution of amino acid attributes, such as relative hydrophobicity, predicted secondary structure, and predicted solvent exposure. Cross-validation testing was performed on 15 of the largest classes. The test shows that proteins were assigned to the correct class (correct positive prediction) with an average accuracy of 71.7%, whereas the inverse prediction of proteins as not belonging to a particular class (correct negative prediction) was 90-95% accurate. When tested on 254 structures used in this study, the top two predictions contained the correct class in 91% of the cases.