Gaining and losing the thermophilic adaptation in prokaryotes

Structural features, temperature adaptation and industrial applications of microbial lipases from psychrophilic, mesophilic and thermophilic origins

Microbial lipases are very prominent biocatalysts because of their ability to catalyze a wide variety of reactions in aqueous and non-aqueous media. Here microbial lipases from different origins (psychrophiles, mesophiles, and thermophiles) have been reviewed. This review emphasizes an update of structural diversity in temperature adaptation and industrial applications, of psychrophilic, mesophilic, and thermophilic lipases. The microbial origins of lipases are logically dynamic, proficient, and also have an extensive range of industrial uses with the manufacturing of altered molecules. It is therefore of interest to understand the molecular mechanisms of adaptation to temperature in occurring lipases. However, lipases from extremophiles (psychrophiles, and thermophiles) are widely used to design biotransformation reactions with higher yields, fewer byproducts, or useful side products and have been predicted to catalyze those reactions also, which otherwise are not possible with the mesophilic lipases. Lipases as a multipurpose biological catalyst have given a favorable vision in meeting the needs of several industries such as biodiesel, foods, and drinks, leather, textile, detergents, pharmaceuticals, and medicals.

Short-Term Exposure to Thermophilic Temperatures Facilitates CO Uptake by Thermophiles Maintained under Predominantly Mesophilic Conditions

Three phylogenetically and phenotypically distinct CO-oxidizing thermophiles (Alicyclobacillus macrosporangiidus CPP55 (Firmicutes), Meiothermus ruber PS4 (Deinococcus-Thermus) and Thermogemmatispora carboxidovorans PM5T (Chloroflexi)) and one CO-oxidizing mesophile (Paraburkholderia paradisi WAT (Betaproteobacteria)) isolated from volcanic soils were used to assess growth responses and CO uptake rates during incubations with constant temperatures (25 °C and 55 °C) and during multi-day incubations with a temperature regime that cycled between 20 °C and 55 °C on a diurnal basis (alternating mesophilic and thermophilic temperatures, AMTT). The results were used to test a conjecture that some thermophiles can survive in mesothermal habitats that experience occasional thermophilic temperatures. Meiothermus ruber PS4, which does not form spores, was able to grow and oxidize CO under all conditions, while the spore-forming Alicyclobacillus macrosporangiidus CPP55 grew and oxidized CO during the AMTT regime and at 55 °C, but was not active at 25 °C. Thermogemmatispora carboxidovorans PM5T, also a spore former, only grew at 55 °C but oxidized CO during AMTT and 55 °C incubations. In contrast, the non-sporing mesophile, Paraburkholderia paradisi WAT, was only able to grow and oxidize CO at 25 °C; growth and CO uptake ceased during the AMTT incubations after exposure to the initial round of thermophilic temperatures. Collectively, these results suggest that temporary, periodic exposure to permissive growth temperatures could help maintain populations of thermophiles in mesothermal habitats after deposition from the atmosphere or other sources.

The Molecular Determinants of Thermoadaptation: Methanococcales as a Case Study

Previous reports have shown that environmental temperature impacts proteome evolution in Bacteria and Archaea. However, it is unknown whether thermoadaptation mainly occurs via the sequential accumulation of substitutions, massive horizontal gene transfers, or both. Measuring the real contribution of amino acid substitution to thermoadaptation is challenging, because of confounding environmental and genetic factors (e.g., pH, salinity, genomic G + C content) that also affect proteome evolution. Here, using Methanococcales, a major archaeal lineage, as a study model, we show that optimal growth temperature is the major factor affecting variations in amino acid frequencies of proteomes. By combining phylogenomic and ancestral sequence reconstruction approaches, we disclose a sequential substitutional scheme in which lysine plays a central role by fine tuning the pool of arginine, serine, threonine, glutamine, and asparagine, whose frequencies are strongly correlated with optimal growth temperature. Finally, we show that colonization to new thermal niches is not associated with high amounts of horizontal gene transfers. Altogether, although the acquisition of a few key proteins through horizontal gene transfer may have favored thermoadaptation in Methanococcales, our findings support sequential amino acid substitutions as the main factor driving thermoadaptation.

Taxonomic Novelty and Distinctive Genomic Features of Hot Spring Cyanobacteria

Several cyanobacterial species are dominant primary producers in hot spring microbial mats. To date, hot spring cyanobacterial taxonomy, as well as the evolution of their genomic adaptations to high temperatures, are poorly understood, with genomic information currently available for only a few dominant genera, including Fischerella and Synechococcus. To address this knowledge gap, the present study expands the genomic landscape of hot spring cyanobacteria and traces the phylum-wide genomic consequences of evolution in high temperature environments. From 21 globally distributed hot spring metagenomes, with temperatures between 32 and 75°C, 57 medium- and high-quality cyanobacterial metagenome-assembled genomes were recovered, representing taxonomic novelty for 1 order, 3 families, 15 genera and 36 species. Comparative genomics of 93 hot spring genomes (including the 57 metagenome-assembled genomes) and 66 non-thermal genomes, showed that the former have smaller genomes and a higher GC content, as well as shorter proteins that are more hydrophilic and basic, when compared to the non-thermal genomes. Additionally, the core accessory orthogroups from the hot spring genomes of some genera had a greater abundance of functional categories, such as inorganic ion metabolism, translation and post-translational modifications. Moreover, hot spring genomes showed increased abundances of inorganic ion transport and amino acid metabolism, as well as less replication and transcription functions in the protein coding sequences. Furthermore, they showed a higher dependence on the CRISPR-Cas defense system against exogenous nucleic acids, and a reduction in secondary metabolism biosynthetic gene clusters. This suggests differences in the cyanobacterial response to environment-specific microbial communities. This phylum-wide study provides new insights into cyanobacterial genomic adaptations to a specific niche where they are dominant, which could be essential to trace bacterial evolution pathways in a warmer world, such as the current global warming scenario.

The tale of a versatile enzyme: Alpha-amylase evolution, structure, and potential biotechnological applications for the bioremediation of n-alkanes

As the primary source of a wide range of industrial products, the study of petroleum-derived compounds is of pivotal importance. However, the process of oil extraction and refinement is among the most environmentally hazardous practices, impacting almost all levels of the ecological chain. So far, the most appropriate strategy to overcome such an issue is through bioremediation, which revolves around the employment of different microorganisms to degrade hazardous compounds, generating less environmental impact and lower monetary costs. In this sense, a myriad of organisms and enzymes are considered possible candidates for the bioremediation process. Amidst the potential candidates is α-amylase, an evolutionary conserved starch-degrading enzyme. Notably, α-amylase was not only seen to degrade n-alkanes, a subclass of alkanes considered the most abundant petroleum-derived compounds but also low-density polyethylene, a dangerous pollutant produced from petroleum. Thus, due to its high conservation in both eukaryotic and prokaryotic lineages, in addition to the capability to degrade different types of hazardous compounds, the study of α-amylase becomes a rising interest. Nevertheless, there are no studies that review all biotechnological applications of α-amylase for bioremediation. In this work, we critically review the potential biotechnological applications of α-amylase, focusing on the biodegradation of petroleum-derived compounds. Evolutionary aspects are discussed, as well for all structural information and all features that could impact on the employment of this protein in the biotechnological industry, such as pH, temperature, and medium conditions. New perspectives and critical assessments are conducted regarding the application of α-amylase in the bioremediation of n-alkanes.

Predicting the optimal growth temperatures of prokaryotes using only genome derived features

MOTIVATION

Optimal growth temperature is a fundamental characteristic of all living organisms. Knowledge of this temperature is central to the study of a prokaryote, the thermal stability and temperature dependent activity of its genes, and the bioprospecting of its genome for thermally adapted proteins. While high throughput sequencing methods have dramatically increased the availability of genomic information, the growth temperatures of the source organisms are often unknown. This limits the study and technological application of these species and their genomes. Here, we present a novel method for the prediction of growth temperatures of prokaryotes using only genomic sequences.

RESULTS

By applying the reverse ecology principle that an organism's genome includes identifiable adaptations to its native environment, we can predict a species' optimal growth temperature with an accuracy of 5.17°C root-mean-square error and a coefficient of determination of 0.835. The accuracy can be further improved for specific taxonomic clades or by excluding psychrophiles. This method provides a valuable tool for the rapid calculation of organism growth temperature when only the genome sequence is known.

AVAILABILITY AND IMPLEMENTATION

Source code, genomes analyzed and features calculated are available at: https://github.com/DavidBSauer/OGT_prediction.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

An “overlooked” habitat for thermophilic bacteria: the phyllosphere

Thermophilic microbes are present everywhere around us and their only known natural biotope is far away and most usually associated with geothermal energy. To answer this paradox, we explore the hypothesis that the phyllosphere (surface of leaves), due to its exposition to the sun, could well be a thermophilic habitat for microbes and thus a source of thermophilic microbes growing around 50°C – 60°C. To support this hypothesis, we reviewed the heat sources on earth and associated microbial habitats, as well as the difficult identification of thermophilic microbes. We further present an experiment to show the presence and activity of thermophilic bacteria in the phyllosphere. Leaves were collected from eleven tree species from five locations on three continents belonging to three different biomes. On fresh leaves, 16S rDNA sequencing reveals the presence of 0.2 to 7% of clearly identified thermophilic bacteria. Moreover, after incubation at 55°C under aerobic and anaerobic conditions, 16S rDNA sequencing reveals the presence of 4 to 99% of clearly identified thermophilic bacteria. The accumulation of observations provides coherence to our hypothesis and allows the emergence of a new vision of leaves as a thermophilic biotope. We then propose a life cycle of microbes belonging to the thermophilic biotope associated with leaf surfaces.

Horizontal Gene Transfer as an Indispensable Driver for Evolution of Neocallimastigomycota into a Distinct Gut-Dwelling Fungal Lineage

Survival and growth of the anaerobic gut fungi (AGF; Neocallimastigomycota) in the herbivorous gut necessitate the possession of multiple abilities absent in other fungal lineages. We hypothesized that horizontal gene transfer (HGT) was instrumental in forging the evolution of AGF into a phylogenetically distinct gut-dwelling fungal lineage. The patterns of HGT were evaluated in the transcriptomes of 27 AGF strains, 22 of which were isolated and sequenced in this study, and 4 AGF genomes broadly covering the breadth of AGF diversity. We identified 277 distinct incidents of HGT in AGF transcriptomes, with subsequent gene duplication resulting in an HGT frequency of 2 to 3.5% in AGF genomes. The majority of HGT events were AGF specific (91.7%) and wide (70.8%), indicating their occurrence at early stages of AGF evolution. The acquired genes allowed AGF to expand their substrate utilization range, provided new venues for electron disposal, augmented their biosynthetic capabilities, and facilitated their adaptation to anaerobiosis. The majority of donors were anaerobic fermentative bacteria prevalent in the herbivorous gut. This study strongly indicates that HGT indispensably forged the evolution of AGF as a distinct fungal phylum and provides a unique example of the role of HGT in shaping the evolution of a high-rank taxonomic eukaryotic lineage.IMPORTANCE The anaerobic gut fungi (AGF) represent a distinct basal phylum lineage (Neocallimastigomycota) commonly encountered in the rumen and alimentary tracts of herbivores. Survival and growth of anaerobic gut fungi in these anaerobic, eutrophic, and prokaryote-dominated habitats necessitates the acquisition of several traits absent in other fungal lineages. We assess here the role of horizontal gene transfer as a relatively fast mechanism for trait acquisition by the Neocallimastigomycota postsequestration in the herbivorous gut. Analysis of 27 transcriptomes that represent the broad diversity of Neocallimastigomycota identified 277 distinct HGT events, with subsequent gene duplication resulting in an HGT frequency of 2 to 3.5% in AGF genomes. These HGT events have allowed AGF to survive in the herbivorous gut by expanding their substrate utilization range, augmenting their biosynthetic pathway, providing new routes for electron disposal by expanding fermentative capacities, and facilitating their adaptation to anaerobiosis. HGT in the AGF is also shown to be mainly a cross-kingdom affair, with the majority of donors belonging to the bacteria. This study represents a unique example of the role of HGT in shaping the evolution of a high-rank taxonomic eukaryotic lineage.

Low Temperature Adaptation Is Not the Opposite Process of High Temperature Adaptation in Terms of Changes in Amino Acid Composition

Previous studies focused on psychrophilic adaptation generally have demonstrated that multiple mechanisms work together to increase protein flexibility and activity, as well as to decrease the thermostability of proteins. However, the relationship between high and low temperature adaptations remains unclear. To investigate this issue, we collected the available predicted whole proteome sequences of species with different optimal growth temperatures, and analyzed amino acid variations and substitutional asymmetry in pairs of homologous proteins from related species. We found that changes in amino acid composition associated with low temperature adaptation did not exhibit a coherent opposite trend when compared with changes in amino acid composition associated with high temperature adaptation. This result indicates that during their evolutionary histories the proteome-scale evolutionary patterns associated with prokaryotes exposed to low temperature environments were distinct from the proteome-scale evolutionary patterns associated with prokaryotes exposed to high temperature environments in terms of changes in amino acid composition of the proteins.

Insights into thermoadaptation and the evolution of mesophily from the bacterial phylum Thermotogae

Thermophiles are extremophiles that grow optimally at temperatures >45 °C. To survive and maintain function of their biological molecules, they have a suite of characteristics not found in organisms that grow at moderate temperature (mesophiles). At the cellular level, thermophiles have mechanisms for maintaining their membranes, nucleic acids, and other cellular structures. At the protein level, each of their proteins remains stable and retains activity at temperatures that would denature their mesophilic homologs. Conversely, cellular structures and proteins from thermophiles may not function optimally at moderate temperatures. These differences between thermophiles and mesophiles presumably present a barrier for evolutionary transitioning between the 2 lifestyles. Therefore, studying closely related thermophiles and mesophiles can help us determine how such lifestyle transitions may happen. The bacterial phylum Thermotogae contains hyperthermophiles, thermophiles, mesophiles, and organisms with temperature ranges wide enough to span both thermophilic and mesophilic temperatures. Genomic, proteomic, and physiological differences noted between other bacterial thermophiles and mesophiles are evident within the Thermotogae. We argue that the Thermotogae is an ideal group of organisms for understanding of the response to fluctuating temperature and of long-term evolutionary adaptation to a different growth temperature range.

Evolution of complete proteomes: guanine-cytosine pressure, phylogeny and environmental influences blend the proteomic architecture

Background

Guanine-cytosine (GC) composition is an important feature of genomes. Likewise, amino acid composition is a distinct, but less valued, feature of proteomes. A major concern is that it is not clear what valuable information can be acquired from amino acid composition data. To address this concern, in-depth analyses of the amino acid composition of the complete proteomes from 63 archaea, 270 bacteria, and 128 eukaryotes were performed.

Results

Principal component analysis of the amino acid matrices showed that the main contributors to proteomic architecture were genomic GC variation, phylogeny, and environmental influences. GC pressure drove positive selection on Ala, Arg, Gly, Pro, Trp, and Val, and adverse selection on Asn, Lys, Ile, Phe, and Tyr. The physico-chemical framework of the complete proteomes withstood GC pressure by frequency complementation of GC-dependent amino acid pairs with similar physico-chemical properties. Gln, His, Ser, and Val were responsible for phylogeny and their constituted components could differentiate archaea, bacteria, and eukaryotes. Environmental niche was also a significant factor in determining proteomic architecture, especially for archaea for which the main amino acids were Cys, Leu, and Thr. In archaea, hyperthermophiles, acidophiles, mesophiles, psychrophiles, and halophiles gathered successively along the environment-based principal component. Concordance between proteomic architecture and the genetic code was also related closely to genomic GC content, phylogeny, and lifestyles.

Conclusions

Large-scale analyses of the complete proteomes of a wide range of organisms suggested that amino acid composition retained the trace of GC variation, phylogeny, and environmental influences during evolution. The findings from this study will help in the development of a global understanding of proteome evolution, and even biological evolution.

Horizontal gene transfer of a chloroplast DnaJ-Fer protein to Thaumarchaeota and the evolutionary history of the DnaK chaperone system in Archaea

Background

In 2004, we discovered an atypical protein in metagenomic data from marine thaumarchaeotal species. This protein, referred as DnaJ-Fer, is composed of a J domain fused to a Ferredoxin (Fer) domain. Surprisingly, the same protein was also found in Viridiplantae (green algae and land plants). Because J domain-containing proteins are known to interact with the major chaperone DnaK/Hsp70, this suggested that a DnaK protein was present in Thaumarchaeota. DnaK/Hsp70, its co-chaperone DnaJ and the nucleotide exchange factor GrpE are involved, among others, in heat shocks and heavy metal cellular stress responses.

Results

Using phylogenomic approaches we have investigated the evolutionary history of the DnaJ-Fer protein and of interacting proteins DnaK, DnaJ and GrpE in Thaumarchaeota. These proteins have very complex histories, involving several inter-domain horizontal gene transfers (HGTs) to explain the contemporary distribution of these proteins in archaea. These transfers include one from Cyanobacteria to Viridiplantae and one from Viridiplantae to Thaumarchaeota for the DnaJ-Fer protein, as well as independent HGTs from Bacteria to mesophilic archaea for the DnaK/DnaJ/GrpE system, followed by HGTs among mesophilic and thermophilic archaea.

Conclusions

We highlight the chimerical origin of the set of proteins DnaK, DnaJ, GrpE and DnaJ-Fer in Thaumarchaeota and suggest that the HGT of these proteins has played an important role in the adaptation of several archaeal groups to mesophilic and thermophilic environments from hyperthermophilic ancestors. Finally, the evolutionary history of DnaJ-Fer provides information useful for the relative dating of the diversification of Archaeplastida and Thaumarchaeota.

Effect of codon message on xylanase thermal activity

Because the genetic codon is known for degeneracy, its effect on enzyme thermal property is seldom investigated. A dataset was constructed for GH10 xylanase coding sequences and optimal temperatures for activity (T(opt)). Codon contents and relative synonymous codon usages were calculated and respectively correlated with the enzyme T(opt) values, which were used to describe the xylanase thermophilic tendencies without dividing them into two thermophilic and mesophilic groups. After analyses of codon content and relative synonymous codon usages were checked by the Bonferroni correction, we found five codons, with three (AUA, AGA, and AGG) correlating positively and two (CGU and AGC) correlating negatively with the T(opt) value. The three positive codons are purine-rich codons, and the two negative codons have A-ends. The two negative codons are pyridine-rich codons, and one has a C-end. Comparable with the codon C- and A-ending features, C- and A-content within mRNA correlated negatively and positively with the T(opt) value, respectively. Thereby, codons have effects on enzyme thermal property. When the issue is analyzed at the residual level, the effect of codon message is lost. The codons relating to enzyme thermal property are selected by thermophilic force at nucleotide level.

The tree and net components of prokaryote evolution

Phylogenetic trees of individual genes of prokaryotes (archaea and bacteria) generally have different topologies, largely owing to extensive horizontal gene transfer (HGT), suggesting that the Tree of Life (TOL) should be replaced by a "net of life" as the paradigm of prokaryote evolution. However, trees remain the natural representation of the histories of individual genes given the fundamentally bifurcating process of gene replication. Therefore, although no single tree can fully represent the evolution of prokaryote genomes, the complete picture of evolution will necessarily combine trees and nets. A quantitative measure of the signals of tree and net evolution is derived from an analysis of all quartets of species in all trees of the "Forest of Life" (FOL), which consists of approximately 7,000 phylogenetic trees for prokaryote genes including approximately 100 nearly universal trees (NUTs). Although diverse routes of net-like evolution collectively dominate the FOL, the pattern of tree-like evolution that reflects the consistent topologies of the NUTs is the most prominent coherent trend. We show that the contributions of tree-like and net-like evolutionary processes substantially differ across bacterial and archaeal lineages and between functional classes of genes. Evolutionary simulations indicate that the central tree-like signal cannot be realistically explained by a self-reinforcing pattern of biased HGT.

Reduction in structural disorder and functional complexity in the thermal adaptation of prokaryotes

Genomic correlates of evolutionary adaptation to very low or very high optimal growth temperature (OGT) values have been the subject of many studies. Whereas these provided a protein-structural rationale of the activity and stability of globular proteins/enzymes, the point has been neglected that adaptation to extreme temperatures could also have resulted from an increased use of intrinsically disordered proteins (IDPs), which are resistant to these conditions in vitro. Contrary to these expectations, we found a conspicuously low level of structural disorder in bacteria of very high (and very low) OGT values. This paucity of disorder does not reflect phylogenetic relatedness, i.e. it is a result of genuine adaptation to extreme conditions. Because intrinsic disorder correlates with important regulatory functions, we asked how these bacteria could exist without IDPs by studying transcription factors, known to harbor a lot of function-related intrinsic disorder. Hyperthermophiles have much less transcription factors, which have reduced disorder compared to their mesophilic counterparts. On the other hand, we found by systematic categorization of proteins with long disordered regions that there are certain functions, such as translation and ribosome biogenesis that depend on structural disorder even in hyperthermophiles. In all, our observations suggest that adaptation to extreme conditions is achieved by a significant functional simplification, apparent at both the level of the genome and individual genes/proteins.

Temperature adaptation at homologous sites in proteins from nine thermophile-mesophile species pairs

Whether particular amino acids are favored by selection at high temperatures over others has long been an open question in protein evolution. One way to approach this question is to compare homologous sites in proteins from one thermophile and a closely related mesophile; asymmetrical substitution patterns have been taken as evidence for selection favoring certain amino acids over others. However, most pairs of prokaryotic species that differ in optimum temperature also differ in genome-wide GC content, and amino acid content is known to be associated with GC content. Here, I compare homologous sites in nine thermophilic prokaryotes and their mesophilic relatives, all with complete published genome sequences. After adjusting for the effects of differing GC content with logistic regression, 139 of the 190 pairs of amino acids show significant substitutional asymmetry, evidence of widespread adaptive amino acid substitution. The patterns are fairly consistent across the nine pairs of species (after taking the effects of differing GC content into account), suggesting that much of the asymmetry results from adaptation to temperature. Some amino acids in some species pairs deviate from the overall pattern in ways indicating that adaptation to other environmental or physiological differences between the species may also play a role. The property that is best correlated with the patterns of substitutional asymmetry is transfer free energy, a measure of hydrophobicity, with more hydrophobic amino acids favored at higher temperatures. The correlation of asymmetry and hydrophobicity is fairly weak, suggesting that other properties may also be important.

Distance-dependent statistical potentials for discriminating thermophilic and mesophilic proteins

Identification of the characteristic structural patterns responsible for protein thermostability is theoretically important and practically useful but largely remains an open problem. These patterns may be revealed through comparative study on thermophilic and mesophilic proteins that have distinct thermostability. In this study, we constructed several distance-dependant potentials from thermophilic and mesophilic proteins. These potentials were then used to evaluate the structural difference between thermophilic and mesophilic proteins. We found that using the subtraction or division of the potentials derived from thermophilic and mesophilic proteins can dramatically increase the discriminatory ability. This approach revealed that the ability to distinct the subtle structural features responsible for protein thermostability may be effectively enhanced through rationally designed comparative study.

Shape and evolution of thermostable protein structure

Organisms evolved at high temperatures must maintain their proteins' structures in the face of increased thermal disorder. This challenge results in differences in residue utilization and overall structure. Focusing on thermostable/mesostable pairs of homologous structures, we have examined these differences using novel geometric measures: specifically burial depth (distance from the molecular surface to each atom) and travel depth (distance from the convex hull to the molecular surface that avoids the protein interior). These along with common metrics like packing and Wadell Sphericity are used to gain insight into the constraints experienced by thermophiles. Mean travel depth of hyperthermostable proteins is significantly less than that of their mesostable counterparts, indicating smaller, less numerous and less deep pockets. The mean burial depth of hyperthermostable proteins is significantly higher than that of mesostable proteins indicating that they bury more atoms further from the surface. The burial depth can also be tracked on the individual residue level, adding a finer level of detail to the standard exposed surface area analysis. Hyperthermostable proteins for the first time are shown to be more spherical than their mesostable homologues, regardless of when and how they adapted to extreme temperature. Additionally, residue specific burial depth examinations reveal that charged residues stay unburied, most other residues are slightly more buried and Alanine is more significantly buried in hyperthermostable proteins.

Evolution: like any other science it is predictable

Evolutionary biology rejoices in the diversity of life, but this comes at a cost: other than working in the common framework of neo-Darwinian evolution, specialists in, for example, diatoms and mammals have little to say to each other. Accordingly, their research tends to track the particularities and peculiarities of a given group and seldom enquires whether there are any wider or deeper sets of explanations. Here, I present evidence in support of the heterodox idea that evolution might look to a general theory that does more than serve as a tautology ('evolution explains evolution'). Specifically, I argue that far from its myriad of products being fortuitous and accidental, evolution is remarkably predictable. Thus, I urge a move away from the continuing obsession with Darwinian mechanisms, which are entirely uncontroversial. Rather, I emphasize why we should seek explanations for ubiquitous evolutionary convergence, as well as the emergence of complex integrated systems. At present, evolutionary theory seems to be akin to nineteenth-century physics, blissfully unaware of the imminent arrival of quantum mechanics and general relativity. Physics had its Newton, biology its Darwin: evolutionary biology now awaits its Einstein.

Trade-off between positive and negative design of protein stability: from lattice models to real proteins

Two different strategies for stabilizing proteins are (i) positive design in which the native state is stabilized and (ii) negative design in which competing non-native conformations are destabilized. Here, the circumstances under which one strategy might be favored over the other are explored in the case of lattice models of proteins and then generalized and discussed with regard to real proteins. The balance between positive and negative design of proteins is found to be determined by their average "contact-frequency", a property that corresponds to the fraction of states in the conformational ensemble of the sequence in which a pair of residues is in contact. Lattice model proteins with a high average contact-frequency are found to use negative design more than model proteins with a low average contact-frequency. A mathematical derivation of this result indicates that it is general and likely to hold also for real proteins. Comparison of the results of correlated mutation analysis for real proteins with typical contact-frequencies to those of proteins likely to have high contact-frequencies (such as disordered proteins and proteins that are dependent on chaperonins for their folding) indicates that the latter tend to have stronger interactions between residues that are not in contact in their native conformation. Hence, our work indicates that negative design is employed when insufficient stabilization is achieved via positive design owing to high contact-frequencies.