A mechanistic view of enzyme evolution

Unlocking the potential of stimuli-responsive biomaterials for bone regeneration

Bone defects caused by tumors, osteoarthritis, and osteoporosis attract great attention. Because of outstanding biocompatibility, osteogenesis promotion, and less secondary infection incidence ratio, stimuli-responsive biomaterials are increasingly used to manage this issue. These biomaterials respond to certain stimuli, changing their mechanical properties, shape, or drug release rate accordingly. Thereafter, the activated materials exert instructive or triggering effects on cells and tissues, match the properties of the original bone tissues, establish tight connection with ambient hard tissue, and provide suitable mechanical strength. In this review, basic definitions of different categories of stimuli-responsive biomaterials are presented. Moreover, possible mechanisms, advanced studies, and pros and cons of each classification are discussed and analyzed. This review aims to provide an outlook on the future developments in stimuli-responsive biomaterials.

Potential involvement of beta-lactamase homologous proteins in resistance to beta-lactam antibiotics in gram-negative bacteria of the ESKAPEE group

BACKGROUND

Enzymatic degradation mediated by beta-lactamases constitutes one of the primary mechanisms of resistance to beta-lactam antibiotics in gram-negative bacteria. This enzyme family comprises four molecular classes, categorized into serine beta-lactamases (Classes A, C, and D) and zinc-dependent metallo-beta-lactamases (Class B). Gram-negative bacteria producing beta-lactamase are of significant concern, particularly due to their prevalence in nosocomial infections. A comprehensive understanding of the evolution and dissemination of this enzyme family is essential for effective control of these pathogens. In this study, we conducted the prospecting, phylogenetic analysis, and in silico analysis of beta-lactamases and homologous proteins identified in 1827 bacterial genomes with phenotypic data on beta-lactam resistance. These genomes were distributed among Klebsiella pneumoniae (45%), Acinetobacter baumannii (31%), Pseudomonas aeruginosa (14%), Escherichia coli (6%), and Enterobacter spp. (4%). Using an HMM profile and searching for conserved domains, we mined 2514, 8733, 5424, and 2957 proteins for molecular classes A, B, C, and D, respectively. This set of proteins encompasses canonical subfamilies of beta-lactamases as well as hypothetical proteins and other functional groups. Canonical beta-lactamases were found to be phylogenetically distant from hypothetical proteins, which, in turn, are closer to other representatives of the penicillin-binding-protein (PBP-like) and metallo-beta-lactamase (MBL) families. The catalytic amino acid residues characteristic of beta-lactamases were identified from the sequence alignment and revealed that motifs are less conserved in homologous groups than in beta-lactamases. After comparing the frequency of protein groups in genomes of resistant strains with those of sensitive ones applying Fisher's exact test and relative risk, it was observed that some groups of homologous proteins to classes B and C are more common in the genomes of resistant strains, particularly to carbapenems. We identified the beta-lactamase-like domain widely distributed in gram-negative species of the ESKAPEE group, which highlights its importance in the context of beta-lactam resistance. Some hypothetical homologous proteins have been shown to potentially possess promiscuous activity against beta-lactam antibiotics, however, they do not appear to expressly determine the resistance phenotype. The selective pressure due to the widespread use of antibiotics may favor the optimization of these functions for specialized resistance enzymes.

Modulating the pH profile of vanillin dehydrogenase enzyme from extremophile Bacillus ligniniphilus L1 through computational guided site-directed mutagenesis

Vanillin dehydrogenase (VDH) has recently come forward as an important enzyme for the commercial production of vanillic acid from vanillin in a one-step enzymatic process. However, VDH with high alkaline tolerance and efficiency is desirable to meet the biorefinery requirements. In this study, computationally guided site-directed mutagenesis was performed by increasing the positive and negative charges on the surface and near the active site of the VDH from the alkaliphilic marine bacterium Bacillus ligniniphilus L1, respectively. In total, 20 residues including 15 from surface amino acids and 5 near active sites were selected based on computational analysis and were subjected to site-directed mutations. The optimum pH of the two screened mutants including I132R, and T235E from surface residue and near active site mutant was shifted to 9, and 8.6, with a 2.82- and 2.95-fold increase in their activity compared to wild enzyme at pH 9, respectively. A double mutant containing both these mutations i.e., I132R/T235E was produced which showed a shift in optimum pH of VDH from 7.4 to 9, with an increase of 74.91 % in enzyme activity. Therefore, the double mutant of VDH from the L1 strain (I132R/T235E) produced in this study represents a potential candidate for industrial applications.

Peptide and Enzyme Catalysts Work in Concert in Stereoselective Cascade Reactions-Oxidation followed by Conjugate Addition

Enzymes and peptide catalysts consist of the same building blocks but require vastly different environments to operate best. Herein, we show that an enzyme and a peptide catalyst can work together in a single reaction vessel to catalyze a two-step cascade reaction with high chemo- and stereoselectivity. Abundant linear alcohols, nitroolefins, an alcohol oxidase, and a tripeptide catalyst provided chiral γ-nitroaldehydes in aqueous buffer. High yields (up to 92 %) and stereoselectivities (up to 98 % ee) were achieved for the cascade through the rational design of the peptide catalyst and the identification of common reaction conditions.

Exploring the Specificity of Nanozymes

Nanozymes, nanomaterials exhibiting enzyme-like activities, have emerged as a prominent interdisciplinary field over the past decade. To date, over 1200 different nanomaterials have been identified as nanozymes, covering four catalytic categories: oxidoreductases, hydrolases, isomerases, and lyases. Catalytic activity and specificity are two pivotal benchmarks for evaluating enzymatic performance. Despite substantial progress being made in quantifying and optimizing the catalytic activity of nanozymes, there is still a lack of in-depth research on the catalytic specificity of nanozymes, preventing the formation of consensual knowledge and impeding a more refined and systematic classification of nanozymes. Recently, debates have emerged regarding whether nanozymes could possess catalytic specificity similar to that of enzymes. This Perspective discusses the specificity of nanozymes by referring to the catalytic specificity of enzymes, highlights the specificity gap between nanozymes and enzymes, and concludes by offering our perspective on future research on the specificity of nanozymes.

Pervasive epistasis exposes intramolecular networks in adaptive enzyme evolution

Enzyme evolution is characterized by constant alterations of the intramolecular residue networks supporting their functions. The rewiring of these network interactions can give rise to epistasis. As mutations accumulate, the epistasis observed across diverse genotypes may appear idiosyncratic, that is, exhibit unique effects in different genetic backgrounds. Here, we unveil a quantitative picture of the prevalence and patterns of epistasis in enzyme evolution by analyzing 41 fitness landscapes generated from seven enzymes. We show that >94% of all mutational and epistatic effects appear highly idiosyncratic, which greatly distorted the functional prediction of the evolved enzymes. By examining seemingly idiosyncratic changes in epistasis along adaptive trajectories, we expose several instances of higher-order, intramolecular rewiring. Using complementary structural data, we outline putative molecular mechanisms explaining higher-order epistasis along two enzyme trajectories. Our work emphasizes the prevalence of epistasis and provides an approach to exploring this phenomenon through a molecular lens.

From random to rational: improving enzyme design through electric fields, second coordination sphere interactions, and conformational dynamics

Enzymes are versatile and efficient biological catalysts that drive numerous cellular processes, motivating the development of enzyme design approaches to tailor catalysts for diverse applications. In this perspective, we investigate the unique properties of natural, evolved, and designed enzymes, recognizing their strengths and shortcomings. We highlight the challenges and limitations of current enzyme design protocols, with a particular focus on their limited consideration of long-range electrostatic and dynamic effects. We then delve deeper into the impact of the protein environment on enzyme catalysis and explore the roles of preorganized electric fields, second coordination sphere interactions, and protein dynamics for enzyme function. Furthermore, we present several case studies illustrating successful enzyme-design efforts incorporating enzyme strategies mentioned above to achieve improved catalytic properties. Finally, we envision the future of enzyme design research, spotlighting the challenges yet to be overcome and the synergy of intrinsic electric fields, second coordination sphere interactions, and conformational dynamics to push the state-of-the-art boundaries.

Conformational Tuning Shapes the Balance between Functional Promiscuity and Specialization in Paralogous Plasmodium Acyl-CoA Binding Proteins

Paralogous proteins confer enhanced fitness to organisms via complex sequence-conformation codes that shape functional divergence, specialization, or promiscuity. Here, we dissect the underlying mechanism of promiscuous binding versus partial subfunctionalization in paralogues by studying structurally identical acyl-CoA binding proteins (ACBPs) from Plasmodium falciparum that serve as promising drug targets due to their high expression during the protozoan proliferative phase. Combining spectroscopic measurements, solution NMR, SPR, and simulations on two of the paralogues, A16 and A749, we show that minor sequence differences shape nearly every local and global conformational feature. A749 displays a broader and heterogeneous native ensemble, weaker thermodynamic coupling and cooperativity, enhanced fluctuations, and a larger binding pocket volume compared to A16. Site-specific tryptophan probes signal a graded reduction in the sampling of substates in the holo form, which is particularly apparent in A749. The paralogues exhibit a spectrum of binding affinities to different acyl-CoAs with A749, the more promiscuous and hence the likely ancestor, binding 1000-fold stronger to lauroyl-CoA under physiological conditions. We thus demonstrate how minor sequence changes modulate the extent of long-range interactions and dynamics, effectively contributing to the molecular evolution of contrasting functional repertoires in paralogues.

Engineering and Expression Strategies for Optimization of L-Asparaginase Development and Production

Genetic engineering for heterologous expression has advanced in recent years. Model systems such as Escherichia coli, Bacillus subtilis and Pichia pastoris are often used as host microorganisms for the enzymatic production of L-asparaginase, an enzyme widely used in the clinic for the treatment of leukemia and in bakeries for the reduction of acrylamide. Newly developed recombinant L-asparaginase (L-ASNase) may have a low affinity for asparagine, reduced catalytic activity, low stability, and increased glutaminase activity or immunogenicity. Some successful commercial preparations of L-ASNase are now available. Therefore, obtaining novel L-ASNases with improved properties suitable for food or clinical applications remains a challenge. The combination of rational design and/or directed evolution and heterologous expression has been used to create enzymes with desired characteristics. Computer design, combined with other methods, could make it possible to generate mutant libraries of novel L-ASNases without costly and time-consuming efforts. In this review, we summarize the strategies and approaches for obtaining and developing L-ASNase with improved properties.

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.

Spinocerebellar ataxia 38: structure-function analysis shows ELOVL5 G230V is proteotoxic, conformationally altered and a mutational hotspot

Fatty acid elongase ELOVL5 is part of a protein family of multipass transmembrane proteins that reside in the endoplasmic reticulum where they regulate long-chain fatty acid elongation. A missense variant (c.689G>T p.Gly230Val) in ELOVL5 causes Spinocerebellar Ataxia subtype 38 (SCA38), a neurodegenerative disorder characterized by autosomal dominant inheritance, cerebellar Purkinje cell demise and adult-onset ataxia. Having previously showed aberrant accumulation of p.G230V in the Golgi complex, here we further investigated the pathogenic mechanisms triggered by p.G230V, integrating functional studies with bioinformatic analyses of protein sequence and structure. Biochemical analysis showed that p.G230V enzymatic activity was normal. In contrast, SCA38-derived fibroblasts showed reduced expression of ELOVL5, Golgi complex enlargement and increased proteasomal degradation with respect to controls. By heterologous overexpression, p.G230V was significantly more active than wild-type ELOVL5 in triggering the unfolded protein response and in decreasing viability in mouse cortical neurons. By homology modelling, we generated native and p.G230V protein structures whose superposition revealed a shift in Loop 6 in p.G230V that altered a highly conserved intramolecular disulphide bond. The conformation of this bond, connecting Loop 2 and Loop 6, appears to be elongase-specific. Alteration of this intramolecular interaction was also observed when comparing wild-type ELOVL4 and the p.W246G variant which causes SCA34. We demonstrate by sequence and structure analyses that ELOVL5 p.G230V and ELOVL4 p.W246G are position-equivalent missense variants. We conclude that SCA38 is a conformational disease and propose combined loss of function by mislocalization and gain of toxic function by ER/Golgi stress as early events in SCA38 pathogenesis.

Biological evolution requires an emergent, self-organizing principle

In this perspective review, we assess fundamental flaws in Darwinian evolution, including its modern versions. Fixed mutations 'explain' microevolution but not macroevolution including speciation events and the origination of all the major body plans of the Cambrian explosion. Complex, multifactorial change is required for speciation events and inevitably requires self-organization beyond what is accomplished by known mechanisms. The assembly of ribosomes and ATP synthase are specific examples. We propose their origin is a model for what is unexplained in biological evolution. Probability of evolution is modeled in Section 9 and values are absurdly improbable. Speciation and higher taxonomic changes become exponentially less probable as the number of required, genetically-based events increase. Also, the power required of the proposed selection mechanism (survival of the fittest) is nil for any biological advance requiring multiple changes, because they regularly occur in multiple generations (different genomes) and would not be selectively conserved by the concept survival of the fittest (a concept ultimately centered on the individual). Thus, survival of the fittest cannot 'explain' the origin of the millions of current and extinct species. We also focus on the inadequacies of laboratory chemistry to explain the complex, required biological self-organization seen in cells. We propose that a 'bioelectromagnetic' field/principle emerges in living cells. Synthesis by self-organization of massive molecular complexes involves biochemical responses to this emergent field/principle. There are ramifications for philosophy, science, and religion. Physics and mathematics must be more strongly integrated with biology and integration should receive dedicated funding with special emphasis for medical applications; treatment of cancer and genetic diseases are examples.

Second-Shell Residues Contribute to Catalysis by Predominately Preorganizing the Apo State in PafA

Residues beyond the first coordination shell are often observed to make considerable cumulative contributions in enzymes. Due to typically indirect perturbations of multiple physicochemical properties of the active site, however, their individual and specific roles in enzyme catalysis and disease-causing mutations remain difficult to predict and understand at the molecular level. Here we analyze the contributions of several second-shell residues in phosphate-irrepressible alkaline phosphatase of flavobacterium (PafA), a representative system as one of the most efficient enzymes. By adopting a multifaceted approach that integrates quantum-mechanical/molecular-mechanical free energy computations, molecular-mechanical molecular dynamics simulations, and density functional theory cluster model calculations, we probe the rate-limiting phosphoryl transfer step and structural properties of all relevant enzyme states. In combination with available experimental data, our computational results show that mutations of the studied second-shell residues impact catalytic efficiency mainly by perturbation of the apo state and therefore substrate binding, while they do not affect the ground state or alter the nature of phosphoryl transfer transition state significantly. Several second-shell mutations also modulate the active site hydration level, which in turn influences the energetics of phosphoryl transfer. These mechanistic insights also help inform strategies that may improve the efficiency of enzyme design and engineering by going beyond the current focus on the first coordination shell.

Leveraging a Structural Blueprint to Rationally Engineer the Rieske Oxygenase TsaM

Rieske nonheme iron oxygenases use two metallocenters, a Rieske-type [2Fe-2S] cluster and a mononuclear iron center, to catalyze oxidation reactions on a broad range of substrates. These enzymes are widely used by microorganisms to degrade environmental pollutants and to build complexity in a myriad of biosynthetic pathways that are industrially interesting. However, despite the value of this chemistry, there is a dearth of understanding regarding the structure-function relationships in this enzyme class, which limits our ability to rationally redesign, optimize, and ultimately exploit the chemistry of these enzymes. Therefore, in this work, by leveraging a combination of available structural information and state-of-the-art protein modeling tools, we show that three "hotspot" regions can be targeted to alter the site selectivity, substrate preference, and substrate scope of the Rieske oxygenase p-toluenesulfonate methyl monooxygenase (TsaM). Through mutation of six to 10 residues distributed between three protein regions, TsaM was engineered to behave as either vanillate monooxygenase (VanA) or dicamba monooxygenase (DdmC). This engineering feat means that TsaM was rationally engineered to catalyze an oxidation reaction at the meta and ortho positions of an aromatic substrate, rather than its favored native para position, and that TsaM was redesigned to perform chemistry on dicamba, a substrate that is not natively accepted by the enzyme. This work thus contributes to unlocking our understanding of structure-function relationships in the Rieske oxygenase enzyme class and expands foundational principles for future engineering of these metalloenzymes.

Improving the catalytic activity of β-glucosidase from Coniophora puteana via semi-rational design for efficient biomass cellulose degradation

In order to improve the degradation activity of β-glucosidase (CpBgl) from Coniophora puteana, the structural modification was conducted. The enzyme activity of mutants CpBgl-Q20C and CpBgl-A240S was increased by 65.75% and 58.58%, respectively. These mutants exhibited maximum activity under the same conditions as wild-type CpBgl (65 ℃ and pH 5.0), slightly improved stabilities compared that of the wild-type, and remarkably enhanced activities in the presence of Mn²⁺ or Fe²⁺. The V_max of CpBgl-Q20C and CpBgl-A240S was increased to 138.18 and 125.14 μmol/mg/min, respectively, from 81.34 μmol/mg/min of the wild-type, and the catalysis efficiency (k_cat/K_m) of CpBgl-Q20C (335.79 min^-1/mM) and CpBgl-A240S (281.51 min^-1/mM) was significantly improved compared with that of the wild-type (149.12 min^-1/mM). When the mutant CpBgl-Q20C were used in the practical degradation of different biomasses, the glucose yields of filter paper, corncob residue, and fungi mycelia residue were increased by 17.68%, 25.10%, and 20.37%, respectively. The spatial locations of the mutation residues in the architecture of CpBgl and their unique roles in the enzyme-substrate binding and catalytic efficiency were probed in this work. These results laid a foundation for evolution of other glycoside hydrolases and the industrial bio-degradation of cellulosic biomass in nature.

Experimental and In Silico Analysis of TEM β-Lactamase Adaptive Evolution

Multiple mutations often have non-additive (epistatic) phenotypic effects. Epistasis is of fundamental biological relevance but is not well understood mechanistically. Adaptive evolution, i.e., the evolution of new biochemical activities, is rich in epistatic interactions. To better understand the principles underlying epistasis during genetic adaptation, we studied the evolution of TEM-1 β-lactamase variants exhibiting cefotaxime resistance. We report the collection of a library of 487 observed evolutionary trajectories for TEM-1 and determine the epistasis status based on cefotaxime resistance phenotype for 206 combinations of 2-3 TEM-1 mutations involving 17 positions under adaptive selective pressure. Gain-of-function (GOF) mutations are gatekeepers for adaptation. To see if GOF phenotypes can be inferred based solely on sequence data, we calculated the enrichment of GOF mutations in the different categories of epistatic pairs. Our results suggest that this is possible because GOF mutations are particularly enriched in sign and reciprocal sign epistasis, which leave a major imprint on the sequence space accessible to evolution. We also used FoldX to explore the relationship between thermodynamic stability and epistasis. We found that mutations in observed evolutionary trajectories tend to destabilize the folded structure of the protein, albeit their cumulative effects are consistently below the protein's free energy of folding. The destabilizing effect is stronger for epistatic pairs, suggesting that modest or local alterations in folding stability can modulate catalysis. Finally, we report a significant relationship between epistasis and the degree to which two protein positions are structurally and dynamically coupled, even in the absence of ligand.

Synthetic Biology: Bottom-Up Assembly of Molecular Systems

The bottom-up assembly of biological and chemical components opens exciting opportunities to engineer artificial vesicular systems for applications with previously unmet requirements. The modular combination of scaffolds and functional building blocks enables the engineering of complex systems with biomimetic or new-to-nature functionalities. Inspired by the compartmentalized organization of cells and organelles, lipid or polymer vesicles are widely used as model membrane systems to investigate the translocation of solutes and the transduction of signals by membrane proteins. The bottom-up assembly and functionalization of such artificial compartments enables full control over their composition and can thus provide specifically optimized environments for synthetic biological processes. This review aims to inspire future endeavors by providing a diverse toolbox of molecular modules, engineering methodologies, and different approaches to assemble artificial vesicular systems. Important technical and practical aspects are addressed and selected applications are presented, highlighting particular achievements and limitations of the bottom-up approach. Complementing the cutting-edge technological achievements, fundamental aspects are also discussed to cater to the inherently diverse background of the target audience, which results from the interdisciplinary nature of synthetic biology. The engineering of proteins as functional modules and the use of lipids and block copolymers as scaffold modules for the assembly of functionalized vesicular systems are explored in detail. Particular emphasis is placed on ensuring the controlled assembly of these components into increasingly complex vesicular systems. Finally, all descriptions are presented in the greater context of engineering valuable synthetic biological systems for applications in biocatalysis, biosensing, bioremediation, or targeted drug delivery.

How gene duplication diversifies the landscape of protein oligomeric state and function

Oligomeric proteins are central to cellular life and the duplication and divergence of their genes is a key driver of evolutionary innovations. The duplication of a gene coding for an oligomeric protein has numerous possible outcomes, which motivates questions on the relationship between structural and functional divergence. How do protein oligomeric states diversify after gene duplication? In the simple case of duplication of a homo-oligomeric protein gene, what properties can influence the fate of descendant paralogs toward forming independent homomers or maintaining their interaction as a complex? Furthermore, how are functional innovations associated with the diversification of oligomeric states? Here, we review recent literature and present specific examples in an attempt to illustrate and answer these questions.

The sphingolipid anteome: implications for evolution of the sphingolipid metabolic pathway

Modern cell membranes contain a bewildering complexity of lipids, among them sphingolipids (SLs). Advances in mass spectrometry have led to the realization that the number and combinatorial complexity of lipids, including SLs, is much greater than previously appreciated. SLs are generated de novo by four enzymes, namely serine palmitoyltransferase, 3-ketodihydrosphingosine reductase, ceramide synthase and dihydroceramide Δ4-desaturase 1. Some of these enzymes depend on the availability of specific substrates and cofactors, which are themselves supplied by other complex metabolic pathways. The evolution of these four enzymes is poorly understood and likely depends on the co-evolution of the metabolic pathways that supply the other essential reaction components. Here, we introduce the concept of the 'anteome', from the Latin ante ('before') to describe the network of metabolic ('omic') pathways that must have converged in order for these pathways to co-evolve and permit SL synthesis. We also suggest that current origin of life and evolutionary models lack appropriate experimental support to explain the appearance of this complex metabolic pathway and its anteome.

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.

Resurrecting Enzymes by Ancestral Sequence Reconstruction

Ancestral Sequence Reconstruction (ASR) allows one to infer the sequences of extinct proteins using the phylogeny of extant proteins. It consists of disclosing the evolutionary history-i.e., the phylogeny-of a protein family of interest and then inferring the sequences of its ancestors-i.e., the nodes in the phylogeny. Assisted by gene synthesis, the selected ancestors can be resurrected in the lab and experimentally characterized. The crucial step to succeed with ASR is starting from a reliable phylogeny. At the same time, it is of the utmost importance to have a clear idea on the evolutionary history of the family under study and the events that influenced it. This allows us to implement ASR with well-defined hypotheses and to apply the appropriate experimental methods. In the last years, ASR has become popular to test hypotheses about the origin of functionalities, changes in activities, understanding physicochemical properties of proteins, among others. In this context, the aim of this chapter is to present the ASR approach applied to the reconstruction of enzymes-i.e., proteins with catalytic roles. The spirit of this contribution is to provide a basic, hands-to-work guide for biochemists and biologists who are unfamiliar with molecular phylogenetics.

Learning Strategies in Protein Directed Evolution

Synthetic biology is a fast-evolving research field that combines biology and engineering principles to develop new biological systems for medical, pharmacological, and industrial applications. Synthetic biologists use iterative "design, build, test, and learn" cycles to efficiently engineer genetic systems that are reliable, reproducible, and predictable. Protein engineering by directed evolution can benefit from such a systematic engineering approach for various reasons. Learning can be carried out before starting, throughout or after finalizing a directed evolution project. Computational tools, bioinformatics, and scanning mutagenesis methods can be excellent starting points, while molecular dynamics simulations and other strategies can guide engineering efforts. Similarly, studying protein intermediates along evolutionary pathways offers fascinating insights into the molecular mechanisms shaped by evolution. The learning step of the cycle is not only crucial for proteins or enzymes that are not suitable for high-throughput screening or selection systems, but it is also valuable for any platform that can generate a large amount of data that can be aided by machine learning algorithms. The main challenge in protein engineering is to predict the effect of a single mutation on one functional parameter-to say nothing of several mutations on multiple parameters. This is largely due to nonadditive mutational interactions, known as epistatic effects-beneficial mutations present in a genetic background may not be beneficial in another genetic background. In this work, we provide an overview of experimental and computational strategies that can guide the user to learn protein function at different stages in a directed evolution project. We also discuss how epistatic effects can influence the success of directed evolution projects. Since machine learning is gaining momentum in protein engineering and the field is becoming more interdisciplinary thanks to collaboration between mathematicians, computational scientists, engineers, molecular biologists, and chemists, we provide a general workflow that familiarizes nonexperts with the basic concepts, dataset requirements, learning approaches, model capabilities and performance metrics of this intriguing area. Finally, we also provide some practical recommendations on how machine learning can harness epistatic effects for engineering proteins in an "outside-the-box" way.

Integrative Role of Albumin: Evolutionary, Biochemical and Pathophysiological Aspects

Being one of the main proteins in the human body and many animal species, albumin plays a crucial role in the transport of various ions, electrically neutral molecules and in maintaining the colloidal osmotic pressure of the blood. Albumin is able to bind almost all known drugs, many nutraceuticals and toxic substances, determining their pharmaco- and toxicokinetics. However, albumin is not only the passive but also the active participant of the pharmacokinetic and toxicokinetic processes possessing a number of enzymatic activities. Due to the thiol group of Cys34, albumin can serve as a trap for reactive oxygen and nitrogen species, thus participating in redox processes. The interaction of the protein with blood cells, blood vessels, and also with tissue cells outside the vascular bed is of great importance. The interaction of albumin with endothelial glycocalyx and vascular endothelial cells largely determines its integrative role. This review provides information of a historical nature, information on evolutionary changes, inflammatory and antioxidant properties of albumin, on its structural and functional modifications and their significance in the pathogenesis of some diseases.

Identification of New Specificity Determinants in Bacterial Purine Nucleobase Transporters based on an Ancestral Sequence Reconstruction Approach

The relation of sequence with specificity in membrane transporters is challenging to explore. Most relevant studies until now rely on comparisons of present-day homologs. In this work, we study a set of closely related transporters by employing an evolutionary, ancestral-reconstruction approach and reveal unexpected new specificity determinants. We analyze a monophyletic group represented by the xanthine-specific XanQ of Escherichia coli in the Nucleobase-Ascorbate Transporter/Nucleobase-Cation Symporter-2 (NAT/NCS2) family. We reconstructed AncXanQ, the putative common ancestor of this clade, expressed it in E. coli K-12, and found that, in contrast to XanQ, it encodes a high-affinity permease for both xanthine and guanine, which also recognizes adenine, hypoxanthine, and a range of analogs. AncXanQ conserves all binding-site residues of XanQ and differs substantially in only five intramembrane residues outside the binding site. We subjected both homologs to rationally designed mutagenesis and present evidence that these five residues are linked with the specificity change. In particular, we reveal Ser377 of XanQ (Gly in AncXanQ) as a major determinant. Replacement of this Ser with Gly enlarges the specificity of XanQ towards an AncXanQ-phenotype. The ortholog from Neisseria meningitidis retaining Gly at this position is also a xanthine/guanine transporter with extended substrate profile like AncXanQ. Molecular Dynamics shows that the S377G replacement tilts transmembrane helix 12 resulting in rearrangement of Phe376 relative to Phe94 in the XanQ binding pocket. This effect may rationalize the enlarged specificity. On the other hand, the specificity effect of S377G can be masked by G27S or other mutations through epistatic interactions.

A shared mechanistic pathway for pyridoxal phosphate-dependent arginine oxidases

The mechanism by which molecular oxygen is activated by the organic cofactor pyridoxal phosphate (PLP) for oxidation reactions remains poorly understood. Recent work has identified arginine oxidases that catalyze desaturation or hydroxylation reactions. Here, we investigate a desaturase from the Pseudoalteromonas luteoviolacea indolmycin pathway. Our work, combining X-ray crystallographic, biochemical, spectroscopic, and computational studies, supports a shared mechanism with arginine hydroxylases, involving two rounds of single-electron transfer to oxygen and superoxide rebound at the 4' carbon of the PLP cofactor. The precise positioning of a water molecule in the active site is proposed to control the final reaction outcome. This proposed mechanism provides a unified framework to understand how oxygen can be activated by PLP-dependent enzymes for oxidation of arginine and elucidates a shared mechanistic pathway and intertwined evolutionary history for arginine desaturases and hydroxylases.

Rate-Perturbing Single Amino Acid Mutation for Hydrolases: A Statistical Profiling

Hydrolases are a critical component for modern chemical, pharmaceutical, and environmental sciences. Identifying mutations that enhance catalytic efficiency presents a roadblock to design and to discover new hydrolases for broad academic and industrial uses. Here, we report the statistical profiling for rate-perturbing mutant hydrolases with a single amino acid substitution. We constructed an integrated structure-kinetics database for hydrolases, IntEnzyDB, which contains 3907 k_cats, 4175 K_Ms, and 2715 Protein Data Bank IDs. IntEnzyDB adopts a relational architecture with a flattened data structure, enabling facile and efficient access to clean and tabulated data for machine learning uses. We conducted statistical analyses on how single amino acids mutations influence the turnover number (i.e., k_cat) and efficiency (i.e., k_cat/K_M), with a particular emphasis on profiling the features for rate-enhancing mutations. The results show that mutation to bulky nonpolar residues with a hydrocarbon chain involves a higher likelihood for rate acceleration than to other types of residues. Linear regression models reveal geometric descriptors of substrate and mutation residues that mediate rate-perturbing outcomes for hydrolases with bulky nonpolar mutations. On the basis of the analyses of the structure-kinetics relationship, we observe that the propensity for rate enhancement is independent of protein sizes. In addition, we observe that distal mutations (i.e., >10 Å from the active site) in hydrolases are significantly more prone to induce efficiency neutrality and avoid efficiency deletion but involve similar propensity for rate enhancement. The studies reveal the statistical features for identifying rate-enhancing mutations in hydrolases, which will potentially guide hydrolase discovery in biocatalysis.

Serum Albumin in Health and Disease: Esterase, Antioxidant, Transporting and Signaling Properties

Being one of the main proteins in the human body and many animal species, albumin plays a decisive role in the transport of various ions-electrically neutral and charged molecules-and in maintaining the colloidal osmotic pressure of the blood. Albumin is able to bind to almost all known drugs, as well as many nutraceuticals and toxic substances, largely determining their pharmaco- and toxicokinetics. Albumin of humans and respective representatives in cattle and rodents have their own structural features that determine species differences in functional properties. However, albumin is not only passive, but also an active participant of pharmacokinetic and toxicokinetic processes, possessing a number of enzymatic activities. Numerous experiments have shown esterase or pseudoesterase activity of albumin towards a number of endogeneous and exogeneous esters. Due to the free thiol group of Cys34, albumin can serve as a trap for reactive oxygen and nitrogen species, thus participating in redox processes. Glycated albumin makes a significant contribution to the pathogenesis of diabetes and other diseases. The interaction of albumin with blood cells, blood vessels and tissue cells outside the vascular bed is of great importance. Interactions with endothelial glycocalyx and vascular endothelial cells largely determine the integrative role of albumin. This review considers the esterase, antioxidant, transporting and signaling properties of albumin, as well as its structural and functional modifications and their significance in the pathogenesis of certain diseases.

Combined Theoretical and Experimental Study to Unravel the Differences in Promiscuous Amidase Activity of Two Nonhomologous Enzymes

Convergent evolution has resulted in nonhomologous enzymes that contain similar active sites that catalyze the same primary and secondary reactions. Comparing how these enzymes achieve their reaction promiscuity can yield valuable insights to develop functions from the optimization of latent activities. In this work, we have focused on the promiscuous amidase activity in the esterase from Bacillus subtilis (Bs2) and compared with the same activity in the promiscuous lipase B from Candida antarctica (CALB). The study, combining multiscale quantum mechanics/molecular mechanics (QM/MM) simulations, deep machine learning approaches, and experimental characterization of Bs2 kinetics, confirms the amidase activity of Bs2 and CALB. The computational results indicate that both enzymes offer a slightly different reaction environment reflected by electrostatic effects within the active site, thus resulting in a different reaction mechanism during the acylation step. A convolutional neural network (CNN) has been used to understand the conserved amino acids among the evolved protein family and suggest that Bs2 provides a more robust protein scaffold to perform future mutagenesis studies. Results derived from this work will help reveal the origin of enzyme promiscuity, which will find applications in enzyme (re)design, particularly in creating a highly active amidase.

Evolution of β-lactamases and enzyme promiscuity

β-Lactamases represent one of the most prevalent resistance mechanisms against β-lactam antibiotics. Beyond their clinical importance, they have also become key models in enzymology and evolutionary biochemistry. A global understanding of their evolution and sequence and functional diversity can therefore aid a wide set of different disciplines. Interestingly, β-lactamases have evolved multiple times from distinct evolutionary origins, with ancestries that reach back billions of years. It is therefore no surprise that these enzymes exhibit diverse structural features and enzymatic mechanisms. In this review, we provide a bird's eye view on the evolution of β-lactamases within the two enzyme superfamilies-i.e. the penicillin-binding protein-like and metallo-β-lactamase superfamily-through phylogenetics. We further discuss potential evolutionary origins of each β-lactamase class by highlighting signs of evolutionary connections in protein functions between β-lactamases and other enzymes, especially cases of enzyme promiscuity.

Epistasis and intramolecular networks in protein evolution

Proteins are molecular machines composed of complex, highly connected amino acid networks. Their functional optimization requires the reorganization of these intramolecular networks by evolution. In this review, we discuss the mechanisms by which epistasis, that is, the dependence of the effect of a mutation on the genetic background, rewires intramolecular interactions to alter protein function. Deciphering the biophysical basis of epistasis is crucial to our understanding of evolutionary dynamics and the elucidation of sequence-structure-function relationships. We featured recent studies that provide insights into the molecular mechanisms giving rise to epistasis, particularly at the structural level. These studies illustrate the convoluted and fascinating nature of the intramolecular networks co-opted by epistasis during the evolution of protein function.

Stability of Ligand-induced Protein Conformation Influences Affinity in Maltose-binding Protein

Our understanding of what determines ligand affinity of proteins is poor, even with high-resolution structures available. Both the non-covalent ligand-protein interactions and the relative free energies of available conformations contribute to the affinity of a protein for a ligand. Distant, non-binding site residues can influence the ligand affinity by altering the free energy difference between a ligand-free and ligand-bound conformation. Our hypothesis is that when different ligands induce distinct ligand-bound conformations, it should be possible to tweak their affinities by changing the free energies of the available conformations. We tested this idea for the maltose-binding protein (MBP) from Escherichia coli. We used single-molecule Förster resonance energy transfer (smFRET) to distinguish several unique ligand-bound conformations of MBP. We engineered mutations, distant from the binding site, to affect the stabilities of different ligand-bound conformations. We show that ligand affinity can indeed be altered in a conformation-dependent manner. Our studies provide a framework for the tuning of ligand affinity, apart from modifying binding site residues.

Cryptic β-Lactamase Evolution Is Driven by Low β-Lactam Concentrations

Very low antibiotic concentrations have been shown to drive the evolution of antimicrobial resistance. While substantial progress has been made to understand the driving role of low concentrations during resistance development for different antimicrobial classes, the importance of β-lactams, the most commonly used antibiotics, is still poorly studied.

Our current understanding of how low antibiotic concentrations shape the evolution of contemporary β-lactamases is limited. Using the widespread carbapenemase OXA-48, we tested the long-standing hypothesis that selective compartments with low antibiotic concentrations cause standing genetic diversity that could act as a gateway to developing clinical resistance. Here, we subjected Escherichia coli expressing bla_OXA-48, on a clinical plasmid, to experimental evolution at sub-MICs of ceftazidime. We identified and characterized seven single variants of OXA-48. Susceptibility profiles and dose-response curves showed that they increased resistance only marginally. However, in competition experiments at sub-MICs of ceftazidime, they demonstrated strong selectable fitness benefits. Increased resistance was also reflected in elevated catalytic efficiencies toward ceftazidime. These changes are likely caused by enhanced flexibility of the Ω- and β5-β6 loops and fine-tuning of preexisting active site residues. In conclusion, low-level concentrations of β-lactams can drive the evolution of β-lactamases through cryptic phenotypes which may act as stepping-stones toward clinical resistance.

IMPORTANCE Very low antibiotic concentrations have been shown to drive the evolution of antimicrobial resistance. While substantial progress has been made to understand the driving role of low concentrations during resistance development for different antimicrobial classes, the importance of β-lactams, the most commonly used antibiotics, is still poorly studied. Here, we shed light on the evolutionary impact of low β-lactam concentrations on the widespread β-lactamase OXA-48. Our data indicate that the exposure to β-lactams at very low concentrations enhances β-lactamase diversity and drives the evolution of β-lactamases by significantly influencing their substrate specificity. Thus, in contrast to high concentrations, low levels of these drugs may substantially contribute to the diversification and divergent evolution of these enzymes, providing a standing genetic diversity that can be selected and mobilized when antibiotic pressure increases.

Protein engineering for natural product biosynthesis and synthetic biology applications

As protein engineering grows more salient, many strategies have emerged to alter protein structure and function, with the goal of redesigning and optimizing natural product biosynthesis. Computational tools, including machine learning and molecular dynamics simulations, have enabled the rational mutagenesis of key catalytic residues for enhanced or altered biocatalysis. Semi-rational, directed evolution and microenvironment engineering strategies have optimized catalysis for native substrates and increased enzyme promiscuity beyond the scope of traditional rational approaches. These advances are made possible using novel high-throughput screens, including designer protein-based biosensors with engineered ligand specificity. Herein, we detail the most recent of these advances, focusing on polyketides, non-ribosomal peptides and isoprenoids, including their native biosynthetic logic to provide clarity for future applications of these technologies for natural product synthetic biology.

Designing better enzymes: Insights from directed evolution

De novo enzymes can be created by computational design and directed evolution. Here, we review recent insights into the origins of catalytic power in evolved designer enzymes to pinpoint opportunities for next-generation designs: Evolution precisely organizes active sites, introduces catalytic H-bonding networks, invokes electrostatic catalysis, and creates dynamical networks embedding the active site in a reactive protein scaffold. Such insights foster our fundamental knowledge of enzyme catalysis and fuel the future design of tailor-made enzymes.

Exploring mutable conserved sites and fatal non-conserved sites by random mutation of esterase from Sulfolobus tokodaii and subtilisin from Thermococcus kodakarensis

Homologous proteins differ in their amino acid sequences at several positions. Generally, conserved sites are recognized as not suitable for amino acid substitution, and thus in evolutionary protein engineering, non-conserved sites are often selected as mutation sites. However, there have also been reports of possible mutations in conserved sites. In this study, we explored mutable conserved sites and immutable non-conserved sites by testing random mutations of two thermostable proteins, an esterase from Sulfolobus tokodaii (Sto-Est) and a subtilisin from Thermococcus kodakarensis (Tko-Sub). The subtilisin domain of Tko-Sub needs Ca²⁺ ions and the propeptide domain for stability, folding and maturation. The results from the two proteins showed that about one-third of the mutable sites were detected in conserved sites and some non-conserved sites lost enzymatic activity at high temperatures due to mutation. Of the conserved sites in Sto-Est, the sites on the loop, on the surface, and far from the active site are more resistant to mutation. In Tko-Sub, the sites flanking Ca²⁺-binding sites and propeptide were undesirable for mutation. The results presented here serve as an index for selecting mutation sites and contribute to the expansion of available sequence range by introducing mutations at conserved sites.

Known Evolutionary Paths Are Accessible to Engineered ß-Lactamases Having Altered Protein Motions at the Timescale of Catalytic Turnover

The evolution of new protein functions is dependent upon inherent biophysical features of proteins. Whereas, it has been shown that changes in protein dynamics can occur in the course of directed molecular evolution trajectories and contribute to new function, it is not known whether varying protein dynamics modify the course of evolution. We investigate this question using three related ß-lactamases displaying dynamics that differ broadly at the slow timescale that corresponds to catalytic turnover yet have similar fast dynamics, thermal stability, catalytic, and substrate recognition profiles. Introduction of substitutions E104K and G238S, that are known to have a synergistic effect on function in the parent ß-lactamase, showed similar increases in catalytic efficiency toward cefotaxime in the related ß-lactamases. Molecular simulations using Protein Energy Landscape Exploration reveal that this results from stabilizing the catalytically-productive conformations, demonstrating the dominance of the synergistic effect of the E014K and G238S substitutions in vitro in contexts that vary in terms of sequence and dynamics. Furthermore, three rounds of directed molecular evolution demonstrated that known cefotaximase-enhancing mutations were accessible regardless of the differences in dynamics. Interestingly, specific sequence differences between the related ß-lactamases were shown to have a higher effect in evolutionary outcomes than did differences in dynamics. Overall, these ß-lactamase models show tolerance to protein dynamics at the timescale of catalytic turnover in the evolution of a new function.

A mechanistic view of enzyme evolution

New enzyme functions often evolve through the recruitment and optimization of latent promiscuous activities. How do mutations alter the molecular architecture of enzymes to enhance their activities? Can we infer general mechanisms that are common to most enzymes, or does each enzyme require a unique optimization process? The ability to predict the location and type of mutations necessary to enhance an enzyme's activity is critical to protein engineering and rational design. In this review, via the detailed examination of recent studies that have shed new light on the molecular changes underlying the optimization of enzyme function, we provide a mechanistic perspective of enzyme evolution. We first present a global survey of the prevalence of activity-enhancing mutations and their distribution within protein structures. We then delve into the molecular solutions that mediate functional optimization, specifically highlighting several common mechanisms that have been observed across multiple examples. As distinct protein sequences encounter different evolutionary bottlenecks, different mechanisms are likely to emerge along evolutionary trajectories toward improved function. Identifying the specific mechanism(s) that need to be improved upon, and tailoring our engineering efforts to each sequence, may considerably improve our chances to succeed in generating highly efficient catalysts in the future.