How protein stability and new functions trade off

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.

Mutation-driven evolution of antibacterial function in an ancestral antifungal scaffold: Significance for peptide engineering

Mutation-driven evolution of novel function on an old gene has been documented in many development- and adaptive immunity-related genes but is poorly understood in immune effector molecules. Drosomycin-type antifungal peptides (DTAFPs) are a family of defensin-type effectors found in plants and ecdysozoans. Their primitive function was to control fungal infection and then co-opted for fighting against bacterial infection in plants, insects, and nematodes. This provides a model to study the structural and evolutionary mechanisms behind such functional diversification. In the present study, we determined the solution structure of mehamycin, a DTAFP from the Northern root-knot nematode Meloidogyne hapla with antibacterial activity and an 18-mer insert, and studied the mutational effect through using a mutant with the insert deleted. Mehamycin adopts an expected cysteine-stabilized α-helix and β-sheet fold in its core scaffold and the inserted region, called single Disulfide Bridge-linked Domain (abbreviated as sDBD), forms an extended loop protruding from the scaffold. The latter folds into an amphipathic architecture stabilized by one disulfide bridge, which likely confers mehamycin a bacterial membrane permeability. Deletion of the sDBD remarkably decreased the ability but accompanying an increase in thermostability, indicative of a structure-function trade-off in the mehamycin evolution. Allosteric analysis revealed an interior interaction between the two domains, which might promote point mutations at some key sites of the core domain and ultimately give rise to the emergence of antibacterial function. Our work may be valuable in guiding protein engineering of mehamycin to improve its activity and stability.

Bacterial Hsp90 predominantly buffers but does not potentiate the phenotypic effects of deleterious mutations during fluorescent protein evolution

Chaperones facilitate the folding of other ("client") proteins and can thus affect the adaptive evolution of these clients. Specifically, chaperones affect the phenotype of proteins via two opposing mechanisms. On the one hand, they can buffer the effects of mutations in proteins and thus help preserve an ancestral, premutation phenotype. On the other hand, they can potentiate the effects of mutations and thus enhance the phenotypic changes caused by a mutation. We study that how the bacterial Hsp90 chaperone (HtpG) affects the evolution of green fluorescent protein. To this end, we performed directed evolution of green fluorescent protein under low and high cellular concentrations of Hsp90. Specifically, we evolved green fluorescent protein under both stabilizing selection for its ancestral (green) phenotype and directional selection toward a new (cyan) phenotype. While Hsp90 did only affect the rate of adaptive evolution transiently, it did affect the phenotypic effects of mutations that occurred during adaptive evolution. Specifically, Hsp90 allowed strongly deleterious mutations to accumulate in evolving populations by buffering their effects. Our observations show that the role of a chaperone for adaptive evolution depends on the organism and the trait being studied.

Optimizing Antibody Affinity and Developability Using a Framework-CDR Shuffling Approach-Application to an Anti-SARS-CoV-2 Antibody

The computational methods used for engineering antibodies for clinical development have undergone a transformation from three-dimensional structure-guided approaches to artificial-intelligence- and machine-learning-based approaches that leverage the large sequence data space of hundreds of millions of antibodies generated by next-generation sequencing (NGS) studies. Building on the wealth of available sequence data, we implemented a computational shuffling approach to antibody components, using the complementarity-determining region (CDR) and the framework region (FWR) to optimize an antibody for improved affinity and developability. This approach uses a set of rules to suitably combine the CDRs and FWRs derived from naturally occurring antibody sequences to engineer an antibody with high affinity and specificity. To illustrate this approach, we selected a representative SARS-CoV-2-neutralizing antibody, H4, which was identified and isolated previously based on the predominant germlines that were employed in a human host to target the SARS-CoV-2-human ACE2 receptor interaction. Compared to screening vast CDR libraries for affinity enhancements, our approach identified fewer than 100 antibody framework-CDR combinations, from which we screened and selected an antibody (CB79) that showed a reduced dissociation rate and improved affinity against the SARS-CoV-2 spike protein (7-fold) when compared to H4. The improved affinity also translated into improved neutralization (>75-fold improvement) of SARS-CoV-2. Our rapid and robust approach for optimizing antibodies from parts without the need for tedious structure-guided CDR optimization will have broad utility for biotechnological applications.

Directed Evolution of Soluble α-1,2-Fucosyltransferase Using Kanamycin Resistance Protein as a Phenotypic Reporter for Efficient Production of 2'-Fucosyllactose

2'-Fucosyllactose (2'-FL), the most abundant fucosylated oligosaccharide in human milk, has multiple beneficial effects on human health. However, its biosynthesis by metabolically engineered Escherichia coli is often hampered owing to the insolubility and instability of α-1,2-fucosyltransferase (the rate-limiting enzyme). In this study, we aimed to enhance 2'-FL production by increasing the expression of soluble α-1,2-fucosyltransferase from Helicobacter pylori (FucT2). Because structural information regarding FucT2 has not been unveiled, we decided to improve the expression of soluble FucT2 in E. coli via directed evolution using a protein solubility biosensor that links protein solubility to antimicrobial resistance. For such a system to be viable, the activity of kanamycin resistance protein (Kan^R) should be dependent on FucT2 solubility. Kan^R was fused to the C-terminus of mutant libraries of FucT2, which were generated using a combination of error-prone PCR and DNA shuffling. Notably, one round of the directed evolution process, which consisted of mutant library generation and selection based on kanamycin resistance, resulted in a significant increase in the expression level of soluble FucT2. As a result, a batch fermentation with the ΔL M15 pBCGW strain, expressing the FucT2 mutant (F#1-5) isolated from the first round of the directed evolution process, resulted in the production of 0.31 g/l 2'-FL with a yield of 0.22 g 2'-FL/g lactose, showing 1.72- and 1.51-fold increase in the titer and yield, respectively, compared to those of the control strain. The simple and powerful method developed in this study could be applied to enhance the solubility of other unstable enzymes.

Orthogonal glycolytic pathway enables directed evolution of noncanonical cofactor oxidase

Noncanonical cofactor biomimetics (NCBs) such as nicotinamide mononucleotide (NMN+) provide enhanced scalability for biomanufacturing. However, engineering enzymes to accept NCBs is difficult. Here, we establish a growth selection platform to evolve enzymes to utilize NMN+-based reducing power. This is based on an orthogonal, NMN+-dependent glycolytic pathway in Escherichia coli which can be coupled to any reciprocal enzyme to recycle the ensuing reduced NMN+. With a throughput of >106 variants per iteration, the growth selection discovers a Lactobacillus pentosus NADH oxidase variant with ~10-fold increase in NMNH catalytic efficiency and enhanced activity for other NCBs. Molecular modeling and experimental validation suggest that instead of directly contacting NCBs, the mutations optimize the enzyme's global conformational dynamics to resemble the WT with the native cofactor bound. Restoring the enzyme's access to catalytically competent conformation states via deep navigation of protein sequence space with high-throughput evolution provides a universal route to engineer NCB-dependent enzymes.

Design and optimization of enzymatic activity in a de novo β-barrel scaffold

While native scaffolds offer a large diversity of shapes and topologies for enzyme engineering, their often unpredictable behavior in response to sequence modification makes de novo generated scaffolds an exciting alternative. Here we explore the customization of the backbone and sequence of a de novo designed eight stranded β-barrel protein to create catalysts for a retro-aldolase model reaction. We show that active and specific catalysts can be designed in this fold and use directed evolution to further optimize activity and stereoselectivity. Our results support previous suggestions that different folds have different inherent amenability to evolution and this property could account, in part, for the distribution of natural enzymes among different folds.

Stability and expression of SARS-CoV-2 spike-protein mutations

Protein fold stability likely plays a role in SARS-CoV-2 S-protein evolution, together with ACE2 binding and antibody evasion. While few thermodynamic stability data are available for S-protein mutants, many systematic experimental data exist for their expression. In this paper, we explore whether such expression levels relate to the thermodynamic stability of the mutants. We studied mutation-induced SARS-CoV-2 S-protein fold stability, as computed by three very distinct methods and eight different protein structures to account for method- and structure-dependencies. For all methods and structures used (24 comparisons), computed stability changes correlate significantly (99% confidence level) with experimental yeast expression from the literature, such that higher expression is associated with relatively higher fold stability. Also significant, albeit weaker, correlations were seen between stability and ACE2 binding effects. The effect of thermodynamic fold stability may be direct or a correlate of amino acid or site properties, notably the solvent exposure of the site. Correlation between computed stability and experimental expression and ACE2 binding suggests that functional properties of the SARS-CoV-2 S-protein mutant space are largely determined by a few simple features, due to underlying correlations. Our study lends promise to the development of computational tools that may ideally aid in understanding and predicting SARS-CoV-2 S-protein evolution.

Evolution avoids a pathological stabilizing interaction in the immune protein S100A9

Stability constrains evolution. While much is known about constraints on destabilizing mutations, less is known about the constraints on stabilizing mutations. We recently identified a mutation in the innate immune protein S100A9 that provides insight into such constraints. When introduced into human S100A9, M63F simultaneously increases the stability of the protein and disrupts its natural ability to activate Toll-like receptor 4. Using chemical denaturation, we found that M63F stabilizes a calcium-bound conformation of hS100A9. We then used NMR to solve the structure of the mutant protein, revealing that the mutation distorts the hydrophobic binding surface of hS100A9, explaining its deleterious effect on function. Hydrogen-deuterium exchange (HDX) experiments revealed stabilization of the region around M63F in the structure, notably Phe37. In the structure of the M63F mutant, the Phe37 and Phe63 sidechains are in contact, plausibly forming an edge-face π-stack. Mutating Phe37 to Leu abolished the stabilizing effect of M63F as probed by both chemical denaturation and HDX. It also restored the biological activity of S100A9 disrupted by M63F. These findings reveal that Phe63 creates a molecular staple with Phe37 that stabilizes a nonfunctional conformation of the protein, thus disrupting function. Using a bioinformatic analysis, we found that S100A9 proteins from different organisms rarely have Phe at both positions 37 and 63, suggesting that avoiding a pathological stabilizing interaction indeed constrains S100A9 evolution. This work highlights an important evolutionary constraint on stabilizing mutations, namely, that they must avoid inappropriately stabilizing nonfunctional protein conformations.

De Novo Design of a Highly Stable Ovoid TIM Barrel: Unlocking Pocket Shape towards Functional Design

The ability to finely control the structure of protein folds is an important prerequisite to functional protein design. The TIM barrel fold is an important target for these efforts as it is highly enriched for diverse functions in nature. Although a TIM barrel protein has been designed de novo, the ability to finely alter the curvature of the central beta barrel and the overall architecture of the fold remains elusive, limiting its utility for functional design. Here, we report the de novo design of a TIM barrel with ovoid (twofold) symmetry, drawing inspiration from natural beta and TIM barrels with ovoid curvature. We use an autoregressive backbone sampling strategy to implement our hypothesis for elongated barrel curvature, followed by an iterative enrichment sequence design protocol to obtain sequences which yield a high proportion of successfully folding designs. Designed sequences are highly stable and fold to the designed barrel curvature as determined by a 2.1 Å resolution crystal structure. The designs show robustness to drastic mutations, retaining high melting temperatures even when multiple charged residues are buried in the hydrophobic core or when the hydrophobic core is ablated to alanine. As a scaffold with a greater capacity for hosting diverse hydrogen bonding networks and installation of binding pockets or active sites, the ovoid TIM barrel represents a major step towards the de novo design of functional TIM barrels.

The Statistical Trends of Protein Evolution: A Lesson from AlphaFold Database

The recent development of artificial intelligence provides us with new and powerful tools for studying the mysterious relationship between organism evolution and protein evolution. In this work, based on the AlphaFold Protein Structure Database (AlphaFold DB), we perform comparative analyses of the proteins of different organisms. The statistics of AlphaFold-predicted structures show that, for organisms with higher complexity, their constituent proteins will have larger radii of gyration, higher coil fractions, and slower vibrations, statistically. By conducting normal mode analysis and scaling analyses, we demonstrate that higher organismal complexity correlates with lower fractal dimensions in both the structure and dynamics of the constituent proteins, suggesting that higher functional specialization is associated with higher organismal complexity. We also uncover the topology and sequence bases of these correlations. As the organismal complexity increases, the residue contact networks of the constituent proteins will be more assortative, and these proteins will have a higher degree of hydrophilic-hydrophobic segregation in the sequences. Furthermore, by comparing the statistical structural proximity across the proteomes with the phylogenetic tree of homologous proteins, we show that, statistical structural proximity across the proteomes may indirectly reflect the phylogenetic proximity, indicating a statistical trend of protein evolution in parallel with organism evolution. This study provides new insights into how the diversity in the functionality of proteins increases and how the dimensionality of the manifold of protein dynamics reduces during evolution, contributing to the understanding of the origin and evolution of lives.

L’impact des mutations neutres sur l’évolvabilité et l’évolution des génomes

Les mutations bénéfiques à forts effets sont rares et les mutations délétères sont éliminées par la sélection naturelle. La majorité des mutations qui s’accumulent dans les génomes ont donc des effets sélectifs très faibles, voire nuls ; elles sont alors appelées mutations neutres. Au cours des deux dernières décennies, il a été montré que les mutations, même en l’absence d’effet sur la valeur sélective des organismes, affectent leur évolvabilité, en donnant accès à de nouveaux phénotypes par le biais de mutations apparaissant ultérieurement, et qui n’auraient pas été disponibles autrement. En plus de cet effet, de nombreuses mutations neutres – indépendamment de leurs effets sélectifs – peuvent affecter la mutabilité de séquences d’ADN voisines, et moduler l’efficacité de la recombinaison homologue. De telles mutations ne modifient pas le spectre des phénotypes accessibles, mais plutôt la vitesse à laquelle de nouveaux phénotypes seront produits, un processus qui a des conséquences à long terme mais aussi potentiellement à court terme, en lien avec l’émergence de cancers.

Intragenic compensation through the lens of deep mutational scanning

A significant fraction of mutations in proteins are deleterious and result in adverse consequences for protein function, stability, or interaction with other molecules. Intragenic compensation is a specific case of positive epistasis when a neutral missense mutation cancels effect of a deleterious mutation in the same protein. Permissive compensatory mutations facilitate protein evolution, since without them all sequences would be extremely conserved. Understanding compensatory mechanisms is an important scientific challenge at the intersection of protein biophysics and evolution. In human genetics, intragenic compensatory interactions are important since they may result in variable penetrance of pathogenic mutations or fixation of pathogenic human alleles in orthologous proteins from related species. The latter phenomenon complicates computational and clinical inference of an allele's pathogenicity. Deep mutational scanning is a relatively new technique that enables experimental studies of functional effects of thousands of mutations in proteins. We review the important aspects of the field and discuss existing limitations of current datasets. We reviewed ten published DMS datasets with quantified functional effects of single and double mutations and described rates and patterns of intragenic compensation in eight of them.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12551-022-01005-w.

A Trait-Based Approach to Predicting Viral Host-Range Evolvability

Predicting the evolution of virus host range has proven to be extremely difficult, in part because of the sheer diversity of viruses, each with unique biology and ecological interactions. We have not solved this problem, but to make the problem more tractable, we narrowed our focus to three traits intrinsic to all viruses that may play a role in host-range evolvability: mutation rate, recombination rate, and phenotypic heterogeneity. Although each trait should increase evolvability, they cannot do so unbounded because fitness trade-offs limit the ability of all three traits to maximize evolvability. By examining these constraints, we can begin to identify groups of viruses with suites of traits that make them especially concerning, as well as ecological and environmental conditions that might push evolution toward accelerating host-range expansion.

Active site center redesign increases protein stability preserving catalysis in thioredoxin

The stabilization of natural proteins is a long-standing desired goal in protein engineering. Optimizing the hydrophobicity of the protein core often results in extensive stability enhancements. However, the presence of totally or partially buried catalytic charged residues, essential for protein function, has limited the applicability of this strategy. Here, focusing on the thioredoxin, we aimed to augment protein stability by removing buried charged residues in the active site without loss of catalytic activity. To this end, we performed a charged-to-hydrophobic substitution of a buried and functional group, resulting in a significant stability increase yet abolishing catalytic activity. Then, to simulate the catalytic role of the buried ionizable group, we designed a combinatorial library of variants targeting a set of seven surface residues adjacent to the active site. Notably, more than 50% of the library variants restored, to some extent, the catalytic activity. The combination of experimental study of 2% of the library with the prediction of the whole mutational space by partial least squares regression revealed that a single point mutation at the protein surface is sufficient to fully restore the catalytic activity without thermostability cost. As a result, we engineered one of the highest thermal stabilities reported for a protein with a natural occurring fold (137°C). Further, our hyperstable variant preserves the catalytic activity both in vitro and in vivo.

A Simulation Analysis and Screening of Deleterious Nonsynonymous Single Nucleotide Polymorphisms (nsSNPs) in Sheep LEP Gene

Leptin is a polypeptide hormone produced in the adipose tissue and governs many processes in the body. Recently, polymorphisms in the LEP gene revealed a significant change in body weight regulation, energy balance, food intake, and reproductive hormone secretion. This study considers its crucial role in the regulation of the economically important traits of sheep. Several computational tools, including SIFT, Predict SNP2, SNAP2, and PROVEAN, have been used to screen out the deleterious nsSNPs. Following the screening of 11 nsSNPs in the sheep genome, 5 nsSNPs, T86M (C → T), D98N (G → A), N136T (A → C), R142Q (G → A), and P157Q (C → A), were predicted to have a significant deleterious effect on the LEP protein function, leading to phenotypic difference. The analysis of proteins’ stability change due to amino acid substitution using the I-stable, SDM, and DynaMut consistently confirmed that three nsSNPs (T86M (C → T), D98N (G → A), and P157Q (C → A)) increased protein stability. It is suggested that these three nsSNPs may enhance the evolvability of LEP protein, which is vital for the evolutionary adaptation of sheep. Our findings demonstrate that the five nsSNPs reported in this study might be responsible for sheep’s structural and functional modifications of LEP protein. This is the first comprehensive report on the sheep LEP gene. It narrow downs the candidate nsSNPs for in vitro experiments to facilitate the development of reliable molecular markers for associated traits.

The biogenesis of β-lactamase enzymes

The discovery of penicillin by Alexander Fleming marked a new era for modern medicine, allowing not only the treatment of infectious diseases, but also the safe performance of life-saving interventions, like surgery and chemotherapy. Unfortunately, resistance against penicillin, as well as more complex β-lactam antibiotics, has rapidly emerged since the introduction of these drugs in the clinic, and is largely driven by a single type of extra-cytoplasmic proteins, hydrolytic enzymes called β-lactamases. While the structures, biochemistry and epidemiology of these resistance determinants have been extensively characterized, their biogenesis, a complex process including multiple steps and involving several fundamental biochemical pathways, is rarely discussed. In this review, we provide a comprehensive overview of the journey of β-lactamases, from the moment they exit the ribosomal channel until they reach their final cellular destination as folded and active enzymes.

Low protein expression enhances phenotypic evolvability by intensifying selection on folding stability

Protein abundance affects the evolution of protein genotypes, but we do not know how it affects the evolution of protein phenotypes. Here we investigate the role of protein abundance in the evolvability of green fluorescent protein (GFP) towards the novel phenotype of cyan fluorescence. We evolve GFP in E. coli through multiple cycles of mutation and selection and show that low GFP expression facilitates the evolution of cyan fluorescence. A computational model whose predictions we test experimentally helps explain why: lowly expressed proteins are under stronger selection for proper folding, which facilitates their evolvability on short evolutionary time scales. The reason is that high fluorescence can be achieved by either few proteins that fold well or by many proteins that fold less well. In other words, we observe a synergy between a protein's scarcity and its stability. Because many proteins meet the essential requirements for this scarcity-stability synergy, it may be a widespread mechanism by which low expression helps proteins evolve new phenotypes and functions.

In silico Prediction of Deleterious Single Nucleotide Polymorphism in S100A4 Metastatic Gene: Potential Early Diagnostic Marker

S100A4 protein overexpression has been reported in different types of cancer and plays a key role by interacting with the tumor suppressor protein Tp53. Single nucleotide polymorphisms (SNP) in S100A4 could directly influence the biomolecular interaction with the tumor suppressor protein Tp53 due to their aberrant conformations. Hence, the study was designed to predict the deleterious SNP and its effect on the S100A4 protein structure and function. Twenty-one SNP data sets were screened for nonsynonymous mutations and subsequently subjected to deleterious mutation prediction using different computational tools. The screened deleterious mutations were analyzed for their changes in functionality and their interaction with the tumor suppressor protein Tp53 by protein-protein docking analysis. The structural effects were studied using the 3DMissense mutation tool to estimate the solvation energy and torsion angle of the screened mutations on the predicted structures. In our study, 21 deleterious nonsynonymous mutations were screened, including F72V, E74G, L5P, D25E, N65S, A28V, A8D, S20L, L58P, and K26N were found to be remarkably conserved by exhibiting the interaction either with the EF-hand 1 or EF-hand 2 domain. The solvation and torsion values significantly deviated for the mutant-type structures with S20L, N65S, and F72L mutations and showed a marked reduction in their binding affinity with the Tp53 protein. Hence, these deleterious mutations might serve as prospective targets for diagnosing and developing personalized treatments for cancer and other related diseases.

Genetics of adaptation

The rediscovery of Mendel’s work showing that the heredity of phenotypes is controlled by discrete genes was followed by the reconciliation of Mendelian genetics with evolution by natural selection in the middle of the last century with the Modern Synthesis. In the past two decades, dramatic advances in genomic methods have facilitated the identification of the loci, genes, and even individual mutations that underlie phenotypic variants that are the putative targets of natural selection. Moreover, these methods have also changed how we can study adaptation by flipping the problem around, allowing us to first examine what loci show evidence of having been under selection, and then connecting these genetic variants to phenotypic variation. As a result, we now have an expanding list of actual genetic changes that underlie potentially adaptive phenotypic variation. Here, we synthesize how considering the effects of these adaptive loci in the context of cellular environments, genomes, organisms, and populations has provided new insights to the genetic architecture of adaptation.

Genome sequence diversity of SARS-CoV-2 obtained from clinical samples in Uzbekistan

Tracking temporal and spatial genomic changes and evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) are among the most urgent research topics worldwide, which help to elucidate the coronavirus disease 2019 (COVID-19) pathogenesis and the effect of deleterious variants. Our current study concentrates genetic diversity of SARS-CoV-2 variants in Uzbekistan and their associations with COVID-19 severity. Thirty-nine whole genome sequences (WGS) of SARS-CoV-2 isolated from PCR-positive patients from Tashkent, Uzbekistan for the period of July-August 2021, were generated and further subjected to further genomic analysis. Genome-wide annotations of clinical isolates from our study have revealed a total of 223 nucleotide-level variations including SNPs and 34 deletions at different positions throughout the entire genome of SARS-CoV-2. These changes included two novel mutations at the Nonstructural protein (Nsp) 13: A85P and Nsp12: Y479N, which were unreported previously. There were two groups of co-occurred substitution patterns: the missense mutations in the Spike (S): D614G, Open Reading Frame (ORF) 1b: P314L, Nsp3: F924, 5`UTR:C241T; Nsp3:P2046L and Nsp3:P2287S, and the synonymous mutations in the Nsp4:D2907 (C8986T), Nsp6:T3646A and Nsp14:A1918V regions, respectively. The “Nextstrain” clustered the largest number of SARS-CoV-2 strains into the Delta clade (n = 32; 82%), followed by two Alpha-originated (n = 4; 10,3%) and 20A (n = 3; 7,7%) clades. Geographically the Delta clade sample sequences were grouped into several clusters with the SARS-CoV genotypes from Russia, Denmark, USA, Egypt and Bangladesh. Phylogenetically, the Delta isolates in our study belong to the two main subclades 21A (56%) and 21J (44%). We found that females were more affected by 21A, whereas males by 21J variant (χ2 = 4.57; p ≤ 0.05, n = 32). The amino acid substitution ORF7a:P45L in the Delta isolates found to be significantly associated with disease severity. In conclusion, this study evidenced that Identified novel substitutions Nsp13: A85P and Nsp12: Y479N, have a destabilizing effect, while missense substitution ORF7a: P45L significantly associated with disease severity.

DnaK response to expression of protein mutants is dependent on translation rate and stability

Chaperones play a central part in the quality control system in cells by clearing misfolded and aggregated proteins. The chaperone DnaK acts as a sensor for molecular stress by recognising short hydrophobic stretches of misfolded proteins. As the level of unfolded protein is a function of protein stability, we hypothesised that the level of DnaK response upon overexpression of recombinant proteins would be correlated to stability. Using a set of mutants of the λ-repressor with varying thermal stabilities and a fluorescent reporter system, the effect of stability on DnaK response and protein abundance was investigated. Our results demonstrate that the initial DnaK response is largely dependent on protein synthesis rate but as the recombinantly expressed protein accumulates and homeostasis is approached the response correlates strongly with stability. Furthermore, we observe a large degree of cell-cell variation in protein abundance and DnaK response in more stable proteins.

Tuning the Catalytic Activity of Synthetic Enzyme KE15 with DNA

Efficiency improvement of synthetic enzymes through scaffold modifications suffers from limitations in terms of effectiveness, cost, and potential devastating consequences for protein structural stability. Here, we propose an alternative to scaffold modification, within electrostatic preorganization theory, where the enzyme's greater environment is designed to support the evolution of the reaction in the active site. We demonstrate the feasibility of such an approach by placing a (polar) DNA fragment in the surroundings of the Kemp eliminase enzyme KE15 (structure from Houk's group) and computing the resulting change in catalytic activity. We find that the introduction of a DNA fragment magnifies the contribution of protein residues to the stabilization of the transition state, estimated from electric field calculations with polarizable molecular dynamics. Our randomly generated test systems reveal a 2.0 kcal/mol reduction in activation energy, suggesting that even more significant catalytic improvements could be made by optimizing DNA size, sequence, and orientation with respect to the enzyme, validating our approach.

Functionalization of a symmetric protein scaffold: Redundant folding nuclei and alternative oligomeric folding pathways

Successful de novo protein design ideally targets specific folding kinetics, stability thermodynamics, and biochemical functionality, and the simultaneous achievement of all these criteria in a single step design is challenging. Protein design is potentially simplified by separating the problem into two steps: (a) an initial design of a protein "scaffold" having appropriate folding kinetics and stability thermodynamics, followed by (b) appropriate functional mutation-possibly involving insertion of a peptide functional "cassette." This stepwise approach can also separate the orthogonal effects of the "stability/function" and "foldability/function" tradeoffs commonly observed in protein design. If the scaffold is a protein architecture having an exact rotational symmetry, then there is the potential for redundant folding nuclei and multiple equivalent sites of functionalization; thereby enabling broader functional adaptation. We describe such a "scaffold" and functional "cassette" design strategy applied to a β-trefoil threefold symmetric architecture and a heparin ligand functionality. The results support the availability of redundant folding nuclei within this symmetric architecture, and also identify a minimal peptide cassette conferring heparin affinity. The results also identify an energy barrier of destabilization that switches the protein folding pathway from monomeric to trimeric, thereby identifying another potential advantage of symmetric protein architecture in de novo design.

Enhanced thermal and alkaline stability of L-lysine decarboxylase CadA by combining directed evolution and computation-guided virtual screening

As an important monomer for bio-based nylons PA5X, cadaverine is mainly produced by enzymatic decarboxylation of L-lysine. A key issue with this process is the instability of L-lysine decarboxylase (CadA) during the reaction due to the dissociation of CadA subunits with the accumulation of alkaline cadaverine. In this work, we attempted to improve the thermal and alkaline stability of CadA by combining directed evolution and computation-guided virtual screening. Interestingly, site 477 residue located at the protein surface and not the decamer interface was found as a hotspot in directed evolution. By combinatorial mutagenesis of the positive mutations obtained by directed evolution and virtual screening with the previously reported T88S mutation, K477R/E445Q/T88S/F102V was generated as the best mutant, delivering 37% improvement of cadaverine yield at 50 ºC and pH 8.0. Molecular dynamics simulations suggested the improved rigidity of regional structures, increased number of salt bridges, and enhancement of hydrogen bonds at the multimeric interface as possible origins of the improved stability of the mutant. Using this four-point mutant, 160.7 g/L of cadaverine was produced from 2.0 M Lysine hydrochloride at 50 °C without pH regulation, with a conversion of 78.5%, whereas the wild type produced 143.7 g/L cadaverine, corresponding to 70% conversion. This work shows the combination of directed evolution and virtual screening as an efficient protein engineering strategy.

Engineering Strategies to Overcome the Stability-Function Trade-Off in Proteins

In addition to its biological function, the stability of a protein is a major determinant for its applicability. Unfortunately, engineering proteins for improved functionality usually results in destabilization of the protein. This so-called stability–function trade-off can be explained by the simple fact that the generation of a novel protein function—or the improvement of an existing one—necessitates the insertion of mutations, i.e., deviations from the evolutionarily optimized wild-type sequence. In fact, it was demonstrated that gain-of-function mutations are not more destabilizing than other random mutations. The stability–function trade-off is a universal phenomenon during protein evolution that has been observed with completely different types of proteins, including enzymes, antibodies, and engineered binding scaffolds. In this review, we discuss three types of strategies that have been successfully deployed to overcome this omnipresent obstacle in protein engineering approaches: (i) using highly stable parental proteins, (ii) minimizing the extent of destabilization during functional engineering (by library optimization and/or coselection for stability and function), and (iii) repairing damaged mutants through stability engineering. The implementation of these strategies in protein engineering campaigns will facilitate the efficient generation of protein variants that are not only functional but also stable and therefore better-suited for subsequent applications.

Efficient Unnatural Protein Production by Pyrrolysyl-tRNA Synthetase With Genetically Fused Solubility Tags

Introducing non-canonical amino acids (ncAAs) by engineered orthogonal pairs of aminoacyl-tRNA synthetases and tRNAs has proven to be a highly useful tool for the expansion of the genetic code. Pyrrolysyl-tRNA synthetase (PylRS) from methanogenic archaeal and bacterial species is particularly attractive due to its natural orthogonal reactivity in bacterial and eukaryotic cells. However, the scope of such a reprogrammed translation is often limited, due to low yields of chemically modified target protein. This can be the result of substrate specificity engineering, which decreases the aminoacyl-tRNA synthetase stability and reduces the abundance of active enzyme. We show that the solubility and folding of these engineered enzymes can become a bottleneck for the production of ncAA-containing proteins in vivo. Solubility tags derived from various species provide a strategy to remedy this issue. We find the N-terminal fusion of the small metal binding protein from Nitrosomonas europaea to the PylRS sequence to improve enzyme solubility and to boost orthogonal translation efficiency. Our strategy enhances the production of site-specifically labelled proteins with a variety of engineered PylRS variants by 200–540%, and further allows triple labeling. Even the wild-type enzyme gains up to 245% efficiency for established ncAA substrates.

GroEL/S overexpression helps to purge deleterious mutations and reduce genetic diversity during adaptive protein evolution

Chaperones are proteins that help other proteins fold. They also affect the adaptive evolution of their client proteins by buffering the effect of deleterious mutations and increasing the genetic diversity of evolving proteins. We study how the bacterial chaperone GroE (GroEL + GroES) affects the evolution of green fluorescent protein (GFP). To this end we subjected GFP to multiple rounds of mutation and selection for its color phenotype in four replicate E. coli populations, and studied its evolutionary dynamics through high-throughput sequencing and mutant engineering. We evolved GFP both under stabilizing selection for its ancestral (green) phenotype, and to directional selection for a new (cyan) phenotype. We did so both under low and high expression of the chaperone GroE. In contrast to previous work, we observe that GroE does not just buffer but also helps purge deleterious (fluorescence reducing) mutations from evolving populations. In doing so, GroE helps reduce the genetic diversity of evolving populations. In addition, it causes phenotypic heterogeneity in mutants with the same genotype, helping to enhance their fluorescence in some cells, and reducing it in others. Our observations show that chaperones can affect adaptive evolution in more than one way.

Data set and fitting dependencies when estimating protein mutant stability: Toward simple, balanced, and interpretable models

Accurate prediction of protein stability changes upon mutation (ΔΔG) is increasingly important to evolution studies, protein engineering, and screening of disease-causing gene variants but is challenged by biases in training data. We investigated 45 linear regression models trained on data sets that account systematically for destabilization bias and mutation-type bias B_M . The models were externally validated on three test data sets probing different pathologies and for internal consistency (symmetry and neutrality). Model structure and performance substantially depended on training data and even fitting method. We developed two final models: SimBa-IB for typical natural mutations and SimBa-SYM for situations where stabilizing and destabilizing mutations occur to a similar extent. SimBa-SYM, despite is simplicity, is essentially non-biased (vs. the S^sym data set) while still performing well for all data sets (R ~ 0.46-0.54, MAE = 1.16-1.24 kcal/mol). The simple models provide advantage in terms of interpretability, use and future improvement, and are freely available on GitHub.

Structure and Mutations of SARS-CoV-2 Spike Protein: A Focused Overview

The spike protein (S-protein) of SARS-CoV-2, the protein that enables the virus to infect human cells, is the basis for many vaccines and a hotspot of concerning virus evolution. Here, we discuss the outstanding progress in structural characterization of the S-protein and how these structures facilitate analysis of virus function and evolution. We emphasize the differences in reported structures and that analysis of structure-function relationships is sensitive to the structure used. We show that the average residue solvent exposure in nearly complete structures is a good descriptor of open vs closed conformation states. Because of structural heterogeneity of functionally important surface-exposed residues, we recommend using averages of a group of high-quality protein structures rather than a single structure before reaching conclusions on specific structure-function relationships. To illustrate these points, we analyze some significant chemical tendencies of prominent S-protein mutations in the context of the available structures. In the discussion of new variants, we emphasize the selectivity of binding to ACE2 vs prominent antibodies rather than simply the antibody escape or ACE2 affinity separately. We note that larger chemical changes, in particular increased electrostatic charge or side-chain volume of exposed surface residues, are recurring in mutations of concern, plausibly related to adaptation to the negative surface potential of human ACE2. We also find indications that the fixated mutations of the S-protein in the main variants are less destabilizing than would be expected on average, possibly pointing toward a selection pressure on the S-protein. The richness of available structures for all of these situations provides an enormously valuable basis for future research into these structure-function relationships.

Site-wise Diversification of Combinatorial Libraries Using Insights from Structure-guided Stability Calculations

Many auspicious clinical and industrial accomplishments have improved the human condition by means of protein engineering. Despite these achievements, our incomplete understanding of the sequence-structure-function relationship prevents rapid innovation. To tackle this problem, we must develop and integrate new and existing technologies. To date, directed evolution and rational design have dominated as protein engineering principles. Even so, prior to screening for novel or improved functions, a large collection of variants, within a protein library, exist along an ambiguous mutational terrain. Complicating things further, the choice of where to initialize investigation along a vast sequence space becomes even more difficult given that the majority of any sequence lacks function entirely. Unfortunately, even when considering functionally relevant positions, random substitutions can prove to be destabilizing, causing a hindrance to an otherwise function-inducing, stability-reliant folding process. To enhance productivity in the field, we seek to address this issue of destabilization, and subsequent disfunction, at protein-protein and protein-ligand interacting regions. Herein, the process of choosing amenable positions - and amino acids at those positions - allows for a refined, knowledge-based approach to combinatorial library design. Using structural data, we perform computational stability prediction with FoldX's PositionScan and Rosetta's ddG_monomer in tandem, allowing for the refinement of our thermodynamic stability data through the comparison of results. In turn, we provide a process for selecting in silico predicted mutually stabilizing positions and avoiding overly destabilizing ones that guides the site-wise diversification of combinatorial libraries.

Recombination of Single Beneficial Substitutions Obtained from Protein Engineering by Computer-Assisted Recombination (CompassR)

A large number of beneficial substitutions can be obtained from a successful directed enzyme evolution campaign and/or (semi)rational design. It is expected that the recombination of some beneficial substitutions leads to a much higher degree of performance through synergistic effect. However, systematic recombination studies show that poorly performing variants are often obtained after recombination of three to four individual beneficial substitutions and this limits protein engineers to exploit nature's potential in generating better performing enzymes. Computer-assisted Recombination (CompassR) strategy allows the recombination of identified beneficial substitutions in an effective and efficient manner in order to generate active enzymes with improved performance. Here, we describe in detail the CompassR procedure with an example of recombining four substitutions and discuss some important practical issues that should be considered (such as the selection of protein structures, number of FoldX runs, evaluation of calculations) for application of the CompassR rule. The core part of this protocol (system setup, ΔΔG_fold calculation, and CompassR application) is transferable to other enzymes and any recombination of single beneficial substitutions.

Conformational dynamics promotes disordered regions from function-dispensable to essential in evolved site-specific DNA recombinases

•

New functional regions emerging in evolution of DNA site-specific recombinase tails.

•

Transient structural nucleation promotes function-dispensable regions to essential.

•

Molecular dynamics reveals conformational diversity and its functional implications.

•

Evolved disordered molecular mechanisms of N-term tails for protein stability.

•

Structural disorder-based link between protein evolution, stability and function.

Protein intrinsically disordered regions (IDRs) play pivotal roles in molecular recognition and regulatory processes through structural disorder-to-order transitions. To understand and exploit the distinctive functional implications of IDRs and to unravel the underlying molecular mechanisms, structural disorder-to-function relationships need to be deciphered. The DNA site-specific recombinase system Cre/loxP represents an attractive model to investigate functional molecular mechanisms of IDRs. Cre contains a functionally dispensable disordered N-terminal tail, which becomes indispensable in the evolved Tre/loxLTR recombinase system. The difficulty to experimentally obtain structural information about this tail has so far precluded any mechanistic study on its involvement in DNA recombination. Here, we use in vitro and in silico evolution data, conformational dynamics, AI-based folding simulations, thermodynamic stability calculations, mutagenesis and DNA recombination assays to investigate how evolution and the dynamic behavior of this IDR may determine distinct functional properties. Our studies suggest that partial conformational order in the N-terminal tail of Tre recombinase and its packing to a conserved hydrophobic surface on the protein provide thermodynamic stability. Based on our results, we propose a link between protein stability and function, offering new plausible atom-detailed mechanistic insights into disorder-function relationships. Our work highlights the potential of N-terminal tails to be exploited for regulation of the activity of Cre-like tyrosine-type SSRs, which merits future investigations and could be of relevance in future rational engineering for their use in biotechnology and genomic medicine.

Enabling biocatalysis in high-concentration organic cosolvent by enzyme gate engineering

Biocatalysis in high-concentration organic solvents (OSs) offers many advantages, but realizing this process remains a huge challenge. An R-selective ω-amine transaminase variant (AcATA_M2 ) exhibited high activity toward 50 g/L pro-sitagliptin ketone 1-[1-piperidinyl]-4-[2,4,5-trifluorophenyl]-1,3-butanedione (PTfpB). However, AcATA_M2 displayed unsatisfactory organic-cosolvent resistance against high-concentration dimethyl sulfoxide (DMSO), which is required to enhance the solubility of the hydrophobic substrate PTfpB. Located in the substrate-binding tunnel, enzyme gates are structural elements that undergo reversible conformational transitions, thus affecting the accessibility of the binding pocket to solvent molecules. Depending on the conformation of the enzyme gates, one can define an open or closed conformation on which the enzyme activity in OSs may depend. To enhance the DMSO resistance of AcATA_M2 , we identified the beneficial residues at the "enzyme gate" region via computational analysis, alanine scanning, and site-saturation mutagenesis. Two beneficial variants, namely, AcATA_M2 ^F56D and AcATA_M2 ^F56V , not only displayed improved enzyme activity but also exhibited enhanced DMSO resistance (the half-life value increased from 25.71 to 42.49 h under 60% DMSO). Molecular dynamic simulations revealed that the increase in DMSO resistance was mainly caused by the decrease in the number of DMSO molecules in the substrate-binding pocket. Moreover, in the kilogram-scale experiment, the conversion of 80 g/L substrate was increased from 50% (AcATA_M2 ) to 85% (M2^F56D in 40% DMSO) with a high e.e. of >99% within 24 h.

Climbing Up and Down Binding Landscapes through Deep Mutational Scanning of Three Homologous Protein-Protein Complexes

Protein-protein interactions (PPIs) have evolved to display binding affinities that can support their function. As such, cognate and noncognate PPIs could be highly similar structurally but exhibit huge differences in binding affinities. To understand this phenomenon, we study three homologous protease-inhibitor PPIs that span 9 orders of magnitude in binding affinity. Using state-of-the-art methodology that combines protein randomization, affinity sorting, deep sequencing, and data normalization, we report quantitative binding landscapes consisting of ΔΔGbind values for the three PPIs, gleaned from tens of thousands of single and double mutations. We show that binding landscapes of the three complexes are strikingly different and depend on the PPI evolutionary optimality. We observe different patterns of couplings between mutations for the three PPIs with negative and positive epistasis appearing most frequently at hot-spot and cold-spot positions, respectively. The evolutionary trends observed here are likely to be universal to other biological complexes in the cell.

Enzyme engineering and in vivo testing of a formate reduction pathway

Formate is an attractive feedstock for sustainable microbial production of fuels and chemicals, but its potential is limited by the lack of efficient assimilation pathways. The reduction of formate to formaldehyde would allow efficient downstream assimilation, but no efficient enzymes are known for this transformation. To develop a 2-step formate reduction pathway, we screened natural variants of acyl-CoA synthetase (ACS) and acylating aldehyde dehydrogenase (ACDH) for activity on one-carbon substrates and identified active and highly expressed homologs of both enzymes. We then performed directed evolution, increasing ACDH-specific activity by 2.5-fold and ACS lysate activity by 5-fold. To test for the in vivo activity of our pathway, we expressed it in a methylotroph which can natively assimilate formaldehyde. Although the enzymes were active in cell extracts, we could not detect formate assimilation into biomass, indicating that further improvement will be required for formatotrophy. Our work provides a foundation for further development of a versatile pathway for formate assimilation.

Shifts in amino acid preferences as proteins evolve: A synthesis of experimental and theoretical work

Amino acid preferences vary across sites and time. While variation across sites is widely accepted, the extent and frequency of temporal shifts are contentious. Our understanding of the drivers of amino acid preference change is incomplete: To what extent are temporal shifts driven by adaptive versus nonadaptive evolutionary processes? We review phenomena that cause preferences to vary (e.g., evolutionary Stokes shift, contingency, and entrenchment) and clarify how they differ. To determine the extent and prevalence of shifted preferences, we review experimental and theoretical studies. Analyses of natural sequence alignments often detect decreases in homoplasy (convergence and reversions) rates, and variation in replacement rates with time-signals that are consistent with temporally changing preferences. While approaches inferring shifts in preferences from patterns in natural alignments are valuable, they are indirect since multiple mechanisms (both adaptive and nonadaptive) could lead to the observed signal. Alternatively, site-directed mutagenesis experiments allow for a more direct assessment of shifted preferences. They corroborate evidence from multiple sequence alignments, revealing that the preference for an amino acid at a site varies depending on the background sequence. However, shifts in preferences are usually minor in magnitude and sites with significantly shifted preferences are low in frequency. The small yet consistent perturbations in preferences could, nevertheless, jeopardize the accuracy of inference procedures, which assume constant preferences. We conclude by discussing if and how such shifts in preferences might influence widely used time-homogenous inference procedures and potential ways to mitigate such effects.

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.

SCONES: Self-Consistent Neural Network for Protein Stability Prediction Upon Mutation

Engineering proteins to have desired properties by mutating amino acids at specific sites is commonplace. Such engineered proteins must be stable to function. Experimental methods used to determine stability at throughputs required to scan the protein sequence space thoroughly are laborious. To this end, many machine learning based methods have been developed to predict thermodynamic stability changes upon mutation. These methods have been evaluated for symmetric consistency by testing with hypothetical reverse mutations. In this work, we propose transitive data augmentation, evaluating transitive consistency with our new S^transitive data set, and a new machine learning based method, the first of its kind, that incorporates both symmetric and transitive properties into the architecture. Our method, called SCONES, is an interpretable neural network that predicts small relative protein stability changes for missense mutations that do not significantly alter the structure. It estimates a residue's contributions toward protein stability (ΔG) in its local structural environment, and the difference between independently predicted contributions of the reference and mutant residues is reported as ΔΔG. We show that this self-consistent machine learning architecture is immune to many common biases in data sets, relies less on data than existing methods, is robust to overfitting, and can explain a substantial portion of the variance in experimental data.

Evidence that toxin resistance in poison birds and frogs is not rooted in sodium channel mutations and may rely on "toxin sponge" proteins

Many poisonous organisms carry small-molecule toxins that alter voltage-gated sodium channel (NaV) function. Among these, batrachotoxin (BTX) from Pitohui poison birds and Phyllobates poison frogs stands out because of its lethality and unusual effects on NaV function. How these toxin-bearing organisms avoid autointoxication remains poorly understood. In poison frogs, a NaV DIVS6 pore-forming helix N-to-T mutation has been proposed as the BTX resistance mechanism. Here, we show that this variant is absent from Pitohui and poison frog NaVs, incurs a strong cost compromising channel function, and fails to produce BTX-resistant channels in poison frog NaVs. We also show that captivity-raised poison frogs are resistant to two NaV-directed toxins, BTX and saxitoxin (STX), even though they bear NaVs sensitive to both. Moreover, we demonstrate that the amphibian STX "toxin sponge" protein saxiphilin is able to protect and rescue NaVs from block by STX. Taken together, our data contradict the hypothesis that BTX autoresistance is rooted in the DIVS6 N→T mutation, challenge the idea that ion channel mutations are a primary driver of toxin resistance, and suggest the possibility that toxin sequestration mechanisms may be key for protecting poisonous species from the action of small-molecule toxins.

Structural determinants and distribution of phosphate specificity in ribonucleotide reductases

Ribonucleotide reductases (RNRs) catalyze the reduction of ribonucleotides to the corresponding deoxyribonucleotides, the building blocks of DNA. RNRs are specific for either ribonucleoside diphosphates or triphosphates as substrates. As far as is known, oxygen-dependent class I RNRs (NrdAB) all reduce ribonucleoside diphosphates, and oxygen-sensitive class III RNRs (NrdD) are all ribonucleoside triphosphate reducers, whereas the adenosylcobalamin-dependent class II (NrdJ) contains both ribonucleoside diphosphate and triphosphate reducers. However, it is unknown how this specificity is conveyed by the active site of the enzymes and how this feature developed in RNR evolution. By structural comparison of the active sites in different RNRs, we identified the apical loop of the phosphate-binding site as a potential structural determinant of substrate specificity. Grafting two residues from this loop from a diphosphate- to a triphosphate-specific RNR caused a change in preference from ribonucleoside triphosphate to diphosphate substrates in a class II model enzyme, confirming them as the structural determinants of phosphate specificity. The investigation of the phylogenetic distribution of this motif in class II RNRs yielded a likely monophyletic clade with the diphosphate-defining motif. This indicates a single evolutionary-split event early in NrdJ evolution in which diphosphate specificity developed from the earlier triphosphate specificity. For those interesting cases where organisms contain more than one nrdJ gene, we observed a preference for encoding enzymes with diverse phosphate specificities, suggesting that this varying phosphate specificity confers a selective advantage.

Evolution-Based Protein Engineering for Antifungal Peptide Improvement

Antimicrobial peptides (AMPs) have been considered as the alternatives to antibiotics because of their less susceptibility to microbial resistance. However, compared with conventional antibiotics they show relatively low activity and the consequent high cost and nonspecific cytotoxicity, hindering their clinical application. What’s more, engineering of AMPs is a great challenge due to the inherent complexity in their sequence, structure, and function relationships. Here, we report an evolution-based strategy for improving the antifungal activity of a nematode-sourced defensin (Cremycin-5). This strategy utilizes a sequence-activity comparison between Cremycin-5 and its functionally diverged paralogs to identify sites associated with antifungal activity for screening of enhanceable activity-modulating sites for subsequent saturation mutagenesis. Using this strategy, we identified a site (Glu-15) whose mutations with nearly all other types of amino acids resulted in a universally enhanced activity against multiple fungal species, which is thereby defined as a Universally Enhanceable Activity-Modulating Site (UEAMS). Especially, Glu15Lys even exhibited >9-fold increased fungicidal potency against several clinical isolates of Candida albicans through inhibiting cytokinesis. This mutant showed high thermal and serum stability and quicker killing kinetics than clotrimazole without detectable hemolysis. Molecular dynamic simulations suggest that the mutations at the UEAMS likely limit the conformational flexibility of a distant functional residue via allostery, enabling a better peptide–fungus interaction. Further sequence, structural, and mutational analyses of the Cremycin-5 ortholog uncover an epistatic interaction between the UEAMS and another site that may constrain its evolution. Our work lights one new road to success of engineering AMP drug leads.

Chaperonin Abundance Enhances Bacterial Fitness

The ability of chaperonins to buffer mutations that affect protein folding pathways suggests that their abundance should be evolutionarily advantageous. Here, we investigate the effect of chaperonin overproduction on cellular fitness in Escherichia coli. We demonstrate that chaperonin abundance confers 1) an ability to tolerate higher temperatures, 2) improved cellular fitness, and 3) enhanced folding of metabolic enzymes, which is expected to lead to enhanced energy harvesting potential.

In silico analysis of non-synonymous missense SNPs (nsSNPs) in CPE, GNAS genes and experimental validation in type II diabetes mellitus through Next Generation Sequencing

Non-synonymous missense SNPs (nsSNPs) in CPE and GNAS genes were investigated computationally. In silico identified nsSNPs were experimentally validated in type II diabetes mellitus (T2DM) in Pakistani Pathan population using next generation sequencing (NGS). Sixty two high-risk nsSNPs in CPE and 44 in GNAS were identified. Only 12 in GNAS were clinically significant. Thirty six high-risk nsSNPs in CPE and 08 clinically significant nsSNPs in GNAS lies in the most conserved regions. I-mutant predicted that nsSNPs decrease the proteins stability and ModPred predicted 20 and 12 post-translational modification sites in CPE and GNAS proteins respectively. Ramachandran plot showed 88.7% residues are in the most favored region of protein models. By experimentation, none of the nsSNPs were found to be associated with T2DM. In conclusion, this study differentiates the deleterious nsSNPs from the neutral ones. Although nsSNPs are not associated with T2DM, they can be targeted in other CPE and GNAS genes related disorders.

Resurrection of Ancestral Malate Dehydrogenases Reveals the Evolutionary History of Halobacterial Proteins : Deciphering Gene Trajectories and Changes in Biochemical Properties

Extreme halophilic Archaea thrive in high salt, where, through proteomic adaptation, they cope with the strong osmolarity and extreme ionic conditions of their environment. In spite of wide fundamental interest, however, studies providing insights into this adaptation are scarce, because of practical difficulties inherent to the purification and characterization of halophilic enzymes. In this work, we describe the evolutionary history of malate dehydrogenases (MalDH) within Halobacteria (a class of the Euryarchaeota phylum). We resurrected nine ancestors along the inferred halobacterial MalDH phylogeny, including the Last Common Ancestral MalDH of Halobacteria (LCAHa) and compared their biochemical properties with those of five modern halobacterial MalDHs. We monitored the stability of these various MalDHs, their oligomeric states and enzymatic properties, as a function of concentration for different salts in the solvent. We found that a variety of evolutionary processes such as amino acid replacement, gene duplication, loss of MalDH gene and replacement owing to horizontal transfer resulted in significant differences in solubility, stability and catalytic properties between these enzymes in the three Halobacteriales, Haloferacales and Natrialbales orders since the LCAHa MalDH.We also showed how a stability trade-off might favor the emergence of new properties during adaptation to diverse environmental conditions. Altogether, our results suggest a new view of halophilic protein adaptation in Archaea.

Resource uptake and the evolution of moderately efficient enzymes

Enzymes speed up reactions that would otherwise be too slow to sustain the metabolism of selfreplicators. Yet, most enzymes seem only moderately efficient, exhibiting kinetic parameters orders of magnitude lower than their expected physically achievable maxima and spanning over surprisingly large ranges of values. Here, we question how these parameters evolve using a mechanistic model where enzyme efficiency is a key component of individual competition for resources. We show that kinetic parameters are under strong directional selection only up to a point, above which enzymes appear to evolve under near-neutrality, thereby confirming the qualitative observation of other modeling approaches. While the existence of a large fitness plateau could potentially explain the extensive variation in enzyme features reported, we show using a population genetics model that such a widespread distribution is an unlikely outcome of evolution on a common landscape, as mutation–selection–drift balance occupy a narrow area even when very moderate biases towards lower efficiency are considered. Instead, differences in the evolutionary context encountered by each enzyme should be involved, such that each evolves on an individual, unique landscape. Our results point to drift and effective population size playing an important role, along with the kinetics of nutrient transporters, the tolerance to high concentrations of intermediate metabolites, and the reversibility of reactions. Enzyme concentration also shapes selection on kinetic parameters, but we show that the joint evolution of concentration and efficiency does not yield extensive variance in evolutionary outcomes when documented costs to protein expression are applied.

Extra disulfide and ionic salt bridge improves the thermostability of lignin peroxidase H8 under acidic condition

The development of a lignin peroxidase (LiP) that is thermostable even under acidic pH conditions is a main issue for efficient enzymatic lignin degradation due to reduced repolymerization of free phenolic products at acidic pH (< 3). Native LiP under mild conditions (half-life (t_1/2) of 8.2 days at pH 6) exhibits a marked decline in thermostability under acidic conditions (t_1/2 of only 14 min at pH 2.5). Thus, improving the thermostability of LiP in acidic environments is required for effective lignin depolymerization in practical applications. Here, we show the improved thermostability of a synthetic LiPH8 variant (S49C/A67C/H239E, PDB: 6ISS) capable of strengthening the helix-loop interactions under acidic conditions. This variant retained excellent thermostability at pH 2.5 with a 10-fold increase in t_1/2 (2.52 h at 25 °C) compared with that of the native enzyme. X-ray crystallography analysis showed that the recombinant LiPH8 variant is the only unique lignin peroxidase containing five disulfide bridges, and the helix-loop interactions of the synthetic disulfide bridge and ionic salt bridge in its structure are responsible for stabilizing the Ca²⁺-binding region and heme environment, resulting in an increase in overall structural resistance against acidic conditions. Our work will allow the design of biocatalysts for ligninolytic enzyme engineering and for efficient biocatalytic degradation of plant biomass in lignocellulose biorefineries.