On constructing folding heteropolymers.

Simplified models of the protein-folding process have led to valuable insights into the generic properties of the folding of heteropolymers. On the basis of theoretical arguments, Shakhnovich and Gutin [(1993) Proc. Natl. Acad. Sci. USA 90, 7195-7199] have proposed a specific method to generate folding sequences for one of these. Here we present a model of folding in heteropolymers that is comparable in simplicity but different in spirit to the one studied by Shakhnovich and Gutin. In our model, the proposed recipe for constructing folding sequence fails. We find that, as a rule, the construction of folding sequences is impossible to achieve by looking at the native conformation only. Rather, competing conformations have to be taken into account too. An evolutionary algorithm that generates folding sequences by optimizing both stability of the native state and folding time is described. Remarkably, this algorithm produces, among others, sequences that fold reproducibly to metastable states.

Calibrating bacterial evolution

Attempts to calibrate bacterial evolution have relied on the assumption that rates of molecular sequence divergence in bacteria are similar to those of higher eukaryotes, or to those of the few bacterial taxa for which ancestors can be reliably dated from ecological or geological evidence. Despite similarities in the substitution rates estimated for some lineages, comparisons of the relative rates of evolution at different classes of nucleotide sites indicate no basis for their universal application to all bacteria. However, there is evidence that bacteria have a constant genome-wide mutation rate on an evolutionary time scale but that this rate differs dramatically from the rate estimated by experimental methods.

A mammalian mitochondrial drug receptor functions as a bacterial “oxygen” sensor

The rat mitochondrial outer membrane-localized benzodiazepine receptor (MBR) was expressed in wild-type and TspO− (tryptophan-rich sensory protein) strains of the facultative photoheterotroph, Rhodobacter sphaeroides 2.4.1, and was shown to retain its structure within the bacterial outer membrane as assayed by its binding properties with a variety of MBR ligands. Functionally, it was able to substitute for TspO by negatively regulating the expression of photosynthesis genes in response to oxygen. This effect was reversed pharmacologically with the MBR ligand PK11195. These results suggest a close evolutionary and functional relationship between the bacterial TspO and the MBR. This relationship provides further support for the origin of the mammalian mitochondrion from a “photosynthetic” precursor. Finally, these findings provide novel insights into the physiological role that has been obscure for the MBR in situ.

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.

Plants and animals share functionally common bacterial virulence factors

By exploiting the ability of Pseudomonas aeruginosa to infect a variety of vertebrate and nonvertebrate hosts, we have developed model systems that use plants and nematodes as adjuncts to mammalian models to help elucidate the molecular basis of P. aeruginosa pathogenesis. Our studies reveal a remarkable degree of conservation in the virulence mechanisms used by P. aeruginosa to infect hosts of divergent evolutionary origins.

Evolutionary genomics of Staphylococcus aureus: Insights into the origin of methicillin-resistant strains and the toxic shock syndrome epidemic

An emerging theme in medical microbiology is that extensive variation exists in gene content among strains of many pathogenic bacterial species. However, this topic has not been investigated on a genome scale with strains recovered from patients with well-defined clinical conditions. Staphylococcus aureus is a major human pathogen and also causes economically important infections in cows and sheep. A DNA microarray representing >90% of the S. aureus genome was used to characterize genomic diversity, evolutionary relationships, and virulence gene distribution among 36 strains of divergent clonal lineages, including methicillin-resistant strains and organisms causing toxic shock syndrome. Genetic variation in S. aureus is very extensive, with ≈22% of the genome comprised of dispensable genetic material. Eighteen large regions of difference were identified, and 10 of these regions have genes that encode putative virulence factors or proteins mediating antibiotic resistance. We find that lateral gene transfer has played a fundamental role in the evolution of S. aureus. The mec gene has been horizontally transferred into distinct S. aureus chromosomal backgrounds at least five times, demonstrating that methicillin-resistant strains have evolved multiple independent times, rather than from a single ancestral strain. This finding resolves a long-standing controversy in S. aureus research. The epidemic of toxic shock syndrome that occurred in the 1970s was caused by a change in the host environment, rather than rapid geographic dissemination of a new hypervirulent strain. DNA microarray analysis of large samples of clinically characterized strains provides broad insights into evolution, pathogenesis, and disease emergence.

A common evolutionary origin for mitochondria and hydrogenosomes.

Trichomonads are among the earliest eukaryotes to diverge from the main line of eukaryotic descent. Keeping with their ancient nature, these facultative anaerobic protists lack two "hallmark" organelles found in most eukaryotes: mitochondria and peroxisomes. Trichomonads do, however, contain an unusual organelle involved in carbohydrate metabolism called the hydrogenosome. Like mitochondria, hydrogenosomes are double-membrane bounded organelles that produce ATP using pyruvate as the primary substrate. Hydrogenosomes are, however, markedly different from mitochondria as they lack DNA, cytochromes and the citric acid cycle. Instead, they contain enzymes typically found in anaerobic bacteria and are capable of producing molecular hydrogen. We show here that hydrogenosomes contain heat shock proteins, Hsp70, Hsp60, and Hsp10, with signature sequences that are conserved only in mitochondrial and alpha-Gram-negative purple bacterial Hsps. Biochemical analysis of hydrogenosomal Hsp60 shows that the mature protein isolated from the organelle lacks a short, N-terminal sequence, similar to that observed for most nuclear-encoded mitochondrial matrix proteins. Moreover, phylogenetic analyses of hydrogenosomal Hsp70, Hsp60, and Hsp10 show that these proteins branch within a monophyletic group composed exclusively of mitochondrial homologues. These data establish that mitochondria and hydrogenosomes have a common eubacterial ancestor and imply that the earliest-branching eukaryotes contained the endosymbiont that gave rise to mitochondria in higher eukaryotes.

cAMP-dependent SOS induction and mutagenesis in resting bacterial populations.

The inducible SOS system increases the survival of bacteria exposed to DNA-damaging agents by increasing the capacity of error-free and error-prone DNA repair systems. The inducible mutator effect is expected to contribute to the adaptation of bacterial populations to these adverse life conditions by increasing their genetic variability. The evolutionary impact of the SOS system would be even greater if it was also induced under conditions common in nature, such as in resting bacterial populations. The results presented here show that SOS induction and mutagenesis do occur in bacteria in aging colonies on agar plates. The observed SOS induction and mutagenesis are controlled by the LexA repressor and are RecA- and cAMP-dependent.

Sequence similarity analysis of Escherichia coli proteins: functional and evolutionary implications.

A computer analysis of 2328 protein sequences comprising about 60% of the Escherichia coli gene products was performed using methods for database screening with individual sequences and alignment blocks. A high fraction of E. coli proteins--86%--shows significant sequence similarity to other proteins in current databases; about 70% show conservation at least at the level of distantly related bacteria, and about 40% contain ancient conserved regions (ACRs) shared with eukaryotic or Archaeal proteins. For > 90% of the E. coli proteins, either functional information or sequence similarity, or both, are available. Forty-six percent of the E. coli proteins belong to 299 clusters of paralogs (intraspecies homologs) defined on the basis of pairwise similarity. Another 10% could be included in 70 superclusters using motif detection methods. The majority of the clusters contain only two to four members. In contrast, nearly 25% of all E. coli proteins belong to the four largest superclusters--namely, permeases, ATPases and GTPases with the conserved "Walker-type" motif, helix-turn-helix regulatory proteins, and NAD(FAD)-binding proteins. We conclude that bacterial protein sequences generally are highly conserved in evolution, with about 50% of all ACR-containing protein families represented among the E. coli gene products. With the current sequence databases and methods of their screening, computer analysis yields useful information on the functions and evolutionary relationships of the vast majority of genes in a bacterial genome. Sequence similarity with E. coli proteins allows the prediction of functions for a number of important eukaryotic genes, including several whose products are implicated in human diseases.

Calculating the probability of multitaxon evolutionary trees: bootstrappers Gambit.

The reconstruction of multitaxon trees from molecular sequences is confounded by the variety of algorithms and criteria used to evaluate trees, making it difficult to compare the results of different analyses. A global method of multitaxon phylogenetic reconstruction described here, Bootstrappers Gambit, can be used with any four-taxon algorithm, including distance, maximum likelihood, and parsimony methods. It incorporates a Bayesian-Jeffreys'-bootstrap analysis to provide a uniform probability-based criterion for comparing the results from diverse algorithms. To examine the usefulness of the method, the origin of the eukaryotes has been investigated by the analysis of ribosomal small subunit RNA sequences. Three common algorithms (paralinear distances, Jukes-Cantor distances, and Kimura distances) support the eocyte topology, whereas one (maximum parsimony) supports the archaebacterial topology, suggesting that the eocyte prokaryotes are the closest prokaryotic relatives of the eukaryotes.

Considerations on the folding topology and evolutionary origin of cadherin domains.

Cell-cell adhesion in zonula adherens and desmosomal junctions is mediated by cadherins, and recent crystal structures of the first domain from murine N-cadherin provide a plausible molecular basis for this adhesive action. A structure-based sequence analysis of this adhesive domain indicates that its fold is common to all extracellular cadherin domains. The cadherin folding topology is also shown to be similar to immunoglobulin-like domains and to other Greek-key beta-sandwich structures, as diverse as domains from plant cytochromes, bacterial cellulases, and eukaryotic transcription factors. Sequence similarities between cadherins and these other molecules are very low, however, and intron patterns are also different. On balance, independent origins for a favorable folding topology seem more likely than evolutionary divergence from an ancestor common to cadherins and immunoglobulins.