Ancestral Protein Reconstruction Yields Insights into Adaptive Evolution of Binding Specificity in Solute-Binding Proteins

Rugged fitness landscapes minimize promiscuity in the evolution of transcriptional repressors

How a protein's function influences the shape of its fitness landscape, smooth or rugged, is a fundamental question in evolutionary biochemistry. Smooth landscapes arise when incremental mutational steps lead to a progressive change in function, as commonly seen in enzymes and binding proteins. On the other hand, rugged landscapes are poorly understood because of the inherent unpredictability of how sequence changes affect function. Here, we experimentally characterize the entire sequence phylogeny, comprising 1,158 extant and ancestral sequences, of the DNA-binding domain (DBD) of the LacI/GalR transcriptional repressor family. Our analysis revealed an extremely rugged landscape with rapid switching of specificity, even between adjacent nodes. Further, the ruggedness arises due to the necessity of the repressor to simultaneously evolve specificity for asymmetric operators and disfavors potentially adverse regulatory crosstalk. Our study provides fundamental insight into evolutionary, molecular, and biophysical rules of genetic regulation through the lens of fitness landscapes.

Flexible active-site loops fine-tune substrate specificity of hyperthermophilic metallo-oxidases

Hyperthermophilic ('superheat-loving') archaea found in high-temperature environments such as Pyrobaculum aerophilum contain multicopper oxidases (MCOs) with remarkable efficiency for oxidizing cuprous and ferrous ions. In this work, directed evolution was used to expand the substrate specificity of P. aerophilum McoP for organic substrates. Six rounds of error-prone PCR and DNA shuffling followed by high-throughput screening lead to the identification of a hit variant with a 220-fold increased efficiency (kcat/Km) than the wild-type for 2,2'-azino-bis(3-ethylbenzothiazoline-6-sulphonic acid) (ABTS) without compromising its intrinsic activity for metal ions. The analysis of the X-ray crystal structure reveals four proximal mutations close to the T1Cu active site. One of these mutations is within the 23-residues loop that occludes this site, a distinctive feature of prokaryotic MCOs. The increased flexibility of this loop results in an enlarged tunnel and one additional pocket that facilitates bulky substrate-enzyme interactions. These findings underscore the synergy between mutations that modulate the dynamics of the active-site loop enabling enhanced catalytic function. This study highlights the potential of targeting loops close to the T1Cu for engineering improvements suitable for biotechnological applications.

Engineering the thermal stability of a polyphosphate kinase by ancestral sequence reconstruction to expand the temperature boundary for an industrially applicable ATP regeneration system

ATP-dependent energy-consuming enzymatic reactions are widely used in cell-free biocatalysis. However, the direct addition of large amounts of expensive ATP can greatly increase cost, and enzymatic production is often difficult to achieve as a result. Although a polyphosphate kinase (PPK)-polyphosphate-based ATP regeneration system has the potential to solve this challenge, the generally poor thermal stability of PPKs limits the widespread use of this method. In this paper, we evaluated the thermal stability of a PPK from Sulfurovum lithotrophicum (SlPPK2). After directed evolution and computation-supported design, we found that SlPPK2 is very recalcitrant and cannot acquire beneficial mutations. Inspired by the usually outstanding stability of ancestral enzymes, we reconstructed the ancestral sequence of the PPK family and used it as a guide to construct three heat-stable variants of SlPPK2, of which the L35F/T144S variant has a half-life of more than 14 h at 60°C. Molecular dynamics simulations were performed on all enzymes to analyze the reasons for the increased thermal stability. The results showed that mutations at these two positions act synergistically from the interior and surface of the protein, leading to a more compact structure. Finally, the robustness of the L35F/T144S variant was verified in the synthesis of nucleotides at high temperature. In practice, the use of this high-temperature ATP regeneration system can effectively avoid byproduct accumulation. Our work extends the temperature boundary of ATP regeneration and has great potential for industrial applications.IMPORTANCEATP regeneration is an important basic applied study in the field of cell-free biocatalysis. Polyphosphate kinase (PPK) is an enzyme tool widely used for energy regeneration during enzymatic reactions. However, the thermal stability of the PPKs reported to date that can efficiently regenerate ATP is usually poor, which greatly limits their application. In this study, the thermal stability of a difficult-to-engineer PPK from Sulfurovum lithotrophicum was improved, guided by an ancestral sequence reconstruction strategy. The optimal variant has a 4.5-fold longer half-life at 60°C than the wild-type enzyme, thus enabling the extension of the temperature boundary for ATP regeneration. The ability of this variant to regenerate ATP was well demonstrated during high-temperature enzymatic production of nucleotides.

A journey from molecule to physiology and in silico tools for drug discovery targeting the transient receptor potential vanilloid type 1 (TRPV1) channel

The heat and capsaicin receptor TRPV1 channel is widely expressed in nerve terminals of dorsal root ganglia (DRGs) and trigeminal ganglia innervating the body and face, respectively, as well as in other tissues and organs including central nervous system. The TRPV1 channel is a versatile receptor that detects harmful heat, pain, and various internal and external ligands. Hence, it operates as a polymodal sensory channel. Many pathological conditions including neuroinflammation, cancer, psychiatric disorders, and pathological pain, are linked to the abnormal functioning of the TRPV1 in peripheral tissues. Intense biomedical research is underway to discover compounds that can modulate the channel and provide pain relief. The molecular mechanisms underlying temperature sensing remain largely unknown, although they are closely linked to pain transduction. Prolonged exposure to capsaicin generates analgesia, hence numerous capsaicin analogs have been developed to discover efficient analgesics for pain relief. The emergence of in silico tools offered significant techniques for molecular modeling and machine learning algorithms to indentify druggable sites in the channel and for repositioning of current drugs aimed at TRPV1. Here we recapitulate the physiological and pathophysiological functions of the TRPV1 channel, including structural models obtained through cryo-EM, pharmacological compounds tested on TRPV1, and the in silico tools for drug discovery and repositioning.

Retracing the evolution of a modern periplasmic binding protein

Investigating the evolution of structural features in modern multidomain proteins helps to understand their immense diversity and functional versatility. The class of periplasmic binding proteins (PBPs) offers an opportunity to interrogate one of the main processes driving diversification: the duplication and fusion of protein sequences to generate new architectures. The symmetry of their two-lobed topology, their mechanism of binding, and the organization of their operon structure led to the hypothesis that PBPs arose through a duplication and fusion event of a single common ancestor. To investigate this claim, we set out to reverse the evolutionary process and recreate the structural equivalent of a single-lobed progenitor using ribose-binding protein (RBP) as our model. We found that this modern PBP can be deconstructed into its lobes, producing two proteins that represent possible progenitor halves. The isolated halves of RBP are well folded and monomeric proteins, albeit with a lower thermostability, and do not retain the original binding function. However, the two entities readily form a heterodimer in vitro and in-cell. The x-ray structure of the heterodimer closely resembles the parental protein. Moreover, the binding function is fully regained upon formation of the heterodimer with a ligand affinity similar to that observed in the modern RBP. This highlights how a duplication event could have given rise to a stable and functional PBP-like fold and provides insights into how more complex functional structures can evolve from simpler molecular components.

Computational remodeling of an enzyme conformational landscape for altered substrate selectivity

Structural plasticity of enzymes dictates their function. Yet, our ability to rationally remodel enzyme conformational landscapes to tailor catalytic properties remains limited. Here, we report a computational procedure for tuning conformational landscapes that is based on multistate design of hinge-mediated domain motions. Using this method, we redesign the conformational landscape of a natural aminotransferase to preferentially stabilize a less populated but reactive conformation and thereby increase catalytic efficiency with a non-native substrate, resulting in altered substrate selectivity. Steady-state kinetics of designed variants reveals activity increases with the non-native substrate of approximately 100-fold and selectivity switches of up to 1900-fold. Structural analyses by room-temperature X-ray crystallography and multitemperature nuclear magnetic resonance spectroscopy confirm that conformational equilibria favor the target conformation. Our computational approach opens the door to targeted alterations of conformational states and equilibria, which should facilitate the design of biocatalysts with customized activity and selectivity.

Identification and evaluation of potential inhibitor molecules against TcyA from Candidatus Liberibacter asiaticus

Of the two putative amino acid binding periplasmic receptors of ABC transporter family in Candidatus Liberibacter asiaticus (CLas), cystine binding receptor (CLasTcyA) has been shown to mainly express in phloem of citrus plant and is a target for inhibitor development. The crystal structure of CLasTcyA in complex with substrates has been reported earlier. The present work reports the identification and evaluation of potential candidates for their inhibitory potential against CLasTcyA. Among many compounds, selected through virtual screening, and MD simulation, pimozide, clidinium, sulfasalazine and folic acid showed significantly higher affinities and stability in complex with CLasTcyA. The SPR studies with CLasTcyA revealed significantly higher binding affinities for pimozide and clidinium (Kd, 2.73 nM and 70 nM, respectively) as compared to cystine (Kd, 1.26 µM). The higher binding affinities could be attributed to significantly increased number of interactions in the binding pocket as evident from the crystal structures of CLasTcyA in complex with pimozide and clidinium as compared to cystine. The CLasTcyA possess relatively large binding pocket where bulkier inhibitors fit quite well. In planta studies, carried out to assess the effect of inhibitors on HLB infected Mosambi plants, showed significant reduction in CLas titre in plants treated with inhibitors as compared to control plants. The results showed that pimozide exhibited higher efficiency as compared to clidinium in reducing CLas titre in treated plants. Our results showed that the inhibitor development against critical proteins like CLasTcyA can be an important strategy in management of HLB.

Discovery of novel amino acid production traits by evolution of synthetic co-cultures

BACKGROUND

Amino acid production features of Corynebacterium glutamicum were extensively studied in the last two decades. Many metabolic pathways, regulatory and transport principles are known, but purely rational approaches often provide only limited progress in production optimization. We recently generated stable synthetic co-cultures, termed Communities of Niche-optimized Strains (CoNoS), that rely on cross-feeding of amino acids for growth. This setup has the potential to evolve strains with improved production by selection of faster growing communities.

RESULTS

Here we performed adaptive laboratory evolution (ALE) with a CoNoS to identify mutations that are relevant for amino acid production both in mono- and co-cultures. During ALE with the CoNoS composed of strains auxotrophic for either L-leucine or L-arginine, we obtained a 23% growth rate increase. Via whole-genome sequencing and reverse engineering, we identified several mutations involved in amino acid transport that are beneficial for CoNoS growth. The L-leucine auxotrophic strain carried an expression-promoting mutation in the promoter region of brnQ (cg2537), encoding a branched-chain amino acid transporter in combination with mutations in the genes for the Na⁺/H⁺-antiporter Mrp1 (cg0326-cg0321). This suggested an unexpected link of Mrp1 to L-leucine transport. The L-arginine auxotrophic partner evolved expression-promoting mutations near the transcriptional start site of the yet uncharacterized operon argTUV (cg1504-02). By mutation studies and ITC, we characterized ArgTUV as the only L-arginine uptake system of C. glutamicum with an affinity of K_D = 30 nM. Finally, deletion of argTUV in an L-arginine producer strain resulted in a faster and 24% higher L-arginine production in comparison to the parental strain.

CONCLUSION

Our work demonstrates the power of the CoNoS-approach for evolution-guided identification of non-obvious production traits, which can also advance amino acid production in monocultures. Further rounds of evolution with import-optimized strains can potentially reveal beneficial mutations also in metabolic pathway enzymes. The approach can easily be extended to all kinds of metabolite cross-feeding pairings of different organisms or different strains of the same organism, thereby enabling the identification of relevant transport systems and other favorable mutations.

Ancestral Sequence Reconstruction Identifies Structural Changes Underlying the Evolution of Ideonella sakaiensis PETase and Variants with Improved Stability and Activity

The improved production, recycling, and removal of plastic waste, such as polyethylene terephthalate (PET), are pressing environmental and economic issues for society. Biocatalytic (enzymatic) PET depolymerization is potentially a sustainable, low-energy solution to PET recycling, especially when compared with current disposal methods such as landfills, incineration, or gasification. IsPETase has been extensively studied for its use in PET depolymerization; however, its evolution from cutinases is not fully understood, and most engineering studies have neglected the majority of the available sequence space remote from the active site. In this study, ancestral protein reconstruction (ASR) has been used to trace the evolutionary trajectory from ancient serine hydrolases to IsPETase, while ASR and the related design approach, protein repair one-stop shop, were used to identify enzyme variants with improved activity and stability. Kinetic and structural characterization of these variants reveals new insights into the evolution of PETase activity and the role of second-shell mutations around the active site. Among the designed and reconstructed variants, we identified several with melting points 20 °C higher than that of IsPETase and two variants with significantly higher catalytic activity.

The present and the future of protein biosensor engineering

Protein biosensors play increasingly important roles in cell and neurobiology and have the potential to revolutionise the way clinical and industrial analytics are performed. The gradual transition from multicomponent biosensors to fully integrated single chain allosteric biosensors has brought the field closer to commercial applications. We evaluate various approaches for converting constitutively active protein reporter domains into analyte operated switches. We discuss the paucity of the natural receptors that undergo conformational changes sufficiently large to control the activity of allosteric reporter domains. This problem can be overcome by constructing artificial versions of such receptors. The design path to such receptors involves the construction of Chemically Induced Dimerisation systems (CIDs) that can be configured to operate single and two-component biosensors.

Engineering and screening of novel β-1,3-xylanases with desired hydrolysate type by optimized ancestor sequence reconstruction and data mining

•

A novel integrative strategy for engineering β-1,3-xylanases with desired products.

•

AncXyl10 is the first successful example of ASR to shift the hydrolysate types.

•

The hydrolysates of AncXyl10 was only β-1,3-xylobiose and β-1,3-xylotriose.

•

The underlying mechanism laid a new groundwork towards hydrolase engineering.

Engineering of hydrolases to shift their hydrolysate types has not been attempted so far, though computer-assisted enzyme design has been successful. A novel integrative strategy for engineering and screening the β-1,3-xylanase with desired hydrolysate types was proposed, with the purpose to solve problems that the separation and preparation of β-1,3-xylo-oligosaccharides was in high cost yet in low yield as monosaccharides existed in the hydrolysates. By classifying the hydrolysate types and coding them into numerical values, two robust mathematical models with five selected attributes from molecular docking were established based on LogitBoost and partial least squares regression with overall accuracy of 83.3% and 100%, respectively. Then, they were adopted for efficient screening the potential mutagenesis library of β-1,3-xylanases that only product oligosaccharides. The virtually designed AncXyl10 was selected and experimentally verified to produce only β-1,3-xylobiose (60.38%) and β-1,3-xylotriose (39.62%), which facilitated the preparation of oligosaccharides with high purity. The underlying mechanism of AncXyl10 may associated with the gap processing and ancestral amino acid substitution in the process of ancestral sequence reconstruction. Since many carbohydrate-active enzymes have highly conserved active sites, the strategy and their biomolecular basis will shield a new light for engineering carbohydrates hydrolase to produce specific oligosaccharides.

Distal Mutations Shape Substrate-Binding Sites during Evolution of a Metallo-Oxidase into a Laccase

Laccases are in increasing demand as innovative solutions in the biorefinery fields. Here, we combine mutagenesis with structural, kinetic, and in silico analyses to characterize the molecular features that cause the evolution of a hyperthermostable metallo-oxidase from the multicopper oxidase family into a laccase (k _cat 273 s^-1 for a bulky aromatic substrate). We show that six mutations scattered across the enzyme collectively modulate dynamics to improve the binding and catalysis of a bulky aromatic substrate. The replacement of residues during the early stages of evolution is a stepping stone for altering the shape and size of substrate-binding sites. Binding sites are then fine-tuned through high-order epistasis interactions by inserting distal mutations during later stages of evolution. Allosterically coupled, long-range dynamic networks favor catalytically competent conformational states that are more suitable for recognizing and stabilizing the aromatic substrate. This work provides mechanistic insight into enzymatic and evolutionary molecular mechanisms and spots the importance of iterative experimental and computational analyses to understand local-to-global changes.

Structural insights into choline-O-sulfatase reveal the molecular determinants for ligand binding

Choline-O-sulfatase (COSe; EC 3.1.6.6) is a member of the alkaline phosphatase (AP) superfamily, and its natural function is to hydrolyze choline-O-sulfate into choline and sulfate. Despite its natural function, the major interest in this enzyme resides in the landmark catalytic/substrate promiscuity of sulfatases, which has led to attention in the biotechnological field due to their potential in protein engineering. In this work, an in-depth structural analysis of wild-type Sinorhizobium (Ensifer) meliloti COSe (SmeCOSe) and its C54S active-site mutant is reported. The binding mode of this AP superfamily member to both products of the reaction (sulfate and choline) and to a substrate-like compound are shown for the first time. The structures further confirm the importance of the C-terminal extension of the enzyme in becoming part of the active site and participating in enzyme activity through dynamic intra-subunit and inter-subunit hydrogen bonds (Asn146A-Asp500B-Asn498B). These residues act as the `gatekeeper' responsible for the open/closed conformations of the enzyme, in addition to assisting in ligand binding through the rearrangement of Leu499 (with a movement of approximately 5 Å). Trp129 and His145 clamp the quaternary ammonium moiety of choline and also connect the catalytic cleft to the C-terminus of an adjacent protomer. The structural information reported here contrasts with the proposed role of conformational dynamics in promoting the enzymatic catalytic proficiency of an enzyme.

Substrate Dynamics Contribute to Enzymatic Specificity in Human and Bacterial Methionine Adenosyltransferases

Protein conformational changes can facilitate the binding of noncognate substrates and underlying promiscuous activities. However, the contribution of substrate conformational dynamics to this process is comparatively poorly understood. Here, we analyze human (hMAT2A) and Escherichia coli (eMAT) methionine adenosyltransferases that have identical active sites but different substrate specificity. In the promiscuous hMAT2A, noncognate substrates bind in a stable conformation to allow catalysis. In contrast, noncognate substrates sample stable productive binding modes less frequently in eMAT owing to altered mobility in the enzyme active site. Different cellular concentrations of substrates likely drove the evolutionary divergence of substrate specificity in these orthologues. The observation of catalytic promiscuity in hMAT2A led to the detection of a new human metabolite, methyl thioguanosine, that is produced at elevated levels in a cancer cell line. This work establishes that identical active sites can result in different substrate specificity owing to the effects of substrate and enzyme dynamics.

Functional characterisation of substrate-binding proteins to address nutrient uptake in marine picocyanobacteria

Marine cyanobacteria are key primary producers, contributing significantly to the microbial food web and biogeochemical cycles by releasing and importing many essential nutrients cycled through the environment. A subgroup of these, the picocyanobacteria (Synechococcus and Prochlorococcus), have colonised almost all marine ecosystems, covering a range of distinct light and temperature conditions, and nutrient profiles. The intra-clade diversities displayed by this monophyletic branch of cyanobacteria is indicative of their success across a broad range of environments. Part of this diversity is due to nutrient acquisition mechanisms, such as the use of high-affinity ATP-binding cassette (ABC) transporters to competitively acquire nutrients, particularly in oligotrophic (nutrient scarce) marine environments. The specificity of nutrient uptake in ABC transporters is primarily determined by the peripheral substrate-binding protein (SBP), a receptor protein that mediates ligand recognition and initiates translocation into the cell. The recent availability of large numbers of sequenced picocyanobacterial genomes indicates both Synechococcus and Prochlorococcus apportion >50% of their transport capacity to ABC transport systems. However, the low degree of sequence homology among the SBP family limits the reliability of functional assignments using sequence annotation and prediction tools. This review highlights the use of known SBP structural representatives for the uptake of key nutrient classes by cyanobacteria to compare with predicted SBP functionalities within sequenced marine picocyanobacteria genomes. This review shows the broad range of conserved biochemical functions of picocyanobacteria and the range of novel and hypothetical ABC transport systems that require further functional characterisation.

Structural dynamics in the evolution of a bilobed protein scaffold

Proteins conduct numerous complex biological functions by use of tailored structural dynamics. The molecular details of how these emerged from ancestral peptides remains mysterious. How does nature utilize the same repertoire of folds to diversify function? To shed light on this, we analyzed bilobed proteins with a common structural core, which is spread throughout the tree of life and is involved in diverse biological functions such as transcription, enzymatic catalysis, membrane transport, and signaling. We show here that the structural dynamics of the structural core differentiate predominantly via terminal additions during a long-period evolution. This diversifies substrate specificity and, ultimately, biological function.

Novel biophysical tools allow the structural dynamics of proteins and the regulation of such dynamics by binding partners to be explored in unprecedented detail. Although this has provided critical insights into protein function, the means by which structural dynamics direct protein evolution remain poorly understood. Here, we investigated how proteins with a bilobed structure, composed of two related domains from the periplasmic-binding protein–like II domain family, have undergone divergent evolution, leading to adaptation of their structural dynamics. We performed a structural analysis on ∼600 bilobed proteins with a common primordial structural core, which we complemented with biophysical studies to explore the structural dynamics of selected examples by single-molecule Förster resonance energy transfer and Hydrogen–Deuterium exchange mass spectrometry. We show that evolutionary modifications of the structural core, largely at its termini, enable distinct structural dynamics, allowing the diversification of these proteins into transcription factors, enzymes, and extracytoplasmic transport-related proteins. Structural embellishments of the core created interdomain interactions that stabilized structural states, reshaping the active site geometry, and ultimately altered substrate specificity. Our findings reveal an as-yet-unrecognized mechanism for the emergence of functional promiscuity during long periods of evolution and are applicable to a large number of domain architectures.

Through-Space Scalar 19F-19F Couplings between Fluorinated Noncanonical Amino Acids for the Detection of Specific Contacts in Proteins

Fluorine atoms are known to display scalar ¹⁹F-¹⁹F couplings in nuclear magnetic resonance (NMR) spectra when they are sufficiently close in space for nonbonding orbitals to overlap. We show that fluorinated noncanonical amino acids positioned in the hydrophobic core or on the surface of a protein can be linked by scalar through-space ¹⁹F-¹⁹F (^TSJ_FF) couplings even if the ¹⁹F spins are in the time average separated by more than the van der Waals distance. Using two different aromatic amino acids featuring CF₃ groups, O-trifluoromethyl-tyrosine and 4-trifluoromethyl-phenylalanine, we show that ¹⁹F-¹⁹F TOCSY experiments are sufficiently sensitive to detect ^TSJ_FF couplings between 2.5 and 5 Hz in the 19 kDa protein PpiB measured on a two-channel 400 MHz NMR spectrometer with a regular room temperature probe. A quantitative J evolution experiment enables the measurement of ^TSJ_FF coupling constants that are up to five times smaller than the ¹⁹F NMR line width. In addition, a new aminoacyl-tRNA synthetase was identified for genetic encoding of N⁶-(trifluoroacetyl)-l-lysine (TFA-Lys) and ¹⁹F-¹⁹F TOCSY peaks were observed between two TFA-Lys residues incorporated into the proteins AncCDT-1 and mRFP despite high solvent exposure and flexibility of the TFA-Lys side chains. With the ready availability of systems for site-specific incorporation of fluorinated amino acids into proteins by genetic encoding, ¹⁹F-¹⁹F interactions offer a straightforward way to probe the spatial proximity of selected sites without any assignments of ¹H NMR resonances.

Sensory perception in bacterial cyclic diguanylate signal transduction

Cyclic diguanylate (c-di-GMP) signal transduction systems provide bacteria with the ability to sense changing cell status or environmental conditions and then execute suitable physiological and social behaviors in response. In this review, we provide a comprehensive census of the stimuli and receptors that are linked to the modulation of intracellular c-di-GMP. Emerging evidence indicates that c-di-GMP networks sense light, surfaces, energy, redox potential, respiratory electron acceptors, temperature, and structurally diverse biotic and abiotic chemicals. Bioinformatic analysis of sensory domains in diguanylate cyclases and c-di-GMP-specific phosphodiesterases as well as the receptor complexes associated with them reveals that these functions are linked to a diverse repertoire of protein domain families. We describe the principles of stimulus perception learned from studying these modular sensory devices, illustrate how they are assembled in varied combinations with output domains, and summarize a system for classifying these sensor proteins based on their complexity. Biological information processing via c-di-GMP signal transduction not only is fundamental to bacterial survival in dynamic environments but also is being used to engineer gene expression circuitry and synthetic proteins with à la carte biochemical functionalities.

Catalytic mechanism of ancestral L-lysine oxidase assigned by sequence data mining

A large number of protein sequences are registered in public databases such as PubMed. Functionally uncharacterized enzymes are included in these databases, some of which likely have potential for industrial applications. However, assignment of the enzymes remained difficult tasks for now. In this study, we assigned a total of 28 original sequences to uncharacterized enzymes in the FAD-dependent oxidase family expressed in some species of bacteria including Chryseobacterium, Flavobacterium, and Pedobactor. Progenitor sequence of the assigned 28 sequences was generated by ancestral sequence reconstruction, and the generated sequence exhibited L-lysine oxidase activity; thus, we named the enzyme AncLLysO. Crystal structures of ligand-free and ligand-bound forms of AncLLysO were determined, indicating that the enzyme recognizes L-Lys by hydrogen bond formation with R76 and E383. The binding of L-Lys to AncLLysO induced dynamic structural change at a plug loop formed by residues 251 to 254. Biochemical assays of AncLLysO variants revealed the functional importance of these substrate recognition residues and the plug loop. R76A and E383D variants were also observed to lose their activity, and the kcat/Km value of G251P and Y253A mutations were approximately 800- to 1800-fold lower than that of AncLLysO, despite the indirect interaction of the substrates with the mutated residues. Taken together, our data demonstrate that combinational approaches to sequence classification from database and ancestral sequence reconstruction may be effective not only to find new enzymes using databases of unknown sequences but also to elucidate their functions.

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.

Ancestral sequence reconstruction for protein engineers

In addition to its value in the study of molecular evolution, ancestral sequence reconstruction (ASR) has emerged as a useful methodology for engineering proteins with enhanced properties. Proteins generated by ASR often exhibit unique or improved activity, stability, and/or promiscuity, all of which are properties that are valued by protein engineers. Comparison between extant proteins and evolutionary intermediates generated by ASR also allows protein engineers to identify substitutions that have contributed to functional innovation or diversification within protein families. As ASR becomes more widely adopted as a protein engineering approach, it is important to understand the applications, limitations, and recent developments of this technique. This review highlights recent exemplifications of ASR, as well as technical aspects of the reconstruction process that are relevant to protein engineering.

Ancestral sequence reconstruction - An underused approach to understand the evolution of gene function in plants?

Whilst substantial research effort has been placed on understanding the interactions of plant proteins with their molecular partners, relatively few studies in plants - by contrast to work in other organisms - address how these interactions evolve. It is thought that ancestral proteins were more promiscuous than modern proteins and that specificity often evolved following gene duplication and subsequent functional refining. However, ancestral protein resurrection studies have found that some modern proteins have evolved de novo from ancestors lacking those functions. Intriguingly, the new interactions evolved as a consequence of just a few mutations and, as such, acquisition of new functions appears to be neither difficult nor rare, however, only a few of them are incorporated into biological processes before they are lost to subsequent mutations. Here, we detail the approach of ancestral sequence reconstruction (ASR), providing a primer to reconstruct the sequence of an ancestral gene. We will present case studies from a range of different eukaryotes before discussing the few instances where ancestral reconstructions have been used in plants. As ASR is used to dig into the remote evolutionary past, we will also present some alternative genetic approaches to investigate molecular evolution on shorter timescales. We argue that the study of plant secondary metabolism is particularly well suited for ancestral reconstruction studies. Indeed, its ancient evolutionary roots and highly diverse landscape provide an ideal context in which to address the focal issue around the emergence of evolutionary novelties and how this affects the chemical diversification of plant metabolism.

Bioinformatics analysis and biochemical characterisation of ABC transporter-associated periplasmic substrate-binding proteins ModA and MetQ from Helicobacter pylori strain SS1

The human gastric pathogen Helicobacter pylori relies on the uptake of host-provided nutrients for its proliferation and pathogenicity. ABC transporters that mediate import of small molecules into the cytoplasm of H. pylori employ their cognate periplasmic substrate-binding proteins (SBPs) for ligand capture in the periplasm. The genome of the mouse-adapted strain SS1 of H. pylori encodes eight ABC transporter-associated SBPs, but little is known about their specificity or structure. In this study, we demonstrated that the SBP annotated as ModA binds molybdate (MoO₄^2-, K_D = 3.8 nM) and tungstate (WO₄^2-, K_D = 7.8 nM). In addition, we showed that MetQ binds D-methionine (K_D = 9.5 μM), but not L-methionine, which suggests the existence of as yet unknown pathway for L-methionine uptake. Homology modelling has led to identification of the ligand-binding residues.

A comprehensive binding study illustrates ligand recognition in the periplasmic binding protein PotF

Periplasmic binding proteins (PBPs) are ubiquitous receptors in gram-negative bacteria. They sense solutes and play key roles in nutrient uptake. Escherichia coli's putrescine receptor PotF has been reported to bind putrescine and spermidine. We reveal that several similar biogenic polyamines are recognized by PotF. Using isothermal titration calorimetry paired with X-ray crystallography of the different complexes, we unveil PotF's binding modes in detail. The binding site for PBPs is located between two lobes that undergo a large conformational change upon ligand recognition. Hence, analyzing the influence of ligands on complex formation is crucial. Therefore, we solved crystal structures of an open and closed apo state and used them as a basis for molecular dynamics simulations. In addition, we accessed structural behavior in solution for all complexes by ¹H-¹⁵N HSQC NMR spectroscopy. This combined analysis provides a robust framework for understanding ligand binding for future developments in drug design and protein engineering.

Altered conformational sampling along an evolutionary trajectory changes the catalytic activity of an enzyme

Several enzymes are known to have evolved from non-catalytic proteins such as solute-binding proteins (SBPs). Although attention has been focused on how a binding site can evolve to become catalytic, an equally important question is: how do the structural dynamics of a binding protein change as it becomes an efficient enzyme? Here we performed a variety of experiments, including propargyl-DO3A-Gd(III) tagging and double electron-electron resonance (DEER) to study the rigid body protein dynamics of reconstructed evolutionary intermediates to determine how the conformational sampling of a protein changes along an evolutionary trajectory linking an arginine SBP to a cyclohexadienyl dehydratase (CDT). We observed that primitive dehydratases predominantly populate catalytically unproductive conformations that are vestiges of their ancestral SBP function. Non-productive conformational states, including a wide-open state, are frozen out of the conformational landscape via remote mutations, eventually leading to extant CDT that exclusively samples catalytically relevant compact states. These results show that remote mutations can reshape the global conformational landscape of an enzyme as a mechanism for increasing catalytic activity.

Beta class amino methyltransferases from bacteria to humans: evolution and structural consequences

S-adenosyl-l-methionine dependent methyltransferases catalyze methyl transfers onto a wide variety of target molecules, including DNA and RNA. We discuss a family of methyltransferases, those that act on the amino groups of adenine or cytosine in DNA, have conserved motifs in a particular order in their amino acid sequence, and are referred to as class beta MTases. Members of this class include M.EcoGII and M.EcoP15I from Escherichia coli, Caulobacter crescentus cell cycle-regulated DNA methyltransferase (CcrM), the MTA1-MTA9 complex from the ciliate Oxytricha, and the mammalian MettL3-MettL14 complex. These methyltransferases all generate N6-methyladenine in DNA, with some members having activity on single-stranded DNA as well as RNA. The beta class of methyltransferases has a unique multimeric feature, forming either homo- or hetero-dimers, allowing the enzyme to use division of labor between two subunits in terms of substrate recognition and methylation. We suggest that M.EcoGII may represent an ancestral form of these enzymes, as its activity is independent of the nucleic acid type (RNA or DNA), its strandedness (single or double), and its sequence (aside from the target adenine).

Structure-based screening of binding affinities via small-angle X-ray scattering

Protein-protein and protein-ligand interactions often involve conformational changes or structural rearrangements that can be quantified by solution small-angle X-ray scattering (SAXS). These scattering intensity measurements reveal structural details of the bound complex, the number of species involved and, additionally, the strength of interactions if carried out as a titration. Although a core part of structural biology workflows, SAXS-based titrations are not commonly used in drug discovery contexts. This is because prior knowledge of expected sample requirements, throughput and prediction accuracy is needed to develop reliable ligand screens. This study presents the use of the histidine-binding protein (26 kDa) and other periplasmic binding proteins to benchmark ligand screen performance. Sample concentrations and exposure times were varied across multiple screening trials at four beamlines to investigate the accuracy and precision of affinity prediction. The volatility ratio between titrated scattering curves and a common apo reference is found to most reliably capture the extent of structural and population changes. This obviates the need to explicitly model scattering intensities of bound complexes, which can be strongly ligand-dependent. Where the dissociation constant is within 102 of the protein concentration and the total exposure times exceed 20 s, the titration protocol presented at 0.5 mg ml-1 yields affinities comparable to isothermal titration calorimetry measurements. Estimated throughput ranges between 20 and 100 ligand titrations per day at current synchrotron beamlines, with the limiting step imposed by sample handling and cleaning procedures.

Evolutionary model of protein secondary structure capable of revealing new biological relationships

Ancestral sequence reconstruction has had recent success in decoding the origins and the determinants of complex protein functions. However, phylogenetic analyses of remote homologues must handle extreme amino acid sequence diversity resulting from extended periods of evolutionary change. We exploited the wealth of protein structures to develop an evolutionary model based on protein secondary structure. The approach follows the differences between discrete secondary structure states observed in modern proteins and those hypothesized in their immediate ancestors. We implemented maximum likelihood-based phylogenetic inference to reconstruct ancestral secondary structure. The predictive accuracy from the use of the evolutionary model surpasses that of comparative modeling and sequence-based prediction; the reconstruction extracts information not available from modern structures or the ancestral sequences alone. Based on a phylogenetic analysis of a sequence-diverse protein family, we showed that the model can highlight relationships that are evolutionarily rooted in structure and not evident in amino acid-based analysis.

The physical basis and practical consequences of biological promiscuity

Proteins interact with metabolites, nucleic acids, and other proteins to orchestrate the myriad catalytic, structural and regulatory functions that support life from the simplest microbes to the most complex multicellular organisms. These molecular interactions are often exquisitely specific, but never perfectly so. Adventitious "promiscuous" interactions are ubiquitous due to the thousands of macromolecules and small molecules crowded together in cells. Such interactions may perturb protein function at the molecular level, but as long as they do not compromise organismal fitness, they will not be removed by natural selection. Although promiscuous interactions are physiologically irrelevant, they are important because they can provide a vast reservoir of potential functions that can provide the starting point for evolution of new functions, both in nature and in the laboratory.

Water in Ras Superfamily Evolution

The Ras GTPase superfamily of proteins coordinates a diverse set of cellular outcomes, including cell morphology, vesicle transport, and cell proliferation. Primary amino acid sequence analysis has identified Specificity determinant positions (SDPs) that drive diversified functions specific to the Ras, Rho, Rab, and Arf subfamilies (Rojas et al. 2012, J Cell Biol 196:189-201). The inclusion of water molecules in structural and functional adaptation is likely to be a major response to the selection pressures that drive evolution, yet hydration patterns are not included in phylogenetic analysis. This article shows that conserved crystallographic water molecules coevolved with SDP residues in the differentiation of proteins within the Ras superfamily of small GTPases. The patterns of water conservation between protein subfamilies parallel those of sequence-based evolutionary trees. Thus, hydration patterns have the potential to help elucidate functional significance in the evolution of amino acid residues observed in phylogenetic analysis of homologous proteins. © 2019 Wiley Periodicals, Inc.

Strategies for Increasing Protein Stability

The stability of wild-type proteins is often a hurdle to their practical use in research, industry, and medicine. The route to engineering stability of a protein of interest lies largely with the available data. Where high-resolution structural data is available, rational design, based on fundamental principles of protein chemistry, can improve protein stability. Recent advances in computational biology and the use of nonnatural amino acids have also provided novel rational methods for improving protein stability. Likewise, the explosion of sequence and structural data available in public databases, in combination with improvements in freely available computational tools, has produced accessible phylogenetic approaches. Trawling modern sequence databases can identify the thermostable homologs of a target protein, and evolutionary data can be quickly generated using available phylogenetic tools. Grafting features from those thermostable homologs or ancestors provides stability improvement through a semi-rational approach. Further, molecular techniques such as directed evolution have shown great promise in delivering designer proteins. These strategies are well documented and newly accessible to the molecular biologist, allowing for rapid enhancements of protein stability.

Origin of conformational dynamics in a globular protein

Protein structures are dynamic, undergoing motions that can play a vital role in function. However, the link between primary sequence and conformational dynamics remains poorly understood. Here, we studied how conformational dynamics can arise in a globular protein by evaluating the impact of individual core-residue substitutions in DANCER-3, a streptococcal protein G domain β1 variant that we previously designed to undergo a specific mode of conformational exchange that has never been observed in the wild-type protein. Using a combination of solution NMR experiments and molecular dynamics simulations, we demonstrate that only two mutations are necessary to create this conformational exchange, and that these mutations work synergistically, with one destabilizing the native structure and the other allowing two new conformational states to be accessed on the energy landscape. Overall, our results show how dynamics can appear in a stable globular fold, a critical step in the molecular evolution of dynamics-linked functions.

Highly thermostable carboxylic acid reductases generated by ancestral sequence reconstruction

Carboxylic acid reductases (CARs) are biocatalysts of industrial importance. Their properties, especially their poor stability, render them sub-optimal for use in a bioindustrial pipeline. Here, we employed ancestral sequence reconstruction (ASR) - a burgeoning engineering tool that can identify stabilizing but enzymatically neutral mutations throughout a protein. We used a three-algorithm approach to reconstruct functional ancestors of the Mycobacterial and Nocardial CAR1 orthologues. Ancestral CARs (AncCARs) were confirmed to be CAR enzymes with a preference for aromatic carboxylic acids. Ancestors also showed varied tolerances to solvents, pH and in vivo-like salt concentrations. Compared to well-studied extant CARs, AncCARs had a Tm up to 35 °C higher, with half-lives up to nine times longer than the greatest previously observed. Using ancestral reconstruction we have expanded the existing CAR toolbox with three new thermostable CAR enzymes, providing access to the high temperature biosynthesis of aldehydes to drive new applications in biocatalysis.

Towards descriptor of elementary functions for protein design

We review studies of the protein evolution that help to formulate rules for protein design. Acknowledging the fundamental importance of Dayhoff's provision on the emergence of functional proteins from short peptides, we discuss multiple evidences of the omnipresent partitioning of protein globules into structural/functional units, using which greatly facilitates the engineering and design efforts. Closed loops and elementary functional loops, which are descendants of ancient ring-like peptides that formed fist protein domains in agreement with Dayhoff's hypothesis, can be considered as basic units of protein structure and function. We argue that future developments in protein design approaches should consider descriptors of the elementary functions, which will help to complement designed scaffolds with functional signatures and flexibility necessary for their functions.

Following the Evolutionary Track of a Highly Specific l-Arginine Oxidase by Reconstruction and Biochemical Analysis of Ancestral and Native Enzymes

Following the evolutionary track of enzymes can help elucidate how enzymes attain their characteristic functions, such as thermal adaptation and substrate selectivity, during the evolutionary process. Ancestral sequence reconstruction (ASR) is effective for following evolutionary processes if sufficient sequence data are available. Selecting sequences from the data to generate a curated sequence library is necessary for the successful design of artificial proteins by ASR. In this study, we tried to follow the evolutionary track of l-arginine oxidase (AROD), a flavin adenine dinucleotide (FAD)-dependent amino acid oxidase (LAAO) that exhibits high specificity for l-arginine. The library was generated by selecting sequences in which the 15th, 50th, 332nd, and 580th residues are Gly, Ser, Trp, and Thr, respectively. We excluded sequences that are either extremely short or long and those with a low degree of sequence identity. Three ancestral ARODs (AncARODn0, AncARODn1, and AncARODn2) were designed using the library. Subsequently, we expressed the ancestral ARODs as well as native Oceanobacter kriegii AROD (OkAROD) in bacteria. AncARODn0 is phylogenetically most remote from OkAROD, whereas AncARODn2 is most similar to OkAROD. Thermal stability was gradually increased by extending AROD sequences back to the progenitor, while the temperature at which the residual activity is half of the maximum measured activity (T _1/2) of AncARODn0 was >20°C higher than that of OkAROD. Remarkably, only AncARODn0 exhibited broad substrate selectivity similar to that of conventional promiscuous LAAO. Taken together, our findings led us to infer that AROD may have evolved from a highly thermostable and promiscuous LAAO.IMPORTANCE In this study, we attempted to infer the molecular evolution of a recently isolated FAD-dependent l-arginine oxidase (AROD) that oxidizes l-arginine to 2-ketoarginine. Utilizing 10 candidate AROD sequences, we obtained a total of three ancestral ARODs. In addition, one native AROD was obtained by cloning one of the candidate ARODs. The candidate sequences were selected utilizing a curation method defined in this study. All the ARODs were successfully expressed in Escherichia coli for analysis of their biochemical functions. The catalytic activity of our bacterially expressed ancestral ARODs suggests that our ASR was successful. The ancestral AROD that is phylogenetically most remote from a native AROD has the highest thermal stability and substrate promiscuity. Our findings led us to infer that AROD evolved from a highly thermostable and promiscuous LAAO. As an application, we can design artificial ARODs with improved functions compared with those of native ones.

Protein engineers turned evolutionists-the quest for the optimal starting point

The advent of laboratory directed evolution yielded a fruitful crosstalk between the disciplines of molecular evolution and bio-engineering. Here, we outline recent developments in both disciplines with respect to how one can identify the best starting points for directed evolution, such that highly efficient and robust tailor-made enzymes can be obtained with minimal optimization. Directed evolution studies have highlighted essential features of engineer-able enzymes: highly stable, mutationally robust enzymes with the capacity to accept a broad range of substrates. Robust, evolvable enzymes can be inferred from the natural sequence record. Broad substrate spectrum relates to conformational plasticity and can also be predicted by phylogenetic analyses and/or by computational design. Overall, an increasingly powerful toolkit is becoming available for identifying optimal starting points including network analyses of enzyme superfamilies and other bioinformatics methods.

High-Throughput Reconstruction of Ancestral Protein Sequence, Structure, and Molecular Function

Ancestral protein sequence reconstruction is a powerful technique for explicitly testing hypotheses about the evolution of molecular function, allowing researchers to meticulously dissect how historical changes in protein sequence impacted functional repertoire by altering the protein's 3D structure. These techniques have provided concrete, experimentally validated insights into ancient evolutionary processes and help illuminate the complex relationship between protein sequence, structure, and function. Inferring the protein family phylogenies on which ancestral sequence reconstruction depends and reconstructing the sequences, themselves, are amenable to high-throughput computational analysis. However, determining the structures of ancestral-reconstructed proteins and characterizing their functions typically rely on time-consuming and expensive laboratory analyses, limiting most current studies to examining a relatively small number of specific hypotheses. For this reason, we have little detailed, unbiased information about how molecular function evolves across large protein family phylogenies. Here we describe a generalized protocol that integrates ancestral sequence reconstruction with structural homology modeling and structure-based molecular affinity prediction to characterize historical changes in protein function across families with thousands of individual sequences. We highlight key steps in the analysis protocol requiring particularly careful attention to avoid introducing potential errors as well as steps for which computationally efficient subroutines can be substituted for more intensive approaches, allowing researchers to scale the analysis up or down, depending on available resources and requirements for reproducibility and scientific rigor. In our view, this approach provides a compelling compliment to more laboratory-intensive procedures, generating important contextual information that can help guide detailed experiments.

Heavy Metal Removal by Bioaccumulation Using Genetically Engineered Microorganisms

Wastewater effluents from mines and metal refineries are often contaminated with heavy metal ions, so they pose hazards to human and environmental health. Conventional technologies to remove heavy metal ions are well-established, but the most popular methods have drawbacks: chemical precipitation generates sludge waste, and activated carbon and ion exchange resins are made from unsustainable non-renewable resources. Using microbial biomass as the platform for heavy metal ion removal is an alternative method. Specifically, bioaccumulation is a natural biological phenomenon where microorganisms use proteins to uptake and sequester metal ions in the intracellular space to utilize in cellular processes (e.g., enzyme catalysis, signaling, stabilizing charges on biomolecules). Recombinant expression of these import-storage systems in genetically engineered microorganisms allows for enhanced uptake and sequestration of heavy metal ions. This has been studied for over two decades for bioremediative applications, but successful translation to industrial-scale processes is virtually non-existent. Meanwhile, demands for metal resources are increasing while discovery rates to supply primary grade ores are not. This review re-thinks how bioaccumulation can be used and proposes that it can be developed for bioextractive applications-the removal and recovery of heavy metal ions for downstream purification and refining, rather than disposal. This review consolidates previously tested import-storage systems into a biochemical framework and highlights efforts to overcome obstacles that limit industrial feasibility, thereby identifying gaps in knowledge and potential avenues of research in bioaccumulation.

A novel bacterial l-arginine sensor controlling c-di-GMP levels in Pseudomonas aeruginosa

Nutrients such as amino acids play key roles in shaping the metabolism of microorganisms in natural environments and in host-pathogen interactions. Beyond taking part to cellular metabolism and to protein synthesis, amino acids are also signaling molecules able to influence group behavior in microorganisms, such as biofilm formation. This lifestyle switch involves complex metabolic reprogramming controlled by local variation of the second messenger 3', 5'-cyclic diguanylic acid (c-di-GMP). The intracellular levels of this dinucleotide are finely tuned by the opposite activity of dedicated diguanylate cyclases (GGDEF signature) and phosphodiesterases (EAL and HD-GYP signatures), which are usually allosterically controlled by a plethora of environmental and metabolic clues. Among the genes putatively involved in controlling c-di-GMP levels in P. aeruginosa, we found that the multidomain transmembrane protein PA0575, bearing the tandem signature GGDEF-EAL, is an l-arginine sensor able to hydrolyse c-di-GMP. Here, we investigate the basis of arginine recognition by integrating bioinformatics, molecular biophysics and microbiology. Although the role of nutrients such as l-arginine in controlling the cellular fate in P. aeruginosa (including biofilm, pathogenicity and virulence) is already well established, we identified the first l-arginine sensor able to link environment sensing, c-di-GMP signaling and biofilm formation in this bacterium.

The burst-phase folding intermediate of ribonuclease H changes conformation over evolutionary history

The amino acid sequence encodes the energy landscape of a protein. Therefore, we expect evolutionary mutations to change features of the protein energy landscape, including the conformations adopted by a polypeptide as it folds to its native state. Ribonucleases H (RNase H) from Escherichia coli and Thermus thermophilus both fold via a partially folded intermediate in which the core region of the protein (helices A-D and strands 4-5) is structured. Strand 1, however, uniquely contributes to the T. thermophilus RNase H folding intermediate (Icore+1 ), but not the E. coli RNase H intermediate (Icore ) (Rosen & Marqusee, PLoS One 2015). We explore the origin of this difference by characterizing the folding intermediate of seven ancestral RNases H spanning the evolutionary history of these two homologs. Using fragment models with or without strand 1 and FRET probes to characterize the folding intermediate of each ancestor, we find a distinct evolutionary trend across the family-the involvement of strand 1 in the folding intermediate is an ancestral feature that is maintained in the thermophilic lineage and is gradually lost in the mesophilic lineage. Evolutionary sequence changes indeed modulate the conformations present on the folding landscape and altered the folding trajectory of RNase H.

Role of conformational dynamics in the evolution of novel enzyme function

The free energy landscape concept that describes enzymes as an ensemble of differently populated conformational sub-states in dynamic equilibrium is key for evaluating enzyme activity, enantioselectivity, and specificity. Mutations introduced in the enzyme sequence can alter the populations of the pre-existing conformational states, thus strongly modifying the enzyme ability to accommodate alternative substrates, revert its enantiopreferences, and even increase the activity for some residual promiscuous reactions. In this feature article, we present an overview of the current experimental and computational strategies to explore the conformational free energy landscape of enzymes. We provide a series of recent publications that highlight the key role of conformational dynamics for the enzyme evolution towards new functions and substrates, and provide some perspectives on how conformational dynamism should be considered in future computational enzyme design protocols.

Benchmark Analysis of Native and Artificial NAD+-Dependent Enzymes Generated by a Sequence-Based Design Method with or without Phylogenetic Data

The expansion of protein sequence databases has enabled us to design artificial proteins by sequence-based design methods, such as full-consensus design (FCD) and ancestral-sequence reconstruction (ASR). Artificial proteins with enhanced activity levels compared with native ones can potentially be generated by such methods, but successful design is rare because preparing a sequence library by curating the database and selecting a method is difficult. Utilizing a curated library prepared by reducing conservation energies, we successfully designed two artificial l-threonine 3-dehydrogenases (SDR-TDH) with higher activity levels than native SDR-TDH, FcTDH-N1, and AncTDH, using FCD and ASR, respectively. The artificial SDR-TDHs had excellent thermal stability and NAD⁺ recognition compared to native SDR-TDH from Cupriavidus necator (CnTDH); the melting temperatures of FcTDH-N1 and AncTDH were about 10 and 5 °C higher than that of CnTDH, respectively, and the dissociation constants toward NAD⁺ of FcTDH-N1 and AncTDH were 2- and 7-fold lower than that of CnTDH, respectively. Enzymatic efficiency of the artificial SDR-TDHs were comparable to that of CnTDH. Crystal structures of FcTDH-N1 and AncTDH were determined at 2.8 and 2.1 Å resolution, respectively. Structural and MD simulation analysis of the SDR-TDHs indicated that only the flexibility at specific regions was changed, suggesting that multiple mutations introduced in the artificial SDR-TDHs altered their flexibility and thereby affected their enzymatic properties. Benchmark analysis of the SDR-TDHs indicated that both FCD and ASR can generate highly functional proteins if a curated library is prepared appropriately.

Evolution of cyclohexadienyl dehydratase from an ancestral solute-binding protein

The emergence of enzymes through the neofunctionalization of noncatalytic proteins is ultimately responsible for the extraordinary range of biological catalysts observed in nature. Although the evolution of some enzymes from binding proteins can be inferred by homology, we have a limited understanding of the nature of the biochemical and biophysical adaptations along these evolutionary trajectories and the sequence in which they occurred. Here we reconstructed and characterized evolutionary intermediate states linking an ancestral solute-binding protein to the extant enzyme cyclohexadienyl dehydratase. We show how the intrinsic reactivity of a desolvated general acid was harnessed by a series of mutations radiating from the active site, which optimized enzyme-substrate complementarity and transition-state stabilization and minimized sampling of noncatalytic conformations. Our work reveals the molecular evolutionary processes that underlie the emergence of enzymes de novo, which are notably mirrored by recent examples of computational enzyme design and directed evolution.

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.