On roads not taken in the evolution of protein catalysts: Antibody steroid isomerases that use an enamine mechanism

Reactive immunization has emerged as a new tool for the study of biological catalysis. A powerful application resulted in catalytic antibodies that use an enamine mechanism akin to that used by the class I aldolases. With regard to the evolution of enzyme mechanisms, we investigated the utility of an enamine pathway for the allylic rearrangement exemplified by Δ5-3-ketosteroid isomerase (KSI; EC 5.3.3.1). Our aldolase antibodies were found to catalyze the isomerization of both steroid model compounds and steroids. The kinetic and chemical studies showed that the antibodies afforded rate accelerations up to a factor of 104 by means of an enamine mechanism in which imine formation was the rate-determining step. In light of our observations and the enzyme studies by other workers, we suggest that an enamine pathway could have been an early, viable KSI mechanism. Although this pathway is amenable to optimization for increased catalytic power, it appears that certain factors precluded its evolution in known KSI enzymes.

In vitro evolution of horse heart myoglobin to increase peroxidase activity

Random mutagenesis and screening for enzymatic activity has been used to engineer horse heart myoglobin to enhance its intrinsic peroxidase activity. A chemically synthesized gene encoding horse heart myoglobin was subjected to successive cycles of PCR random mutagenesis. The mutated myoglobin gene was expressed in Escherichia coli LE392, and the variants were screened for peroxidase activity with a plate assay. Four cycles of mutagenesis and screening produced a series of single, double, triple, and quadruple variants with enhanced peroxidase activity. Steady-state kinetics analysis demonstrated that the quadruple variant T39I/K45D/F46L/I107F exhibits peroxidase activity significantly greater than that of the wild-type protein with k1 (for H2O2 oxidation of metmyoglobin) of 1.34 × 104 M−1 s−1 (≈25-fold that of wild-type myoglobin) and k3 [for reducing the substrate (2, 2′-azino-di-(3-ethyl)benzthiazoline-6-sulfonic acid] of 1.4 × 106 M−1 s−1 (1.6-fold that of wild-type myoglobin). Thermal stability of these variants as measured with circular dichroism spectroscopy demonstrated that the Tm of the quadruple variant is decreased only slightly compared with wild-type (74.1°C vs. 76.5°C). The rate constants for binding of dioxygen exhibited by the quadruple variant are identical to the those observed for wild-type myoglobin (kon, 22.2 × 10−6 M−1 s−1 vs. 22.3 × 10−6 M−1 s−1; koff, 24.3 s−1 vs. 24.2 s−1; KO2, 0.91 × 10−6 M−1 vs. 0.92 × 10−6 M−1). The affinity of the quadruple variant for CO is increased slightly (kon, 0.90 × 10−6 M−1s−1 vs. 0.51 × 10−6 M−1s−1; koff, 5.08 s−1 vs. 3.51 s−1; KCO, 1.77 × 10−7 M−1 vs. 1.45 × 10−7 M−1). All four substitutions are in the heme pocket and within 5 Å of the heme group.

Transgenic tobacco plants with reduced capability to detoxify reactive oxygen intermediates are hyperresponsive to pathogen infection

Reactive oxygen intermediates (ROI) play a critical role in the defense of plants against invading pathogens. Produced during the “oxidative burst,” they are thought to activate programmed cell death (PCD) and induce antimicrobial defenses such as pathogenesis-related proteins. It was shown recently that during the interaction of plants with pathogens, the expression of ROI-detoxifying enzymes such as ascorbate peroxidase (APX) and catalase (CAT) is suppressed. It was suggested that this suppression, occurring upon pathogen recognition and coinciding with an enhanced rate of ROI production, plays a key role in elevating cellular ROI levels, thereby potentiating the induction of PCD and other defenses. To examine the relationship between the suppression of antioxidative mechanisms and the induction of PCD and other defenses during pathogen attack, we studied the interaction between transgenic antisense tobacco plants with reduced APX or CAT and a bacterial pathogen that triggers the hypersensitive response. Transgenic plants with reduced capability to detoxify ROI (i.e., antisense APX or CAT) were found to be hyperresponsive to pathogen attack. They activated PCD in response to low amounts of pathogens that did not trigger the activation of PCD in control plants. Our findings support the hypothesis that suppression of ROI-scavenging enzymes during the hypersensitive response plays an important role in enhancing pathogen-induced PCD.

The carbamate kinase-like carbamoyl phosphate synthetase of the hyperthermophilic archaeon Pyrococcus furiosus, a missing link in the evolution of carbamoyl phosphate biosynthesis

Microbial carbamoyl phosphate synthetases (CPS) use glutamine as nitrogen donor and are composed of two subunits (or domains), one exhibiting glutaminase activity, the other able to synthesize carbamoyl phosphate (CP) from bicarbonate, ATP, and ammonia. The pseudodimeric organization of this synthetase suggested that it has evolved by duplication of a smaller kinase, possibly a carbamate kinase (CK). In contrast to other prokaryotes the hyperthermophilic archaeon Pyrococcus furiosus was found to synthesize CP by using ammonia and not glutamine. We have purified the cognate enzyme and found it to be a dimer of two identical subunits of Mr 32,000. Its thermostability is considerable, 50% activity being retained after 1 h at 100°C or 3 h at 95°C. The corresponding gene was cloned by PCR and found to present about 50% amino acid identity with known CKs. The stoichiometry of the reaction (two ATP consumed per CP synthesized) and the ability of the enzyme to catalyze at high rate a bicarbonate-dependent ATPase reaction however clearly distinguish P. furiosus CPS from ordinary CKs. Thus the CPS of P. furiosus could represent a primeval step in the evolution of CPS from CK. Our results suggest that the first event in this evolution was the emergence of a primeval synthetase composed of subunits able to synthesize both carboxyphosphate and CP; this step would have preceded the duplication assumed to have generated the two subdomains of modern CPSs. The gene coding for this CK-like CPS was called cpkA.

An epidemic of tuberculosis with a high rate of tuberculin anergy among a population previously unexposed to tuberculosis, the Yanomami Indians of the Brazilian Amazon

A survey of an emerging tuberculosis epidemic among the Yanomami Indians of the Amazonian rain forest provided a unique opportunity to study the impact of tuberculosis on a population isolated from contact with the tubercle bacillus for millennia until the mid-1960s. Within the Yanomami population, an extraordinary high prevalence of active tuberculosis (6.4% of 625 individuals clinically examined) was observed, indicating a high susceptibility to disease, even among bacille Calmette–Guérin-vaccinated individuals. Observational studies on cell-mediated and humoral immune responses of the Yanomami Indians compared with contemporary residents of the region suggest profound differences in immunological responsiveness to Mycobacterium tuberculosis infection. Among the Yanomami, a very high prevalence of tuberculin skin test anergy was found. Of patients with active tuberculosis, 46% had purified protein derivative of tuberculosis reactions <10 mm; similarly 58% of recent bacillus Calmette–Guérin vaccines exhibited skin test reactions <5 mm. The Yanomami also had higher titers of antibodies against M. tuberculosis glycolipid antigens (>70%) than the control subjects comprised of Brazilians of European descent (14%). The antibodies were mostly of the IgM isotype. Among the tuberculosis patients who also produced IgG antibodies, the titers of IgG4 were significantly higher among the Yanomami than in the control population. Although it was not possible to analyze T-cell responses or patterns of lymphokine production in vitro because of the remoteness of the villages from laboratory facilities, the results suggest that the first encounter of the Yanomami Indian population with tuberculosis engenders a diminished cell-mediated immune response and an increased production antibody responses, relative to other populations with extensive previous contact with the pathogen. These findings suggest that tuberculosis may represent a powerful selective pressure on human evolution that over centuries has shaped the nature of human immune responses to infection.

The mutation rate and cancer

The selection of advantageous mutations underlies tumorigenesis. The growth of a tumor is therefore a form of evolution at the somatic level, in which the population is comprised of individual cells within the tumor. Models of tumorigenesis have considered the relative importance of mutation and selection. We show that selection is more important than an increased mutation rate in the growth of a tumor. Some cancers may acquire a “mutator phenotype,” probably leading to faster growth, but mutator phenotypes are not necessary for carcinogenesis.

Directed evolution of antibody fragments with monovalent femtomolar antigen-binding affinity

Single-chain antibody mutants have been evolved in vitro with antigen-binding equilibrium dissociation constant Kd = 48 fM and slower dissociation kinetics (half-time > 5 days) than those for the streptavidin–biotin complex. These mutants possess the highest monovalent ligand-binding affinity yet reported for an engineered protein by over two orders of magnitude. Optimal kinetic screening of randomly mutagenized libraries of 105–107 yeast surface-displayed antibodies enabled a >1,000-fold decrease in the rate of dissociation after four cycles of affinity mutagenesis and screening. The consensus mutations are generally nonconservative by comparison with naturally occurring mouse Fv sequences and with residues that do not contact the fluorescein antigen in the wild-type complex. The existence of these mutants demonstrates that the antibody Fv architecture is not intrinsically responsible for an antigen-binding affinity ceiling during in vivo affinity maturation.

The origin of atmospheric oxygen on Earth: The innovation of oxygenic photosynthesis

The evolution of O2-producing cyanobacteria that use water as terminal reductant transformed Earth's atmosphere to one suitable for the evolution of aerobic metabolism and complex life. The innovation of water oxidation freed photosynthesis to invade new environments and visibly changed the face of the Earth. We offer a new hypothesis for how this process evolved, which identifies two critical roles for carbon dioxide in the Archean period. First, we present a thermodynamic analysis showing that bicarbonate (formed by dissolution of CO2) is a more efficient alternative substrate than water for O2 production by oxygenic phototrophs. This analysis clarifies the origin of the long debated “bicarbonate effect” on photosynthetic O2 production. We propose that bicarbonate was the thermodynamically preferred reductant before water in the evolution of oxygenic photosynthesis. Second, we have examined the speciation of manganese(II) and bicarbonate in water, and find that they form Mn-bicarbonate clusters as the major species under conditions that model the chemistry of the Archean sea. These clusters have been found to be highly efficient precursors for the assembly of the tetramanganese-oxide core of the water-oxidizing enzyme during biogenesis. We show that these clusters can be oxidized at electrochemical potentials that are accessible to anoxygenic phototrophs and thus the most likely building blocks for assembly of the first O2 evolving photoreaction center, most likely originating from green nonsulfur bacteria before the evolution of cyanobacteria.

Long term trends in the evolution of H(3) HA1 human influenza type A

We have studied the HA1 domain of 254 human influenza A(H3N2) virus genes for clues that might help identify characteristics of hemagglutinins (HAs) of circulating strains that are predictive of that strain’s epidemic potential. Our preliminary findings include the following. (i) The most parsimonious tree found requires 1,260 substitutions of which 712 are silent and 548 are replacement substitutions. (ii) The HA1 portion of the HA gene is evolving at a rate of 5.7 nucleotide substitutions/year or 5.7 × 10−3 substitutions/site per year. (iii) The replacement substitutions are distributed randomly across the three positions of the codon when allowance is made for the number of ways each codon can change the encoded amino acid. (iv) The replacement substitutions are not distributed randomly over the branches of the tree, there being 2.2 times more changes per tip branch than for non-tip branches. This result is independent of how the virus was amplified (egg grown or kidney cell grown) prior to sequencing or if sequencing was carried out directly on the original clinical specimen by PCR. (v) These excess changes on the tip branches are probably the result of a bias in the choice of strains to sequence and the detection of deleterious mutations that had not yet been removed by negative selection. (vi) There are six hypervariable codons accumulating replacement substitutions at an average rate that is 7.2 times that of the other varied codons. (vii) The number of variable codons in the trunk branches (the winners of the competitive race against the immune system) is 47 ± 5, significantly fewer than in the twigs (90 ± 7), which in turn is significantly fewer variable codons than in tip branches (175 ± 8). (viii) A minimum of one of every 12 branches has nodes at opposite ends representing viruses that reside on different continents. This is, however, no more than would be expected if one were to randomly reassign the continent of origin of the isolates. (ix) Of 99 codons with at least four mutations, 31 have ratios of non-silent to silent changes with probabilities less than 0.05 of occurring by chance, and 14 of those have probabilities <0.005. These observations strongly support positive Darwinian selection. We suggest that the small number of variable positions along the successful trunk lineage, together with knowledge of the codons that have shown positive selection, may provide clues that permit an improved prediction of which strains will cause epidemics and therefore should be used for vaccine production.

The evolution of plant nuclear genes

We analyze the evolutionary dynamics of three of the best-studied plant nuclear multigene families. The data analyzed derive from the genes that encode the small subunit of ribulose-1,5-bisphosphate carboxylase (rbcS), the gene family that encodes the enzyme chalcone synthase (Chs), and the gene family that encodes alcohol dehydrogenases (Adh). In addition, we consider the limited evolutionary data available on plant transposable elements. New Chs and rbcS genes appear to be recruited at about 10 times the rate estimated for Adh genes, and this is correlated with a much smaller average gene family size for Adh genes. In addition, duplication and divergence in function appears to be relatively common for Chs genes in flowering plant evolution. Analyses of synonymous nucleotide substitution rates for Adh genes in monocots reject a linear relationship with clock time. Replacement substitution rates vary with time in a complex fashion, which suggests that adaptive evolution has played an important role in driving divergence following gene duplication events. Molecular population genetic studies of Adh and Chs genes reveal high levels of molecular diversity within species. These studies also reveal that inter- and intralocus recombination are important forces in the generation allelic novelties. Moreover, illegitimate recombination events appear to be an important factor in transposable element loss in plants. When we consider the recruitment and loss of new gene copies, the generation of allelic diversity within plant species, and ectopic exchange among transposable elements, we conclude that recombination is a pervasive force at all levels of plant evolution.

Atmospheric oxygen over Phanerozoic time

It is quite possible that the level of atmospheric oxygen has varied (roughly between 15 and 30% O2) over the past 550 million years. This variation is suggested by modeling of the carbon and sulfur cycles, by the excessive sediment burial of organic matter that accompanied the advent of large vascular land plants, and by recent physiological studies that relate to biological evolution.

Blood flow and oxygen delivery to human brain during functional activity: Theoretical modeling and experimental data

Coupling of cerebral blood flow (CBF) and cerebral metabolic rate for oxygen (CMRO2) in physiologically activated brain states remains the subject of debates. Recently it was suggested that CBF is tightly coupled to oxidative metabolism in a nonlinear fashion. As part of this hypothesis, mathematical models of oxygen delivery to the brain have been described in which disproportionately large increases in CBF are necessary to sustain even small increases in CMRO2 during activation. We have explored the coupling of CBF and oxygen delivery by using two complementary methods. First, a more complex mathematical model was tested that differs from those recently described in that no assumptions were made regarding tissue oxygen level. Second, [15O] water CBF positron emission tomography (PET) studies in nine healthy subjects were conducted during states of visual activation and hypoxia to examine the relationship of CBF and oxygen delivery. In contrast to previous reports, our model showed adequate tissue levels of oxygen could be maintained without the need for increased CBF or oxygen delivery. Similarly, the PET studies demonstrated that the regional increase in CBF during visual activation was not affected by hypoxia. These findings strongly indicate that the increase in CBF associated with physiological activation is regulated by factors other than local requirements in oxygen.

Measuring genome evolution

The determination of complete genome sequences provides us with an opportunity to describe and analyze evolution at the comprehensive level of genomes. Here we compare nine genomes with respect to their protein coding genes at two levels: (i) we compare genomes as “bags of genes” and measure the fraction of orthologs shared between genomes and (ii) we quantify correlations between genes with respect to their relative positions in genomes. Distances between the genomes are related to their divergence times, measured as the number of amino acid substitutions per site in a set of 34 orthologous genes that are shared among all the genomes compared. We establish a hierarchy of rates at which genomes have changed during evolution. Protein sequence identity is the most conserved, followed by the complement of genes within the genome. Next is the degree of conservation of the order of genes, whereas gene regulation appears to evolve at the highest rate. Finally, we show that some genomes are more highly organized than others: they show a higher degree of the clustering of genes that have orthologs in other genomes.

Molecular clock or erratic evolution? A tale of two genes.

We have investigated the evolution of glycerol-3-phosphate dehydrogenase (Gpdh). The rate of amino acid replacements is 1 x 10(-10)/site/year when Drosophila species are compared. The rate is 2.7 times greater when Drosophila and Chymomyza species are compared; and about 5 times greater when any of those species are compared with the medfly Ceratitis capitata. This rate of 5 x 10(-10)/site/year is also the rate observed in comparisons between mammals, or between different animal phyla, or between the three multicellular kingdoms. We have also studied the evolution of Cu,Zn superoxide dismutase (Sod). The rate of amino acid replacements is about 17 x 10(-10)/site/year when comparisons are made between dipterans or between mammals, but only 5 x 10(-10) when animal phyla are compared, and only 3 x 10(-10) when the multicellular kingdoms are compared. The apparent decrease by about a factor of 5 in the rate of SOD evolution as the divergence between species increases can be consistent with the molecular clock hypothesis by assuming the covarion hypothesis (namely, that the number of amino acids that can change is constant, but the set of such amino acids changes from time to time and from lineage to lineage). However, we know of no model consistent with the molecular clock hypothesis that would account for the increase in the rate of GPDH evolution as the divergence between species increases.

Molecular phylogeny analysis of fiddler crabs: test of the hypothesis of increasing behavioral complexity in evolution.

The current phylogenetic hypothesis for the evolution and biogeography of fiddler crabs relies on the assumption that complex behavioral traits are assumed to also be evolutionary derived. Indo-west Pacific fiddler crabs have simpler reproductive social behavior and are more marine and were thought to be ancestral to the more behaviorally complex and more terrestrial American species. It was also hypothesized that the evolution of more complex social and reproductive behavior was associated with the colonization of the higher intertidal zones. Our phylogenetic analysis, based upon a set of independent molecular characters, however, demonstrates how widely entrenched ideas about evolution and biogeography led to a reasonable, but apparently incorrect, conclusion about the evolutionary trends within this pantropical group of crustaceans. Species bearing the set of "derived traits" are phylogenetically ancestral, suggesting an alternative evolutionary scenario: the evolution of reproductive behavioral complexity in fiddler crabs may have arisen multiple times during their evolution. The evolution of behavioral complexity may have arisen by coopting of a series of other adaptations for high intertidal living and antipredator escape. A calibration of rates of molecular evolution from populations on either side of the Isthmus of Panama suggest a sequence divergence rate for 16S rRNA of 0.9% per million years. The divergence between the ancestral clade and derived forms is estimated to be approximately 22 million years ago, whereas the divergence between the American and Indo-west Pacific is estimated to be approximately 17 million years ago.