The heme-copper oxidases of Thermus thermophilus catalyze the reduction of nitric oxide: Evolutionary implications

We show that the heme-copper terminal oxidases of Thermus thermophilus (called ba3 and caa3) are able to catalyze the reduction of nitric oxide (NO) to nitrous oxide (N2O) under reducing anaerobic conditions. The rate of NO consumption and N2O production were found to be linearly dependent on enzyme concentration, and activity was abolished by enzyme denaturation. Thus, contrary to the eukaryotic enzyme, both T. thermophilus oxidases display a NO reductase activity (3.0 ± 0.7 mol NO/mol ba3 × min and 32 ± 8 mol NO/mol caa3 × min at [NO] ≈ 50 μM and 20°C) that, though considerably lower than that of bona fide NO reductases (300–4,500 mol NO/mol enzyme × min), is definitely significant. We also show that for ba3 oxidase, NO reduction is associated to oxidation of cytochrome b at a rate compatible with turnover, suggesting a mechanism consistent with the stoichiometry of the overall reaction. We propose that the NO reductase activity of T. thermophilus oxidases may depend on a peculiar CuB+ coordination, which may be revealed by the forthcoming three-dimensional structure. These findings support the hypothesis of a common phylogeny of aerobic respiration and bacterial denitrification, which was proposed on the basis of structural similarities between the Pseudomonas stutzeri NO reductase and the cbb3 terminal oxidases. Our findings represent functional evidence in support of this hypothesis.

Evolutionary instability of the major histocompatibility complex class I loci in New World primates

Homologues of the human major histocompatibility complex (MHC) HLA-A, -B, -E, -F, and -G loci are present in all the Catarrhini (Old World primates, apes, and humans), and some of their allelic lineages have survived several speciation events. Analysis of 26 MHC class I cDNAs from seven different genera of New World primates revealed that the Callitrichinae (tamarins and marmosets) are an exception to these rules of MHC stability. In gene trees of primate MHC class I genes, sequences from the Callitrichinae cluster in a genus-specific fashion, whereas in the other genera of New World primates, as in the Catarrhini, they cluster in a transgeneric way. The genus-specific clustering of the Callitrichinae cDNAs indicates that there is no orthology between MHC class I loci in genera of this phyletic group. Additionally, the Callitrichinae genera exhibit limited variability of their MHC class I genes, in contrast to the high variability displayed by all other primates. Each Callitrichinae genus, therefore, expresses its own set of MHC class I genes, suggesting that an unusually high rate of turnover of loci occurs in this subfamily. The limited variability of MHC class I genes in the Callitrichinae is likely the result of the recent origin of these loci.

DNA polymerase active site is highly mutable: Evolutionary consequences

DNA polymerases contain active sites that are structurally superimposable and highly conserved in sequence. To assess the significance of this preservation and to determine the mutational burden that active sites can tolerate, we randomly mutated a stretch of 13 amino acids within the polymerase catalytic site (motif A) of Thermus aquaticus DNA polymerase I. After selection, by using genetic complementation, we obtained a library of approximately 8,000 active mutant DNA polymerases, of which 350 were sequenced and analyzed. This is the largest collection of physiologically active polymerase mutants. We find that all residues of motif A, except one (Asp-610), are mutable while preserving wild-type activity. A wide variety of amino acid substitutions were obtained at sites that are evolutionarily maintained, and conservative substitutions predominate at regions that stabilize tertiary structures. Several mutants exhibit unique properties, including DNA polymerase activity higher than the wild-type enzyme or the ability to incorporate ribonucleotide analogs. Bacteria dependent on these mutated polymerases for survival are fit to replicate repetitively. The high mutability of the polymerase active site in vivo and the ability to evolve altered enzymes may be required for survival in environments that demand increased mutagenesis. The inherent substitutability of the polymerase active site must be addressed relative to the constancy of nucleotide sequence found in nature.

Frequency of migrants and migratory activity are genetically correlated in a bird population: Evolutionary implications

Most migratory bird populations are composed of individuals that migrate and individuals that remain resident. While the role of ecological factors in maintaining this behavioral dimorphism has received much attention, the importance of genetic constraints on the evolution of avian migration has not yet been considered. Drawing on the recorded migratory activities of 775 blackcaps (Sylvia atricapilla) from a partially migratory population in southern France, we tested two alternative genetic models about the relationship between incidence and amount of migratory activity. The amount of migratory activity could be the continuous variable “underlying” the phenotypic expression of migratory urge, or, alternatively, the expression of both traits could be controlled by two separate genetic systems. The distributions of migratory activities in five different cohorts and the inheritance pattern derived from selective breeding experiments both indicate that incidence and amount of migratory activity are two aspects of one trait. Thus, all birds without measurable activity have activity levels at the low end of a continuous distribution, below the limit of expression or detection. The phenotypic dichotomy “migrant–nonmigrant” is caused by a threshold which may not be fixed but influenced both genetically and environmentally. This finding has profound implications for the evolution of migration: the transition from migratoriness to residency should not only be driven by selection favoring resident birds but also by selection for lower migratory activity. This potential for selection on two aspects, residency and migration distance, of the same trait may enable extremely rapid evolutionary changes to occur in migratory behavior.

Evolutionary theory predicts late-life mortality plateaus

Most demographic data indicate a roughly exponential increase in adult mortality with age, a phenomenon that has been explained in terms of a decline in the force of natural selection acting on age-specific mortality. Scattered demographic findings suggest the existence of a late-life mortality plateau in both humans and dipteran insects, seemingly at odds with both prior data and evolutionary theory. Extensions to the evolutionary theory of aging are developed which indicate that such late-life mortality plateaus are to be expected when enough late-life data are collected. This expanded theory predicts late-life mortality plateaus, with both antagonistic pleiotropy and mutation accumulation as driving population genetic mechanisms.

Evolutionary conservation of apoptosis mechanisms: Lepidopteran and baculoviral inhibitor of apoptosis proteins are inhibitors of mammalian caspase-9

We cloned a new inhibitor of apoptosis protein (IAP) homolog, SfIAP, from Spodoptera frugiperda Sf-21 cells, a host of insect baculoviruses. SfIAP contains two baculovirus IAP repeat domains followed by a RING domain. SfIAP has striking amino acid sequence similarity with baculoviral IAPs, CpIAP and OpIAP, suggesting that baculoviral IAPs may be host-derived genes. SfIAP and baculoviral CpIAP inhibit Bax but not Fas-induced apoptosis in human cells. Their apoptosis-suppressing activity in mammalian cells requires both baculovirus IAP repeat and RING domains. Further biochemical data suggest that SfIAP and CpIAP are specific inhibitors of mammalian caspase-9, the pinnacle caspase in the mitochondria/cytochrome c pathway for apoptosis, but are not inhibitors of downstream caspase-3 and caspase-7. Thus the mechanisms by which insect and baculoviral IAPs suppress apoptosis may involve inhibition of an insect caspase-9 homologue. Peptides representing the IAP-binding domain of the Drosophila cell death protein Grim abrogated human caspase suppression by SfIAP and CpIAP, implying evolutionary conservation of the functions of IAPs and their inhibitors.

Cryptocyanin, a crustacean molting protein: Evolutionary link with arthropod hemocyanins and insect hexamerins

Cryptocyanin, a copper-free hexameric protein in crab (Cancer magister) hemolymph, has been characterized and the amino acid sequence has been deduced from its cDNA. It is markedly similar in sequence, size, and structure to hemocyanin, the copper-containing oxygen-transport protein found in many arthropods. Cryptocyanin does not bind oxygen, however, and lacks three of the six highly conserved copper-binding histidine residues of hemocyanin. Cryptocyanin has no phenoloxidase activity, although a phenoloxidase is present in the hemolymph. The concentration of cryptocyanin in the hemolymph is closely coordinated with the molt cycle and reaches levels higher than hemocyanin during premolt. Cryptocyanin resembles insect hexamerins in the lack of copper, molt cycle patterns of biosynthesis, and potential contributions to the new exoskeleton. Phylogenetic analysis of sequence similarities between cryptocyanin and other members of the hemocyanin gene family shows that cryptocyanin is closely associated with crustacean hemocyanins and suggests that cryptocyanin arose as a result of a hemocyanin gene duplication. The presence of both hemocyanin and cryptocyanin in one animal provides an example of how insect hexamerins might have evolved from hemocyanin. Our results suggest that multiple members of the hemocyanin gene family—hemocyanin, cryptocyanin, phenoloxidase, and hexamerins—may participate in two vital functions of molting animals, oxygen binding and molting. Cryptocyanin may provide important molecular data to further investigate evolutionary relationships among all molting animals.

Evolutionary relationships among diverse bacteriophages and prophages: All the world’s a phage

We report DNA and predicted protein sequence similarities, implying homology, among genes of double-stranded DNA (dsDNA) bacteriophages and prophages spanning a broad phylogenetic range of host bacteria. The sequence matches reported here establish genetic connections, not always direct, among the lambdoid phages of Escherichia coli, phage φC31 of Streptomyces, phages of Mycobacterium, a previously unrecognized cryptic prophage, φflu, in the Haemophilus influenzae genome, and two small prophage-like elements, φRv1 and φRv2, in the genome of Mycobacterium tuberculosis. The results imply that these phage genes, and very possibly all of the dsDNA tailed phages, share common ancestry. We propose a model for the genetic structure and dynamics of the global phage population in which all dsDNA phage genomes are mosaics with access, by horizontal exchange, to a large common genetic pool but in which access to the gene pool is not uniform for all phage.

Toward functional genomics in bacteria: Analysis of gene expression in Escherichia coli from a bacterial artificial chromosome library of Bacillus cereus

As the study of microbes moves into the era of functional genomics, there is an increasing need for molecular tools for analysis of a wide diversity of microorganisms. Currently, biological study of many prokaryotes of agricultural, medical, and fundamental scientific interest is limited by the lack of adequate genetic tools. We report the application of the bacterial artificial chromosome (BAC) vector to prokaryotic biology as a powerful approach to address this need. We constructed a BAC library in Escherichia coli from genomic DNA of the Gram-positive bacterium Bacillus cereus. This library provides 5.75-fold coverage of the B. cereus genome, with an average insert size of 98 kb. To determine the extent of heterologous expression of B. cereus genes in the library, we screened it for expression of several B. cereus activities in the E. coli host. Clones expressing 6 of 10 activities tested were identified in the library, namely, ampicillin resistance, zwittermicin A resistance, esculin hydrolysis, hemolysis, orange pigment production, and lecithinase activity. We analyzed selected BAC clones genetically to identify rapidly specific B. cereus loci. These results suggest that BAC libraries will provide a powerful approach for studying gene expression from diverse prokaryotes.

Identification of eIF2Bγ and eIF2γ as cofactors of hepatitis C virus internal ribosome entry site-mediated translation using a functional genomics approach

The 5′-untranslated region of hepatitis C virus (HCV) is highly conserved, folds into a complex secondary structure, and functions as an internal ribosome entry site (IRES) to initiate translation of HCV proteins. We have developed a selection system based on a randomized hairpin ribozyme gene library to identify cellular factors involved in HCV IRES function. A retroviral vector ribozyme library with randomized target recognition sequences was introduced into HeLa cells, stably expressing a bicistronic construct encoding the hygromycin B phosphotransferase gene and the herpes simplex virus thymidine kinase gene (HSV-tk). Translation of the HSV-tk gene was mediated by the HCV IRES. Cells expressing ribozymes that inhibit HCV IRES-mediated translation of HSV-tk were selected via their resistance to both ganciclovir and hygromycin B. Two ribozymes reproducibly conferred the ganciclovir-resistant phenotype and were shown to inhibit IRES-mediated translation of HCV core protein but did not inhibit cap-dependent protein translation or cell growth. The functional targets of these ribozymes were identified as the gamma subunits of human eukaryotic initiation factors 2B (eIF2Bγ) and 2 (eIF2γ), respectively. The involvement of eIF2Bγ and eIF2γ in HCV IRES-mediated translation was further validated by ribozymes directed against additional sites within the mRNAs of these genes. In addition to leading to the identification of cellular IRES cofactors, ribozymes obtained from this cellular selection system could be directly used to specifically inhibit HCV viral translation, thereby facilitating the development of new antiviral strategies for HCV infection.

Molecular phylogenetic analysis of evolutionary trends in stonefly wing structure and locomotor behavior

Insects in the order Plecoptera (stoneflies) use a form of two-dimensional aerodynamic locomotion called surface skimming to move across water surfaces. Because their weight is supported by water, skimmers can achieve effective aerodynamic locomotion even with small wings and weak flight muscles. These mechanical features stimulated the hypothesis that surface skimming may have been an intermediate stage in the evolution of insect flight, which has perhaps been retained in certain modern stoneflies. Here we present a phylogeny of Plecoptera based on nucleotide sequence data from the small subunit rRNA (18S) gene. By mapping locomotor behavior and wing structural data onto the phylogeny, we distinguish between the competing hypotheses that skimming is a retained ancestral trait or, alternatively, a relatively recent loss of flight. Our results show that basal stoneflies are surface skimmers, and that various forms of surface skimming are distributed widely across the plecopteran phylogeny. Stonefly wings show evolutionary trends in the number of cross veins and the thickness of the cuticle of the longitudinal veins that are consistent with elaboration and diversification of flight-related traits. These data support the hypothesis that the first stoneflies were surface skimmers, and that wing structures important for aerial flight have become elaborated and more diverse during the radiation of modern stoneflies.

A chemical genomics approach toward understanding the global functions of the target of rapamycin protein (TOR)

The target of rapamycin protein (TOR) is a highly conserved ataxia telangiectasia-related protein kinase essential for cell growth. Emerging evidence indicates that TOR signaling is highly complex and is involved in a variety of cellular processes. To understand its general functions, we took a chemical genomics approach to explore the genetic interaction between TOR and other yeast genes on a genomic scale. In this study, the rapamycin sensitivity of individual deletion mutants generated by the Saccharomyces Genome Deletion Project was systematically measured. Our results provide a global view of the rapamycin-sensitive functions of TOR. In contrast to conventional genetic analysis, this approach offers a simple and thorough analysis of genetic interaction on a genomic scale and measures genetic interaction at different possible levels. It can be used to study the functions of other drug targets and to identify novel protein components of a conserved core biological process such as DNA damage checkpoint/repair that is interfered with by a cell-permeable chemical compound.

Assessing sequence comparison methods with reliable structurally identified distant evolutionary relationships

Pairwise sequence comparison methods have been assessed using proteins whose relationships are known reliably from their structures and functions, as described in the scop database [Murzin, A. G., Brenner, S. E., Hubbard, T. & Chothia C. (1995) J. Mol. Biol. 247, 536–540]. The evaluation tested the programs blast [Altschul, S. F., Gish, W., Miller, W., Myers, E. W. & Lipman, D. J. (1990). J. Mol. Biol. 215, 403–410], wu-blast2 [Altschul, S. F. & Gish, W. (1996) Methods Enzymol. 266, 460–480], fasta [Pearson, W. R. & Lipman, D. J. (1988) Proc. Natl. Acad. Sci. USA 85, 2444–2448], and ssearch [Smith, T. F. & Waterman, M. S. (1981) J. Mol. Biol. 147, 195–197] and their scoring schemes. The error rate of all algorithms is greatly reduced by using statistical scores to evaluate matches rather than percentage identity or raw scores. The E-value statistical scores of ssearch and fasta are reliable: the number of false positives found in our tests agrees well with the scores reported. However, the P-values reported by blast and wu-blast2 exaggerate significance by orders of magnitude. ssearch, fasta ktup = 1, and wu-blast2 perform best, and they are capable of detecting almost all relationships between proteins whose sequence identities are >30%. For more distantly related proteins, they do much less well; only one-half of the relationships between proteins with 20–30% identity are found. Because many homologs have low sequence similarity, most distant relationships cannot be detected by any pairwise comparison method; however, those which are identified may be used with confidence.

Proposed mechanism for stability of proteins to evolutionary mutations

It is shown that the sequence-ordering tendencies induced by design into different fast-folding, thermally stable native structures interfere. This interference results in a type of quasiorthogonality between optimal native structures, which divides sequence space into fast-folding, thermally stable families surrounded by slow-folding, low stability shells. A concrete example of this effect is provided by using a simple α carbon type model in which a complete correspondence is established between sequence and structure. It is speculated that gaps can occur in the space of protein-like sequences separating the sequence families and resulting in a mechanism for stability and diversity of protein sequence information.

Structure-based assignment of the biochemical function of a hypothetical protein: A test case of structural genomics

Many small bacterial, archaebacterial, and eukaryotic genomes have been sequenced, and the larger eukaryotic genomes are predicted to be completely sequenced within the next decade. In all genomes sequenced to date, a large portion of these organisms’ predicted protein coding regions encode polypeptides of unknown biochemical, biophysical, and/or cellular functions. Three-dimensional structures of these proteins may suggest biochemical or biophysical functions. Here we report the crystal structure of one such protein, MJ0577, from a hyperthermophile, Methanococcus jannaschii, at 1.7-Å resolution. The structure contains a bound ATP, suggesting MJ0577 is an ATPase or an ATP-mediated molecular switch, which we confirm by biochemical experiments. Furthermore, the structure reveals different ATP binding motifs that are shared among many homologous hypothetical proteins in this family. This result indicates that structure-based assignment of molecular function is a viable approach for the large-scale biochemical assignment of proteins and for discovering new motifs, a basic premise of structural genomics.