Engineering functional thermostable proteins using ancestral sequence reconstruction

Designing cytochrome P450 enzymes for use in cancer gene therapy

Cancer is a significant global socioeconomic burden, as millions of new cases and deaths occur annually. In 2020, almost 10 million cancer deaths were recorded worldwide. Advancements in cancer gene therapy have revolutionized the landscape of cancer treatment. An approach with promising potential for cancer gene therapy is introducing genes to cancer cells that encode for chemotherapy prodrug metabolizing enzymes, such as Cytochrome P450 (CYP) enzymes, which can contribute to the effective elimination of cancer cells. This can be achieved through gene-directed enzyme prodrug therapy (GDEPT). CYP enzymes can be genetically engineered to improve anticancer prodrug conversion to its active metabolites and to minimize chemotherapy side effects by reducing the prodrug dosage. Rational design, directed evolution, and phylogenetic methods are some approaches to developing tailored CYP enzymes for cancer therapy. Here, we provide a compilation of genetic modifications performed on CYP enzymes aiming to build highly efficient therapeutic genes capable of bio-activating different chemotherapeutic prodrugs. Additionally, this review summarizes promising preclinical and clinical trials highlighting engineered CYP enzymes' potential in GDEPT. Finally, the challenges, limitations, and future directions of using CYP enzymes for GDEPT in cancer gene therapy are discussed.

Self-controlled in silico gene knockdown strategies to enhance the sustainable production of heterologous terpenoid by Saccharomyces cerevisiae

Microbial bioengineering is a growing field for producing plant natural products (PNPs) in recent decades, using heterologous metabolic pathways in host cells. Once heterologous metabolic pathways have been introduced into host cells, traditional metabolic engineering techniques are employed to enhance the productivity and yield of PNP biosynthetic routes, as well as to manage competing pathways. The advent of computational biology has marked the beginning of a novel epoch in strain design through in silico methods. These methods utilize genome-scale metabolic models (GEMs) and flux optimization algorithms to facilitate rational design across the entire cellular metabolic network. However, the implementation of in silico strategies can often result in an uneven distribution of metabolic fluxes due to the rigid knocking out of endogenous genes, which can impede cell growth and ultimately impact the accumulation of target products. In this study, we creatively utilized synthetic biology to refine in silico strain design for efficient PNPs production. OptKnock simulation was performed on the GEM of Saccharomyces cerevisiae OA07, an engineered strain for oleanolic acid (OA) bioproduction that has been reported previously. The simulation predicted that the single deletion of fol1, fol2, fol3, abz1, and abz2, or a combined knockout of hfd1, ald2 and ald3 could improve its OA production. Consequently, strains EK1∼EK7 were constructed and cultivated. EK3 (OA07△fol3), EK5 (OA07△abz1), and EK6 (OA07△abz2) had significantly higher OA titers in a batch cultivation compared to the original strain OA07. However, these increases were less pronounced in the fed-batch mode, indicating that gene deletion did not support sustainable OA production. To address this, we designed a negative feedback circuit regulated by malonyl-CoA, a growth-associated intermediate whose synthesis served as a bypass to OA synthesis, at fol3, abz1, abz2, and at acetyl-CoA carboxylase-encoding gene acc1, to dynamically and autonomously regulate the expression of these genes in OA07. The constructed strains R_3A, R_5A and R_6A had significantly higher OA titers than the initial strain and the responding gene-knockout mutants in either batch or fed-batch culture modes. Among them, strain R_3A stand out with the highest OA titer reported to date. Its OA titer doubled that of the initial strain in the flask-level fed-batch cultivation, and achieved at 1.23 ± 0.04 g L-1 in 96 h in the fermenter-level fed-batch mode. This indicated that the integration of optimization algorithm and synthetic biology approaches was efficiently rational for PNP-producing strain design.

Engineering Biocatalysts for the C-H Activation of Fatty Acids by Ancestral Sequence Reconstruction

Selective, one-step C-H activation of fatty acids from biomass is an attractive concept in sustainable chemistry. Biocatalysis has shown promise for generating high-value hydroxy acids, but to date enzyme discovery has relied on laborious screening and produced limited hits, which predominantly oxidise the subterminal positions of fatty acids. Herein we show that ancestral sequence reconstruction (ASR) is an effective tool to explore the sequence-activity landscape of a family of multidomain, self-sufficient P450 monooxygenases. We resurrected 11 catalytically active CYP116B ancestors, each with a unique regioselectivity fingerprint that varied from subterminal in the older ancestors to mid-chain in the lineage leading to the extant, P450-TT. In lineages leading to extant enzymes in thermophiles, thermostability increased from ancestral to extant forms, as expected if thermophily had arisen de novo. Our studies show that ASR can be applied to multidomain enzymes to develop active, self-sufficient monooxygenases as regioselective biocatalysts for fatty acid hydroxylation.

Click, Compute, Create: A Review of Web-based Tools for Enzyme Engineering

Enzyme engineering, though pivotal across various biotechnological domains, is often plagued by its time-consuming and labor-intensive nature. This review aims to offer an overview of supportive in silico methodologies for this demanding endeavor. Starting from methods to predict protein structures, to classification of their activity and even the discovery of new enzymes we continue with describing tools used to increase thermostability and production yields of selected targets. Subsequently, we discuss computational methods to modulate both, the activity as well as selectivity of enzymes. Last, we present recent approaches based on cutting-edge machine learning methods to redesign enzymes. With exception of the last chapter, there is a strong focus on methods easily accessible via web-interfaces or simple Python-scripts, therefore readily useable for a diverse and broad community.

Structure and Activity of Reconstructed Pseudo-Ancestral Cyclotides

Cyclotides are cyclic peptides that are promising scaffolds for the design of drug candidates and chemical tools. However, despite there being hundreds of reported cyclotides, drug design studies have commonly focussed on a select few prototypic examples. Here, we explored whether ancestral sequence reconstruction could be used to generate new cyclotides for further optimization. We show that the reconstructed 'pseudo-ancestral' sequences, named Ancy-m (for the ancestral cyclotide of the Möbius sub-family) and Ancy-b (for the bracelet sub-family), have well-defined structures like their extant members, comprising the core structural feature of a cyclic cystine knot. This motif underpins efforts to re-engineer cyclotides for agrochemical and therapeutic applications. We further show that the reconstructed sequences are resistant to temperatures approaching boiling, bind to phosphatidyl-ethanolamine lipid bilayers at micromolar affinity, and inhibit the growth of insect cells at inhibitory concentrations in the micromolar range. Interestingly, the Ancy-b cyclotide had a higher oxidative folding yield than its comparator cyclotide cyO2, which belongs to the bracelet cyclotide subfamily known to be notoriously difficult to fold. Overall, this study provides new cyclotide sequences not yet found naturally that could be valuable starting points for the understanding of cyclotide evolution and for further optimization as drug leads.

Killer yeasts: expanding frontiers in the age of synthetic biology

Killer yeasts secrete protein toxins that are selectively lethal to other yeast and filamentous fungi. These exhibit exceptional genetic and functional diversity, and have several biotechnological applications. However, despite decades of research, several limitations hinder their widespread adoption. In this perspective we contend that technical advances in synthetic biology present an unprecedented opportunity to unlock the full potential of yeast killer systems across a spectrum of applications. By leveraging these new technologies, engineered killer toxins may emerge as a pivotal new tool to address antifungal resistance and food security. Finally, we speculate on the biotechnological potential of re-engineering host double-stranded (ds) RNA mycoviruses, from which many toxins derive, as a safe and noninfectious system to produce designer RNA.

Extant Sequence Reconstruction: The Accuracy of Ancestral Sequence Reconstructions Evaluated by Extant Sequence Cross-Validation

Ancestral sequence reconstruction (ASR) is a phylogenetic method widely used to analyze the properties of ancient biomolecules and to elucidate mechanisms of molecular evolution. Despite its increasingly widespread application, the accuracy of ASR is currently unknown, as it is generally impossible to compare resurrected proteins to the true ancestors. Which evolutionary models are best for ASR? How accurate are the resulting inferences? Here we answer these questions using a cross-validation method to reconstruct each extant sequence in an alignment with ASR methodology, a method we term "extant sequence reconstruction" (ESR). We thus can evaluate the accuracy of ASR methodology by comparing ESR reconstructions to the corresponding known true sequences. We find that a common measure of the quality of a reconstructed sequence, the average probability, is indeed a good estimate of the fraction of correct amino acids when the evolutionary model is accurate or overparameterized. However, the average probability is a poor measure for comparing reconstructions from different models, because, surprisingly, a more accurate phylogenetic model often results in reconstructions with lower probability. While better (more predictive) models may produce reconstructions with lower sequence identity to the true sequences, better models nevertheless produce reconstructions that are more biophysically similar to true ancestors. In addition, we find that a large fraction of sequences sampled from the reconstruction distribution may have fewer errors than the single most probable (SMP) sequence reconstruction, despite the fact that the SMP has the lowest expected error of all possible sequences. Our results emphasize the importance of model selection for ASR and the usefulness of sampling sequence reconstructions for analyzing ancestral protein properties. ESR is a powerful method for validating the evolutionary models used for ASR and can be applied in practice to any phylogenetic analysis of real biological sequences. Most significantly, ESR uses ASR methodology to provide a general method by which the biophysical properties of resurrected proteins can be compared to the properties of the true protein.

Evolution of Phosphorylase Activity in an Ancestral Glycosyltransferase

The reconstruction of ancestral sequences can offer a glimpse into the fascinating process of molecular evolution by exposing the adaptive pathways that shape the proteins found in nature today. Here, we track the evolution of the carbohydrate-active enzymes responsible for the synthesis and turnover of mannogen, a critical carbohydrate reserve in Leishmania parasites. Biochemical characterization of resurrected enzymes demonstrated that mannoside phosphorylase activity emerged in an ancestral bacterial mannosyltransferase, and later disappeared in the process of horizontal gene transfer and gene duplication in Leishmania. By shuffling through plausible historical sequence space in an ancestral mannosyltransferase, we found that mannoside phosphorylase activity could be toggled on through various combinations of mutations at positions outside of the active site. Molecular dynamics simulations showed that such mutations can affect loop rigidity and shield the active site from water molecules that disrupt key interactions, allowing α-mannose 1-phosphate to adopt a catalytically productive conformation. These findings highlight the importance of subtle distal mutations in protein evolution and suggest that the vast collection of natural glycosyltransferases may be a promising source of engineering templates for the design of tailored phosphorylases.

Computation-Aided Phylogeny-Oriented Engineering of β-Xylosidase: Modification of "Blades" to Enhance Stability and Activity for the Bioconversion of Hemicellulose to Produce Xylose

Hemicellulose is a highly abundant, ubiquitous, and renewable natural polysaccharide, widely present in agricultural and forestry residues. The enzymatic hydrolysis of hemicellulose has generally been accomplished using β-xylosidases, but concomitantly increasing the stability and activity of these enzymes remains challenging. Here, we rationally engineered a β-xylosidase from Bacillus clausii to enhance its stability by computation-aided design combining ancestral sequence reconstruction and structural analysis. The resulting combinatorial mutant rXYLOM25I/S51L/S79E exhibited highly improved robustness, with a 6.9-fold increase of the half-life at 60 °C, while also exhibiting improved pH stability, catalytic efficiency, and hydrolytic activity. Structural analysis demonstrated that additional interactions among the propeller blades in the catalytic module resulted in a much more compact protein structure and induced the rearrangement of the opposing catalytic pocket to mediate the observed improvement of activity. Our work provides a robust biocatalyst for the hydrolysis of agricultural waste to produce various high-value-added chemicals and biofuels.

Uncovering the role of leucine 59 in Renilla luciferase stability and activity with error-prone PCR: Implications for protein engineering

A new variant of Renilla luciferase, named Met C-SRLuc 8, was obtained from a random mutagenesis library and expressed in Escherichia coli BL21 (DE3) plys and purified. The results of the enzyme's binding affinity, kinetic stability, and bioinformatic studies demonstrated that leucine 59, located within the hot-spot foldon in the N-terminal domain of the protein, plays a significant role in the stability and activity of Renilla luciferase. These findings may facilitate the engineering of different variants of this enzyme to achieve thermally stable versions for various biotechnological applications.

Improving Protein Expression, Stability, and Function with ProteinMPNN

Natural proteins are highly optimized for function but are often difficult to produce at a scale suitable for biotechnological applications due to poor expression in heterologous systems, limited solubility, and sensitivity to temperature. Thus, a general method that improves the physical properties of native proteins while maintaining function could have wide utility for protein-based technologies. Here, we show that the deep neural network ProteinMPNN, together with evolutionary and structural information, provides a route to increasing protein expression, stability, and function. For both myoglobin and tobacco etch virus (TEV) protease, we generated designs with improved expression, elevated melting temperatures, and improved function. For TEV protease, we identified multiple designs with improved catalytic activity as compared to the parent sequence and previously reported TEV variants. Our approach should be broadly useful for improving the expression, stability, and function of biotechnologically important proteins.

Evolution and separation of actinobacterial pyranose and C-glycoside-3-oxidases

FAD-dependent pyranose oxidase (POx) and C-glycoside-3-oxidase (CGOx) are both members of the glucose-methanol-choline superfamily of oxidoreductases and belong to the same sequence space. Pyranose oxidases had been studied for their oxidation of monosaccharides such as D-glucose, but recently, a bacterial C-glycoside-3-oxidase that is phylogenetically related to POx and that reacts with C-glycosides such as carminic acid, mangiferin or puerarin has been described. Since these actinobacterial CGOx enzymes belong to the same sequence space as bacterial POx, they must have evolved from the same ancestor. Here, we performed a phylogenetic analysis of actinobacterial sequences and resurrected seven ancestral enzymes of the POx/CGOx sequence space to study the evolutionary trajectory of substrate preferences for monosaccharides and C-glycosides. Clade I, with its dimeric member POx from Kitasatospora aureofaciens, shows strict preference for monosaccharides (D-glucose and D-xylose) and does not react with any of the glycosides tested. No extant member of clade II has been studied to date. The two extant members of clades III and IV, monomeric POx/CGOx from Pseudoarthrobacter siccitolerans and Streptomyces canus, oxidized both monosaccharides as well as various C-glycosides (homoorientin, isovitexin, mangiferin, and puerarin). Steady-state kinetic parameters of several clades III and IV ancestral enzymes indicate that the generalist ancestor N35 slowly evolved to present-day enzymes with a much higher preference for C-glycosides than monosaccharides. Based on structural predictions of ancestors, we hypothesize that the strict specificity of bacterial clade I POx (and also fungal POx) is the result of oligomerization, which in turn results from the evolution of protein segments that were shown to be important for oligomerization, the arm, and the head domain.IMPORTANCEC-Glycosides often form active compounds in various plants. Breakage of the C-C bond in these glycosides to release the aglycone is challenging and proceeds via a two-step reaction, the oxidation of the sugar and subsequent cleavage of the C-C bond. Recently, an enzyme from a soil bacterium, FAD-dependent C-glycoside-3-oxidase (CGOx), was shown to catalyze the initial oxidation reaction. Here, we show that CGOx belongs to the same sequence space as pyranose oxidase (POx), and that an actinobacterial ancestor of the POx/CGOx family evolved into four clades, two of which show a high preference for C-glycosides.

Comparative analysis of reconstructed ancestral proteins with their extant counterparts suggests primitive life had an alkaline habitat

To understand the origin and early evolution of life it is crucial to establish characteristics of the primordial environment that facilitated the emergence and evolution of life. One important environmental factor is the pH of the primordial environment. Here, we assessed the pH-dependent thermal stabilities of previously reconstructed ancestral nucleoside diphosphate kinases and ribosomal protein uS8s. The selected proteins were likely to be present in ancient organisms such as the last common ancestor of bacteria and that of archaea. We also assessed the thermal stability of homologous proteins from extant acidophilic, neutralophilic, and alkaliphilic microorganisms as a function of pH. Our results indicate that the reconstructed ancestral proteins are more akin to those of extant alkaliphilic bacteria, which display greater stability under alkaline conditions. These findings suggest that the common ancestors of bacterial and archaeal species thrived in an alkaline environment. Moreover, we demonstrate the reconstruction method employed in this study is a valuable technique for generating alkali-tolerant proteins that can be used in a variety of biotechnological and environmental applications.

Study on CAT activity of tomato leaf cells under salt stress based on microhyperspectral imaging and transfer learning algorithm

Salt stress easily leads to oxidative stress and promotes the catalase (CAT) response in tomato leaves. For the changes in catalase activity in leaf subcells, there is a need for a visual in situ detection method and mechanism analysis. This paper, taking catalase in leaf subcells under salt stress as the starting point, describes the use of microscopic hyperspectral imaging technology to dynamically detect and study catalase activity from a microscopic perspective, and lay the theoretical foundation for exploring the detection limit of catalase activity under salt stress. In this study, a total of 298 microscopic images were obtained under different concentrations of salt stress (0 g/L, 1 g/L, 2 g/L, 3 g/L) in the spectral range of 400-1000 nm. With the increase in salt solution concentration and the advancement of the growth period, the CAT activity value increased. Regions of interest were extracted according to the reflectance of the samples, and the model was established by combining CAT activity. The characteristic wavelength was extracted by five methods (SPA, IVISSA, IRFJ, GAPLSR and CARS), and four models (PLSR, PCR, CNN and LSSVM) were established according to the characteristic wavelengths. The results show that the random sampling (RS) method was better for the selection samples of the correction set and prediction set. Raw wavelengths are optimized as the pretreatment method. The partial least-squares regression model based on the IRFJ method is the best, and the coefficient of correlation (Rp) and root mean square error of the prediction set (RMSEP) are 0.81 and 58.03 U/g, respectively. According to the ratio of microarea area to the area of the macroscopic tomato leaf slice, the Rp and RMSEP of the prediction model for the detection of microarea cells are 0.71 and 23.00 U/g, respectively. Finally, the optimal model was used for quantitative visualization analysis of CAT activity in tomato leaves, and the distribution of CAT activity was consistent with its color trend. The results show that it is feasible to detect the CAT activity in tomato leaves by microhyperspectral imaging combined with stoichiometry.

Impact of ancestral sequence reconstruction on mechanistic and structural enzymology

Ancestral sequence reconstruction (ASR) provides insight into the changes within a protein sequence across evolution. More specifically, it can illustrate how specific amino acid changes give rise to different phenotypes within a protein family. Over the last few decades it has established itself as a powerful technique for revealing molecular common denominators that govern enzyme function. Here, we describe the strength of ASR in unveiling catalytic mechanisms and emerging phenotypes for a range of different proteins, also highlighting biotechnological applications the methodology can provide.

Computer-Aided Lipase Engineering for Improving Their Stability and Activity in the Food Industry: State of the Art

As some of the most widely used biocatalysts, lipases have exhibited extreme advantages in many processes, such as esterification, amidation, and transesterification reactions, which causes them to be widely used in food industrial production. However, natural lipases have drawbacks in terms of organic solvent resistance, thermostability, selectivity, etc., which limits some of their applications in the field of foods. In this systematic review, the application of lipases in various food processes was summarized. Moreover, the general structure of lipases is discussed in-depth, and the engineering strategies that can be used in lipase engineering are also summarized. The protocols of some classical methods are compared and discussed, which can provide some information about how to choose methods of lipase engineering. Thermostability engineering and solvent tolerance engineering are highlighted in this review, and the basic principles for improving thermostability and solvent tolerance are summarized. In the future, comput er-aided technology should be more emphasized in the investigation of the mechanisms of reactions catalyzed by lipases and guide the engineering of lipases. The engineering of lipase tunnels to improve the diffusion of substrates is also a promising prospect for further enhanced lipase activity and selectivity.

Ion-Powered Rotary Motors: Where Did They Come from and Where They Are Going?

Molecular motors are found in many living organisms. One such molecular machine, the ion-powered rotary motor (IRM), requires the movement of ions across a membrane against a concentration gradient to drive rotational movement. The bacterial flagellar motor (BFM) is an example of an IRM which relies on ion movement through the stator proteins to generate the rotation of the flagella. There are many ions which can be used by the BFM stators to power motility and different ions can be used by a single bacterium expressing multiple stator variants. The use of ancestral sequence reconstruction (ASR) and functional analysis of reconstructed stators shows promise for understanding how these proteins evolved and when the divergence in ion use may have occurred. In this review, we discuss extant BFM stators and the ions that power them as well as recent examples of the use of ASR to study ion-channel selectivity and how this might be applied to further study of the BFM stator complex.

Recombinant pseudorabies virus (PRV) expressing stabilized E2 of classical swine fever virus (CSFV) protects against both PRV and CSFV

Pseudorabies (PR) and classical swine fever (CSF) are economically important infectious diseases of pigs. Most pig farms in China are immunized against these two diseases. Here, we describe a stabilized E2 protein as an immunogen inserted into the PRV genome as a bivalent live virus-vectored vaccine. The E2 protein has 48 variant sites, there are 2-5 candidate amino acids per variant site, and the relative energy contribution of each amino acid to E2 energy was calculated. Combined substitutions of amino acids at the neighbor variant site (neighbor substitution) were performed to obtain the E2 protein sequence with the lowest energy (stabilized E2). Multiple amino acid substitutions at 48 variant sites were performed, and the results were consistent with neighbor substitutions. The stabilized E2 sequence was obtained, and its energy decreased by 22 Rosetta Energy Units (REUs) compared with the original sequence. After the recombinant PRV expressing stabilized E2 of CSFV was constructed, the secretion efficiency of stabilized E2 was increased by 2.97 times, and the thermal stability was increased by 10.5 times. Immunization of mice resulted in a 2-fold increase in antibody production, and a balanced antibody level against subtype 1.1 and subtype 2.1d E2 was achieved. In rabbits immunized, the lethal challenge of PRV-ZJ and the fever response induced by CSFV could be prevented simultaneously. These findings suggest that rPRV-muta/287aaE2 is a promising bivalent vaccine against CSFV and PRV infections.