Oxidative Turnover of Soybean Root Glutamine Synthetase. In Vitro and in Vivo Studies1

Glutamine synthetase (GS) is the key enzyme in ammonia assimilation and catalyzes the ATP-dependent condensation of NH3 with glutamate to produce glutamine. GS in plants is an octameric enzyme. Recent work from our laboratory suggests that GS activity in plants may be regulated at the level of protein turnover (S.J. Temple, T.J. Knight, P.J. Unkefer, C. Sengupta-Gopalan [1993] Mol Gen Genet 236: 315–325; S.J. Temple, S. Kunjibettu, D. Roche, C. Sengupta-Gopalan [1996] Plant Physiol 112: 1723–1733; S.J. Temple, C. Sengupta-Gopalan [1997] In C.H. Foyer, W.P. Quick, eds, A Molecular Approach to Primary Metabolism in Higher Plants. Taylor & Francis, London, pp 155–177). Oxidative modification of GS has been implicated as the first step in the turnover of GS in bacteria. By incubating soybean (Glycine max) root extract enriched in GS in a metal-catalyzed oxidation system to produce the ·OH radical, we have shown that GS is oxidized and that oxidized GS is inactive and more susceptible to degradation than nonoxidized GS. Histidine and cysteine protect GS from metal-catalyzed inactivation, indicating that oxidation modifies the GS active site and that cysteine and histidine residues are the site of modification. Similarly, ATP and particularly ATP/glutamate give the enzyme the greatest protection against oxidative inactivation. The roots of plants fed ammonium nitrate showed a 3-fold increase in the level of GS polypeptides and activity compared with plants not fed ammonium nitrate but without a corresponding increase in the GS transcript level. This would suggest either translational or posttranslational control of GS levels.

High-resolution minisatellite-based typing as a portable approach to global analysis of Mycobacterium tuberculosis molecular epidemiology

The worldwide threat of tuberculosis to human health emphasizes the need to develop novel approaches to a global epidemiological surveillance. The current standard for Mycobacterium tuberculosis typing based on IS6110 restriction fragment length polymorphism (RFLP) suffers from the difficulty of comparing data between independent laboratories. Here, we propose a high-resolution typing method based on variable number tandem repeats (VNTRs) of genetic elements named mycobacterial interspersed repetitive units (MIRUs) in 12 human minisatellite-like regions of the M. tuberculosis genome. MIRU-VNTR profiles of 72 different M. tuberculosis isolates were established by PCR analysis of all 12 loci. From 2 to 8 MIRU-VNTR alleles were identified in the 12 regions in these strains, which corresponds to a potential of over 16 million different combinations, yielding a resolution power close to that of IS6110-RFLP. All epidemiologically related isolates tested were perfectly clustered by MIRU-VNTR typing, indicating that the stability of these MIRU-VNTRs is adequate to track outbreak episodes. The correlation between genetic relationships inferred from MIRU-VNTR and IS6110-RFLP typing was highly significant. Compared with IS6110-RFLP, high-resolution MIRU-VNTR typing has the considerable advantages of being fast, appropriate for all M. tuberculosis isolates, including strains that have a few IS6110 copies, and permitting easy and rapid comparison of results from independent laboratories. This typing method opens the way to the construction of digital global databases for molecular epidemiology studies of M. tuberculosis.

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

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

Detection of 400-year-old Yersinia pestis DNA in human dental pulp: An approach to the diagnosis of ancient septicemia

Ancient septicemic plague epidemics were reported to have killed millions of people for 2 millenniums. However, confident diagnosis of ancient septicemia solely on the basis of historical clinical observations is not possible. The lack of suitable infected material has prevented direct demonstration of ancient septicemia; thus, the history of most infections such as plague remains hypothetical. The durability of dental pulp, together with its natural sterility, makes it a suitable material on which to base such research. We hypothesized that it would be a lasting refuge for Yersinia pestis, the plague agent. DNA extracts were made from the dental pulp of 12 unerupted teeth extracted from skeletons excavated from 16th and 18th century French graves of persons thought to have died of plague (“plague teeth”) and from 7 ancient negative control teeth. PCRs incorporating ancient DNA extracts and primers specific for the human β-globin gene demonstrated the absence of inhibitors in these preparations. The incorporation of primers specific for Y. pestis rpoB (the RNA polymerase β-subunit-encoding gene) and the recognized virulence-associated pla (the plasminogen activator-encoding gene) repeatedly yielded products that had a nucleotide sequence indistinguishable from that of modern day isolates of the bacterium. The specific pla sequence was obtained from 6 of 12 plague skeleton teeth but 0 of 7 negative controls (P < 0.034, Fisher exact test). A nucleic acid-based confirmation of ancient plague was achieved for historically identified victims, and we have confirmed the presence of the disease at the end of 16th century in France. Dental pulp is an attractive target in the quest to determine the etiology of septicemic illnesses detected in ancient corpses. Molecular techniques could be applied to this material to resolve historical outbreaks.

Direct molecular level measurements of the electrostatic properties of a protein surface

In this work, we used direct measurements with the surface force apparatus to determine the pH-dependent electrostatic charge density of a single binding face of streptavidin. Mean field calculations have been used with considerable success to model electrostatic potential fields near protein surfaces, but these models and their inherent assumptions have not been tested directly at the molecular level. Using the force apparatus and immobilized, oriented monolayers of streptavidin, we measured a pI of 5–5.5 for the biotin-binding face of the protein. This differs from the pI of 6.3 for the soluble protein and confirms that we probed the local electrostatic features of the macromolecule. With finite difference solutions of the linearized Poisson–Boltzmann equation, we then calculated the pH-dependent charge densities adjacent to the same face of the protein. These calculated values agreed quantitatively with those obtained by direct force measurements. Although our study focuses on the pH-dependence of surface electrostatics, this direct approach to probing the electrostatic features of proteins is applicable to investigations of any perturbations that alter the charge distribution of the surfaces of immobilized molecules.

Cell surface expression of the epithelial Na channel and a mutant causing Liddle syndrome: A quantitative approach

The epithelial amiloride-sensitive sodium channel (ENaC) controls transepithelial Na+ movement in Na+-transporting epithelia and is associated with Liddle syndrome, an autosomal dominant form of salt-sensitive hypertension. Detailed analysis of ENaC channel properties and the functional consequences of mutations causing Liddle syndrome has been, so far, limited by lack of a method allowing specific and quantitative detection of cell-surface-expressed ENaC. We have developed a quantitative assay based on the binding of 125I-labeled M2 anti-FLAG monoclonal antibody (M2Ab*) directed against a FLAG reporter epitope introduced in the extracellular loop of each of the α, β, and γ ENaC subunits. Insertion of the FLAG epitope into ENaC sequences did not change its functional and pharmacological properties. The binding specificity and affinity (Kd = 3 nM) allowed us to correlate in individual Xenopus oocytes the macroscopic amiloride-sensitive sodium current (INa) with the number of ENaC wild-type and mutant subunits expressed at the cell surface. These experiments demonstrate that: (i) only heteromultimeric channels made of α, β, and γ ENaC subunits are maximally and efficiently expressed at the cell surface; (ii) the overall ENaC open probability is one order of magnitude lower than previously observed in single-channel recordings; (iii) the mutation causing Liddle syndrome (β R564stop) enhances channel activity by two mechanisms, i.e., by increasing ENaC cell surface expression and by changing channel open probability. This quantitative approach provides new insights on the molecular mechanisms underlying one form of salt-sensitive hypertension.

Molecular computation: RNA solutions to chess problems

We have expanded the field of “DNA computers” to RNA and present a general approach for the solution of satisfiability problems. As an example, we consider a variant of the “Knight problem,” which asks generally what configurations of knights can one place on an n × n chess board such that no knight is attacking any other knight on the board. Using specific ribonuclease digestion to manipulate strands of a 10-bit binary RNA library, we developed a molecular algorithm and applied it to a 3 × 3 chessboard as a 9-bit instance of this problem. Here, the nine spaces on the board correspond to nine “bits” or placeholders in a combinatorial RNA library. We recovered a set of “winning” molecules that describe solutions to this problem.

Molecular mechanisms underlying differential odor responses of a mouse olfactory receptor

The prevailing paradigm for G protein-coupled receptors is that each receptor is narrowly tuned to its ligand and closely related agonists. An outstanding problem is whether this paradigm applies to olfactory receptor (ORs), which is the largest gene family in the genome, in which each of 1,000 different G protein-coupled receptors is believed to interact with a range of different odor molecules from the many thousands that comprise “odor space.” Insights into how these interactions occur are essential for understanding the sense of smell. Key questions are: (i) Is there a binding pocket? (ii) Which amino acid residues in the binding pocket contribute to peak affinities? (iii) How do affinities change with changes in agonist structure? To approach these questions, we have combined single-cell PCR results [Malnic, B., Hirono, J., Sato, T. & Buck, L. B. (1999) Cell 96, 713–723] and well-established molecular dynamics methods to model the structure of a specific OR (OR S25) and its interactions with 24 odor compounds. This receptor structure not only points to a likely odor-binding site but also independently predicts the two compounds that experimentally best activate OR S25. The results provide a mechanistic model for olfactory transduction at the molecular level and show how the basic G protein-coupled receptor template is adapted for encoding the enormous odor space. This combined approach can significantly enhance the identification of ligands for the many members of the OR family and also may shed light on other protein families that exhibit broad specificities, such as chemokine receptors and P450 oxidases.

Molecular evolution of two vertebrate aryl hydrocarbon (dioxin) receptors (AHR1 and AHR2) and the PAS family

The aryl hydrocarbon receptor (AHR) is a ligand-activated transcription factor through which halogenated aromatic hydrocarbons such as 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) cause altered gene expression and toxicity. The AHR belongs to the basic helix–loop–helix/Per-ARNT-Sim (bHLH-PAS) family of transcriptional regulatory proteins, whose members play key roles in development, circadian rhythmicity, and environmental homeostasis; however, the normal cellular function of the AHR is not yet known. As part of a phylogenetic approach to understanding the function and evolutionary origin of the AHR, we sequenced the PAS homology domain of AHRs from several species of early vertebrates and performed phylogenetic analyses of these AHR amino acid sequences in relation to mammalian AHRs and 24 other members of the PAS family. AHR sequences were identified in a teleost (the killifish Fundulus heteroclitus), two elasmobranch species (the skate Raja erinacea and the dogfish Mustelus canis), and a jawless fish (the lamprey Petromyzon marinus). Two putative AHR genes, designated AHR1 and AHR2, were found both in Fundulus and Mustelus. Phylogenetic analyses indicate that the AHR2 genes in these two species are orthologous, suggesting that an AHR gene duplication occurred early in vertebrate evolution and that multiple AHR genes may be present in other vertebrates. Database searches and phylogenetic analyses identified four putative PAS proteins in the nematode Caenorhabditis elegans, including possible AHR and ARNT homologs. Phylogenetic analysis of the PAS gene family reveals distinct clades containing both invertebrate and vertebrate PAS family members; the latter include paralogous sequences that we propose have arisen by gene duplication early in vertebrate evolution. Overall, our analyses indicate that the AHR is a phylogenetically ancient protein present in all living vertebrate groups (with a possible invertebrate homolog), thus providing an evolutionary perspective to the study of dioxin toxicity and AHR function.

Molecular cloning of a cytochrome P450 taxane 10β-hydroxylase cDNA from Taxus and functional expression in yeast

The early steps in the biosynthesis of Taxol involve the cyclization of geranylgeranyl diphosphate to taxa-4(5),11(12)-diene followed by cytochrome P450-mediated hydroxylation at C5, acetylation of this intermediate, and a second cytochrome P450-dependent hydroxylation at C10 to yield taxadien-5α-acetoxy-10β-ol. Subsequent steps of the pathway involve additional cytochrome P450 catalyzed oxygenations and CoA-dependent acylations. The limited feasibility of reverse genetic cloning of cytochrome P450 oxygenases led to the use of Taxus cell cultures induced for Taxol production and the development of an approach based on differential display of mRNA-reverse transcription-PCR, which ultimately provided full-length forms of 13 unique but closely related cytochrome P450 sequences. Functional expression of these enzymes in yeast was monitored by in situ spectrophotometry coupled to in vivo screening of oxygenase activity by feeding taxoid substrates. This strategy yielded a family of taxoid-metabolizing enzymes and revealed the taxane 10β-hydroxylase as a 1494-bp cDNA that encodes a 498-residue cytochrome P450 capable of transforming taxadienyl acetate to the 10β-hydroxy derivative; the identity of this latter pathway intermediate was confirmed by chromatographic and spectrometric means. The 10β-hydroxylase represents the initial cytochrome P450 gene of Taxol biosynthesis to be isolated by an approach that should provide access to the remaining oxygenases of the pathway.

Integrated fossil and molecular data reconstruct bat echolocation

Molecular and morphological data have important roles in illuminating evolutionary history. DNA data often yield well resolved phylogenies for living taxa, but are generally unattainable for fossils. A distinct advantage of morphology is that some types of morphological data may be collected for extinct and extant taxa. Fossils provide a unique window on evolutionary history and may preserve combinations of primitive and derived characters that are not found in extant taxa. Given their unique character complexes, fossils are critical in documenting sequences of character transformation over geologic time and may elucidate otherwise ambiguous patterns of evolution that are not revealed by molecular data alone. Here, we employ a methodological approach that allows for the integration of molecular and paleontological data in deciphering one of the most innovative features in the evolutionary history of mammals—laryngeal echolocation in bats. Molecular data alone, including an expanded data set that includes new sequences for the A2AB gene, suggest that microbats are paraphyletic but do not resolve whether laryngeal echolocation evolved independently in different microbat lineages or evolved in the common ancestor of bats and was subsequently lost in megabats. When scaffolds from molecular phylogenies are incorporated into parsimony analyses of morphological characters, including morphological characters for the Eocene taxa Icaronycteris, Archaeonycteris, Hassianycteris, and Palaeochiropteryx, the resulting trees suggest that laryngeal echolocation evolved in the common ancestor of fossil and extant bats and was subsequently lost in megabats. Molecular dating suggests that crown-group bats last shared a common ancestor 52 to 54 million years ago.

An approach to three-dimensional structures of biomolecules by using single-molecule diffraction images

We describe an approach to the high-resolution three-dimensional structural determination of macromolecules that utilizes ultrashort, intense x-ray pulses to record diffraction data in combination with direct phase retrieval by the oversampling technique. It is shown that a simulated molecular diffraction pattern at 2.5-Å resolution accumulated from multiple copies of single rubisco biomolecules, each generated by a femtosecond-level x-ray free electron laser pulse, can be successfully phased and transformed into an accurate electron density map comparable to that obtained by more conventional methods. The phase problem is solved by using an iterative algorithm with a random phase set as an initial input. The convergence speed of the algorithm is reasonably fast, typically around a few hundred iterations. This approach and phasing method do not require any ab initio information about the molecule, do not require an extended ordered lattice array, and can tolerate high noise and some missing intensity data at the center of the diffraction pattern. With the prospects of the x-ray free electron lasers, this approach could provide a major new opportunity for the high-resolution three-dimensional structure determination of single biomolecules.

New approach for inhibiting Rev function and HIV-1 production using the influenza virus NS1 protein.

The Rev protein of HIV-1, which facilitates the nuclear export of HIV-1 pre-mRNAs, has been a target for antiviral therapy. Here we describe a new strategy for inhibiting Rev function and HIV-1 replication. In contrast to previous approaches, we use a wild-type rather than a mutant Rev protein and covalently link this Rev sequence to the NS1 protein of influenza A virus, a protein that inhibits the nuclear export of mRNAs. The NS1 protein contains an RNA-binding domain mutation (RM), so that the only functional RNA-binding domain in the chimeric protein (NS1RM-Rev) is in the Rev protein sequence. In the presence of the NS1RM-Rev chimeric protein, HIV-1 pre-mRNAs were retained in, rather than exported from, the nucleus. In addition, this chimeric protein effectively inhibited Rev function in trans in transfection experiments and effectively inhibited the production of HIV-1 in tissue culture cells transfected with an infectious molecular clone of HIV-1 DNA. The inhibitory activities of the NS1RM-Rev chimera were at least equivalent to those of the Rev M10 mutant protein, which has been considered to be the prototype trans inhibitor of Rev function and is currently in phase I clinical trials for the treatment of AIDS patients. We discuss (i) the potential for increasing the inhibitory activity of NS1-Rev chimeras against HIV-1 and (ii) the need for additional studies to evaluate these chimeras for the treatment of AIDS.

Human immunodeficiency virus (HIV) type 2-mediated inhibition of HIV type 1: a new approach to gene therapy of HIV-infection.

Human immunodeficiency virus (HIV) type 2, the second AIDS-associated human retrovirus, differs from HIV-1 in its natural history, infectivity, and pathogenicity, as well as in details of its genomic structure and molecular behavior. We report here that HIV-2 inhibits the replication of HIV-1 at the molecular level. This inhibition was selective, dose-dependent, and nonreciprocal. The closely related simian immunodeficiency provirus also inhibited HIV-1. The selectivity of inhibition was shown by the observation that HIV-2 did not significantly downmodulate the expression of the unrelated murine leukemia virus; neither did the murine leukemia virus markedly affect HIV-1 or HIV-2 expression. Moreover, while HIV-2 potently inhibited HIV-1, the reverse did not happen, thus identifying yet another and remarkable difference between HIV-1 and HIV-2. Mutational analysis of the HIV-2 genome suggested that the inhibition follows a complex pathway, possibly involving multiple genes and redundant mechanisms. Introduction of inactivating mutations into the structural and regulatory/accessory genes did not render the HIV-2 provirus ineffective. Some of the HIV-2 gene defects, such as that of tat and rev genes, were phenotypically transcomplemented by HIV-1. The HIV-2 proviruses with deletions in the putative packaging signal and defective for virus replication were effective in inducing the suppressive phenotype. Though the exact mechanism remains to be defined, the inhibition appeared to be mainly due to an intracellular molecular event because it could not be explained solely on the basis of cell surface receptor mediated interference. The results support the notion that the inhibition likely occurred at the level of viral RNA, possibly involving competition between viral RNAs for some transcriptional factor essential for virus replication. Induction of a cytokine is another possibility. These findings might be relevant to the clinical-epidemiological data suggesting that infection with HIV-2 may offer some protection against HIV-1 infection.

Molecular biology of retinal ganglion cells.

Retinal ganglion cells are the output neurons that encode and transmit information from the eye to the brain. Their diverse physiologic and anatomic properties have been intensively studied and appear to account well for a number of psychophysical phenomena such as lateral inhibition and chromatic opponency. In this paper, we summarize our current view of retinal ganglion cell properties and pose a number of questions regarding underlying molecular mechanisms. As an example of one approach to understanding molecular mechanisms, we describe recent work on several POU domain transcription factors that are expressed in subsets of retinal ganglion cells and that appear to be involved in ganglion cell development.