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

Deciphering Structural Traits for Thermal and Kinetic Stability across Protein Family Evolution through Ancestral Sequence Reconstruction

Natural proteins are frequently marginally stable, and an increase in environmental temperature can easily lead to unfolding. As a result, protein engineering to improve protein stability is an area of intensive research. Nonetheless, since there is usually a high degree of structural homology between proteins from thermophilic organisms and their mesophilic counterparts, the identification of structural determinants for thermoadaptation is challenging. Moreover, in many cases, it has become clear that the success of stabilization strategies is often dependent on the evolutionary history of a protein family. In the last few years, the use of ancestral sequence reconstruction (ASR) as a tool for elucidation of the evolutionary history of functional traits of a protein family has gained strength. Here, we used ASR to trace the evolutionary pathways between mesophilic and thermophilic kinases that participate in the biosynthetic pathway of vitamin B1 in bacteria. By combining biophysics approaches, X-ray crystallography, and molecular dynamics simulations, we found that the thermal stability of these enzymes correlates with their kinetic stability, where the highest thermal/kinetic stability is given by an increase in small hydrophobic amino acids that allow a higher number of interatomic hydrophobic contacts, making this type of interaction the main support for stability in this protein architecture. The results highlight the potential benefits of using ASR to explore the evolutionary history of protein sequence and structure to identify traits responsible for the kinetic and thermal stability of any protein architecture.

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.

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.

Crystal structure of chloroplastic thioredoxin z defines a type-specific target recognition

Thioredoxins (TRXs) are ubiquitous disulfide oxidoreductases structured according to a highly conserved fold. TRXs are involved in a myriad of different processes through a common chemical mechanism. Plant TRXs evolved into seven types with diverse subcellular localization and distinct protein target selectivity. Five TRX types coexist in the chloroplast, with yet scarcely described specificities. We solved the crystal structure of a chloroplastic z-type TRX, revealing a conserved TRX fold with an original electrostatic surface potential surrounding the redox site. This recognition surface is distinct from all other known TRX types from plant and non-plant sources and is exclusively conserved in plant z-type TRXs. We show that this electronegative surface endows thioredoxin z (TRXz) with a capacity to activate the photosynthetic Calvin-Benson cycle enzyme phosphoribulokinase. The distinct electronegative surface of TRXz thereby extends the repertoire of TRX-target recognitions.

Evolution, folding, and design of TIM barrels and related proteins

Proteins are chief actors in life that perform a myriad of exquisite functions. This diversity has been enabled through the evolution and diversification of protein folds. Analysis of sequences and structures strongly suggest that numerous protein pieces have been reused as building blocks and propagated to many modern folds. This information can be traced to understand how the protein world has diversified. In this review, we discuss the latest advances in the analysis of protein evolutionary units, and we use as a model system one of the most abundant and versatile topologies, the TIM-barrel fold, to highlight the existing common principles that interconnect protein evolution, structure, folding, function, and design.

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.

•

Phage adapts to a host modified to hinder the essential recruitment of a host protein

•

No genetic mutations are fixed at the engineered virus-host interaction interface

•

Adaptation is likely linked to phenotypic mutations caused by transcription errors

•

Sub-genomic RNAs may enable this kind of adaptation mechanism in coronaviruses

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.

Family 1 glycosidases (GH1) are present in the three domains of life and share classical TIM-barrel fold. Structural and biochemical analyses of a resurrected ancestral GH1 enzyme reveal heme binding, not known in its modern descendants. Heme rigidifies the TIM-barrel and allosterically enhances catalysis.

Protein folding stability and binding interactions through the lens of evolution: a dynamical perspective

While the function of a protein depends heavily on its ability to fold into a correct 3D structure, billions of years of evolution have tailored proteins from highly stable objects to flexible molecules as they adapted to environmental changes. Nature maintains the fine balance of protein folding and stability while still evolving towards new function through generations of fine-tuning necessary interactions with other proteins and small molecules. Here we focus on recent computational and experimental studies that shed light onto how evolution molds protein folding and the functional landscape from a conformational dynamics' perspective. Particularly, we explore the importance of dynamic allostery throughout protein evolution and discuss how the protein anisotropic network can give rise to allosteric and epistatic interactions.

Molecular origins of folding rate differences in the thioredoxin family

Thioredoxins are a family of conserved oxidoreductases responsible for maintaining redox balance within cells. They have also served as excellent model systems for protein design and engineering studies particularly through ancestral sequence reconstruction methods. The recent work by Gamiz-Arco et al. [Biochem J (2019) 476, 3631–3647] answers fundamental questions on how specific sequence differences can contribute to differences in folding rates between modern and ancient thioredoxins but also among a selected subset of modern thioredoxins. They surprisingly find that rapid unassisted folding, a feature of ancient thioredoxins, is not conserved in the modern descendants suggestive of co-evolution of better folding machinery that likely enabled the accumulation of mutations that slow-down folding. The work thus provides an interesting take on the expected folding-stability-function constraint while arguing for additional factors that contribute to sequence evolution and hence impact folding efficiency.