In silico detection of control signals: mRNA 3′-end-processing sequences in diverse species

We have investigated mRNA 3′-end-processing signals in each of six eukaryotic species (yeast, rice, arabidopsis, fruitfly, mouse, and human) through the analysis of more than 20,000 3′-expressed sequence tags. The use and conservation of the canonical AAUAAA element vary widely among the six species and are especially weak in plants and yeast. Even in the animal species, the AAUAAA signal does not appear to be as universal as indicated by previous studies. The abundance of single-base variants of AAUAAA correlates with their measured processing efficiencies. As found previously, the plant polyadenylation signals are more similar to those of yeast than to those of animals, with both common content and arrangement of the signal elements. In all species examined, the complete polyadenylation signal appears to consist of an aggregate of multiple elements. In light of these and previous results, we present a broadened concept of 3′-end-processing signals in which no single exact sequence element is universally required for processing. Rather, the total efficiency is a function of all elements and, importantly, an inefficient word in one element can be compensated for by strong words in other elements. These complex patterns indicate that effective tools to identify 3′-end-processing signals will require more than consensus sequence identification.

Identification of protein–protein interfaces by decreased amide proton solvent accessibility

Matrix-assisted laser desorption ionization–time-of-flight mass spectrometry was used to identify peptic fragments from protein complexes that retained deuterium under hydrogen exchange conditions due to decreased solvent accessibility at the interface of the complex. Short deuteration times allowed preferential labeling of rapidly exchanging surface amides so that primarily solvent accessibility changes and not conformational changes were detected. A single mass spectrum of the peptic digest mixture was analyzed to determine the deuterium content of all proteolytic fragments of the protein. The protein–protein interface was reliably indicated by those peptides that retained more deuterons in the complex compared with control experiments in which only one protein was present. The method was used to identify the kinase inhibitor [PKI(5–24)] and ATP-binding sites in the cyclic-AMP-dependent protein kinase. Three overlapping peptides identified the ATP-binding site, three overlapping peptides identified the glycine-rich loop, and two peptides identified the PKI(5–24)-binding site. A complex of unknown structure also was analyzed, human α-thrombin bound to an 83-aa fragment of human thrombomodulin [TMEGF(4–5)]. Five peptides from thrombin showed significantly decreased solvent accessibility in the complex. Three peptides identified the anion-binding exosite I, confirming ligand competition experiments. Two peptides identified a new region of thrombin near the active site providing a potential mechanism of how thrombomodulin alters thrombin substrate specificity.

Identification of Potential Regulatory Elements for the Transport of Emp24p

To examine the possibility of active recycling of Emp24p between the endoplasmic reticulum (ER) and the Golgi, we sought to identify transport signal(s) in the carboxyl-terminal region of Emp24p. Reporter molecules were constructed by replacing parts of a control invertase-Wbp1p chimera with those of Emp24p, and their transport rates were assessed. The transport of the reporter was found to be accelerated by the presence of the cytoplasmic domain of Emp24p. Mutational analyses revealed that the two carboxyl-terminal residues, leucine and valine (LV), were necessary and sufficient to accelerate the transport. The acceleration was sequence specific, and the terminal valine appeared to be more important. The LV residues accelerated not only the overall transport to the vacuole but also the ER to cis-Golgi transport, suggesting its function in the ER export. Hence the LV residues are a novel anterograde transport signal. The double-phenylalanine residues did not affect the transport by itself but attenuated the effect of the anterograde transport signal. On the other hand, the transmembrane domain significantly slowed down the ER to cis-Golgi transport and effectively counteracted the anterograde transport signal at this step. It may also take part in the retrieval of the protein, because the overall transport to the vacuole was more evidently slowed down. Consistently, the mutation of a conserved glutamine residue in the transmembrane domain further slowed down the transport in a step after arriving at the cis-Golgi. Taken together, the existence of the anterograde transport signal and the elements that regulate its function support the active recycling of Emp24p.

Epithelial and endothelial expression of the green fluorescent protein reporter gene under the control of bovine prion protein (PrP) gene regulatory sequences in transgenic mice

The expression of the cellular form of the prion protein (PrPc) gene is required for prion replication and neuroinvasion in transmissible spongiform encephalopathies. The identification of the cell types expressing PrPc is necessary to understanding how the agent replicates and spreads from peripheral sites to the central nervous system. To determine the nature of the cell types expressing PrPc, a green fluorescent protein reporter gene was expressed in transgenic mice under the control of 6.9 kb of the bovine PrP gene regulatory sequences. It was shown that the bovine PrP gene is expressed as two populations of mRNA differing by alternative splicing of one 115-bp 5′ untranslated exon in 17 different bovine tissues. The analysis of transgenic mice showed reporter gene expression in some cells that have been identified as expressing PrP, such as cerebellar Purkinje cells, lymphocytes, and keratinocytes. In addition, expression of green fluorescent protein was observed in the plexus of the enteric nervous system and in a restricted subset of cells not yet clearly identified as expressing PrP: the epithelial cells of the thymic medullary and the endothelial cells of both the mucosal capillaries of the intestine and the renal capillaries. These data provide valuable information on the distribution of PrPc at the cellular level and argue for roles of the epithelial and endothelial cells in the spread of infection from the periphery to the brain. Moreover, the transgenic mice described in this paper provide a model that will allow for the study of the transcriptional activity of the PrP gene promoter in response to scrapie infection.

Identification of the vitamin K-dependent carboxylase active site: Cys-99 and Cys-450 are required for both epoxidation and carboxylation

The vitamin K-dependent carboxylase modifies and renders active vitamin K-dependent proteins involved in hemostasis, cell growth control, and calcium homeostasis. Using a novel mechanism, the carboxylase transduces the free energy of vitamin K hydroquinone (KH2) oxygenation to convert glutamate into a carbanion intermediate, which subsequently attacks CO2, generating the γ-carboxylated glutamate product. How the carboxylase effects this conversion is poorly understood because the active site has not been identified. Dowd and colleagues [Dowd, P., Hershline, R., Ham, S. W. & Naganathan, S. (1995) Science 269, 1684–1691] have proposed that a weak base (cysteine) produces a strong base (oxygenated KH2) capable of generating the carbanion. To define the active site and test this model, we identified the amino acids that participate in these reactions. N-ethyl maleimide inhibited epoxidation and carboxylation, and both activities were equally protected by KH2 preincubation. Amino acid analysis of 14C- N-ethyl maleimide-modified human carboxylase revealed 1.8–2.3 reactive residues and a specific activity of 7 × 108 cpm/hr per mg. Tryptic digestion and liquid chromatography electrospray mass spectrometry identified Cys-99 and Cys-450 as active site residues. Mutation to serine reduced both epoxidation and carboxylation, to 0.2% (Cys-99) or 1% (Cys-450), and increased the Kms for a glutamyl substrate 6- to 8-fold. Retention of some activity indicates a mechanism for enhancing cysteine/serine nucleophilicity, a property shared by many active site thiol enzymes. These studies, which represent a breakthrough in defining the carboxylase active site, suggest a revised model in which the glutamyl substrate indirectly coordinates at least one thiol, forming a catalytic complex that ionizes a thiol to initiate KH2 oxygenation.

Building a dictionary for genomes: Identification of presumptive regulatory sites by statistical analysis

The availability of complete genome sequences and mRNA expression data for all genes creates new opportunities and challenges for identifying DNA sequence motifs that control gene expression. An algorithm, “MobyDick,” is presented that decomposes a set of DNA sequences into the most probable dictionary of motifs or words. This method is applicable to any set of DNA sequences: for example, all upstream regions in a genome or all genes expressed under certain conditions. Identification of words is based on a probabilistic segmentation model in which the significance of longer words is deduced from the frequency of shorter ones of various lengths, eliminating the need for a separate set of reference data to define probabilities. We have built a dictionary with 1,200 words for the 6,000 upstream regulatory regions in the yeast genome; the 500 most significant words (some with as few as 10 copies in all of the upstream regions) match 114 of 443 experimentally determined sites (a significance level of 18 standard deviations). When analyzing all of the genes up-regulated during sporulation as a group, we find many motifs in addition to the few previously identified by analyzing the subclusters individually to the expression subclusters. Applying MobyDick to the genes derepressed when the general repressor Tup1 is deleted, we find known as well as putative binding sites for its regulatory partners.

Identification of rat susceptibility loci for adjuvant-oil-induced arthritis

One intradermal injection of incomplete Freund’s adjuvant-oil induces a T cell-mediated inflammatory joint disease in DA rats. Susceptibility genes for oil-induced arthritis (OIA) are located both within and outside the major histocompatibility complex (MHC, Oia1). We have searched for disease-linked non-MHC loci in an F2 intercross between DA rats and MHC-identical but arthritis-resistant LEW.1AV1 rats. A genome-wide scan with microsatellite markers revealed two major chromosome regions that control disease incidence and severity: Oia2 on chromosome 4 (P = 4 × 10−13) and Oia3 on chromosome 10 (P = 1 × 10−6). All animals homozygous for DA alleles at both loci developed severe arthritis, whereas all those homozygous for LEW.1AV1 alleles were resistant. These results have general implications for situations where nonspecific activation of the immune system (e.g., incomplete Freund’s adjuvant-oil) causes inflammation and disease, either alone or in conjunction with specific antigens. They may also provide clues to the etiology of inflammatory diseases in humans, including rheumatoid arthritis.

Identification of two short internal ribosome entry sites selected from libraries of random oligonucleotides

Sequences that control translation of mRNA may play critical roles in regulating protein levels. One such element is the internal ribosome entry site (IRES). We previously showed that a 9-nt segment in the 5′ leader sequence of the mRNA encoding Gtx homeodomain protein could function as an IRES. To identify other short sequences with similar properties, we designed a selection procedure that uses a retroviral vector to express dicistronic mRNAs encoding enhanced green and cyan fluorescent proteins as the first and second cistrons, respectively. Expression of the second cistron was dependent upon the intercistronic sequences and was indicative of IRES activity. B104 cells were infected with two retroviral libraries that contained random sequences of 9 or 18 nt in the intercistronic region. Cells expressing both cistrons were sorted, and sequences recovered from selected cells were reassayed for IRES activity in a dual luciferase dicistronic mRNA. Two novel IRESes were identified by this procedure, and both contained segments with complementarity to 18S rRNA. When multiple copies of either segment were linked together, IRES activities were dramatically enhanced. Moreover, these synthetic IRESes were differentially active in various cell types. These properties are similar to those of the previously identified 9-nt IRES module from Gtx mRNA. These results provide further evidence that short nucleotide sequences can function as IRESes and support the idea that some cellular IRESes may be composed of shorter functional modules. The ability to identify IRES modules with specific expression properties may be useful in the design of vectors for biotechnology and gene therapy.

A novel approach for the identification of protein–protein interaction with integral membrane proteins

Protein–protein interaction plays a major role in all biological processes. The currently available genetic methods such as the two-hybrid system and the protein recruitment system are relatively limited in their ability to identify interactions with integral membrane proteins. Here we describe the development of a reverse Ras recruitment system (reverse RRS), in which the bait used encodes a membrane protein. The bait is expressed in its natural environment, the membrane, whereas the protein partner (the prey) is fused to a cytoplasmic Ras mutant. Protein–protein interaction between the proteins encoded by the prey and the bait results in Ras membrane translocation and activation of a viability pathway in yeast. We devised the expression of the bait and prey proteins under the control of dual distinct inducible promoters, thus enabling a rapid selection of transformants in which growth is attributed solely to specific protein–protein interaction. The reverse RRS approach greatly extends the usefulness of the protein recruitment systems and the use of integral membrane proteins as baits. The system serves as an attractive approach to explore novel protein–protein interactions with high specificity and selectivity, where other methods fail.

Identification and characterization of a second melanin-concentrating hormone receptor, MCH-2R

Melanin-concentrating hormone (MCH) is a 19-aa cyclic neuropeptide originally isolated from chum salmon pituitaries. Besides its effects on the aggregation of melanophores in fish several lines of evidence suggest that in mammals MCH functions as a regulator of energy homeostasis. Recently, several groups reported the identification of an orphan G protein-coupled receptor as a receptor for MCH (MCH-1R). We hereby report the identification of a second human MCH receptor termed MCH-2R, which shares about 38% amino acid identity with MCH-1R. MCH-2R displayed high-affinity MCH binding, resulting in inositol phosphate turnover and release of intracellular calcium in mammalian cells. In contrast to MCH-1R, MCH-2R signaling is not sensitive to pertussis toxin and MCH-2R cannot reduce forskolin-stimulated cAMP production, suggesting an exclusive Gαq coupling of the MCH-2R in cell-based systems. Northern blot and in situ hybridization analysis of human and monkey tissue shows that expression of MCH-2R mRNA is restricted to several regions of the brain, including the arcuate nucleus and the ventral medial hypothalamus, areas implicated in regulation of body weight. In addition, the human MCH-2R gene was mapped to the long arm of chromosome 6 at band 6q16.2–16.3, a region reported to be associated with cytogenetic abnormalities of obese patients. The characterization of a second mammalian G protein-coupled receptor for MCH potentially indicates that the control of energy homeostasis in mammals by the MCH neuropeptide system may be more complex than initially anticipated.

Identification and characterization of a melanin-concentrating hormone receptor

Melanin-concentrating hormone (MCH), a neuropeptide expressed in central and peripheral nervous systems, plays an important role in the control of feeding behaviors and energy metabolism. An orphan G protein-coupled receptor (SLC-1/GPR24) has recently been identified as a receptor for MCH (MCHR1). We report here the identification and characterization of a G protein-coupled receptor as the MCH receptor subtype 2 (MCHR2). MCHR2 has higher protein sequence homology to MCHR1 than any other G protein-coupled receptor. The expression of MCHR2 has been detected in many regions of the brain. In contrast to MCHR1, which is intronless in the coding region and is located at the chromosomal locus 22q13.3, the MCHR2 gene has multiple exons and is mapped to locus 6q21. MCHR2 is specifically activated by nanomolar concentrations of MCH, binds to MCH with high affinity, and signals through Gq protein. This discovery is important for a full understanding of MCH biology and the development of potential therapeutics for diseases involving MCH, including obesity.

Identification of three physically and functionally distinct binding sites for C3b in human complement factor H by deletion mutagenesis.

Human complement factor H controls spontaneous activation of complement in plasma and appears to play a role in distinguishing host cells from activators of the alternative pathway of complement. In both mice and humans, the protein is composed of 20 homologous short consensus repeat (SCR) domains. The size of the protein suggests that portions of the structure outside the known C3b binding site (SCR 1-4) possess a significant biological role. We have expressed the full-length cDNA of factor H in the baculovirus system and have shown the recombinant protein to be fully active. Mutants of this full-length protein have now been prepared, purified, and examined for cofactor activity and binding to C3b and heparin. The results demonstrate (i) that factor H has at least three sites that bind C3b, (ii) that one of these sites is located in SCR domains 1-4, as has been shown by others, (iii) that a second site exists in the domain 6-10 region, (iv) that a third site resides in the SCR 16-20 region, and (v) that two heparin binding sites exist in factor H, one near SCR 13 and another in the SCR 6-10 region. Functional assays demonstrated that only the first C3b site located in SCR 1-4 expresses factor I cofactor activity. Mutant proteins lacking any one of the three C3b binding sites exhibited 6- to 8-fold reductions in affinity for C3b on sheep erythrocytes, indicating that all three sites contribute to the control of complement activation on erythrocytes. The identification of multiple functionally distinct sites on factor H clarifies many of the heretofore unexplainable behaviors of this protein, including the heterogeneous binding of factor H to surface-bound C3b, the effects of trypsin cleavage, and the differential control of complement activation on activators and nonactivators of the alternative pathway of complement.

Identification of a minimal sequence of the mouse pro-alpha 1(I) collagen promoter that confers high-level osteoblast expression in transgenic mice and that binds a protein selectively present in osteoblasts.

Based on our previous transgenic mice results, which strongly suggested that separate cell-specific cis-acting elements of the mouse pro-alpha 1(I) collagen promoter control the activity of the gene in different type I collagen-producing cells, we attempted to delineate a short segment in this promoter that could direct high-level expression selectively in osteoblasts. By generating transgenic mice harboring various fragments of the promoter, we identified a 117-bp segment (-1656 to -1540) that is a minimal sequence able to confer high-level expression of a lacZ reporter gene selectively in osteoblasts when cloned upstream of the proximal 220-bp pro-alpha 1(I) promoter. This 220-bp promoter by itself was inactive in transgenic mice and unable to direct osteoblast-specific expression. The 117-bp enhancer segment contained two sequences that appeared to have different functions. The A sequence (-1656 to -1628) was required to obtain expression of the lacZ gene in osteoblasts, whereas the C sequence (-1575 to -1540) was essential to obtain consistent and high-level expression of the lacZ gene in osteoblasts. Gel shift assays showed that the A sequence bound a nuclear protein present only in osteoblastic cells. A mutation in the A segment that abolished the binding of this osteoblast-specific protein also abolished lacZ expression in osteoblasts of transgenic mice.

Fate of biological control introductions: monitoring an Australian fungal pathogen of grasshoppers in North America.

In North America there are two generally recognized pathotypes (pathotypes 1 and 2) of the fungus Entomophaga grylli which show host-preferential infection of grasshopper subfamilies. Pathotype 3, discovered in Australia, has a broader grasshopper host range and was considered to be a good biocontrol agent. Between 1989 and 1991 pathotype 3 was introduced at two field sites in North Dakota. Since resting spores are morphologically indistinguishable among pathotypes, we used pathotype-specific DNA probes to confirm pathotype identification in E. grylli-infected grasshoppers collected at the release sites in 1992, 1993, and 1994. In 1992, up to 23% of E. grylli-infected grasshoppers of the subfamilies Melanoplinae, Oedipodinae, and Gomphocerinae were infected by pathotype 3, with no infections > 1 km from the release sites. In 1993, pathotype 3 infections declined to 1.7%. In 1994 grasshopper populations were low and no pathotype 3 infections were found. The frequency of pathotype 3 infection has declined to levels where its long-term survival in North America is questionable. Analyses of biocontrol releases are critical to evaluating the environmental risks associated with these ecological manipulations, and molecular probes are powerful tools for monitoring biocontrol releases.

Identification of a cytoplasmic, phorbol ester-inducible isoform of protein tyrosine phosphatase epsilon.

The protein-tyrosine phosphatase epsilon (PTP epsilon) is a transmembranal, receptor-type protein that possesses two phosphatase catalytic domains characteristic of transmembranal phosphatases. Here we demonstrate the existence of a nontransmembranal isoform of PTP epsilon, PTP epsilon-cytoplasmic. PTP epsilon-cytoplasmic and the transmembranal isoform of PTP epsilon have separate, nonoverlapping expression patterns. Further, the data clearly indicate that control of which of the two isoforms is to be expressed is initiated at the transcriptional level, suggesting that they have distinct physiological roles. PTP epsilon-cytoplasmic mRNA is the product of a delayed early response gene in NIH 3T3 fibroblasts, and its transcription is regulated through a pathway that requires protein kinase C. The human homologue of PTP epsilon-cytoplasmic has also been cloned and is strongly up-regulated in the early stages of phorbol 12-tetradecanoate 13-acetate-induced differentiation of HL-60 cells. Sequence analysis indicates and cellular fractionation experiments confirm that this isoform is a cytoplasmic molecule. PTP epsilon-cytoplasmic is therefore the initial example to our knowledge of a nontransmembranal protein-tyrosine phosphatase that contains two tandem of catalytic domains.