Peptide design by artificial neural networks and computer-based evolutionary search

A technique for systematic peptide variation by a combination of rational and evolutionary approaches is presented. The design scheme consists of five consecutive steps: (i) identification of a “seed peptide” with a desired activity, (ii) generation of variants selected from a physicochemical space around the seed peptide, (iii) synthesis and testing of this biased library, (iv) modeling of a quantitative sequence-activity relationship by an artificial neural network, and (v) de novo design by a computer-based evolutionary search in sequence space using the trained neural network as the fitness function. This strategy was successfully applied to the identification of novel peptides that fully prevent the positive chronotropic effect of anti-β1-adrenoreceptor autoantibodies from the serum of patients with dilated cardiomyopathy. The seed peptide, comprising 10 residues, was derived by epitope mapping from an extracellular loop of human β1-adrenoreceptor. A set of 90 peptides was synthesized and tested to provide training data for neural network development. De novo design revealed peptides with desired activities that do not match the seed peptide sequence. These results demonstrate that computer-based evolutionary searches can generate novel peptides with substantial biological activity.

A structural view of evolutionary divergence

Two directed evolution experiments on p-nitrobenzyl esterase yielded one enzyme with a 100-fold increased activity in aqueous-organic solvents and another with a 17°C increase in thermostability. Structures of the wild type and its organophilic and thermophilic counterparts are presented at resolutions of 1.5 Å, 1.6 Å, and 2.0 Å, respectively. These structures identify groups of interacting mutations and demonstrate how directed evolution can traverse complex fitness landscapes. Early-generation mutations stabilize flexible loops not visible in the wild-type structure and set the stage for further beneficial mutations in later generations. The mutations exert their influence on the esterase structure over large distances, in a manner that would be difficult to predict. The loops with the largest structural changes generally are not the sites of mutations. Similarly, none of the seven amino acid substitutions in the organophile are in the active site, even though the enzyme experiences significant changes in the organization of this site. In addition to reduction of surface loop flexibility, thermostability in the evolved esterase results from altered core packing, helix stabilization, and the acquisition of surface salt bridges, in agreement with other comparative studies of mesophilic and thermophilic enzymes. Crystallographic analysis of the wild type and its evolved counterparts reveals networks of mutations that collectively reorganize the active site. Interestingly, the changes that led to diversity within the α/β hydrolase enzyme family and the reorganization seen in this study result from main-chain movements.

Evidence of evolutionary up-regulation of the single active X chromosome in mammals based on Clc4 expression levels in Mus spretus and Mus musculus

Previous studies have shown that the chloride channel gene Clc4 is X-linked and subject to X inactivation in Mus spretus, but that the same gene is autosomal in laboratory strains of mice. This exception to the conservation of linkage of the X chromosome in one of two interfertile mouse species was exploited to compare expression of Clc4 from the X chromosome to that from the autosome. Clc4 was found to be highly expressed in brain tissues of both mouse species. Quantitative analyses of species-specific expression of Clc4 in brain tissues from mice resulting from M. spretus × laboratory strain crosses, demonstrate that each autosomal locus has half the level of Clc4 expression as compared with the single active X-linked locus. In contrast expression of another chloride channel gene, Clc3, which is autosomal in both mouse species is equal between alleles in F1 animals. There is no evidence of imprinting of the Clc4 autosomal locus. These results are consistent with Ohno’s hypothesis of an evolutionary requirement for a higher expression of genes on the single active X chromosome to maintain balance with autosomal gene expression [Ohno, S. (1967) Sex Chromosomes and Sex-Linked Genes (Springer, Berlin)].

Lactone-ring-cleaving enzyme: Genetic analysis, novel RNA editing, and evolutionary implications

A lactonohydrolase from Fusarium oxysporum AKU 3702 is an enzyme catalyzing the hydrolysis of aldonate lactones to the corresponding aldonic acids. The amino acid sequences of the NH2 terminus and internal peptide fragments of the enzyme were determined to prepare synthetic oligonucleotides as primers for the PCR. An approximate 1,000-base genomic DNA fragment thus amplified was used as the probe to clone both genomic DNA and cDNA for the enzyme. The lactonohydrolase genomic gene consists of six exons separated by five short introns. A novel type of RNA editing, in which lactonohydrolase mRNA included the insertion of guanosine and cytidine residues, was observed. The predicted amino acid sequence of the cloned lactonohydrolase cDNA showed significant similarity to those of the gluconolactonase from Zymomonas mobilis, and paraoxonases from human and rabbit, forming a unique superfamily consisting of C-O cleaving enzymes and P-O cleaving enzymes. Lactonohydrolase was expressed under the control of the lac promoter in Escherichia coli.

Identification of active sites in amidase: Evolutionary relationship between amide bond- and peptide bond-cleaving enzymes

Mainly based on various inhibitor studies previously performed, amidases came to be regarded as sulfhydryl enzymes. Not completely satisfied with this generally accepted interpretation, we performed a series of site-directed mutagenesis studies on one particular amidase of Rhodococcus rhodochrous J1 that was involved in its nitrile metabolism. For these experiments, the recombinant amidase was produced as the inclusion body in Escherichia coli to greatly facilitate its recovery and subsequent purification. With regard to the presumptive active site residue Cys203, a Cys203 → Ala mutant enzyme still retained 11.5% of the original specific activity. In sharp contrast, substitutions in certain other positions in the neighborhood of Cys203 had a far more dramatic effect on the amidase. Glutamic acid substitution of Asp191 reduced the specific activity of the mutant enzyme to 1.33% of the wild-type activity. Furthermore, Asp191 → Asn substitution as well as Ser195 → Ala substitution completely abolished the specific activity. It would thus appear that, among various conserved residues residing within the so-called signature sequence common to all amidases, the real active site residues are Asp191 and Ser195 rather than Cys203. Inasmuch as an amide bond (CO-NH2) in the amide substrate is not too far structurally removed from a peptide bond (CO-NH-), the signature sequences of various amidases were compared with the active site sequences of various types of proteases. It was found that aspartic acid and serine residues corresponding to Asp191 and Ser195 of the Rhodococcus amidase are present within the active site sequences of aspartic proteinases, thus suggesting the evolutionary relationship between the two.

Two hemoglobin genes in Arabidopsis thaliana: The evolutionary origins of leghemoglobins

We cloned two hemoglobin genes from Arabidopsis thaliana. One gene, AHB1, is related in sequence to the family of nonsymbiotic hemoglobin genes previously identified in a number of plant species (class 1). The second hemoglobin gene, AHB2, represents a class of nonsymbiotic hemoglobin (class 2) related in sequence to the symbiotic hemoglobin genes of legumes and Casuarina. The properties of these two hemoglobins suggest that the two families of nonsymbiotic hemoglobins may differ in function from each other and from the symbiotic hemoglobins. AHB1 is induced, in both roots and rosette leaves, by low oxygen levels. Recombinant AHB1 has an oxygen affinity so high as to make it unlikely to function as an oxygen transporter. AHB2 is expressed at a low level in rosette leaves and is low temperature-inducible. AHB2 protein has a lower affinity for oxygen than AHB1 but is similar to AHB1 in having an unusually low, pH-sensitive oxygen off-rate.

Origin and evolutionary relationships of giant Galápagos tortoises

Perhaps the most enduring debate in reptile systematics has involved the giant Galápagos tortoises (Geochelone nigra), whose origins and systematic relationships captivated Charles Darwin and remain unresolved to this day. Here we report a phylogenetic reconstruction based on mitochondrial DNA sequences from Galápagos tortoises and Geochelone from mainland South America and Africa. The closest living relative to the Galápagos tortoise is not among the larger-bodied tortoises of South America but is the relatively small-bodied Geochelone chilensis, or Chaco tortoise. The split between G. chilensis and the Galápagos lineage probably occurred 6 to 12 million years ago, before the origin of the oldest extant Galápagos island. Our data suggest that the four named southern subspecies on the largest island, Isabela, are not distinct genetic units, whereas a genetically distinct northernmost Isabela subspecies is probably the result of a separate colonization. Most unexpectedly, the lone survivor of the abingdoni subspecies from Pinta Island (“Lonesome George”) is very closely related to tortoises from San Cristóbal and Española, the islands farthest from the island of Pinta. To rule out a possible recent transplant of Lonesome George, we sequenced DNA from three tortoises collected on Pinta in 1906. They have sequences identical to Lonesome George, consistent with his being the last survivor of his subspecies. This finding may provide guidance in finding a mate for Lonesome George, who so far has failed to reproduce.

Evolutionary analyses of hedgehog and Hoxd-10 genes in fish species closely related to the zebrafish

The study of development has relied primarily on the isolation of mutations in genes with specific functions in development and on the comparison of their expression patterns in normal and mutant phenotypes. Comparative evolutionary analyses can complement these approaches. Phylogenetic analyses of Sonic hedgehog (Shh) and Hoxd-10 genes from 18 cyprinid fish species closely related to the zebrafish provide novel insights into the functional constraints acting on Shh. Our results confirm and extend those gained from expression and crystalline structure analyses of this gene. Unexpectedly, exon 1 of Shh is found to be almost invariant even in third codon positions among these morphologically divergent species suggesting that this exon encodes for a functionally important domain of the hedgehog protein. This is surprising because the main functional domain of Shh had been thought to be that encoded by exon 2. Comparisons of Shh and Hoxd-10 gene sequences and of resulting gene trees document higher evolutionary constraints on the former than on the latter. This might be indicative of more general evolutionary patterns in networks of developmental regulatory genes interacting in a hierarchical fashion. The presence of four members of the hedgehog gene family in cyprinid fishes was documented and their homologies to known hedgehog genes in other vertebrates were established.

The diversity and evolutionary relationships of the pregnancy-associated glycoproteins, an aspartic proteinase subfamily consisting of many trophoblast-expressed genes

The pregnancy-associated glycoproteins (PAGs) are structurally related to the pepsins, thought to be restricted to the hooved (ungulate) mammals and characterized by being expressed specifically in the outer epithelial cell layer (chorion/trophectoderm) of the placenta. At least some PAGs are catalytically inactive as proteinases, although each appears to possess a cleft capable of binding peptides. By cloning expressed genes from ovine and bovine placental cDNA libraries, by Southern genomic blotting, by screening genomic libraries, and by using PCR to amplify portions of PAG genes from genomic DNA, we estimate that cattle, sheep, and most probably all ruminant Artiodactyla possess many, possibly 100 or more, PAG genes, many of which are placentally expressed. The PAGs are highly diverse in sequence, with regions of hypervariability confined largely to surface-exposed loops. Nonsynonymous (replacement) mutations in the regions of the genes coding for these hypervariable loop segments have accumulated at a higher rate than synonymous (silent) mutations. Construction of distance phylograms, based on comparisons of PAG and related aspartic proteinase amino acid sequences, suggests that much diversification of the PAG genes occurred after the divergence of the Artiodactyla and Perissodactyla, but that at least one gene is represented outside the hooved species. The results also suggest that positive selection of duplicated genes has acted to provide considerable functional diversity among the PAGs, whose presence at the interface between the placenta and endometrium and in the maternal circulation indicates involvement in fetal–maternal interactions.

A single response mechanism is responsible for evolutionary adaptive variation in a bird’s laying date

Many proximate factors determine a bird’s laying date, including environmental and social stimuli as well as individual responses to internal and external factors. However, the relative importance of these factors has not been experimentally demonstrated. Here we show that (i) large differences in the onset of first clutches between different populations result from variation in different responses to photoperiod and not from variation in responses to any other proximate factors and (ii) the same response mechanism causes maladaptive laying dates in habitats modified by humans. We present, to our knowledge, the first experimental demonstration that a single response mechanism is responsible for evolutionary adaptive intraspecific variation in a vertebrate life history trait.

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.