High-resolution structure of the conger eel galectin, congerin I, in lactose-liganded and ligand-free forms: emergence of a new structure class by accelerated evolution.

BACKGROUND: Congerin I is a member of the galectin (animal beta-galactoside-binding lectin) family and is found in the skin mucus of conger eel. The galectin family proteins perform a variety of biological activities. Because of its histological localization and activity against marine bacteria and starfish embryos, congerin I is thought to take part in the eels' biological defense system against parasites. RESULTS: The crystal structure of congerin I has been determined in both lactose-liganded and ligand-free forms to 1. 5 A and 1.6 A resolution, respectively. The protein is a homodimer of 15 kDa subunits. Congerin I has a beta-sheet topology that is markedly different from those of known relatives. One of the beta-strands is exchanged between two identical subunits. This strand swap might increase the dimer stability. Of the known galectin complexes, congerin I forms the most extensive interaction with lactose molecules. Most of these interactions are substituted by similar interactions with water molecules, including a pi-electron hydrogen bond, in the ligand-free form. This observation indicates an increased affinity of congerin I for the ligand. CONCLUSIONS: The genes for congerin I and an isoform, congerin II, are known to have evolved under positive selection pressure. The strand swap and the modification in the carbohydrate-binding site might enhance the cross-linking activity, and should be the most apparent consequence of positive selection. The protein has been adapted to functioning in skin mucus that is in direct contact with surrounding environments by an enhancement in cross-linking activity. The structure of congerin I demonstrates the emergence of a new structure class by accelerated evolution under selection pressure.     

Protein engineering of conger eel galectins by tracing of molecular evolution using probable ancestral mutants

Background  

Conger eel galectins, congerin I (ConI) and congerin II (ConII), show the different molecular characteristics resulting from accelerating evolution. We recently reconstructed a probable ancestral form of congerins, Con-anc. It showed properties similar to those of ConII in terms of thermostability and carbohydrate recognition specificity, although it shares a higher sequence similarity with ConI than ConII.     

In vitro evolutionary thermostabilization of congerin II: a limited reproduction of natural protein evolution by artificial selection pressure

The thermostability of the conger eel galectin, congerin II, was improved by in vitro evolutionary protein engineering. Two rounds of random PCR mutagenesis and selection experiments increased the congerin II thermostability to a level comparative to its naturally thermostable isoform, congerin I. The crystal structures of the most thermostable double mutant, Y16S/T88I, and the related single mutants, Y16S and T88I, were determined at 2.0 angstroms, 1.8 angstroms, and 1.6 angstroms resolution, respectively. The exclusion of two interior water molecules by the Thr88Ile mutation, and the relief of adjacent conformational stress by the Tyr16Ser mutation were the major contributions to the thermostability. These features in the congerin II mutants are similar to those observed in congerin I. The natural evolution of congerin genes, with the K(A)/K(S) ratio of 2.6, was accelerated under natural selection pressures. The thermostabilizing selection pressure artificially applied to congerin II mimicked the implied natural pressure on congerin I. The results showed that the artificial pressure made congerin II partially reproduce the natural evolution of congerin I.     

Reconstruction of a probable ancestral form of conger eel galectins revealed their rapid adaptive evolution process for specific carbohydrate recognition

Recently, many cases of rapid adaptive evolution, which is characterized by the higher substitution rates of nonsynonymous substitutions to synonymous ones, have been identified in the various genes of venomous and biodefense proteins, including the conger eel galectins, congerins I and II (ConI and ConII). To understand the evolutionary process of congerins, we prepared a probable ancestral form, Con-anc, corresponding to the putative amino acid sequence at the divergence of ConI and ConII in phylogenetic tree with 76% and 61% sequence identities to the current proteins, respectively. Con-anc and ConII had comparable thermostability and similar carbohydrate specificities in general, whereas ConI was more thermostable and showed different carbohydrate specificities. Con-anc showed decreased specificity to oligosaccharides with alpha 2,3-sialyl galactose moieties. These suggest that ConI and ConII have evolved via accelerated evolution under significant selective pressure to increase the thermostability and to acquire the activity to bind to alpha2,3-sialyl galactose present in pathogenic bacteria, respectively. Furthermore, comparative mutagenesis analyses of Con-anc and congerins revealed the structural basis for specific recognition of ConII to alpha2,3-sialyl galactose moiety, and strong binding ability of ConI to oligosaccharides including lacto-N-biosyl (Galbeta1-3GlcNAc) or lacto-N-neobiosyl (Galbeta1-4GlcNAc) residues, respectively. Thus, protein engineering using a probable ancestral form presented here is a powerful approach not only to determine the evolutionary process but also to investigate the structure-activity relationships of proteins.     

Crystal structure of a conger eel galectin (congerin II) at 1.45A resolution: implication for the accelerated evolution of a new ligand-binding site following gene duplication

The crystal structure of congerin II, a galectin family lectin from conger eel, was determined at 1.45A resolution. The previously determined structure of its isoform, congerin I, had revealed a fold evolution via strand swap; however, the structure of congerin II described here resembles other prototype galectins. A comparison of the two congerin genes with that of several other galectins suggests acceralated evolution of both congerin genes following gene duplication. The presence of a Mes (2-[N-morpholino]ethanesulfonic acid) molecule near the carbohydrate-binding site in the crystal structure points to the possibility of an additional binding site in congerin II. The binding site consists of a group of residues that had been replaced following gene duplication suggesting that the binding site was built under selective pressure. Congerin II may be a protein specialized for biological defense with an affinity for target carbohydrates on parasites' cell surface.     

Allosteric Regulation of the Carbohydrate-binding Ability of a Novel Conger Eel Galectin by D-Mannoside

Conger eel has two galectins, termed congerins I and II (Con I and II), that function in mucus as biodefense molecules. Con I and II have acquired a novel protein fold via domain swapping and a new ligand-binding site by accelerated evolution, which enables recognition of some marine bacteria. In this study, we identified a new congerin isotype, congerin P (Con-P), from the peritoneal cells of conger eel. Although Con-P displayed obvious homology with galectins, we observed substitution of 7 out of 8 amino acid residues in the carbohydrate recognition domain that are conserved in all other known galectins. To understand the structure-function relationships of this unique galectin, recombinant Con-P was successfully expressed in Escherichia coli by using a Con II-tagged fusion protein system and subsequently characterized. In the presence of d-mannose, Con-P displayed 30-fold greater hemagglutinating activity than Con I; however, no activity was observed without mannose, indicating that d-mannoside can act as a modulator of Con-P. Frontal affinity chromatography analysis showed that activated Con-P, allosterically induced by mannose, displayed affinity for oligomannose-type sugars as well as N-acetyllactosamine-type β-galactosides. Thus, Con-P represents a new member of the galectin family with unique properties.     

Robustness of predictions of extremely thermally stable proteins in ancient organisms

A number of studies have addressed the environmental temperatures experienced by ancient life. Computational studies using a nonhomogeneous evolution model have estimated ancestral G + C contents of ribosomal RNAs and the amino acid compositions of ancestral proteins, generating hypotheses regarding the mesophilic last universal common ancestor. In contrast, our previous study computationally reconstructed ancestral amino acid sequences of nucleoside diphosphate kinases using a homogeneous model and then empirically resurrected the ancestral proteins. The thermal stabilities of these ancestral proteins were equivalent to or greater than those of extant homologous thermophilic proteins, supporting the thermophilic universal ancestor theory. In this study, we reinferred ancestral sequences using a dataset from which hyperthermophilic sequences were excluded. We also reinferred ancestral sequences using a nonhomogeneous evolution model. The newly reconstructed ancestral proteins are still thermally stable, further supporting the hypothesis that the ancient organisms contained thermally stable proteins and therefore that they were thermophilic.     

Functional Diversification and Evolution of Antifreeze Proteins in the Antarctic Fish Lycodichthys dearborni

Antifreeze proteins (AFPs) have independently evolved in many organisms. AFPs act by binding to ice crystals, effectively lowering the freezing point. AFPs are often at high copy number in a genome and diversity exists between copies. Type III antifreeze proteins are found in Arctic and Antarctic eel pouts, and have previously been shown to evolve under positive selection. Here we combine molecular and proteomic techniques to understand the molecular evolution and diversity of Type III antifreeze proteins in a single individual Antarctic fish Lycodichthys dearborni. Our expressed sequence tag (EST) screen reveals that at least seven different AFP variants are transcribed, which are ultimately translated into five different protein isoforms. The isoforms have identical 66 base pair signal sequences and different numbers of subsequent ice-binding domains followed by a stop codon. Isoforms with one ice-binding unit (monomer), two units (dimer), and multiple units (multimer) were present in the EST library. We identify a previously uncharacterized protein dimer, providing further evidence that there is diversity between Type III AFP isoforms, perhaps driven by positive selection for greater thermal hysteresis. Proteomic analysis confirms that several of these isoforms are translated and present in the liver. Our molecular evolution study shows that paralogs have diverged under positive selection. We hypothesize that antifreeze protein diversity is an important contributor to depressing the serum freezing point.     

Evolution of structure and function of V-ATPases

Proton pumping ATPases/ATPsynthases are found in all groups of present-day organisms. The structure of V- and F-type ATPases/ATP synthases is very conserved throughout evolution. Sequence analysis shows that the V- and F-type ATPases evolved from the same enzyme already present in the last common ancestor of all known extant life forms. The catalytic and noncatalytic subunits found in the dissociable head groups of the V/F-type ATPases are paralogous subunits, i.e., these two types of subunits evolved from a common ancestral gene. The gene duplication giving rise to these two genes (i.e., encoding the catalytic and noncatalytic subunits) predates the time of the last common ancestor.Mapping of gene duplication events that occurred in the evolution of the proteolipid, the noncatalytic and the catalytic subunits, onto the tree of life leads to a prediction for the likely subunit structure of the encoded ATPases. A correlation between structure and function of V/F-ATPases has been established for present-day organisms. Implications resulting from this correlation for the bioenergetics operative in proto-eukaryotes and in the last common ancestor are presented. The similarities of the V/F-ATPase subunits to an ATPase-like protein that was implicated to play a role in flagellar assembly are evaluated.Different V-ATPase isoforms have been detected in some higher eukaryotes. These data are analyzed with respect to the possible function of the different isoforms (tissue specific, organelle specific) and with respect to the point in their evolution when these gene duplications giving rise to the isoforms had occurred, i.e., how far these isoforms are distributed.     

High-level Expression and Characterization of Fully Active Recombinant Conger Eel Galectins in Eschericia coli

An expression system for recombinant conger eel galectins, congerins I and II, were constructed using the pTV 118N plasmid vector and Escherichia coli. Recombinant congerins I and II could be obtained in the soluble active form with high quantitative yield. Mutation of codons for Val and Leu located in the N-terminal region of Con I increased the expression efficiency. Purification of recombinant proteins were done by only two chromatographical steps from E. coli extract. The purified recombinant congerins were found to be almost the same as the native ones except for the acetyl group at the N-terminus; that is, they showed the same structures and carbohydrate binding activities, suggesting that N-terminal acetyl groups of congerins were not significant for activity.     

High-level expression and characterization of fully active recombinant conger eel galectins in Eschericia coli

An expression system for recombinant conger eel galectins, congerins I and II, were constructed using the pTV 118N plasmid vector and Escherichia coli. Recombinant congerins I and II could be obtained in the soluble active form with high quantitative yield. Mutation of codons for Val and Leu located in the N-terminal region of Con I increased the expression efficiency. Purification of recombinant proteins were done by only two chromatographical steps from E. coli extract. The purified recombinant congerins were found to be almost the same as the native ones except for the acetyl group at the N-terminus; that is, they showed the same structures and carbohydrate binding activities, suggesting that N-terminal acetyl groups of congerins were not significant for activity.     

The evolution of hominoid cranial diversity: A quantitative genetic approach

Hominoid cranial evolution is characterized by substantial phenotypic diversity, yet the cause of this variability has rarely been explored. Quantitative genetic techniques for investigating evolutionary processes underlying morphological divergence are dependent on the availability of good ancestral models, a problem in hominoids where the fossil record is fragmentary and poorly understood. Here, we use a maximum likelihood approach based on a Brownian motion model of evolutionary change to estimate nested hypothetical ancestral forms from 15 extant hominoid taxa. These ancestors were then used to calculate rates of evolution along each branch of a phylogenetic tree using Lande's generalized genetic distance. Our results show that hominoid cranial evolution is characterized by strong stabilizing selection. Only two instances of directional selection were detected; the divergence of Homo from its last common ancestor with Pan, and the divergence of the lesser apes from their last common ancestor with the great apes. In these two cases, selection gradients reconstructed to identify the specific traits undergoing selection indicated that selection on basicranial flexion, cranial vault expansion, and facial retraction characterizes the divergence of Homo, whereas the divergence of the lesser apes was defined by selection on neurocranial size reduction.     

A Universal Trend among Proteomes Indicates an Oily Last Common Ancestor

Despite progresses in ancestral protein sequence reconstruction, much needs to be unraveled about the nature of the putative last common ancestral proteome that served as the prototype of all extant lifeforms. Here, we present data that indicate a steady decline (oil escape) in proteome hydrophobicity over species evolvedness (node number) evident in 272 diverse proteomes, which indicates a highly hydrophobic (oily) last common ancestor (LCA). This trend, obtained from simple considerations (free from sequence reconstruction methods), was corroborated by regression studies within homologous and orthologous protein clusters as well as phylogenetic estimates of the ancestral oil content. While indicating an inherent irreversibility in molecular evolution, oil escape also serves as a rare and universal reaction-coordinate for evolution (reinforcing Darwin''s principle of Common Descent), and may prove important in matters such as (i) explaining the emergence of intrinsically disordered proteins, (ii) developing composition- and speciation-based “global” molecular clocks, and (iii) improving the statistical methods for ancestral sequence reconstruction.     

Modeling Olfactory Bulb Evolution through Primate Phylogeny

Adaptive characterizations of primates have usually included a reduction in olfactory sensitivity. However, this inference of derivation and directionality assumes an ancestral state of olfaction, usually by comparison to a group of extant non-primate mammals. Thus, the accuracy of the inference depends on the assumed ancestral state. Here I present a phylogenetic model of continuous trait evolution that reconstructs olfactory bulb volumes for ancestral nodes of primates and mammal outgroups. Parent-daughter comparisons suggest that, relative to the ancestral euarchontan, the crown-primate node is plesiomorphic and that derived reduction in olfactory sensitivity is an attribute of the haplorhine lineage. The model also suggests a derived increase in olfactory sensitivity at the strepsirrhine node. This oppositional diversification of the strepsirrhine and haplorhine lineages from an intermediate and non-derived ancestor is inconsistent with a characterization of graded reduction through primate evolution.     

The speciation of conger eel galectins by rapid adaptive evolution

Many cases of accelerated evolution driven by positive Darwinian selection are identified in the genes of venomous and reproductive proteins. This evolutional phenomenon might have important consequences in their gene-products' functions, such as multiple specific toxins for quick immobilization of the prey and the establishment of barriers to fertilization that might lead to speciation, and in the molecular evolution of novel genes. Recently, we analyzed the molecular evolution of two galectins isolated from the skin mucus of conger eel (Conger myriaster), named congerins I and II, by cDNA cloning and X-ray structural analysis, and we found that they have evolved in the rapid adaptive manner to emergence of a new structure including strand-swapping and a unique new ligand-binding site. In this review article we summarize and discuss the molecular evolution, especially the rapid adaptive evolution, and the structure-function relationships of conger eel galectins.     

Estimation of the ancestral effective population sizes of African great apes under different selection regimes

Reliable estimates of ancestral effective population sizes are necessary to unveil the population-level phenomena that shaped the phylogeny and molecular evolution of the African great apes. Although several methods have previously been applied to infer ancestral effective population sizes, an analysis of the influence of the selective regime on the estimates of ancestral demography has not been thoroughly conducted. In this study, three independent data sets under different selective regimes were used were composed to tackle this issue. The results showed that selection had a significant impact on the estimates of ancestral effective population sizes of the African great apes. The inference of the ancestral demography of African great apes was affected by the selection regime. The effects, however, were not homogeneous along the ancestral populations of great apes. The effective population size of the ancestor of humans and chimpanzees was more impacted by the selection regime when compared to the same parameter in the ancestor of humans, chimpanzees and gorillas. Because the selection regime influenced the estimates of ancestral effective population size, it is reasonable to assume that a portion of the discrepancy found in previous studies that inferred the ancestral effective population size may be attributable to the differential action of selection on the genes sampled.     

The origin of ancestral bacterial metabolism

A model for metabolism of the last bacterial common ancestor based on biomimetic analysis of the metabolic systems of phylogenetically ancient bacteria is developed. The mechanism of natural selection and evolution of the autocatalytic chemical systems under the effect of natural homeostatic parameters, such as chemical potentials, temperature, and pressure of environment is proposed. Competition between particular parts of the autocatalytic network with positive-plus-negative feedback resulted in the formation of particular systems of primary autotrophic, mixotrophic, and heterotrophic metabolism. The model of the last common ancestor as a combination of coupled metabolic cycles among population of protocells is discussed. Physicochemical features of these metabolic cycles determined the major principles of natural selection towards ancestral bacterial taxa.     

To Investigate Protein Evolution by Detecting Suppressed Epitope Structures

Material remains of ancestor nucleotides and proteins are largely unavailable, thus sequence comparison among homologous genes in present-day organisms forms the core of current knowledge of molecular evolution. Variation in protein three-dimensional structure is a basis for functional diversity. To study the evolution of three-dimensional structures in related proteins would significantly improve our understanding of protein evolution and function. A protein may contain ancestor conformations that have been allosterically suppressed by evolutionarily additive structures. Using monoclonal antibody probes to detect such conformation in proteins after removing the suppressor structure, our study demonstrates three-dimensional structure evidence for the evolutionary relationship between troponin I and troponin T, two subunits of the troponin complex in the Ca²⁺-regulatory system of striated muscle, and among their muscle type-specific isoforms. The experimental data show the feasibility of detecting evolutionarily suppressed history-telling structural states in proteins by removing conformational modulator segments added during evolution. In addition to identifying structural modifications that were critical to the emergence of diverged proteins, investigating this novel mode of evolution will help us to understand the origin and functional potential of protein structures.     

Estimates of natural selection due to protein tertiary structure inform the ancestry of biallelic loci

We consider the inference of which of two alleles is ancestral when the alleles have a single nonsynonymous difference and when natural selection acts via protein tertiary structure. Whereas the probability that an allele is ancestral under neutrality is equal to its frequency, under selection this probability depends on allele frequency and on the magnitude and direction of selection pressure. Although allele frequencies can be well estimated from intraspecific data, small fitness differences have a large evolutionary impact but can be difficult to estimate with only intraspecific data. Methods for predicting aspects of phenotype from genotype can supplement intraspecific sequence data. Recently developed statistical techniques can assess effects of phenotypes, such as protein tertiary structure on molecular evolution. While these techniques were initially designed for comparing protein-coding genes from different species, the resulting interspecific inferences can be assigned population genetic interpretations to assess the effect of selection pressure, and we use them here along with intraspecific allele frequency data to estimate the probability that an allele is ancestral. We focus on 140 nonsynonymous single nucleotide polymorphisms of humans that are in proteins with known tertiary structures. We find that our technique for employing protein tertiary structure information yields some biologically plausible results but that it does not substantially improve the inference of ancestral human allele types.     

The very early evolution of protein translocation across membranes

In this study, we used a computational approach to investigate the early evolutionary history of a system of proteins that, together, embed and translocate other proteins across cell membranes. Cell membranes comprise the basis for cellularity, which is an ancient, fundamental organizing principle shared by all organisms and a key innovation in the evolution of life on Earth. Two related requirements for cellularity are that organisms are able to both embed proteins into membranes and translocate proteins across membranes. One system that accomplishes these tasks is the signal recognition particle (SRP) system, in which the core protein components are the paralogs, FtsY and Ffh. Complementary to the SRP system is the Sec translocation channel, in which the primary channel-forming protein is SecY. We performed phylogenetic analyses that strongly supported prior inferences that FtsY, Ffh, and SecY were all present by the time of the last universal common ancestor of life, the LUCA, and that the ancestor of FtsY and Ffh existed before the LUCA. Further, we combined ancestral sequence reconstruction and protein structure and function prediction to show that the LUCA had an SRP system and Sec translocation channel that were similar to those of extant organisms. We also show that the ancestor of Ffh and FtsY that predated the LUCA was more similar to FtsY than Ffh but could still have comprised a rudimentary protein translocation system on its own. Duplication of the ancestor of FtsY and Ffh facilitated the specialization of FtsY as a membrane bound receptor and Ffh as a cytoplasmic protein that could bind nascent proteins with specific membrane-targeting signal sequences. Finally, we analyzed amino acid frequencies in our ancestral sequence reconstructions to infer that the ancestral Ffh/FtsY protein likely arose prior to or just after the completion of the canonical genetic code. Taken together, our results offer a window into the very early evolutionary history of cellularity.