Comparative analysis of ribonuclease P RNA using gene sequences from natural microbial populations reveals tertiary structural elements.

PCR amplification of template DNAs extracted from mixed, naturally occurring microbial populations, using oligonucleotide primers complementary to highly conserved sequences, was used to obtain a large collection of diverse RNase P RNA-encoding genes. An alignment of these sequences was used in a comparative analysis of RNase P RNA secondary and tertiary structure. The new sequences confirm the secondary structure model based on sequences from cultivated organisms (with minor alterations in helices P12 and P18), providing additional support for nearly every base pair. Analysis of sequence covariation using the entire RNase P RNA data set reveals elements of tertiary structure in the RNA; the third nucleotides (underlined) of the GNRA tetraloops L14 and L18 are seen to interact with adjacent Watson-Crick base pairs in helix P8, forming A:G/C or G:A/U base triples. These experiments demonstrate one way in which the enormous diversity of natural microbial populations can be used to elucidate molecular structure through comparative analysis.

RNA analysis by ion-pair reversed-phase high performance liquid chromatography

Ion-pair reversed-phase high performance liquid chromatography (IP RP HPLC) is presented as a new, superior method for the analysis of RNA. IP RP HPLC provides a fast and reliable alternative to classical methods of RNA analysis, including separation of different RNA species, quantification and purification. RNA is stable under the analysis conditions used; degradation of RNA during the analyses was not observed. The versatility of IP RP HPLC for RNA analysis is demonstrated. Components of an RNA ladder, ranging in size from 155 to 1770 nt, were resolved. RNA transcripts of up to 5219 nt were analyzed, their integrity determined and they were quantified and purified. Purification of mRNA from total RNA is described, separating mouse rRNA from poly(A)+ mRNA. IP RP HPLC is also suitable for the separation and purification of DIG-labeled from unlabeled RNA. RNA purified by IP RP HPLC exhibits improved stability.

Statistical analysis of shape through triangulation of landmarks: A study of sexual dimorphism in hominids

Two objects with homologous landmarks are said to be of the same shape if the configuration of landmarks of one object can be exactly matched with that of the other by translation, rotation/reflection, and scaling. In an earlier paper, the authors proposed statistical analysis of shape by considering logarithmic differences of all possible Euclidean distances between landmarks. Tests of significance for differences in the shape of objects and methods of discrimination between populations were developed with such data. In the present paper, the corresponding statistical methodology is developed by triangulation of the landmarks and by considering the angles as natural measurements of shape. This method is applied to the study of sexual dimorphism in hominids.

A pair of adjacent glucocorticoid response elements regulate expression of two mouse metallothionein genes

Synthesis of mouse metallothionein (MT)-I and MT-II is transcriptionally induced by the synthetic glucocorticoid, dexamethasone (DEX) or both in vivo as well as in numerous cell lines. However, the location(s) of a glucocorticoid response element (GRE) has not been described. The observation that a marked MT-I gene, as well as heterologous genes, when placed in the context of 17 kb of flanking sequence from the MT locus, are inducible by DEX and lipopolysaccharide in transgenic mice renewed the search for the GRE. Analysis of a series of deletion constructs from this 17-kb region in cultured cells identified a single 455-bp region that conferred DEX induction on a reporter gene. This 455-bp region contains two GREs that bind to the glucocorticoid receptor as assessed by gel mobility shift. Deletion of this fragment from the 17-kb flanking region eliminates the DEX responsiveness of reporter genes. The two GREs, which are located ≈1 kb upstream of the MT-II gene and ≈7 kb upstream of the MT-I gene, are necessary for induction of both genes and can function independently of elements within the proximal promoter region of either gene.

The orientation of the AP-1 heterodimer on DNA strongly affects transcriptional potency

Activation of gene transcription in eukaryotes requires the cooperative assembly of an initiation complex containing many protein subunits. The necessity that these components contact each other and the promoter/enhancer in defined ways suggests that their spatial arrangement might influence the activation response. Indeed, growing evidence indicates that DNA architecture can profoundly affect transcriptional potency. Much less is known about the influence of protein architecture on transcriptional activation. Here, we examine the architectural dependence of activator function through the analysis of matched pairs of AP-1•DNA complexes differing only in their orientation. Mutation of a critical Arg residue in the basic-leucine zipper domain of either Fos or Jun yielded single point-mutant heterodimers that bind DNA in a single defined orientation, as determined directly by native chemical ligation/affinity cleavage; by contrast, the corresponding wild-type protein binds DNA as a roughly equal mixture of two isomeric orientations, which are related by subunit interchange. The stereochemistry of the point-mutant heterodimers could be switched by inversion of a C•G base pair in the center of the AP-1 site, thus providing access to both fixed orientational isomers. Yeast reporter gene assays consistently revealed that one orientational isomer activates transcription at least 10-fold more strongly than the other. These results suggest that protein architecture, especially the spatial relationship of the activation domain to the promoter, can exert a powerful influence on activator potency.

A novel pair of immunoglobulin-like receptors expressed by B cells and myeloid cells

An Fcα receptor probe of human origin was used to identify novel members of the Ig gene superfamily in mice. Paired Ig-like receptors, named PIR-A and PIR-B, are predicted from sequence analysis of the cDNAs isolated from a mouse splenic library. Both type I transmembrane proteins possess similar ectodomains with six Ig-like loops, but have different transmembrane and cytoplasmic regions. The predicted PIR-A protein has a short cytoplasmic tail and a charged Arg residue in the transmembrane region that, by analogy with the FcαR relative, suggests the potential for association with an additional transmembrane protein to form a signal transducing unit. In contrast, the PIR-B protein has an uncharged transmembrane region and a long cytoplasmic tail containing four potential immunoreceptor tyrosine-based inhibitory motifs. These features are shared by the related killer inhibitory receptors. PIR-A proteins appear to be highly variable, in that predicted peptide sequences differ for seven randomly selected PIR-A clones, whereas PIR-B cDNA clones are invariant. Southern blot analysis with PIR-B and PIR-A-specific probes suggests only one PIR-B gene and multiple PIR-A genes. The PIR-A and PIR-B genes are expressed in B lymphocytes and myeloid lineage cells, wherein both are expressed simultaneously. The characteristics of the highly-conserved PIR-A and PIR-B genes and their coordinate cellular expression suggest a potential regulatory role in humoral, inflammatory, and allergic responses.

Proper sorting of the cation-dependent mannose 6-phosphate receptor in endosomes depends on a pair of aromatic amino acids in its cytoplasmic tail

The 67-amino acid cytoplasmic tail of the cation-dependent mannose 6-phosphate receptor (CD-MPR) contains a signal(s) that prevents the receptor from entering lysosomes where it would be degraded. To identify the key residues required for proper endosomal sorting, we analyzed the intracellular distribution of mutant forms of the receptor by Percoll density gradients. A receptor with a Trp19 → Ala substitution in the cytoplasmic tail was highly missorted to lysosomes whereas receptors with either Phe18 → Ala or Phe13 → Ala mutations were partially defective in avoiding transport to lysosomes. Analysis of double and triple mutants confirmed the key role of Trp19 for sorting of the CD-MPR in endosomes, with Phe18, Phe13, and several neighboring residues contributing to this function. The addition of the Phe18-Trp19 motif of the CD-MPR to the cytoplasmic tail of the lysosomal membrane protein Lamp1 was sufficient to partially impair its delivery to lysosomes. Replacing Phe18 and Trp19 with other aromatic amino acids did not impair endosomal sorting of the CD-MPR, indicating that two aromatic residues located at these positions are sufficient to prevent the receptor from trafficking to lysosomes. However, alterations in the spacing of the diaromatic amino acid sequence relative to the transmembrane domain resulted in receptor accumulation in lysosomes. These findings indicate that the endosomal sorting of the CD-MPR depends on the correct presentation of a diaromatic amino acid-containing motif in its cytoplasmic tail. Because a diaromatic amino acid sequence is also present in the cytoplasmic tail of other receptors known to be internalized from the plasma membrane, this feature may prove to be a general determinant for endosomal sorting.

Identification of diverse nerve growth factor-regulated genes by serial analysis of gene expression (SAGE) profiling

Neurotrophic factors such as nerve growth factor (NGF) promote a wide variety of responses in neurons, including differentiation, survival, plasticity, and repair. Such actions often require changes in gene expression. To identify the regulated genes and thereby to more fully understand the NGF mechanism, we carried out serial analysis of gene expression (SAGE) profiling of transcripts derived from rat PC12 cells before and after NGF-promoted neuronal differentiation. Multiple criteria supported the reliability of the profile. Approximately 157,000 SAGE tags were analyzed, representing at least 21,000 unique transcripts. Of these, nearly 800 were regulated by 6-fold or more in response to NGF. Approximately 150 of the regulated transcripts have been matched to named genes, the majority of which were not previously known to be NGF-responsive. Functional categorization of the regulated genes provides insight into the complex, integrated mechanism by which NGF promotes its multiple actions. It is anticipated that as genomic sequence information accrues the data derived here will continue to provide information about neurotrophic factor mechanisms.

Catalytic mechanism of the adenylyl and guanylyl cyclases: Modeling and mutational analysis

The adenylyl and guanylyl cyclases catalyze the formation of 3′,5′-cyclic adenosine or guanosine monophosphate from the corresponding nucleoside 5′-triphosphate. The guanylyl cyclases, the mammalian adenylyl cyclases, and their microbial homologues function as pairs of homologous catalytic domains. The crystal structure of the rat type II adenylyl cyclase C2 catalytic domain was used to model by homology a mammalian adenylyl cyclase C1-C2 domain pair, a homodimeric adenylyl cyclase of Dictyostelium discoideum, a heterodimeric soluble guanylyl cyclase, and a homodimeric membrane guanylyl cyclase. Mg2+ATP or Mg2+GTP were docked into the active sites based on known stereochemical constraints on their conformation. The models are consistent with the activities of seven active-site mutants. Asp-310 and Glu-432 of type I adenylyl cyclase coordinate a Mg2+ ion. The D310S and D310A mutants have 10-fold reduced Vmax and altered [Mg2+] dependence. The NTP purine moieties bind in mostly hydrophobic pockets. Specificity is conferred by a Lys and an Asp in adenylyl cyclase, and a Glu, an Arg, and a Cys in guanylyl cyclase. The models predict that an Asp from one domain is a general base in the reaction, and that the transition state is stabilized by a conserved Asn-Arg pair on the other domain.

Global analysis of proteasomal substrate specificity using positional-scanning libraries of covalent inhibitors

The proteasome is a large protease complex consisting of multiple catalytic subunits that function simultaneously to digest protein substrates. This complexity has made deciphering the role each subunit plays in the generation of specific protein fragments difficult. Positional scanning libraries of peptide vinyl sulfones were generated in which the amino acid located directly at the site of hydrolysis (P1 residue) was held constant and sequences distal to that residue (P2, P3, and P4 positions) were varied across all natural amino acids (except cysteine and methionine). Binding information for each of the individual catalytic subunits was obtained for each library under a variety of different conditions. The resulting specificity profiles indicated that substrate positions distal to P1 are critical for directing substrates to active subunits in the complex. Furthermore, specificity profiles of IFN-γ-regulated subunits closely matched those of their noninducible counterparts, suggesting that subunit swapping may modulate substrate processing by a mechanism that does require a change in the primary sequence specificity of individual catalytic subunits in the complex. Finally, specificity profiles were used to design specific inhibitors of a single active site in the complex. These reagents can be used to further establish the role of each subunit in substrate processing by the proteasome.

The Chlamydomonas PF6 Locus Encodes a Large Alanine/Proline-Rich Polypeptide That Is Required for Assembly of a Central Pair Projection and Regulates Flagellar Motility

Efficient motility of the eukaryotic flagellum requires precise temporal and spatial control of its constituent dynein motors. The central pair and its associated structures have been implicated as important members of a signal transduction cascade that ultimately regulates dynein arm activity. To identify central pair components involved in this process, we characterized a Chlamydomonas motility mutant (pf6-2) obtained by insertional mutagenesis. pf6-2 flagella twitch ineffectively and lack the 1a projection on the C1 microtubule of the central pair. Transformation with constructs containing a full-length, wild-type copy of the PF6 gene rescues the functional, structural, and biochemical defects associated with the pf6 mutation. Sequence analysis indicates that the PF6 gene encodes a large polypeptide that contains numerous alanine-rich, proline-rich, and basic domains and has limited homology to an expressed sequence tag derived from a human testis cDNA library. Biochemical analysis of an epitope-tagged PF6 construct demonstrates that the PF6 polypeptide is an axonemal component that cosediments at 12.6S with several other polypeptides. The PF6 protein appears to be an essential component required for assembly of some of these polypeptides into the C1-1a projection.

Molecular Analysis of (R)-(+)-Mandelonitrile Lyase Microheterogeneity in Black Cherry1

The flavoprotein (R)-(+)-mandelonitrile lyase (MDL; EC 4.1.2.10), which plays a key role in cyanogenesis in rosaceous stone fruits, occurs in black cherry (Prunus serotina Ehrh.) homogenates as several closely related isoforms. Biochemical and molecular biological methods were used to investigate MDL microheterogeneity and function in this species. Three novel MDL cDNAs of high sequence identity (designated MDL2, MDL4, and MDL5) were isolated. Like MDL1 and MDL3 cDNAs (Z. Hu, J.E. Poulton [1997] Plant Physiol 115: 1359–1369), they had open reading frames that predicted a flavin adenine dinucleotide-binding site, multiple N-glycosylation sites, and an N-terminal signal sequence. The N terminus of an MDL isoform purified from seedlings matched the derived amino acid sequence of the MDL4 cDNA. Genomic sequences corresponding to the MDL1, MDL2, and MDL4 cDNAs were obtained by polymerase chain reaction amplification of genomic DNA. Like the previously reported mdl3 gene, these genes are interrupted at identical positions by three short, conserved introns. Given their overall similarity, we conclude that the genes mdl1, mdl2, mdl3, mdl4, and mdl5 are derived from a common ancestral gene and constitute members of a gene family. Genomic Southern-blot analysis showed that this family has approximately eight members. Northern-blot analysis using gene-specific probes revealed differential expression of the genes mdl1, mdl2, mdl3, mdl4, and mdl5.

Genetic fidelity under harsh conditions: Analysis of spontaneous mutation in the thermoacidophilic archaeon Sulfolobus acidocaldarius

Microbes whose genomes are encoded by DNA and for which adequate information is available display similar genomic mutation rates (average 0.0034 mutations per chromosome replication, range 0.0025 to 0.0046). However, this value currently is based on only a few well characterized microbes reproducing within a narrow range of environmental conditions. In particular, no genomic mutation rate has been determined either for a microbe whose natural growth conditions may extensively damage DNA or for any member of the archaea, a prokaryotic lineage deeply diverged from both bacteria and eukaryotes. Both of these conditions are met by the extreme thermoacidophile Sulfolobus acidocaldarius. We determined the genomic mutation rate for this species when growing at pH 3.5 and 75°C based on the rate of forward mutation at the pyrE gene and the nucleotide changes identified in 101 independent mutants. The observed value of about 0.0018 extends the range of DNA-based microbes with rates close to the standard rate simultaneously to an archaeon and to an extremophile whose cytoplasmic pH and normal growth temperature greatly accelerate the spontaneous decomposition of DNA. The mutations include base pair substitutions (BPSs) and additions and deletions of various sizes, but the S. acidocaldarius spectrum differs from those of other DNA-based organisms in being relatively poor in BPSs. The paucity of BPSs cannot yet be explained by known properties of DNA replication or repair enzymes of Sulfolobus spp. It suggests, however, that molecular evolution per genome replication may proceed more slowly in S. acidocaldarius than in other DNA-based organisms examined to date.

Statistical analysis of shape of objects based on landmark data.

Two objects with homologous landmarks are said to be of the same shape if the configurations of landmarks of one object can be exactly matched with that of the other by translation, rotation/reflection, and scaling. The observations on an object are coordinates of its landmarks with reference to a set of orthogonal coordinate axes in an appropriate dimensional space. The origin, choice of units, and orientation of the coordinate axes with respect to an object may be different from object to object. In such a case, how do we quantify the shape of an object, find the mean and variation of shape in a population of objects, compare the mean shapes in two or more different populations, and discriminate between objects belonging to two or more different shape distributions. We develop some methods that are invariant to translation, rotation, and scaling of the observations on each object and thereby provide generalizations of multivariate methods for shape analysis.

Accurate reconstruction of a known HIV-1 transmission history by phylogenetic tree analysis.

Phylogenetic analyses are increasingly used in attempts to clarify transmission patterns of human immunodeficiency virus type 1 (HIV-1), but there is a continuing discussion about their validity because convergent evolution and transmission of minor HIV variants may obscure epidemiological patterns. Here we have studied a unique HIV-1 transmission cluster consisting of nine infected individuals, for whom the time and direction of each virus transmission was exactly known. Most of the transmissions occurred between 1981 and 1983, and a total of 13 blood samples were obtained approximately 2-12 years later. The p17 gag and env V3 regions of the HIV-1 genome were directly sequenced from uncultured lymphocytes. A true phylogenetic tree was constructed based on the knowledge about when the transmissions had occurred and when the samples were obtained. This complex, known HIV-1 transmission history was compared with reconstructed molecular trees, which were calculated from the DNA sequences by several commonly used phylogenetic inference methods [Fitch-Margoliash, neighbor-joining, minimum-evolution, maximum-likelihood, maximum-parsimony, unweighted pair group method using arithmetic averages (UPGMA), and a Fitch-Margoliash method assuming a molecular clock (KITSCH)]. A majority of the reconstructed trees were good estimates of the true phylogeny; 12 of 13 taxa were correctly positioned in the most accurate trees. The choice of gene fragment was found to be more important than the choice of phylogenetic method and substitution model. However, methods that are sensitive to unequal rates of change performed more poorly (such as UPGMA and KITSCH, which assume a constant molecular clock). The rapidly evolving V3 fragment gave better reconstructions than p17, but a combined data set of both p17 and V3 performed best. The accuracy of the phylogenetic methods justifies their use in HIV-1 research and argues against convergent evolution and selective transmission of certain virus variants.