Combining ancestral sequence reconstruction with protein design to identify an interface hotspot in a key metabolic enzyme complex

Preparation and Characterization of an Ancient Aminopeptidase Obtained from Ancestral Sequence Reconstruction for L-Carnosine Synthesis

As a biologically active peptide, L-carnosine has been widely used in the pharmaceutical, cosmetic and health care industries due to its various physiological properties. However, relatively little research is available regarding L-carnosine's enzymatic synthesis function. In this study, a potential enzyme sequence with the function of carnosine synthesizing was screened out using the ancestral sequence reconstruction (ASR) technique. Identified with L-carnosine synthesis activity, this enzyme was further confirmed using autoproteolytic phenomenon via Western blot and N-terminal sequencing. After purification, the enzymatic properties of LUCA-DmpA were characterized. The melting temperature (T_m) and denaturation enthalpy (ΔH) of LUCA-DmpA were 60.27 ± 1.24 °C and 1306.00 ± 26.73 kJ·mol^-1, respectively. Circular dichroism (CD) spectroscopy results showed that this ancestral enzyme was composed of α-helix (35.23 ± 0.06%), β-sheet (11.06 ± 0.06%), β-turn (23.67 ± 0.06%) and random coil (32.03 ± 0.06%). The enzyme was characterized with the optimal temperature and pH of 45 °C and 9.0, respectively. Notably, LUCA-DmpA was also characterized with remarkable pH tolerance based on the observation of more than 85% remaining enzymatic activity after incubation at different pH buffers (pH = 6-11) for 12 h. Additionally, rather than being improved or inhibited by metal ions, its enzymatic activity was found to be promoted by introducing organic solvent with a larger log P value. Based on these homology modeling results, the screened LUCA-DmpA is suggested to have further optimization potential, and thereafter to be offered as a promising candidate for real industrial applications.

Consequences of Substitution Model Selection on Protein Ancestral Sequence Reconstruction

The selection of the best-fitting substitution model of molecular evolution is a traditional step for phylogenetic inferences, including ancestral sequence reconstruction (ASR). However, a few recent studies suggested that applying this procedure does not affect the accuracy of phylogenetic tree reconstruction. Here, we revisited this debate topic by analyzing the influence of selection among substitution models of protein evolution, with focus on exchangeability matrices, on the accuracy of ASR using simulated and real data. We found that the selected best-fitting substitution model produces the most accurate ancestral sequences, especially if the data present large genetic diversity. Indeed, ancestral sequences reconstructed under substitution models with similar exchangeability matrices were similar, suggesting that if the selected best-fitting model cannot be used for the reconstruction, applying a model similar to the selected one is preferred. We conclude that selecting among substitution models of protein evolution is recommended for reconstructing accurate ancestral sequences.

Ancestral reconstruction of duplicated signaling proteins reveals the evolution of signaling specificity

Gene duplication is crucial to generating novel signaling pathways during evolution. However, it remains unclear how the redundant proteins produced by gene duplication ultimately acquire new interaction specificities to establish insulated paralogous signaling pathways. Here, we used ancestral sequence reconstruction to resurrect and characterize a bacterial two-component signaling system that duplicated in α-proteobacteria. We determined the interaction specificities of the signaling proteins that existed before and immediately after this duplication event and then identified key mutations responsible for establishing specificity in the two systems. Just three mutations, in only two of the four interacting proteins, were sufficient to establish specificity of the extant systems. Some of these mutations weakened interactions between paralogous systems to limit crosstalk. However, others strengthened interactions within a system, indicating that the ancestral interaction, although functional, had the potential to be strengthened. Our work suggests that protein-protein interactions with such latent potential may be highly amenable to duplication and divergence.

A simple method to determine changes in the affinity between HisF and HisH in the Imidazole Glycerol Phosphate Synthase heterodimer

The bi-enzyme HisF-HisH heterodimer is part of the pathway that produces histidine and purines in bacteria and lower eukaryotes, but it is absent in mammals. This heterodimer has been largely studied probing the basis of the allosteric effects and the structural stability in proteins. It is also a potential target for antibacterial drugs. In this work, we developed a simple method to evaluate changes in the affinity between HisF and HisH in the heterodimer of the bacteria Thermotoga maritima. HisH contains a single tryptophan residue, which is exposed in the free protein, but buried in the heterodimer interface. Hence, the intrinsic fluorescence maximum of this residue changes to shorter wavelengths upon dimerization. Thus, we used the fluorescence intensity at this shorter wavelength to monitor heterodimer accumulation when HisH was combined with sub-stoichiometric HisF. Under conditions where the HisF-HisH heterodimer is in equilibrium with the free states of these enzymes, when [HisH] > [HisF], we deduced a linear function connecting [HisF-HisH] to [HisF], in which the slope depends on the heterodimer dissociation constant (Kd). Based on this equation, taking fluorescence intensities as proxies of the heterodimer and HisF concentrations, we experimentally determined the Kd at four different temperatures. These Kd values were compared to those evaluated using ITC. Both methods revealed an increase in the HisF and HisH binding affinity as the temperature increases. In spite of differences in their absolute values, the Kd determined using these methods presented an evident linear correlation. To demonstrate the effectiveness of the fluorescence method we determined the effect on the Kd caused by 12 single mutations in HisF. Coherently, this test singled out the only mutation in the binding interface. In brief, the method described here effectively probes qualitative effects on the Kd, can be carried out using common laboratory equipment and is scalable.

Engineering of Ancestors as a Tool to Elucidate Structure, Mechanism, and Specificity of Extant Terpene Cyclase

Structural information is crucial for understanding catalytic mechanisms and to guide enzyme engineering efforts of biocatalysts, such as terpene cyclases. However, low sequence similarity can impede homology modeling, and inherent protein instability presents challenges for structural studies. We hypothesized that X-ray crystallography of engineered thermostable ancestral enzymes can enable access to reliable homology models of extant biocatalysts. We have applied this concept in concert with molecular modeling and enzymatic assays to understand the structure activity relationship of spiroviolene synthase, a class I terpene cyclase, aiming to engineer its specificity. Engineering a surface patch in the reconstructed ancestor afforded a template structure for generation of a high-confidence homology model of the extant enzyme. On the basis of structural considerations, we designed and crystallized ancestral variants with single residue exchanges that exhibited tailored substrate specificity and preserved thermostability. We show how the two single amino acid alterations identified in the ancestral scaffold can be transferred to the extant enzyme, conferring a specificity switch that impacts the extant enzyme's specificity for formation of the diterpene spiroviolene over formation of sesquiterpenes hedycaryol and farnesol by up to 25-fold. This study emphasizes the value of ancestral sequence reconstruction combined with enzyme engineering as a versatile tool in chemical biology.

Rosetta design with co-evolutionary information retains protein function

Computational protein design has the ambitious goal of crafting novel proteins that address challenges in biology and medicine. To overcome these challenges, the computational protein modeling suite Rosetta has been tailored to address various protein design tasks. Recently, statistical methods have been developed that identify correlated mutations between residues in a multiple sequence alignment of homologous proteins. These subtle inter-dependencies in the occupancy of residue positions throughout evolution are crucial for protein function, but we found that three current Rosetta design approaches fail to recover these co-evolutionary couplings. Thus, we developed the Rosetta method ResCue (residue-coupling enhanced) that leverages co-evolutionary information to favor sequences which recapitulate correlated mutations, as observed in nature. To assess the protocols via recapitulation designs, we compiled a benchmark of ten proteins each represented by two, structurally diverse states. We could demonstrate that ResCue designed sequences with an average sequence recovery rate of 70%, whereas three other protocols reached not more than 50%, on average. Our approach had higher recovery rates also for functionally important residues, which were studied in detail. This improvement has only a minor negative effect on the fitness of the designed sequences as assessed by Rosetta energy. In conclusion, our findings support the idea that informing protocols with co-evolutionary signals helps to design stable and native-like proteins that are compatible with the different conformational states required for a complex function.

In homologous proteins, functionally or structurally important residues are strongly conserved. Thus, the consideration of conservation signals during protein design protocols can help to create sequences that are more native-like. However, the number of conserved residues is small in many proteins and not all important residues can be captured by conservation analysis. Residues are forming networks whose composition is dictated by protein structure function and thus is visible through the co-evolutionary analysis. Nowadays, advanced methods allow us to deduce these networks from multiple sequence alignments. Thus, we have implemented the novel Rosetta method termed ‘ResCue’ that informs the design protocol with co-evolutionary signals. Recapitulation designs based on ten difficult benchmarks made clear that this protocol creates sequences that are more native-like than three other, state-of-the-art design protocols.

Ancestral sequence reconstruction produces thermally stable enzymes with mesophilic enzyme-like catalytic properties

Enzymes have high catalytic efficiency and low environmental impact, and are therefore potentially useful tools for various industrial processes. Crucially, however, natural enzymes do not always have the properties required for specific processes. It may be necessary, therefore, to design, engineer, and evolve enzymes with properties that are not found in natural enzymes. In particular, the creation of enzymes that are thermally stable and catalytically active at low temperature is desirable for processes involving both high and low temperatures. In the current study, we designed two ancestral sequences of 3-isopropylmalate dehydrogenase by an ancestral sequence reconstruction technique based on a phylogenetic analysis of extant homologous amino acid sequences. Genes encoding the designed sequences were artificially synthesized and expressed in Escherichia coli. The reconstructed enzymes were found to be slightly more thermally stable than the extant thermophilic homologue from Thermus thermophilus. Moreover, they had considerably higher low-temperature catalytic activity as compared with the T. thermophilus enzyme. Detailed analyses of their temperature-dependent specific activities and kinetic properties showed that the reconstructed enzymes have catalytic properties similar to those of mesophilic homologues. Collectively, our study demonstrates that ancestral sequence reconstruction can produce a thermally stable enzyme with catalytic properties adapted to low-temperature reactions.

Analysis of allosteric communication in a multienzyme complex by ancestral sequence reconstruction

Tryptophan synthase (TS) is a heterotetrameric αββα complex. It is characterized by the channeling of the reaction intermediate indole and the mutual activation of the α-subunit TrpA and the β-subunit TrpB via a complex allosteric network. We have analyzed this allosteric network by means of ancestral sequence reconstruction (ASR), which is an in silico method to resurrect extinct ancestors of modern proteins. Previously, the sequences of TrpA and TrpB from the last bacterial common ancestor (LBCA) have been computed by means of ASR and characterized. LBCA-TS is similar to modern TS by forming a αββα complex with indole channeling taking place. However, LBCA-TrpA allosterically decreases the activity of LBCA-TrpB, whereas, for example, the modern ncTrpA from Neptuniibacter caesariensis allosterically increases the activity of ncTrpB. To identify amino acid residues that are responsible for this inversion of the allosteric effect, all 6 evolutionary TrpA and TrpB intermediates that stepwise link LBCA-TS with ncTS were characterized. Remarkably, the switching from TrpB inhibition to TrpB activation by TrpA occurred between 2 successive TS intermediates. Sequence comparison of these 2 intermediates and iterative rounds of site-directed mutagenesis allowed us to identify 4 of 413 residues from TrpB that are crucial for its allosteric activation by TrpA. The effect of our mutational studies was rationalized by a community analysis based on molecular dynamics simulations. Our findings demonstrate that ancestral sequence reconstruction can efficiently identify residues contributing to allosteric signal propagation in multienzyme complexes.

Sequence selection by FitSS4ASR alleviates ancestral sequence reconstruction as exemplified for geranylgeranylglyceryl phosphate synthase

For evolutionary studies, but also for protein engineering, ancestral sequence reconstruction (ASR) has become an indispensable tool. The first step of every ASR protocol is the preparation of a representative sequence set containing at most a few hundred recent homologs whose composition determines decisively the outcome of a reconstruction. A common approach for sequence selection consists of several rounds of manual recompilation that is driven by embedded phylogenetic analyses of the varied sequence sets. For ASR of a geranylgeranylglyceryl phosphate synthase, we additionally utilized FitSS4ASR, which replaces this time-consuming protocol with an efficient and more rational approach. FitSS4ASR applies orthogonal filters to a set of homologs to eliminate outlier sequences and those bearing only a weak phylogenetic signal. To demonstrate the usefulness of FitSS4ASR, we determined experimentally the oligomerization state of eight predecessors, which is a delicate and taxon-specific property. Corresponding ancestors deduced in a manual approach and by means of FitSS4ASR had the same dimeric or hexameric conformation; this concordance testifies to the efficiency of FitSS4ASR for sequence selection. FitSS4ASR-based results of two other ASR experiments were added to the Supporting Information. Program and documentation are available at https://gitlab.bioinf.ur.de/hek61586/FitSS4ASR.

Ancestral Sequence Reconstruction as a Tool for the Elucidation of a Stepwise Evolutionary Adaptation

Ancestral sequence reconstruction (ASR) is a powerful tool to infer primordial sequences from contemporary, i.e., extant ones. An essential element of ASR is the computation of a phylogenetic tree whose leaves are the chosen extant sequences. Most often, the reconstructed sequence related to the root of this tree is of greatest interest: It represents the common ancestor (CA) of the sequences under study. If this sequence encodes a protein, one can "resurrect" the CA by means of gene synthesis technology and study biochemical properties of this extinct predecessor with the help of wet-lab experiments.However, ASR deduces also sequences for all internal nodes of the tree, and the well-considered analysis of these "intermediates" can help to elucidate evolutionary processes. Moreover, one can identify key mutations that alter proteins or protein complexes and are responsible for the differing properties of extant proteins. As an illustrative example, we describe the protocol for the rapid identification of hotspots determining the binding of the two subunits within the heteromeric complex imidazole glycerol phosphate synthase.

Biotechnological and protein-engineering implications of ancestral protein resurrection

Approximations to the sequences of ancestral proteins can be derived from the sequences of their modern descendants. Proteins encoded by such reconstructed sequences can be prepared in the laboratory and subjected to experimental scrutiny. These 'resurrected' ancestral proteins often display remarkable properties, reflecting ancestral adaptations to intra-cellular and extra-cellular environments that differed from the environments hosting modern/extant proteins. Recent experimental and computational work has specifically discussed high stability, substrate and catalytic promiscuity, conformational flexibility/diversity and altered patterns of interaction with other sub-cellular components. In this review, we discuss these remarkable properties as well as recent attempts to explore their biotechnological and protein-engineering potential.

Characterization of Reconstructed Ancestral Proteins Suggests a Change in Temperature of the Ancient Biosphere

Understanding the evolution of ancestral life, and especially the ability of some organisms to flourish in the variable environments experienced in Earth’s early biosphere, requires knowledge of the characteristics and the environment of these ancestral organisms. Information about early life and environmental conditions has been obtained from fossil records and geological surveys. Recent advances in phylogenetic analysis, and an increasing number of protein sequences available in public databases, have made it possible to infer ancestral protein sequences possessed by ancient organisms. However, the in silico studies that assess the ancestral base content of ribosomal RNAs, the frequency of each amino acid in ancestral proteins, and estimate the environmental temperatures of ancient organisms, show conflicting results. The characterization of ancestral proteins reconstructed in vitro suggests that ancient organisms had very thermally stable proteins, and therefore were thermophilic or hyperthermophilic. Experimental data supports the idea that only thermophilic ancestors survived the catastrophic increase in temperature of the biosphere that was likely associated with meteorite impacts during the early history of Earth. In addition, by expanding the timescale and including more ancestral proteins for reconstruction, it appears as though the Earth’s surface temperature gradually decreased over time, from Archean to present.