酶祖先序列重建与定向进化

酶祖先序列重建是指通过计算机算法推导来自灭绝生物的祖先酶的氨基酸序列的技术.通常可分为6个步骤,依次为现代酶的核酸/氨基酸序列收集、多序列比对、系统发育树构建、祖先酶序列的计算机推测、基因克隆、酶学性质表征.该方法广泛应用于研究分子在行星时间尺度上对环境条件不断变化的适应性和进化机制.随着酶在生物催化领域中扮演越来越重...     

Glycation improves the thermostability of trypsin and chymotrypsin

A novel glycation procedure, in vacuo glycation, was used to attach glucose covalently to the lysine residues of trypsin and chymotrypsin. Glycated trypsin and glycated chymotrypsin have greatly increased thermostability compared to the native enzymes. For example, glycated bovine trypsin, incubated at 50 degrees C and pH 8.0 for 3 h, retained more than 50% of its original activity whereas the native enzyme was inactivated under the same conditions. Similarly, after incubation at 50 degrees C and pH 8.0, glycated bovine chymotrypsin retained 45% of its original activity and the native enzyme was inactivated. Glycated porcine trypsin is exceptionally thermostable and could be used to digest native ribonuclease at 70 degrees C without the need for prior denaturation. The apparent increase in the thermal stability of the glycated proteins observed in activity measurements is also reflected by an increase in the T(m) values determined with differential scanning calorimetry (DSC) and circular dichroism (CD). The glycation does not alter the activity or specificity of these enzymes.     

Structural and Dynamics Comparison of Thermostability in Ancient,Modern, and Consensus Elongation Factor Tus

1. Download high-res image (199KB)
2. Download full-size image

     

The interface of protein structure, protein biophysics, and molecular evolution

Abstract The interface of protein structural biology, protein biophysics, molecular evolution, and molecular population genetics forms the foundations for a mechanistic understanding of many aspects of protein biochemistry. Current efforts in interdisciplinary protein modeling are in their infancy and the state-of-the art of such models is described. Beyond the relationship between amino acid substitution and static protein structure, protein function, and corresponding organismal fitness, other considerations are also discussed. More complex mutational processes such as insertion and deletion and domain rearrangements and even circular permutations should be evaluated. The role of intrinsically disordered proteins is still controversial, but may be increasingly important to consider. Protein geometry and protein dynamics as a deviation from static considerations of protein structure are also important. Protein expression level is known to be a major determinant of evolutionary rate and several considerations including selection at the mRNA level and the role of interaction specificity are discussed. Lastly, the relationship between modeling and needed high-throughput experimental data as well as experimental examination of protein evolution using ancestral sequence resurrection and in vitro biochemistry are presented, towards an aim of ultimately generating better models for biological inference and prediction.     

The evolution and evolutionary consequences of marginal thermostability in proteins

When we seek to explain the characteristics of living systems in their evolutionary context, we are often interested in understanding how and why certain properties arose through evolution, and how these properties then affected the continuing evolutionary process. This endeavor has been assisted by the use of simple computational models that have properties characteristic of natural living systems but allow simulations over evolutionary timescales with full transparency. We examine a model of the evolution of a gene under selective pressure to code for a protein that exists in a prespecified folded state at a given growth temperature. We observe the emergence of proteins with modest stabilities far below those possible with the model, with a denaturation temperature tracking the simulation temperature, despite the absence of selective pressure for such marginal stability. This demonstrates that neither observations of marginally stable proteins, nor even instances where increased stability interferes with function, provide evidence that marginal stability is an adaptation. Instead the marginal stability is the result of a balance between predominantly destabilizing mutations and selection that shifts depending on effective population size. Even if marginal stability is not an adaptation, the natural tendency of proteins toward marginal stability, and the range of stabilities that occur during evolution, may have significant effect on the evolutionary process.     

Tracing the evolution of novel features of human Toll‐like receptor 4

Toll‐like receptor 4 (TLR4) is a critical innate immune protein that activates inflammation in response to extracellular cues. Much of the work to understand how the protein works in humans has been done using mouse models. Although human and mouse TLR4 have many shared features, they have also diverged significantly since their last common ancestor, acquiring 277 sequence differences. Functional differences include the extent of ligand‐independent activation, whether lipid IVa acts as an antagonist or agonist, and the relative species cross‐compatibility of their MD‐2 cofactor. We set out to understand the evolutionary origins for these functional differences between human and mouse TLR4. Using a combination of phylogenetics, ancestral sequence reconstruction, and functional characterization, we found that evolutionary changes to the human TLR4, rather than changes to the mouse TLR4, were largely responsible for these functional changes. Human TLR4 repressed ancestral ligand‐independent activity and gained antagonism to lipid IVa. Additionally, mutations to the human TLR4 cofactor MD‐2 led to lineage‐specific incompatibility between human and opossum TLR4 complex members. These results were surprising, as mouse TLR4 has acquired many more mutations than human TLR4 since their last common ancestor. Our work has polarized this set of transitions and sets up work to study the mechanistic underpinnings for the evolution of new functions in TLR4.     

A mechanistic view of enzyme evolution

New enzyme functions often evolve through the recruitment and optimization of latent promiscuous activities. How do mutations alter the molecular architecture of enzymes to enhance their activities? Can we infer general mechanisms that are common to most enzymes, or does each enzyme require a unique optimization process? The ability to predict the location and type of mutations necessary to enhance an enzyme's activity is critical to protein engineering and rational design. In this review, via the detailed examination of recent studies that have shed new light on the molecular changes underlying the optimization of enzyme function, we provide a mechanistic perspective of enzyme evolution. We first present a global survey of the prevalence of activity‐enhancing mutations and their distribution within protein structures. We then delve into the molecular solutions that mediate functional optimization, specifically highlighting several common mechanisms that have been observed across multiple examples. As distinct protein sequences encounter different evolutionary bottlenecks, different mechanisms are likely to emerge along evolutionary trajectories toward improved function. Identifying the specific mechanism(s) that need to be improved upon, and tailoring our engineering efforts to each sequence, may considerably improve our chances to succeed in generating highly efficient catalysts in the future.     

Effects of helix and fingertip mutations on the thermostability of xyn11A investigated by molecular dynamics simulations and enzyme activity assays

Local conformational changes and global unfolding pathways of wildtype xyn11A recombinant and its mutated structures were studied through a series of atomistic molecular dynamics (MD) simulations, along with enzyme activity assays at three incubation temperatures to investigate the effects of mutations at three different sites to the thermostability. The first mutation was to replace an unstable negatively charged residue at a surface beta turn near the active site (D32G) by a hydrophobic residue. The second mutation was to create a disulphide bond (S100C/N147C) establishing a strong connection between an alpha helix and a distal beta hairpin associated with the thermally sensitive Thumb loop, and the third mutation add an extra hydrogen bond (A155S) to the same alpha helix. From the MD simulations performed, MM/PBSA energy calculations of the unfolding energy were in a good agreement with the enzyme activities measured from the experiment, as all mutated structures demonstrated the improved thermostability, especially the S100C/N147C proved to be the most stable mutant both by the simulations and the experiment. Local conformational analysis at the catalytic sites and the xylan access region also suggested that mutated xyn11A structures could accommodate xylan binding. However, the analysis of global unfolding pathways showed that structural disruptions at the beta sheet regions near the N-terminal were still imminent. These findings could provide the insight on the molecular mechanisms underlying the enhanced thermostability due to mutagenesis and changes in the protein unfolding pathways for further protein engineering of the GH11 family xylanase enzymes.     

Phylogenetic Articulation of Uric Acid Evolution in Mammals and How It Informs a Therapeutic Uricase

The role of uric acid during primate evolution has remained elusive ever since it was discovered over 100 years ago that humans have unusually high levels of the small molecule in our serum. It has been difficult to generate a neutral or adaptive explanation in part because the uricase enzyme evolved to become a pseudogene in apes thus masking typical signals of sequence evolution. Adding to the difficulty is a lack of clarity on the functional role of uric acid in apes. One popular hypothesis proposes that uric acid is a potent antioxidant that increased in concentration to compensate for the lack of vitamin C synthesis in primate species ∼65 Ma. Here, we have expanded on our previous work with resurrected ancient uricase proteins to better resolve the reshaping of uricase enzymatic activity prior to ape evolution. Our results suggest that the pivotal death-knell to uricase activity occurred between 20 and 30 Ma despite small sequential modifications to its catalytic efficiency for the tens of millions of years since primates lost their ability to synthesize vitamin C, and thus the two appear uncorrelated. We also use this opportunity to demonstrate how molecular evolution can contribute to biomedicine by presenting ancient uricases to human immune cells that assay for innate reactivity against foreign antigens. A highly stable and highly catalytic ancient uricase is shown to elicit a lower immune response in more human haplotypes than other uricases currently in therapeutic development.     

Were Ancestral Proteins Less Specific?

Some have hypothesized that ancestral proteins were, on average, less specific than their descendants. If true, this would provide a universal axis along which to organize protein evolution and suggests that reconstructed ancestral proteins may be uniquely powerful tools for protein engineering. Ancestral sequence reconstruction studies are one line of evidence used to support this hypothesis. Previously, we performed such a study, investigating the evolution of peptide-binding specificity for the paralogs S100A5 and S100A6. The modern proteins appeared more specific than their last common ancestor (ancA5/A6), as each paralog bound a subset of the peptides bound by ancA5/A6. In this study, we revisit this transition, using quantitative phage display to measure the interactions of 30,533 random peptides with human S100A5, S100A6, and ancA5/A6. This unbiased screen reveals a different picture. While S100A5 and S100A6 do indeed bind to a subset of the peptides recognized by ancA5/A6, they also acquired new peptide partners outside of the set recognized by ancA5/A6. Our previous work showed that ancA5/A6 had lower specificity than its descendants when measured against biological targets; our new work shows that ancA5/A6 has similar specificity to the modern proteins when measured against a random set of peptide targets. This demonstrates that altered biological specificity does not necessarily indicate altered intrinsic specificity, and sounds a cautionary note for using ancestral reconstruction studies with biological targets as a means to infer global evolutionary trends in specificity.     

Ancestry and evolution of seasonal migration in the Parulidae

Seasonal migration in birds is known to be highly labile and subject to rapid change in response to selection, such that researchers have hypothesized that phylogenetic relationships should neither predict nor constrain the migratory behaviour of a species. Many theories on the evolution of bird migration assume a framework that extant migratory species have evolved repeatedly and relatively recently from sedentary tropical or subtropical ancestors. We performed ancestral state reconstructions of migratory behaviour using a comprehensive, well-supported phylogeny of the Parulidae (the 'wood-warblers'), a large family of Neotropical and Nearctic migratory and sedentary songbirds, and examined the rates of gain and loss of migration throughout the Parulidae. Counter to traditional hypotheses, our results suggest that the ancestral wood-warbler was migratory and that losses of migration have been at least as prevalent as gains throughout the history of Parulidae. Therefore, extant sedentary tropical radiations in the Parulidae represent losses of latitudinal migration and colonization of the tropics from temperate regions. We also tested for phylogenetic signal in migratory behaviour, and our results indicate that although migratory behaviour is variable within some wood-warbler species and clades, phylogeny significantly predicts the migratory distance of species in the Parulidae.     

Divergent Evolution of a Protein–Protein Interaction Revealed through Ancestral Sequence Reconstruction and Resurrection

The postsynaptic density extends across the postsynaptic dendritic spine with discs large (DLG) as the most abundant scaffolding protein. DLG dynamically alters the structure of the postsynaptic density, thus controlling the function and distribution of specific receptors at the synapse. DLG contains three PDZ domains and one important interaction governing postsynaptic architecture is that between the PDZ3 domain from DLG and a protein called cysteine-rich interactor of PDZ3 (CRIPT). However, little is known regarding functional evolution of the PDZ3:CRIPT interaction. Here, we subjected PDZ3 and CRIPT to ancestral sequence reconstruction, resurrection, and biophysical experiments. We show that the PDZ3:CRIPT interaction is an ancient interaction, which was likely present in the last common ancestor of Eukaryotes, and that high affinity is maintained in most extant animal phyla. However, affinity is low in nematodes and insects, raising questions about the physiological function of the interaction in species from these animal groups. Our findings demonstrate how an apparently established protein–protein interaction involved in cellular scaffolding in bilaterians can suddenly be subject to dynamic evolution including possible loss of function.     

The molecular signal for the adaptation to cold temperature during early life on Earth

Several lines of evidence such as the basal location of thermophilic lineages in large-scale phylogenetic trees and the ancestral sequence reconstruction of single enzymes or large protein concatenations support the conclusion that the ancestors of the bacterial and archaeal domains were thermophilic organisms which were adapted to hot environments during the early stages of the Earth. A parsimonious reasoning would therefore suggest that the last universal common ancestor (LUCA) was also thermophilic. Various authors have used branch-wise non-homogeneous evolutionary models that better capture the variation of molecular compositions among lineages to accurately reconstruct the ancestral G + C contents of ribosomal RNAs and the ancestral amino acid composition of highly conserved proteins. They confirmed the thermophilic nature of the ancestors of Bacteria and Archaea but concluded that LUCA, their last common ancestor, was a mesophilic organism having a moderate optimal growth temperature. In this letter, we investigate the unknown nature of the phylogenetic signal that informs ancestral sequence reconstruction to support this non-parsimonious scenario. We find that rate variation across sites of molecular sequences provides information at different time scales by recording the oldest adaptation to temperature in slow-evolving regions and subsequent adaptations in fast-evolving ones.     

Fossils and living taxa agree on patterns of body mass evolution: a case study with Afrotheria

Most of life is extinct, so incorporating some fossil evidence into analyses of macroevolution is typically seen as necessary to understand the diversification of life and patterns of morphological evolution. Here we test the effects of inclusion of fossils in a study of the body size evolution of afrotherian mammals, a clade that includes the elephants, sea cows and elephant shrews. We find that the inclusion of fossil tips has little impact on analyses of body mass evolution; from a small ancestral size (approx. 100 g), there is a shift in rate and an increase in mass leading to the larger-bodied Paenungulata and Tubulidentata, regardless of whether fossils are included or excluded from analyses. For Afrotheria, the inclusion of fossils and morphological character data affect phylogenetic topology, but these differences have little impact upon patterns of body mass evolution and these body mass evolutionary patterns are consistent with the fossil record. The largest differences between our analyses result from the evolutionary model, not the addition of fossils. For some clades, extant-only analyses may be reliable to reconstruct body mass evolution, but the addition of fossils and careful model selection is likely to increase confidence and accuracy of reconstructed macroevolutionary patterns.     

Evolution and molecular basis of a novel allosteric property of crocodilian hemoglobin

1. Download : Download high-res image (85KB)
2. Download : Download full-size image

     

Sequence optimization for native state stability determines the evolution and folding kinetics of a small protein

Investigating the relative importance of protein stability, function, and folding kinetics in driving protein evolution has long been hindered by the fact that we can only compare modern natural proteins, the products of the very process we seek to understand, to each other, with no external references or baselines. Through a large-scale all-atom simulation of protein evolution, we have created a large diverse alignment of SH3 domain sequences which have been selected only for native state stability, with no other influencing factors. Although the average pairwise identity between computationally evolved and natural sequences is only 17%, the residue frequency distributions of the computationally evolved sequences are similar to natural SH3 sequences at 86% of the positions in the domain, suggesting that optimization for the native state structure has dominated the evolution of natural SH3 domains. Additionally, the positions which play a consistent role in the transition state of three well-characterized SH3 domains (by phi-value analysis) are structurally optimized for the native state, and vice versa. Indeed, we see a specific and significant correlation between sequence optimization for native state stability and conservation of transition state structure.     

Catalysing new reactions during evolution: economy of residues and mechanism

The diversity of function in some enzyme superfamilies shows that during evolution, enzymes have evolved to catalyse different reactions on the same structure scaffold. In this analysis, we examine in detail how enzymes can modify their chemistry, through a comparison of the catalytic residues and mechanisms in 27 pairs of homologous enzymes of totally different functions. We find that evolution is very economical. Enzymes retain structurally conserved residues to aid catalysis, including residues that bind catalytic metal ions and modulate cofactor chemistry. We examine the conservation of residue type and residue function in these structurally conserved residue pairs. Additionally, enzymes often retain common mechanistic steps catalyzed by structurally conserved residues. We have examined these steps in the context of their overall reactions.     

Diversity and evolution of male secondary sexual characters in African squeakers and long-fingered frogs

The African frog genera Arthroleptis and Cardioglossa are unique among vertebrates in having males with extremely long third fingers. In some species, this sexual dimorphism is impressive, with male third fingers approaching 40% of body length. However, the diversity of this trait has not been documented thoroughly and several species appear to lack this trait. The present study documents the diversity of male secondary sexual traits in Arthroleptis and Cardioglossa , including elongate third fingers and digital and inguinal spines. Furthermore, it explores hypotheses of trait evolution, including explanations for the absence of male traits. Analyses of covariance suggest that the functional relationship between finger length and snout–vent length (SVL), both within and among species, is different for male finger III than for male fingers I, II, and IV, or for female finger III. Ancestral state reconstruction suggests that all male traits were present in the most recent common ancestor of Arthroleptis and Cardioglossa and that reduction or loss of traits occurred later. Across species, independent contrast analyses find no correlation between SVL and either male relative third digit length or dimorphism in relative third digit length. The number of spines on male fingers II and III are positively correlated but spine number is not correlated with SVL and only weakly correlated with relative third digit length. The diversity of male traits is evolutionarily labile and is not explained by simple hypotheses of character evolution. Arthroleptis and Cardioglossa may thus provide an interesting study system for understanding how changes in sexual selection forces produce male trait diversity.  © 2009 The Linnean Society of London, Biological Journal of the Linnean Society , 2009, 96 , 553–573.