Cell Survival Enabled by Leakage of a Labile Metabolic Intermediate.

Many metabolites are generated in one step of a biochemical pathway and consumed in a subsequent step. Such metabolic intermediates are often reactive molecules which, if allowed to freely diffuse in the intracellular milieu, could lead to undesirable side reactions and even become toxic to the cell. Therefore, metabolic intermediates are often protected as protein-bound species and directly transferred between enzyme active sites in multi-functional enzymes, multi-enzyme complexes, and metabolons. Sequestration of reactive metabolic intermediates thus contributes to metabolic efficiency. It is not known, however, whether this evolutionary adaptation can be relaxed in response to challenges to organismal survival. Here, we report evolutionary repair experiments on Escherichia coli cells in which an enzyme crucial for the biosynthesis of proline has been deleted. The deletion makes cells unable to grow in a culture medium lacking proline. Remarkably, however, cell growth is efficiently restored by many single mutations (12 at least) in the gene of glutamine synthetase. The mutations cause the leakage to the intracellular milieu of a highly reactive phosphorylated intermediate common to the biosynthetic pathways of glutamine and proline. This intermediate is generally assumed to exist only as a protein-bound species. Nevertheless, its diffusion upon mutation-induced leakage enables a new route to proline biosynthesis. Our results support that leakage of sequestered metabolic intermediates can readily occur and contribute to organismal adaptation in some scenarios. Enhanced availability of reactive molecules may enable the generation of new biochemical pathways and the potential of mutation-induced leakage in metabolic engineering is noted.

Folding Free Energy Surfaces from Differential Scanning Calorimetry.

Protein folding/unfolding processes involve a large number of weak, non-covalent interactions and are more appropriately described in terms of the movement of a point representing protein conformation in a plot of internal free energy versus conformational degrees of freedom. While these energy landscapes have an astronomically large number of dimensions, it has been shown that many relevant aspects of protein folding can be understood in terms of their projections onto a few relevant coordinates. Remarkably, such low-dimensional free energy surfaces can be obtained from experimental DSC data using suitable analytical models. Here, we describe the experimental procedures to be used to obtain the high-quality DSC data that are required for free-energy surface analysis.

Combining Ancestral Reconstruction with Folding-Landscape Simulations to Engineer Heterologous Protein Expression.

Obligate symbionts typically exhibit high evolutionary rates. Consequently, their proteins may differ considerably from their modern and ancestral homologs in terms of both sequence and properties, thus providing excellent models to study protein evolution. Also, obligate symbionts are challenging to culture in the lab and proteins from uncultured organisms must be produced in heterologous hosts using recombinant DNA technology. Obligate symbionts thus replicate a fundamental scenario of metagenomics studies aimed at the functional characterization and biotechnological exploitation of proteins from the bacteria in soil. Here, we use the thioredoxin from Candidatus Photodesmus katoptron, an uncultured symbiont of flashlight fish, to explore evolutionary and engineering aspects of protein folding in heterologous hosts. The symbiont protein is a standard thioredoxin in terms of 3D-structure, stability and redox activity. However, its folding outside the original host is severely impaired, as shown by a very slow refolding in vitro and an inefficient expression in E. coli that leads mostly to insoluble protein. By contrast, resurrected Precambrian thioredoxins express efficiently in E. coli, plausibly reflecting an ancient adaptation to unassisted folding. We have used a statistical-mechanical model of the folding landscape to guide back-to-ancestor engineering of the symbiont protein. Remarkably, we find that the efficiency of heterologous expression correlates with the in vitro (i.e., unassisted) folding rate and that the ancestral expression efficiency can be achieved with only 1-2 back-to-ancestor replacements. These results demonstrate a minimal-perturbation, sequence-engineering approach to rescue inefficient heterologous expression which may potentially be useful in metagenomics efforts targeting recent adaptations.

A protocol to study bacteriophage adaptation to new hosts.

A general protocol for the experimental assessment of bacteriophage adaptation to new hosts is described. We use as a model system the lytic phage T7 and an engineered E. coli strain modified to hamper the recruitment of a known proviral factor. Our protocol includes steps of phage amplification, plaque and liquid lysis assays, and DNA extraction for next-generation sequencing of the viral genome over several rounds of laboratory evolution thus allowing the investigation of the sequence determinants of viral adaptation. For complete information on the generation and use of this protocol, please refer to Luzon-Hidalgo et al. (2021).

Evidence for a role of phenotypic mutations in virus adaptation.

Viruses interact extensively with the host molecular machinery, but the underlying mechanisms are poorly understood. Bacteriophage T7 recruits the small protein thioredoxin of the Escherichia coli host as an essential processivity factor for the viral DNA polymerase. We challenged the phage to propagate in a host in which thioredoxin had been extensively modified to hamper its recruitment. The virus adapted to the engineered host without losing the capability to propagate in the original host, but no genetic mutations were fixed in the thioredoxin binding domain of the viral DNA polymerase. Virus adaptation correlated with mutations in the viral RNA polymerase, supporting that promiscuous thioredoxin recruitment was enabled by phenotypic mutations caused by transcription errors. These results point to a mechanism of virus adaptation that may play a role in cross-species transmission. We propose that phenotypic mutations may generally contribute to the capability of viruses to evade antiviral strategies.

Heme-binding enables allosteric modulation in an ancient TIM-barrel glycosidase.

Glycosidases are phylogenetically widely distributed enzymes that are crucial for the cleavage of glycosidic bonds. Here, we present the exceptional properties of a putative ancestor of bacterial and eukaryotic family-1 glycosidases. The ancestral protein shares the TIM-barrel fold with its modern descendants but displays large regions with greatly enhanced conformational flexibility. Yet, the barrel core remains comparatively rigid and the ancestral glycosidase activity is stable, with an optimum temperature within the experimental range for thermophilic family-1 glycosidases. None of the ∼5500 reported crystallographic structures of ∼1400 modern glycosidases show a bound porphyrin. Remarkably, the ancestral glycosidase binds heme tightly and stoichiometrically at a well-defined buried site. Heme binding rigidifies this TIM-barrel and allosterically enhances catalysis. Our work demonstrates the capability of ancestral protein reconstructions to reveal valuable but unexpected biomolecular features when sampling distant sequence space. The potential of the ancestral glycosidase as a scaffold for custom catalysis and biosensor engineering is discussed.

Engineering protein assemblies with allosteric control via monomer fold-switching.

The macromolecular machines of life use allosteric control to self-assemble, dissociate and change shape in response to signals. Despite enormous interest, the design of nanoscale allosteric assemblies has proven tremendously challenging. Here we present a proof of concept of allosteric assembly in which an engineered fold switch on the protein monomer triggers or blocks assembly. Our design is based on the hyper-stable, naturally monomeric protein CI2, a paradigm of simple two-state folding, and the toroidal arrangement with 6-fold symmetry that it only adopts in crystalline form. We engineer CI2 to enable a switch between the native and an alternate, latent fold that self-assembles onto hexagonal toroidal particles by exposing a favorable inter-monomer interface. The assembly is controlled on demand via the competing effects of temperature and a designed short peptide. These findings unveil a remarkable potential for structural metamorphosis in proteins and demonstrate key principles for engineering protein-based nanomachinery.

Non-conservation of folding rates in the thioredoxin family reveals degradation of ancestral unassisted-folding.

Evolution involves not only adaptation, but also the degradation of superfluous features. Many examples of degradation at the morphological level are known (vestigial organs, for instance). However, the impact of degradation on molecular evolution has been rarely addressed. Thioredoxins serve as general oxidoreductases in all cells. Here, we report extensive mutational analyses on the folding of modern and resurrected ancestral bacterial thioredoxins. Contrary to claims from recent literature, in vitro folding rates in the thioredoxin family are not evolutionarily conserved, but span at least a ∼100-fold range. Furthermore, modern thioredoxin folding is often substantially slower than ancestral thioredoxin folding. Unassisted folding, as probed in vitro, thus emerges as an ancestral vestigial feature that underwent degradation, plausibly upon the evolutionary emergence of efficient cellular folding assistance. More generally, our results provide evidence that degradation of ancestral features shapes, not only morphological evolution, but also the evolution of individual proteins.

Using Resurrected Ancestral Proviral Proteins to Engineer Virus Resistance.

Proviral factors are host proteins hijacked by viruses for processes essential for virus propagation such as cellular entry and replication. Pathogens and their hosts co-evolve. It follows that replacing a proviral factor with a functional ancestral form of the same protein could prevent viral propagation without fatally compromising organismal fitness. Here, we provide proof of concept of this notion. Thioredoxins serve as general oxidoreductases in all known cells. We report that several laboratory resurrections of Precambrian thioredoxins display substantial levels of functionality within Escherichia coli. Unlike E. coli thioredoxin, however, these ancestral thioredoxins are not efficiently recruited by the bacteriophage T7 for its replisome and therefore prevent phage propagation in E. coli. These results suggest an approach to the engineering of virus resistance. Diseases caused by viruses may have a devastating effect in agriculture. We discuss how the suggested approach could be applied to the engineering of plant virus resistance.

Engineering ancestral protein hyperstability.

Many experimental analyses and proposed scenarios support that ancient life was thermophilic. In congruence with this hypothesis, proteins encoded by reconstructed sequences corresponding to ancient phylogenetic nodes often display very high stability. Here, we show that such 'reconstructed ancestral hyperstability' can be further engineered on the basis of a straightforward approach that uses exclusively information afforded by the ancestral reconstruction process itself. Since evolution does not imply continuous progression, screening of the mutations between two evolutionarily related resurrected ancestral proteins may identify mutations that further stabilize the most stable one. To explore this approach, we have used a resurrected thioredoxin corresponding to the last common ancestor of the cyanobacterial, Deinococcus and Thermus groups (LPBCA thioredoxin), which has a denaturation temperature of ∼123°C. This high value is within the top 0.1% of the denaturation temperatures in the ProTherm database and, therefore, achieving further stabilization appears a priori as a challenging task. Nevertheless, experimental comparison with a resurrected thioredoxin corresponding to the last common ancestor of bacteria (denaturation temperature of ∼115°C) immediately identifies three mutations that increase the denaturation temperature of LPBCA thioredoxin to ∼128°C. Comparison between evolutionarily related resurrected ancestral proteins thus emerges as a simple approach to expand the capability of ancestral reconstruction to search sequence space for extreme protein properties of biotechnological interest. The fact that ancestral sequences for many phylogenetic nodes can be derived from a single alignment of modern sequences should contribute to the general applicability of this approach.

Selection for Protein Kinetic Stability Connects Denaturation Temperatures to Organismal Temperatures and Provides Clues to Archaean Life.

The relationship between the denaturation temperatures of proteins (Tm values) and the living temperatures of their host organisms (environmental temperatures: TENV values) is poorly understood. Since different proteins in the same organism may show widely different Tm's, no simple universal relationship between Tm and TENV should hold, other than Tm≥TENV. Yet, when analyzing a set of homologous proteins from different hosts, Tm's are oftentimes found to correlate with TENV's but this correlation is shifted upward on the Tm axis. Supporting this trend, we recently reported Tm's for resurrected Precambrian thioredoxins that mirror a proposed environmental cooling over long geological time, while remaining a shocking ~50°C above the proposed ancestral ocean temperatures. Here, we show that natural selection for protein kinetic stability (denaturation rate) can produce a Tm↔TENV correlation with a large upward shift in Tm. A model for protein stability evolution suggests a link between the Tm shift and the in vivo lifetime of a protein and, more specifically, allows us to estimate ancestral environmental temperatures from experimental denaturation rates for resurrected Precambrian thioredoxins. The TENV values thus obtained match the proposed ancestral ocean cooling, support comparatively high Archaean temperatures, and are consistent with a recent proposal for the environmental temperature (above 75°C) that hosted the last universal common ancestor. More generally, this work provides a framework for understanding how features of protein stability reflect the environmental temperatures of the host organisms.

Modern Analysis of Protein Folding by Differential Scanning Calorimetry.

Differential scanning calorimetry (DSC) is a very powerful tool for investigating protein folding and stability because its experimental output reflects the energetics of all conformations that become minimally populated during thermal unfolding. Accordingly, analysis of DSC experiments with simple thermodynamic models has been key for developing our understanding of protein stability during the past five decades. The discovery of ultrafast folding proteins, which have naturally broad conformational ensembles and minimally cooperative unfolding, opens the possibility of probing the complete folding free energy landscape, including those conformations at the top of the barrier to folding, via DSC. Exploiting this opportunity requires high-quality experiments and the implementation of novel analytical methods based on statistical mechanics. Here, we cover the recent exciting developments in this front, describing the new analytical procedures in detail as well as providing experimental guidelines for performing such analysis.

Mutational studies on resurrected ancestral proteins reveal conservation of site-specific amino acid preferences throughout evolutionary history.

Local protein interactions ("molecular context" effects) dictate amino acid replacements and can be described in terms of site-specific, energetic preferences for any different amino acid. It has been recently debated whether these preferences remain approximately constant during evolution or whether, due to coevolution of sites, they change strongly. Such research highlights an unresolved and fundamental issue with far-reaching implications for phylogenetic analysis and molecular evolution modeling. Here, we take advantage of the recent availability of phenotypically supported laboratory resurrections of Precambrian thioredoxins and ß-lactamases to experimentally address the change of site-specific amino acid preferences over long geological timescales. Extensive mutational analyses support the notion that evolutionary adjustment to a new amino acid may occur, but to a large extent this is insufficient to erase the primitive preference for amino acid replacements. Generally, site-specific amino acid preferences appear to remain conserved throughout evolutionary history despite local sequence divergence. We show such preference conservation to be readily understandable in molecular terms and we provide crystallographic evidence for an intriguing structural-switch mechanism: Energetic preference for an ancestral amino acid in a modern protein can be linked to reorganization upon mutation to the ancestral local structure around the mutated site. Finally, we point out that site-specific preference conservation naturally leads to one plausible evolutionary explanation for the existence of intragenic global suppressor mutations.

Identification of the structural basis of thermal lability of a virus provides a rationale for improved vaccines.

Virus stability and dynamics play critical roles during infection. Some viruses, including foot-and-mouth disease virus (FMDV), are surprisingly prone to thermal dissociation outside the cell. The structural bases and functional implications of this distinctive trait were essentially unknown. This study (1) uncovers the structural determinants of FMDV thermolability, (2) investigates the relationship between virus thermolability and infectivity, and (3) provides a structure-based rationale for engineering thermostable virus particles to develop improved vaccines and nanocontainers. The results reveal that negatively charged residues close to protein-protein interfaces exert electrostatic repulsions between capsid subunits and mediate the sensitivity of the virion to thermal dissociation, even at neutral pH. Based on these results, a series of fully infectious virions of increased thermostability were engineered by individually removing different carboxylates involved in intersubunit repulsions. The implications for virus biology and the design of thermostable vaccines are discussed.

Conservation of protein structure over four billion years.

Little is known about the evolution of protein structures and the degree of protein structure conservation over planetary time scales. Here, we report the X-ray crystal structures of seven laboratory resurrections of Precambrian thioredoxins dating up to approximately four billion years ago. Despite considerable sequence differences compared with extant enzymes, the ancestral proteins display the canonical thioredoxin fold, whereas only small structural changes have occurred over four billion years. This remarkable degree of structure conservation since a time near the last common ancestor of life supports a punctuated-equilibrium model of structure evolution in which the generation of new folds occurs over comparatively short periods and is followed by long periods of structural stasis.

Probing the mutational interplay between primary and promiscuous protein functions: a computational-experimental approach.

Protein promiscuity is of considerable interest due its role in adaptive metabolic plasticity, its fundamental connection with molecular evolution and also because of its biotechnological applications. Current views on the relation between primary and promiscuous protein activities stem largely from laboratory evolution experiments aimed at increasing promiscuous activity levels. Here, on the other hand, we attempt to assess the main features of the simultaneous modulation of the primary and promiscuous functions during the course of natural evolution. The computational/experimental approach we propose for this task involves the following steps: a function-targeted, statistical coupling analysis of evolutionary data is used to determine a set of positions likely linked to the recruitment of a promiscuous activity for a new function; a combinatorial library of mutations on this set of positions is prepared and screened for both, the primary and the promiscuous activities; a partial-least-squares reconstruction of the full combinatorial space is carried out; finally, an approximation to the Pareto set of variants with optimal primary/promiscuous activities is derived. Application of the approach to the emergence of folding catalysis in thioredoxin scaffolds reveals an unanticipated scenario: diverse patterns of primary/promiscuous activity modulation are possible, including a moderate (but likely significant in a biological context) simultaneous enhancement of both activities. We show that this scenario can be most simply explained on the basis of the conformational diversity hypothesis, although alternative interpretations cannot be ruled out. Overall, the results reported may help clarify the mechanisms of the evolution of new functions. From a different viewpoint, the partial-least-squares-reconstruction/Pareto-set-prediction approach we have introduced provides the computational basis for an efficient directed-evolution protocol aimed at the simultaneous enhancement of several protein features and should therefore open new possibilities in the engineering of multi-functional enzymes.

How many ionizable groups can sit on a protein hydrophobic core?

Full or partial burial of ionizable groups in the hydrophobic interior of proteins underlies the large modulation in group properties (modified pK value, high nucleophilicity, enhanced capability of interaction with chemical moieties of the substrate, etc.) linked to biological function. Indeed, the few internal ionizable residues found in proteins are known to play important functional roles in catalysis and, in general, in energy transduction processes. However, ionizable-group burial is expected to be seriously disruptive and, it is important to note, most functional sites contain not just one, but several ionizable residues. Hence, the adaptations involved in the development of function in proteins (through in vitro engineering or during the course of natural evolution) are not fully understood. Here, we explore experimentally how proteins respond to the accumulation of hydrophobic-to-ionizable residue substitutions. For this purpose, we have constructed a combinatorial library targeting a hydrophobic cluster in a consensus-engineered, stabilized form of a small model protein. Contrary to naïve expectation, half of the variants randomly selected from the library are soluble, folded, and active, despite including up to four mutations. Furthermore, for these variants, the dependence of stability with the number of mutations is not synergistic and catastrophic, but smooth and approximately linear. Clearly, stabilized protein scaffolds may be robust enough to withstand many disruptive hydrophobic-to-ionizable residue mutations, even when they are introduced in the same region of the structure. These results should be relevant for protein engineering and may have implications for the understanding of the early evolution of enzymes.

Highly anomalous energetics of protein cold denaturation linked to folding-unfolding kinetics.

Despite several careful experimental analyses, it is not yet clear whether protein cold-denaturation is just a "mirror image" of heat denaturation or whether it shows unique structural and energetic features. Here we report that, for a well-characterized small protein, heat denaturation and cold denaturation show dramatically different experimental energetic patterns. Specifically, while heat denaturation is endothermic, the cold transition (studied in the folding direction) occurs with negligible heat effect, in a manner seemingly akin to a gradual, second-order-like transition. We show that this highly anomalous energetics is actually an apparent effect associated to a large folding/unfolding free energy barrier and that it ultimately reflects kinetic stability, a naturally-selected trait in many protein systems. Kinetics thus emerges as an important factor linked to differential features of cold denaturation. We speculate that kinetic stabilization against cold denaturation may play a role in cold adaptation of psychrophilic organisms. Furthermore, we suggest that folding-unfolding kinetics should be taken into account when analyzing in vitro cold-denaturation experiments, in particular those carried out in the absence of destabilizing conditions.

Proteolytic scanning calorimetry: a novel methodology that probes the fundamental features of protein kinetic stability.

We introduce proteolytic scanning calorimetry, a modification of the differential scanning calorimetry approach to the determination of protein stability in which a proteolytic enzyme (thermolysin) is used to mimic a harsh environment. This methodology allows the straightforward calculation of the rate of irreversible denaturation as a function of temperature and concentration of proteolytic enzyme and, as a result, has the potential to probe efficiently the fundamental biophysical features of protein kinetic stability. In the particular case of Escherichia coli thioredoxin (used as an illustrative example in this article), we find that the rate of irreversible denaturation is determined by 1), the global unfolding mechanism at low thermolysin concentrations, indicating that thermodynamic stability may contribute directly to the kinetic stability of thioredoxin under moderately harsh conditions and 2), the rate of unfolding at high thermolysin concentrations, indicating that the free-energy barrier for unfolding may act as a safety mechanism that ensures significant kinetic stability, even in very harsh environments. This thioredoxin picture, however, is by no means expected to be general and different proteins may show different patterns of kinetic stabilization. Proteolytic scanning calorimetry is particularly well-suited to probe this diversity at a fundamental biophysical level.