Rapid directed molecular evolution of fluorescent proteins in mammalian cells

In vivo imaging of model organisms is heavily reliant on fluorescent proteins with high intracellular brightness. Here we describe a practical method for rapid optimization of fluorescent proteins via directed molecular evolution in cultured mammalian cells. Using this method, we were able to perform screening of large gene libraries containing up to 2 × 10⁷ independent random genes of fluorescent proteins expressed in HEK cells, completing one iteration of directed evolution in a course of 8 days. We employed this approach to develop a set of green and near‐infrared fluorescent proteins with enhanced intracellular brightness. The developed near‐infrared fluorescent proteins demonstrated high performance for fluorescent labeling of neurons in culture and in vivo in model organisms such as Caenorhabditis elegans, Drosophila, zebrafish, and mice. Spectral properties of the optimized near‐infrared fluorescent proteins enabled crosstalk‐free multicolor imaging in combination with common green and red fluorescent proteins, as well as dual‐color near‐infrared fluorescence imaging. The described method has a great potential to be adopted by protein engineers due to its simplicity and practicality. We also believe that the new enhanced fluorescent proteins will find wide application for in vivo multicolor imaging of small model organisms.     

Phylogeny and evolution of chloroplast tRNAs in Adoxaceae

Chloroplasts are semiautonomous organelles found in photosynthetic plants. The major functions of chloroplasts include photosynthesis and carbon fixation, which are mainly regulated by its circular genomes. In the highly conserved chloroplast genome, the chloroplast transfer RNA genes (cp tRNA) play important roles in protein translation within chloroplasts. However, the evolution of cp tRNAs remains unclear. Thus, in the present study, we investigated the evolutionary characteristics of chloroplast tRNAs in five Adoxaceae species using 185 tRNA gene sequences. In total, 37 tRNAs encoding 28 anticodons are found in the chloroplast genome in Adoxaceae species. Some consensus sequences are found within the Ψ‐stem and anticodon loop of the tRNAs. Some putative novel structures were also identified, including a new stem located in the variable region of tRNA^Tyr in a similar manner to the anticodon stem. Furthermore, phylogenetic and evolutionary analyses indicated that synonymous tRNAs may have evolved from multiple ancestors and frequent tRNA duplications during the evolutionary process may have been primarily caused by positive selection and adaptive evolution. The transition and transversion rates are uneven among different tRNA isotypes. For all tRNAs, the transition rate is greater with a transition/transversion bias of 3.13. Phylogenetic analysis of cp tRNA suggested that the type I introns in different taxa (including eukaryote organisms and cyanobacteria) share the conserved sequences “U‐U‐x2‐C” and “U‐x‐G‐x2‐T,” thereby indicating the diverse cyanobacterial origins of organelles. This detailed study of cp tRNAs in Adoxaceae may facilitate further investigations of the evolution, phylogeny, structure, and related functions of chloroplast tRNAs.     

Molecular basis of cold adaptation

Cold-adapted, or psychrophilic, organisms are able to thrive at low temperatures in permanently cold environments, which in fact characterize the greatest proportion of our planet. Psychrophiles include both prokaryotic and eukaryotic organisms and thus represent a significant proportion of the living world. These organisms produce cold-evolved enzymes that are partially able to cope with the reduction in chemical reaction rates induced by low temperatures. As a rule, cold-active enzymes display a high catalytic efficiency, associated however, with a low thermal stability. In most cases, the adaptation to cold is achieved through a reduction in the activation energy that possibly originates from an increased flexibility of either a selected area or of the overall protein structure. This enhanced plasticity seems in turn to be induced by the weak thermal stability of psychrophilic enzymes. The adaptation strategies are beginning to be understood thanks to recent advances in the elucidation of the molecular characteristics of cold-adapted enzymes derived from X-ray crystallography, protein engineering and biophysical methods. Psychrophilic organisms and their enzymes have, in recent years, increasingly attracted the attention of the scientific community due to their peculiar properties that render them particularly useful in investigating the possible relationship existing between stability, flexibility and specific activity and as valuable tools for biotechnological purposes.     

Hexamerization and thermostability emerged very early during geranylgeranylglyceryl phosphate synthase evolution

A large number of archaea live in hyperthermophilic environments. In consequence, their proteins need to adopt to these harsh conditions, including the enzymes that catalyze the synthesis of their membrane ether lipids. The enzyme that catalyzes the formation of the first ether bond in these lipids, geranylgeranylglyceryl phosphate synthase (GGGPS), exists as a hexamer in many hyperthermophilic archaea, and a recent study suggested that hexamerization serves for a fine‐tuning of the flexibility – stability trade‐off under hyperthermophilic conditions. We have recently reconstructed the sequences of ancestral group II GGGPS enzymes and now present a detailed biochemical characterization of nine of these predecessors, which allowed us to trace back the evolution of hexameric GGGPS and to draw conclusions about the properties of extant GGGPS branches that were not accessible to experiments up to now. Almost all ancestral GGGPS proteins formed hexamers, which demonstrates that hexamerization is even more widespread among the GGGPS family than previously assumed. Furthermore, all experimentally studied ancestral proteins showed high thermostability. Our results indicate that the hexameric oligomerization state and thermostability were present very early during the evolution of group II GGGPS, while the fine tuning of the flexibility – stability trade‐off developed very late, independent of the emergence of hexamerization.     

Evolutionary conservation of complexins: from choanoflagellates to mice

Complexins are synaptic SNARE complex‐binding proteins that cooperate with synaptotagmins in activating Ca²⁺‐stimulated, synaptotagmin‐dependent synaptic vesicle exocytosis and in clamping spontaneous, synaptotagmin‐independent synaptic vesicle exocytosis. Here, we show that complexin sequences are conserved in some non‐metazoan unicellular organisms and in all metazoans, suggesting that complexins are a universal feature of metazoans that predate metazoan evolution. We show that complexin from Nematostella vectensis, a cnidarian sea anemone far separated from mammals in metazoan evolution, functionally replaces mouse complexins in activating Ca²⁺‐triggered exocytosis, but is unable to clamp spontaneous exocytosis. Thus, the activating function of complexins is likely conserved throughout metazoan evolution.     

Extension of a de novo TIM barrel with a rationally designed secondary structure element

The ability to construct novel enzymes is a major aim in de novo protein design. A popular enzyme fold for design attempts is the TIM barrel. This fold is a common topology for enzymes and can harbor many diverse reactions. The recent de novo design of a four‐fold symmetric TIM barrel provides a well understood minimal scaffold for potential enzyme designs. Here we explore opportunities to extend and diversify this scaffold by adding a short de novo helix on top of the barrel. Due to the size of the protein, we developed a design pipeline based on computational ab initio folding that solves a less complex sub‐problem focused around the helix and its vicinity and adapt it to the entire protein. We provide biochemical characterization and a high‐resolution X‐ray structure for one variant and compare it to our design model. The successful extension of this robust TIM‐barrel scaffold opens opportunities to diversify it towards more pocket like arrangements and as such can be considered a building block for future design of binding or catalytic sites.     

An ultra‐high affinity protein–protein interface displaying sequence‐robustness

Protein–protein interactions are crucial in biology and play roles in for example, the immune system, signaling pathways, and enzyme regulation. Ultra‐high affinity interactions (K _d <0.1 nM) occur in these systems, however, structures and energetics behind stability of ultra‐high affinity protein–protein complexes are not well understood. Regulation of the starch debranching barley limit dextrinase (LD) and its endogenous cereal type inhibitor (LDI) exemplifies an ultra‐high affinity complex (K _d of 42 pM). In this study the LD–LDI complex is investigated to unveil how robust the ultra‐high affinity is to LDI sequence variation at the protein–protein interface and whether alternative sequences can retain the ultra‐high binding affinity. The interface of LD–LDI was engineered using computational protein redesign aiming at identifying LDI variants predicted to retain ultra‐high binding affinity. These variants present a very diverse set of mutations going beyond conservative and alanine substitutions typically used to probe interfaces. Surface plasmon resonance analysis of the LDI variants revealed that high affinity of LD–LDI requires interactions of several residues at the rim of the protein interface, unlike the classical hotspot arrangement where key residues are found at the center of the interface. Notably, substitution of interface residues in LDI, including amino acids with functional groups different from the wild‐type, could occur without loss of affinity. This demonstrates that ultra‐high binding affinity can be conferred without hotspot residues, thus making complexes more robust to mutational drift in evolution. The present mutational analysis also demonstrates how energetic coupling can emerge between residues at large distances at the interface.     

Directed evolution of enzyme stability

Modern enzyme development relies to an increasing extent on strategies based on diversity generation followed by screening for variants with optimised properties. In principle, these directed evolution strategies might be used for optimising any enzyme property, which can be screened for in an economically feasible way, even if the molecular basis of that property is not known. Stability is an interesting property of enzymes because (1) it is of great industrial importance, (2) it is relatively easy to screen for, and (3) the molecular basis of stability relates closely to contemporary issues in protein science such as the protein folding problem and protein folding diseases. Thus, engineering enzyme stability is of both commercial and scientific interest. Here, we review how directed evolution has contributed to the development of stable enzymes and to new insight into the principles of protein stability. Several recent examples are described. These examples show that directed evolution is an effective strategy to obtain stable enzymes, especially when used in combination with rational or semi-rational engineering strategies. With respect to the principles of protein stability, some important lessons to learn from recent efforts in directed evolution are (1) that there are many structural ways to stabilize a protein, which are not always easy to rationalize, (2) that proteins may very well be stabilized by optimizing their surfaces, and (3) that high thermal stability may be obtained without forfeiture of catalytic performance at low temperatures.     

The C‐degron pathway eliminates mislocalized proteins and products of deubiquitinating enzymes

Protein termini are determinants of protein stability. Proteins bearing degradation signals, or degrons, at their amino‐ or carboxyl‐termini are eliminated by the N‐ or C‐degron pathways, respectively. We aimed to elucidate the function of C‐degron pathways and to unveil how normal proteomes are exempt from C‐degron pathway‐mediated destruction. Our data reveal that C‐degron pathways remove mislocalized cellular proteins and cleavage products of deubiquitinating enzymes. Furthermore, the C‐degron and N‐degron pathways cooperate in protein removal. Proteome analysis revealed a shortfall in normal proteins targeted by C‐degron pathways, but not of defective proteins, suggesting proteolysis‐based immunity as a constraint for protein evolution/selection. Our work highlights the importance of protein termini for protein quality surveillance, and the relationship between the functional proteome and protein degradation pathways.     

“Protein aggregates” contain RNA and DNA,entrapped by misfolded proteins but largely rescued by slowing translational elongation

All neurodegenerative diseases feature aggregates, which usually contain disease‐specific diagnostic proteins; non‐protein constituents, however, have rarely been explored. Aggregates from SY5Y‐APP_Sw neuroblastoma, a cell model of familial Alzheimer''s disease, were crosslinked and sequences of linked peptides identified. We constructed a normalized “contactome” comprising 11 subnetworks, centered on 24 high‐connectivity hubs. Remarkably, all 24 are nucleic acid‐binding proteins. This led us to isolate and sequence RNA and DNA from Alzheimer''s and control aggregates. RNA fragments were mapped to the human genome by RNA‐seq and DNA by ChIP‐seq. Nearly all aggregate RNA sequences mapped to specific genes, whereas DNA fragments were predominantly intergenic. These nucleic acid mappings are all significantly nonrandom, making an artifactual origin extremely unlikely. RNA (mostly cytoplasmic) exceeded DNA (chiefly nuclear) by twofold to fivefold. RNA fragments recovered from AD tissue were ~1.5‐to 2.5‐fold more abundant than those recovered from control tissue, similar to the increase in protein. Aggregate abundances of specific RNA sequences were strikingly differential between cultured SY5Y‐APP_Sw glioblastoma cells expressing APOE3 vs. APOE4, consistent with APOE4 competition for E‐box/CLEAR motifs. We identified many G‐quadruplex and viral sequences within RNA and DNA of aggregates, suggesting that sequestration of viral genomes may have driven the evolution of disordered nucleic acid‐binding proteins. After RNA‐interference knockdown of the translational‐procession factor EEF2 to suppress translation in SY5Y‐APP_Sw cells, the RNA content of aggregates declined by >90%, while reducing protein content by only 30% and altering DNA content by ≤10%. This implies that cotranslational misfolding of nascent proteins may ensnare polysomes into aggregates, accounting for most of their RNA content.     

Molecular Structural Basis for the Cold Adaptedness of the Psychrophilic β-Glucosidase BglU in Micrococcus antarcticus

Psychrophilic enzymes play crucial roles in cold adaptation of microbes and provide useful models for studies of protein evolution, folding, and dynamic properties. We examined the crystal structure (2.2-Å resolution) of the psychrophilic β-glucosidase BglU, a member of the glycosyl hydrolase 1 (GH1) enzyme family found in the cold-adapted bacterium Micrococcus antarcticus. Structural comparison and sequence alignment between BglU and its mesophilic and thermophilic counterpart enzymes (BglB and GlyTn, respectively) revealed two notable features distinct to BglU: (i) a unique long-loop L3 (35 versus 7 amino acids in others) involved in substrate binding and (ii) a unique amino acid, His299 (Tyr in others), involved in the stabilization of an ordered water molecule chain. Shortening of loop L3 to 25 amino acids reduced low-temperature catalytic activity, substrate-binding ability, the optimal temperature, and the melting temperature (T_m). Mutation of His299 to Tyr increased the optimal temperature, the T_m, and the catalytic activity. Conversely, mutation of Tyr301 to His in BglB caused a reduction in catalytic activity, thermostability, and the optimal temperature (45 to 35°C). Loop L3 shortening and H299Y substitution jointly restored enzyme activity to the level of BglU, but at moderate temperatures. Our findings indicate that loop L3 controls the level of catalytic activity at low temperatures, residue His299 is responsible for thermolability (particularly heat lability of the active center), and long-loop L3 and His299 are jointly responsible for the psychrophilic properties. The described structural basis for the cold adaptedness of BglU will be helpful for structure-based engineering of new cold-adapted enzymes and for the production of mutants useful in a variety of industrial processes at different temperatures.     

Chimeric dihydrofolate reductases display properties of modularity and biophysical diversity

While reverse genetics and functional genomics have long affirmed the role of individual mutations in determining protein function, there have been fewer studies addressing how large‐scale changes in protein sequences, such as in entire modular segments, influence protein function and evolution. Given how recombination can reassort protein sequences, these types of changes may play an underappreciated role in how novel protein functions evolve in nature. Such studies could aid our understanding of whether certain organismal phenotypes related to protein function—such as growth in the presence or absence of an antibiotic—are robust with respect to the identity of certain modular segments. In this study, we combine molecular genetics with biochemical and biophysical methods to gain a better understanding of protein modularity in dihydrofolate reductase (DHFR), an enzyme target of antibiotics also widely used as a model for protein evolution. We replace an integral α‐helical segment of Escherichia coli DHFR with segments from a number of different organisms (many nonmicrobial) and examine how these chimeric enzymes affect organismal phenotypes (e.g., resistance to an antibiotic) as well as biophysical properties of the enzyme (e.g., thermostability). We find that organismal phenotypes and enzyme properties are highly sensitive to the identity of DHFR modules, and that this chimeric approach can create enzymes with diverse biophysical characteristics.     

Hyperthermophilic Enzymes: Sources, Uses, and Molecular Mechanisms for Thermostability

Enzymes synthesized by hyperthermophiles (bacteria and archaea with optimal growth temperatures of >80°C), also called hyperthermophilic enzymes, are typically thermostable (i.e., resistant to irreversible inactivation at high temperatures) and are optimally active at high temperatures. These enzymes share the same catalytic mechanisms with their mesophilic counterparts. When cloned and expressed in mesophilic hosts, hyperthermophilic enzymes usually retain their thermal properties, indicating that these properties are genetically encoded. Sequence alignments, amino acid content comparisons, crystal structure comparisons, and mutagenesis experiments indicate that hyperthermophilic enzymes are, indeed, very similar to their mesophilic homologues. No single mechanism is responsible for the remarkable stability of hyperthermophilic enzymes. Increased thermostability must be found, instead, in a small number of highly specific alterations that often do not obey any obvious traffic rules. After briefly discussing the diversity of hyperthermophilic organisms, this review concentrates on the remarkable thermostability of their enzymes. The biochemical and molecular properties of hyperthermophilic enzymes are described. Mechanisms responsible for protein inactivation are reviewed. The molecular mechanisms involved in protein thermostabilization are discussed, including ion pairs, hydrogen bonds, hydrophobic interactions, disulfide bridges, packing, decrease of the entropy of unfolding, and intersubunit interactions. Finally, current uses and potential applications of thermophilic and hyperthermophilic enzymes as research reagents and as catalysts for industrial processes are described.     

Innovation and tinkering in the evolution of oxidases

Although molecular oxygen is a relative newcomer to the biosphere, it has had a profound impact on metabolism. About 700 oxygen‐dependent enzymatic reactions are known, the vast majority of which emerged only after the appearance of oxygen in the biosphere, circa 3 billion years ago. Oxygen was a major driving force for evolutionary innovation—~60% of all known oxygen‐dependent enzyme families emerged as such; that is, the founding ancestor was an O₂‐dependent enzyme. The other 40% seem to have diverged by tinkering from pre‐existing proteins whose function was not related to oxygen. Here, we focus on the latter. We describe transitions from various enzyme classes, as well as from non‐enzymatic proteins, and we explore these transitions in terms of catalytic chemistry, metabolism, and protein structure. These transitions vary from subtle ones, such as simply repurposing oxidoreductases by replacing an electron acceptor such as NAD by O₂, to drastic changes in reaction mechanism, such as turning carboxylases and hydrolases into oxidases. The latter is more common and can occur with strikingly minor changes, for example, only one mutation in the active site. We further suggest that engineering enzymes to harness the extraordinary reactivity of oxygen may yield higher catabolic power and versatility.     

Using continuous directed evolution to improve enzymes for plant applications

Continuous directed evolution of enzymes and other proteins in microbial hosts is capable of outperforming classical directed evolution by executing hypermutation and selection concurrently in vivo, at scale, with minimal manual input. Provided that a target enzyme’s activity can be coupled to growth of the host cells, the activity can be improved simply by selecting for growth. Like all directed evolution, the continuous version requires no prior mechanistic knowledge of the target. Continuous directed evolution is thus a powerful way to modify plant or non-plant enzymes for use in plant metabolic research and engineering. Here, we first describe the basic features of the yeast (Saccharomyces cerevisiae) OrthoRep system for continuous directed evolution and compare it briefly with other systems. We then give a step-by-step account of three ways in which OrthoRep can be deployed to evolve primary metabolic enzymes, using a THI4 thiazole synthase as an example and illustrating the mutational outcomes obtained. We close by outlining applications of OrthoRep that serve growing demands (i) to change the characteristics of plant enzymes destined for return to plants, and (ii) to adapt (“plantize”) enzymes from prokaryotes—especially exotic prokaryotes—to function well in mild, plant-like conditions.

Continuous directed evolution using the yeast OrthoRep system is a powerful way to improve enzymes for use in plant engineering as illustrated by “plantizing” a bacterial thiamin synthesis enzyme.     

Prognostic and predictive role of a metabolic rate‐limiting enzyme signature in hepatocellular carcinoma

ObjectivesAbnormal expression of metabolic rate‐limiting enzymes drives the occurrence and progression of hepatocellular carcinoma (HCC). This study aimed to elucidate the comprehensive model of metabolic rate‐limiting enzymes associated with the prognosis of HCC.Materials and MethodsHCC animal model and TCGA project were used to screen out differentially expressed metabolic rate‐limiting enzyme. Cox regression, least absolute shrinkage and selection operation (LASSO) and experimentally verification were performed to identify metabolic rate‐limiting enzyme signature. The area under the receiver operating characteristic curve (AUC) and prognostic nomogram were used to assess the efficacy of the signature in the three HCC cohorts (TCGA training cohort, internal cohort and an independent validation cohort).ResultsA classifier based on three rate‐limiting enzymes (RRM1, UCK2 and G6PD) was conducted and serves as independent prognostic factor. This effect was further confirmed in an independent cohort, which indicated that the AUC at year 5 was 0.715 (95% CI: 0.653‐0.777) for clinical risk score, whereas it was significantly increased to 0.852 (95% CI: 0.798‐0.906) when combination of the clinical with signature risk score. Moreover, a comprehensive nomogram including the signature and clinicopathological aspects resulted in significantly predict the individual outcomes.ConclusionsOur results highlighted the prognostic value of rate‐limiting enzymes in HCC, which may be useful for accurate risk assessment in guiding clinical management and treatment decisions.     

GRASP55 regulates intra‐Golgi localization of glycosylation enzymes to control glycosphingolipid biosynthesis

The Golgi apparatus, the main glycosylation station of the cell, consists of a stack of discontinuous cisternae. Glycosylation enzymes are usually concentrated in one or two specific cisternae along the cis‐trans axis of the organelle. How such compartmentalized localization of enzymes is achieved and how it contributes to glycosylation are not clear. Here, we show that the Golgi matrix protein GRASP55 directs the compartmentalized localization of key enzymes involved in glycosphingolipid (GSL) biosynthesis. GRASP55 binds to these enzymes and prevents their entry into COPI‐based retrograde transport vesicles, thus concentrating them in the trans‐Golgi. In genome‐edited cells lacking GRASP55, or in cells expressing mutant enzymes without GRASP55 binding sites, these enzymes relocate to the cis‐Golgi, which affects glycosphingolipid biosynthesis by changing flux across metabolic branch points. These findings reveal a mechanism by which a matrix protein regulates polarized localization of glycosylation enzymes in the Golgi and controls competition in glycan biosynthesis.     

Mycobacterium tuberculosis protein kinase G acts as an unusual ubiquitinating enzyme to impair host immunity

Upon Mycobacterium tuberculosis (Mtb) infection, protein kinase G (PknG), a eukaryotic‐type serine‐threonine protein kinase (STPK), is secreted into host macrophages to promote intracellular survival of the pathogen. However, the mechanisms underlying this PknG–host interaction remain unclear. Here, we demonstrate that PknG serves both as a ubiquitin‐activating enzyme (E1) and a ubiquitin ligase (E3) to trigger the ubiquitination and degradation of tumor necrosis factor receptor‐associated factor 2 (TRAF2) and TGF‐β‐activated kinase 1 (TAK1), thereby inhibiting the activation of NF‐κB signaling and host innate responses. PknG promotes the attachment of ubiquitin (Ub) to the ubiquitin‐conjugating enzyme (E2) UbcH7 via an isopeptide bond (UbcH7 K82‐Ub), rather than the usual C86‐Ub thiol‐ester bond. PknG induces the discharge of Ub from UbcH7 by acting as an isopeptidase, before attaching Ub to its substrates. These results demonstrate that PknG acts as an unusual ubiquitinating enzyme to remove key components of the innate immunity system, thus providing a potential target for tuberculosis treatment.