Robust enzyme design: bioinformatic tools for improved protein stability

Enhanced molecular stability of ApxII antigen during secretion in Corynebacterium glutamicum by rational design

ApxII is a vaccine antigen used to protect against porcine contagious pleuropneumonia, which is a significant threat to the pig industry. Here, we aimed to improve the proteolytic degradation stability of ApxII during its secretion by establishing a complete screening process of stable variants through bioinformatics and site-directed mutagenesis. We employed a combination of semi-rational and rational design strategies to create 34 single-point variants of ApxII. Among them, R114E and T115D variants exhibited better stability without compromising antigen activity. Furthermore, we constructed a multi-site variant, R114E/T115D, which demonstrated the best stability, activity, and yield. Protein stability and molecular dynamic analysis indicated that the greater solubility and lower structural expansion coefficient might explain the increased stability of R114E/T115D. Additionally, site T115 was identified as a key point of truncated ApxII stability. The R114E/T115D variant, with its proven stability and intact antigenic activity, holds promising prospects for industrial-scale applications in the prevention of porcine contagious pleuropneumonia.

Sustainable production of the drug precursor tyramine by engineered Corynebacterium glutamicum

Tyramine has attracted considerable interest due to recent findings that it is an excellent starting material for the production of high-performance thermoplastics and hydrogels. Furthermore, tyramine is a precursor of a diversity of pharmaceutically relevant compounds, contributing to its growing importance. Given the limitations of chemical synthesis, including lack of selectivity and laborious processes with harsh conditions, the biosynthesis of tyramine by decarboxylation of L-tyrosine represents a promising sustainable alternative. In this study, the de novo production of tyramine from simple nitrogen and sustainable carbon sources was successfully established by metabolic engineering of the L-tyrosine overproducing Corynebacterium glutamicum strain AROM3. A phylogenetic analysis of aromatic-L-amino acid decarboxylases (AADCs) revealed potential candidate enzymes for the decarboxylation of tyramine. The heterologous overexpression of the respective AADC genes resulted in successful tyramine production, with the highest tyramine titer of 1.9 g L-1 obtained for AROM3 overexpressing the tyrosine decarboxylase gene of Levilactobacillus brevis. Further metabolic engineering of this tyramine-producing strain enabled tyramine production from the alternative carbon sources ribose and xylose. Additionally, up-scaling of tyramine production from xylose to a 1.5 L bioreactor batch fermentation was demonstrated to be stable, highlighting the potential for sustainable tyramine production. KEY POINTS: • Phylogenetic analysis revealed candidate l-tyrosine decarboxylases • C. glutamicum was engineered for de novo production of tyramine • Tyramine production from alternative carbon substrates was enabled.

Characterization and rational engineering of an alkaline-tolerant azoreductase derived from Roseibium sp. H3510 for enhanced decolorization of azo dyes

rAzoR2326, an azoreductase derived from Roseibium sp. H3510, functions as an FMN-dependent homodimer utilizing NADH as cofactor. It demonstrated maximum activity at 45 °C and retained moderate activity above 50 °C, exhibiting stability from pH 7-10. Evolution and structure guided rational design of wild-type rAzoR2326 (WT) efficiently yielded 6 single-point mutants with improved thermostability and activity from a 22-variant library. Further combinatorial mutation led to mutant M20 with substantially enhanced thermostability (15-fold longer half-life at 50 °C) and activity (3.24-fold higher kcat/Km). M20 exhibited superior catalytic properties for decolorizing Allura Red compared to WT. Specifically, its decolorization capacity at pH 10.0 was 4.26-fold higher than WT. Additionally, M20 demonstrated remarkable thermostability, retaining 76.83 % decolorization activity for Allura Red after 120 min at 50 °C, whereas WT nearly lost all catalytic activity under the same conditions. Molecular dynamics simulations revealed the structural changes in M20, such as improved hydrogen bonding and a new C-H···π interaction, led to a more compact and rigid enzyme structure. This resulted in a more stable FMN-binding pocket and substrate tunnel, thereby improving the catalytic stability and activity of M20. Given its enhanced dye decolorization ability and alkaline tolerance, M20 shows promise as a biocatalyst for treating azo dye effluents.

ThermoLink: Bridging disulfide bonds and enzyme thermostability through database construction and machine learning prediction

Disulfide bonds, covalently formed by sulfur atoms in cysteine residues, play a crucial role in protein folding and structure stability. Considering their significance, artificial disulfide bonds are often introduced to enhance protein thermostability. Although an increasing number of tools can assist with this task, significant amounts of time and resources are often wasted owing to inadequate consideration. To enhance the accuracy and efficiency of designing disulfide bonds for protein thermostability improvement, we initially collected disulfide bond and protein thermostability data from extensive literature sources. Thereafter, we extracted various sequence- and structure-based features and constructed machine-learning models to predict whether disulfide bonds can improve protein thermostability. Among all models, the neighborhood context model based on the Adaboost-DT algorithm performed the best, yielding "area under the receiver operating characteristic curve" and accuracy scores of 0.773 and 0.714, respectively. Furthermore, we also found AlphaFold2 to exhibit high superiority in predicting disulfide bonds, and to some extent, the coevolutionary relationship between residue pairs potentially guided artificial disulfide bond design. Moreover, several mutants of imine reductase 89 (IR89) with artificially designed thermostable disulfide bonds were experimentally proven to be considerably efficient for substrate catalysis. The SS-bond data have been integrated into an online server, namely, ThermoLink, available at guolab.mpu.edu.mo/thermoLink.

Heterologous expression and characterization of an M4 family extracellular metalloprotease for detergent application

Proteolytic enzymes stand out as the most widely employed category utilized in manufacturing industry. A new protease was separated from Planococcus sp.11815 strain and named as nprS-15615 in this research. The gene of this protease has not been reported, and its enzymatic properties have been studied for the first time. To enhance enzyme production, the Planococcus sp. protease gene was expressed in Bacillus licheniformis 2709. The expression level of nprS-15615 was observed under the control of regulatory elements PaprE. nprS-15615 protease activity reached 1186.24±32.87 U/mL after 48 hours of cultivation in shake flasks which was nearly four times the output of the original bacteria (291.38±25.73U/mL). The optimum temperature and pH of the recombinant protease were 30 ℃ and 8.0, respectively.The enzyme exhibited the highest capacity for hydrolyzing casein and demonstrated resilience towards a NaCl concentration of 10.0% (wt/v). Furthermore, in the presence of 0.5% surfactants, the recombinant protease activity can maintain above 75%, and with the existence of 0.5% liquid detergents, there was basically no loss of enzyme activity which indicated that nprS-15615 had good compatibility with surfactants and liquid detergents. In addition, npS-15615 performed well in the washing experiment, and the washing effect at 20 ℃ can be significantly improved by adding crude enzyme solution in the washing process.

Carbohydrate-active enzymes in animal feed

Considering an ever-growing global population, which hit 8 billion people in the fall of 2022, it is essential to find solutions to avoid the competition between human food and animal feed for croplands. Agricultural co-products have become important components of the circular economy with their use in animal feed. Their implementation was made possible by the addition of exogenous enzymes in the diet, especially carbohydrate-active enzymes (CAZymes). In this review, we describe the diversity and versatility of microbial CAZymes targeting non-starch polysaccharides to improve the nutritional potential of diets containing cereals and protein meals. We focused our attention on cellulases, hemicellulases, pectinases which were often found to be crucial in vivo. We also highlight the performance and health benefits brought by the exogenous addition of enzymatic cocktails containing CAZymes in the diets of monogastric animals. Taking the example of the well-studied commercial cocktail Rovabio™, we discuss the evolution, constraints and future challenges faced by feed enzymes suppliers. We hope that this review will promote the use and development of enzyme solutions for industries to sustainably feed humans in the future.

Protein Stability: Enhancement and Measurement

This chapter defines protein stability, emphasizes its importance, and surveys the field of protein stabilization, with summary reference to a selection of 2014-2021 publications. One can enhance stability, particularly by protein engineering strategies but also by chemical modification and by other means. General protocols are set out on how to measure a given protein's (i) kinetic thermal stability and (ii) oxidative stability and (iii) how to undertake chemical modification of a protein in solution.

Bioinformatic characterization of a triacylglycerol lipase produced by Aspergillus flavus isolated from the decaying seed of Cucumeropsis mannii

Lipases are enzymes of industrial importance responsible for the hydrolysis of ester bonds of triglycerides. A lipolytic fungus was isolated and subsequently identified based on the ITS sequence analysis as putative Aspergillus flavus with accession number LC424503. The gene coding for extracellular triacylglycerol lipase was isolated from Aspergillus flavus species, sequenced, and characterised using bioinformatics tools. An open reading frame of 420 amino acid sequence was obtained and designated as Aspergillus flavus lipase (AFL) sequence. Alignment of the amino acid sequence with other lipases revealed the presence GHSLG sequence which is the lipase consensus sequence Gly-X1-Ser-X2-Gly indicating that it a classical lipase. A catalytic active site lid domain composed of TYITDTIIDLS amino acids sequence was also revealed. This lid protects the active site, control the catalytic activity and substrate selectivity in lipases. The 3-Dimensional structural model shared 34.08% sequence identity with a lipase from Yarrowia lipolytica covering 272 amino acid residues of the template model. A search of the lipase engineering database using AFL sequence revealed that it belongs to the class GX-lipase, superfamily abH23 and homologous family abH23.02, molecular weight and isoelectric point values of 46.95 KDa and 5.7, respectively. N-glycosylation sites were predicted at residues 164, 236 and 333, with potentials of 0.7250, 0.7037 and 0.7048, respectively. O-glycosylation sites were predicted at residues 355, 358, 360 and 366. A signal sequence of 37 amino acids was revealed at the N-terminal of the polypeptide. This is a short peptide sequence that marks a protein for transport across the cell membrane and indicates that AFL is an extracellular lipase. The findings on the structural and molecular properties of Aspergillus flavus lipase in this work will be crucial in future studies aiming at engineering the enzyme for biotechnology applications.Communicated by Ramaswamy H. Sarma.

Enhancing Thermostability of Lipase from Pseudomonas alcaligenes for producing L-menthol by the CREATE Strategy

The Lipase from Pseudomonas alcaligenes (PaL) catalyzes the hydrolysis of racemic menthol propionate to produce L-menthol, one of the most important flavoring agents in food, cosmetics and pharmaceuticals industries. However,...

Exogenous Enzymes as Zootechnical Additives in Animal Feed: A Review

Enzymes are widely used in the food industry. Their use as a supplement to the raw material for animal feed is a current research topic. Although there are several studies on the application of enzyme additives in the animal feed industry, it is necessary to search for new enzymes, as well as to utilize bioinformatics tools for the design of specific enzymes that work in certain environmental conditions and substrates. This will allow the improvement of the productive parameters in animals, reducing costs and making the processes more efficient. Technological needs have considered these catalysts as essential in many industrial sectors and research is constantly being carried out to optimize their use in those processes. This review describes the enzymes used in animal nutrition, their mode of action, their production and new sources of production as well as studies on different animal models to evaluate their effect on the productive performance intended for the production of animal feed.

Conjugation of haloalkane dehalogenase DhaA with arabinogalactan to increase its stability

Haloalkane dehalogenase DhaA can catalyze the hydrolytic cleavage of carbonhalogen bonds, along with production of the corresponding alcohol, a proton and a halide. However, DhaA suffers from poor environmental tolerance, such as sensitivity to high temperature, low pH and hypersaline. Arabinogalactan (AG) is a hydrophilic polysaccharide with highly branched long chains. DhaA was conjugated with AG to improve the environmental stability of DhaA in the present study. Each DhaA was averagely conjugated with 4∼5 AG molecules. Conjugation of AG essentially maintained the enzymatic activity of DhaA (91.4 %) without apparent structural alteration. The hydration layer formed by AG could reduce the solvent accessible area of DhaA and slow the protonation process, thereby improving the pH and high salt stability of DhaA. In particular, the remaining activities of the conjugate (AG-DhaA) were 35.3 % after treatment at pH4.0 for 1 h, and 80.8 % in 1 M NaCl after treatment for 16 h. As compared with DhaA, AG-DhaA showed slightly different kinetic parameters (K M of 1.90 μmol/L and k cat of 2.60 s -1).

Prospects of Using Biocatalysis for the Synthesis and Modification of Polymers

Trends in the dynamically developing application of biocatalysis for the synthesis and modification of polymers over the past 5 years are considered, with an emphasis on the production of biodegradable, biocompatible and functional polymeric materials oriented to medical applications. The possibilities of using enzymes not only as catalysts for polymerization but also for the preparation of monomers for polymerization or oligomers for block copolymerization are considered. Special attention is paid to the prospects and existing limitations of biocatalytic production of new synthetic biopolymers based on natural compounds and monomers from biomass, which can lead to a huge variety of functional biomaterials. The existing experience and perspectives for the integration of bio- and chemocatalysis in this area are discussed.

Mustguseal and Sister Web-Methods: A Practical Guide to Bioinformatic Analysis of Protein Superfamilies

Bioinformatic analysis of functionally diverse superfamilies can help to study the structure-function relationship in proteins, but represents a methodological challenge. The Mustguseal web-server can build large structure-guided sequence alignments of thousands of homologs that cover all currently available sequence variants within a common structural fold. The input to the method is a PDB code of the query protein, which represents the protein superfamily of interest. The collection and subsequent alignment of protein sequences and structures is fully automated and driven by the particular choice of parameters. Four integrated sister web-methods-the Zebra, pocketZebra, visualCMAT, and Yosshi-are available to further analyze the resulting superimposition and identify conserved, subfamily-specific, and co-evolving residues, as well as to classify and study disulfide bonds in protein superfamilies. The integration of these web-based bioinformatic tools provides an out-of-the-box easy-to-use solution, first of its kind, to study protein function and regulation and design improved enzyme variants for practical applications and selective ligands to modulate their functional properties. In this chapter, we provide a step-by-step protocol for a comprehensive bioinformatic analysis of a protein superfamily using a web-browser as the main tool and notes on selecting the appropriate values for the key algorithm parameters depending on your research objective. The web-servers are freely available to all users at https://biokinet.belozersky.msu.ru/m-platform with no login requirement.

De novo sequence redesign of a functional Ras-binding domain globally inverted the surface charge distribution and led to extreme thermostability

To acquire extremely thermostable proteins of given functions is challenging for conventional protein engineering. Here we applied ABACUS, a statistical energy function we developed for de novo amino acid sequence design, to globally redesign a Ras-binding domain (RBD), and obtained an extremely thermostable RBD that unfolds reversibly at above 110°C, the redesigned RBD experimentally confirmed to have expected structure and Ras-binding interface. Directed evolution of the redesigned RBD improved its Ras-binding affinity to the native protein level without excessive loss of thermostability. The designed amino acid substitutions were mostly at the protein surface. For many substitutions, strong epistasis or significantly differentiated effects on thermostability in the native sequence context relative to the redesigned sequence context were observed, suggesting the globally redesigned sequence to be unreachable through combining beneficial mutations of the native sequence. Further analyses revealed that by replacing 38 of a total of 48 non-interfacial surface residues at once, ABACUS redesign was able to globally "invert" the protein's charge distribution pattern in an optimized way. Our study demonstrates that computational protein design provides powerful new tools to solve challenging protein engineering problems.

Bioinformatic analysis of subfamily-specific regions in 3D-structures of homologs to study functional diversity and conformational plasticity in protein superfamilies

Local 3D-structural differences in homologous proteins contribute to functional diversity observed in a superfamily, but so far received little attention as bioinformatic analysis was usually carried out at the level of amino acid sequences. We have developed Zebra3D - the first-of-its-kind bioinformatic software for systematic analysis of 3D-alignments of protein families using machine learning. The new tool identifies subfamily-specific regions (SSRs) - patterns of local 3D-structure (i.e. single residues, loops, or secondary structure fragments) that are spatially equivalent within families/subfamilies, but are different among them, and thus can be associated with functional diversity and function-related conformational plasticity. Bioinformatic analysis of protein superfamilies by Zebra3D can be used to study 3D-determinants of catalytic activity and specific accommodation of ligands, help to prepare focused libraries for directed evolution or assist development of chimeric enzymes with novel properties by exchange of equivalent regions between homologs, and to characterize plasticity in binding sites. A companion Mustguseal web-server is available to automatically construct a 3D-alignment of functionally diverse proteins, thus reducing the minimal input required to operate Zebra3D to a single PDB code. The Zebra3D + Mustguseal combined approach provides the opportunity to systematically explore the value of SSRs in superfamilies and to use this information for protein design and drug discovery. The software is available open-access at https://biokinet.belozersky.msu.ru/Zebra3D.

Improved thermostability of proteinase K and recognizing the synergistic effect of Rosetta and FoldX approaches

Proteinase K (PRK) is a proteolytic enzyme that has been widely used in industrial applications. However, poor stability has severely limited the uses of PRK. In this work, we used two structure-guided rational design methods, Rosetta and FoldX, to modify PRK thermostability. Fifty-two single amino acid conversion mutants were constructed based on software predictions of residues that could affect protein stability. Experimental characterization revealed that 46% (21 mutants) exhibited enhanced thermostability. The top four variants, D260V, T4Y, S216Q, and S219Q, showed improved half-lives at 69°C by 12.4-, 2.6-, 2.3-, and 2.2-fold that of the parent enzyme, respectively. We also found that selecting mutations predicted by both methods could increase the predictive accuracy over that of either method alone, with 73% of the shared predicted mutations resulting in higher thermostability. In addition to providing promising new variants of PRK in industrial applications, our findings also show that combining these programs may synergistically improve their predictive accuracy.

Power of Biocatalysis for Organic Synthesis

Biocatalysis, using defined enzymes for organic transformations, has become a common tool in organic synthesis, which is also frequently applied in industry. The generally high activity and outstanding stereo-, regio-, and chemoselectivity observed in many biotransformations are the result of a precise control of the reaction in the active site of the biocatalyst. This control is achieved by exact positioning of the reagents relative to each other in a fine-tuned 3D environment, by specific activating interactions between reagents and the protein, and by subtle movements of the catalyst. Enzyme engineering enables one to adapt the catalyst to the desired reaction and process. A well-filled biocatalytic toolbox is ready to be used for various reactions. Providing nonnatural reagents and conditions and evolving biocatalysts enables one to play with the myriad of options for creating novel transformations and thereby opening new, short pathways to desired target molecules. Combining several biocatalysts in one pot to perform several reactions concurrently increases the efficiency of biocatalysis even further.

Zebra2: advanced and easy-to-use web-server for bioinformatic analysis of subfamily-specific and conserved positions in diverse protein superfamilies

Zebra2 is a highly automated web-tool to search for subfamily-specific and conserved positions (i.e. the determinants of functional diversity as well as the key catalytic and structural residues) in protein superfamilies. The bioinformatic analysis is facilitated by Mustguseal—a companion web-server to automatically collect and superimpose a large representative set of functionally diverse homologs with high structure similarity but low sequence identity to the selected query protein. The results are automatically prioritized and provided at four information levels to facilitate the knowledge-driven expert selection of the most promising positions on-line: as a sequence similarity network; interfaces to sequence-based and 3D-structure-based analysis of conservation and variability; and accompanied by the detailed annotation of proteins accumulated from the integrated databases with links to the external resources. The integration of Zebra2 and Mustguseal web-tools provides the first of its kind out-of-the-box open-access solution to conduct a systematic analysis of evolutionarily related proteins implementing different functions within a shared 3D-structure of the superfamily, determine common and specific patterns of function-associated local structural elements, assist to select hot-spots for rational design and to prepare focused libraries for directed evolution. The web-servers are free and open to all users at https://biokinet.belozersky.msu.ru/zebra2, no login required.

Mapping enzyme-substrate interactions: its potential to study the mechanism of enzymes

With the increase of the need to use more sustainable processes for the industry in our society, the modeling of enzymes has become crucial to fully comprehend their mechanism of action and use this knowledge to enhance and design their properties. A lot of methods to study enzymes computationally exist and they have been classified on sequence-based, structure-based, and the more new artificial intelligence-based ones. Albeit the abundance of methods to help predict the function of an enzyme, molecular modeling is crucial when trying to understand the enzyme mechanism, as they aim to correlate atomistic information with experimental data. Among them, methods that simulate the system dynamics at a molecular mechanics level of theory (classical force fields) have shown to offer a comprehensive study. In this book chapter, we will analyze these techniques, emphasizing the importance of precise modeling of enzyme-substrate interactions. In the end, a brief explanation of the transference of the information from research studies to the industry is given accompanied with two examples of family enzymes where their modeling has helped their exploitation.

EasyAmber: A comprehensive toolbox to automate the molecular dynamics simulation of proteins

Conformational plasticity of the functionally important regions and binding sites in protein/enzyme structures is one of the key factors affecting their function and interaction with substrates/ligands. Molecular dynamics (MD) can address the challenge of accounting for protein flexibility by predicting the time-dependent behavior of a molecular system. It has a potential of becoming a particularly important tool in protein engineering and drug discovery, but requires specialized training and skills, what impedes practical use by many investigators. We have developed the easyAmber - a comprehensive set of programs to automate the molecular dynamics routines implemented in the Amber package. The toolbox can address a wide set of tasks in computational biology struggling to account for protein flexibility. The automated workflow includes a complete set of steps from the initial "static" molecular model to the MD "production run": the full-atom model building, optimization/equilibration of the molecular system, classical/conventional and accelerated molecular dynamics simulations. The easyAmber implements advanced MD protocols, but is highly automated and easy-to-operate to attract a broad audience. The toolbox can be used on a personal desktop station equipped with a compatible gaming GPU-accelerator, as well as help to manage huge workloads on a powerful supercomputer. The software provides an opportunity to operate multiple simulations of different proteins at the same time, thus significantly increasing work efficiency. The easyAmber takes the molecular dynamics to the next level in terms of usability for complex processing of large volumes of data, thus supporting the recent trend away from inefficient "static" approaches in biology toward a deeper understanding of the dynamics in protein structures. The software is freely available for download at https://biokinet.belozersky.msu.ru/easyAmber, no login required.

Comparative bioinformatic and structural analyses of pepsin and renin

Pepsin, the archetypal pepsin-like aspartic protease, is irreversibly denatured when exposed to neutral pH conditions whereas renin, a structural homologue of pepsin, is fully stable and optimally active in the same conditions despite sharing highly similar enzyme architecture. To gain insight into the structural determinants of differential aspartic protease pH stability, the present study used comparative bioinformatic and structural analyses. In pepsin, an abundance of polar and aspartic acid residues were identified, a common trait with other acid-stable enzymes. Conversely, renin was shown to have increased levels of basic amino acids. In both pepsin and renin, the solvent exposure of these charged groups was high. Having similar overall acidic residue content, the solvent-exposed basic residues may allow for extensive salt bridge formation in renin, whereas in pepsin, these residues are protonated and serve to form stabilizing hydrogen bonds at low pH. Relative differences in structure and sequence in the turn and joint regions of the β-barrel and ψ-loop in both the N- and C-terminal lobes were identified as regions of interest in defining divergent pH stability. Compared to the structural rigidity of renin, pepsin has more instability associated with the N-terminus, specifically the B/C connector. By contrast, renin exhibits greater C-terminal instability in turn and connector regions. Overall, flexibility differences in connector regions, and amino acid composition, particularly in turn and joint regions of the β-barrel and ψ-loops, likely play defining roles in determining pH stability for renin and pepsin.

Yosshi: a web-server for disulfide engineering by bioinformatic analysis of diverse protein families

Disulfide bonds play a significant role in protein stability, function or regulation but are poorly conserved among evolutionarily related proteins. The Yosshi can help to understand the role of S–S bonds by comparing sequences and structures of homologs with diverse properties and different disulfide connectivity patterns within a common structural fold of a superfamily, and assist to select the most promising hot-spots to improve stability of proteins/enzymes or modulate their functions by introducing naturally occurring crosslinks. The bioinformatic analysis is supported by the integrated Mustguseal web-server to construct large structure-guided sequence alignments of functionally diverse protein families that can include thousands of proteins based on all available information in public databases. The Yosshi+Mustguseal is a new integrated web-tool for a systematic homology-driven analysis and engineering of S–S bonds that facilitates a broader interpretation of disulfides not just as a factor of structural stability, but rather as a mechanism to implement functional diversity within a superfamily. The results can be downloaded as a content-rich PyMol session file or further studied online using the HTML5-based interactive analysis tools. Both web-servers are free and open to all users at https://biokinet.belozersky.msu.ru/yosshi and there is no login requirement.

Common mechanism of thermostability in small α- and β-proteins studied by molecular dynamics

Protein thermostability is important to evolution, diseases, and industrial applications. Proteins use diverse molecular strategies to achieve stability at high temperature, yet reducing the entropy of unfolding seems required. We investigated five small α-proteins and five β-proteins with known, distinct structures and thermostability (T_m ) using multi-seed molecular dynamics simulations at 300, 350, and 400 K. The proteins displayed diverse changes in hydrogen bonding, solvent exposure, and secondary structure with no simple relationship to T_m . Our dynamics were in good agreement with experimental B-factors at 300 K and insensitive to force-field choice. Despite the very distinct structures, the native-state (300 + 350 K) free-energy landscapes (FELs) were significantly broader for the two most thermostable proteins and smallest for the three least stable proteins in both the α- and β-group and with both force fields studied independently (tailed t-test, 95% confidence level). Our results suggest that entropic ensembles stabilize proteins at high temperature due to reduced entropy of unfolding, viz., ΔG = ΔH - TΔS. Supporting this mechanism, the most thermostable proteins were also the least kinetically stable, consistent with broader FELs, typified by villin headpiece and confirmed by specific comparison to a mesophilic ortholog of Thermus thermophilus apo-pyrophosphate phosphohydrolase. We propose that molecular strategies of protein thermostabilization, although diverse, tend to converge toward highest possible entropy in the native state consistent with the functional requirements. We speculate that this tendency may explain why many proteins are not optimally structured and why molten-globule states resemble native proteins so much.

Highly thermostable carboxylic acid reductases generated by ancestral sequence reconstruction

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

Dextran Aldehyde in Biocatalysis: More Than a Mere Immobilization System

Dextran aldehyde (dexOx), resulting from the periodate oxidative cleavage of 1,2-diol moiety inside dextran, is a polymer that is very useful in many areas, including as a macromolecular carrier for drug delivery and other biomedical applications. In particular, it has been widely used for chemical engineering of enzymes, with the aim of designing better biocatalysts that possess improved catalytic properties, making them more stable and/or active for different catalytic reactions. This polymer possesses a very flexible hydrophilic structure, which becomes inert after chemical reduction; therefore, dexOx comes to be highly versatile in a biocatalyst design. This paper presents an overview of the multiple applications of dexOx in applied biocatalysis, e.g., to modulate the adsorption of biomolecules on carrier surfaces in affinity chromatography and biosensors design, to serve as a spacer arm between a ligand and the support in biomacromolecule immobilization procedures or to generate artificial microenvironments around the enzyme molecules or to stabilize multimeric enzymes by intersubunit crosslinking, among many other applications.

Mustguseal: a server for multiple structure-guided sequence alignment of protein families

Motivation

Comparative analysis of homologous proteins in a functionally diverse superfamily is a valuable tool at studying structure-function relationship, but represents a methodological challenge.

Results

The Mustguseal web-server can automatically build large structure-guided sequence alignments of functionally diverse protein families that include thousands of proteins basing on all available information about their structures and sequences in public databases. Superimposition of protein structures is implemented to compare evolutionarily distant relatives, whereas alignment of sequences is used to compare close homologues. The final alignment can be downloaded for a local use or operated on-line with the built-in interactive tools and further submitted to the integrated sister web-servers of Mustguseal to analyze conserved, subfamily-specific and co-evolving residues at studying a protein function and regulation, designing improved enzyme variants for practical applications and selective ligands to modulate functional properties of proteins.

Availability and implementation

Freely available on the web at https://biokinet.belozersky.msu.ru/mustguseal.

Contact

vytas@belozersky.msu.ru.

Supplementary information

Supplementary data are available at Bioinformatics online.

Separation, biochemical characterization and salt-tolerant mechanisms of alkaline protease from Aspergillus oryzae

BACKGROUND

The salt tolerance of proteases secreted by Aspergillus oryzae 3.042 closely relates to the utilization of raw materials and the quality of soy sauce. However, little is known about the salt-tolerant proteases and their salt-tolerant mechanisms.

RESULTS

In this study, we isolated and identified a salt-tolerant alkaline protease (AP, approximately 29 kDa) produced by A. oryzae 3.042. It was considered as a metal-ion-independent serine protease. The optimum and stable pH values were both pH 9.0 and the optimum temperature was 40 °C. Over 20% relative activity of AP remained in the presence of 3.0 mol L^-1 NaCl after 7 days, but its K_m and V_max were only mildly influenced by the presence of 3.0 mol L^-1 NaCl, indicating its outstanding salt tolerance. Furthermore, AP was more stable than non-salt-tolerant protease at high salinity. The salt-tolerant mechanisms of AP could be due to more salt bridges, higher proportion of ordered secondary structures and stronger hydrophobic amino acid residues in the interior.

CONCLUSIONS

The above results are vital for maintaining, activating and/or modulating the activity of AP in high-salt environments. They would also provide theoretical guidance for the modification of AP and the engineering of A. oryzae 3.042 so as to secrete more AP. © 2018 Society of Chemical Industry.

Coupled molecular dynamics mediate long- and short-range epistasis between mutations that affect stability and aggregation kinetics

Multiple mutations are typically required to significantly improve protein stability or aggregation kinetics. However, when several substitutions are made in a single protein, the mutations can potentially interact in a nonadditive manner, resulting in epistatic effects, which can hamper protein-engineering strategies to improve thermostability or aggregation kinetics. Here, we have examined the role of protein dynamics in mediating epistasis between pairs of mutations. With Escherichia coli transketolase (TK) as a model, we explored the epistatic interactions between two single variants H192P and A282P, and also between the double-mutant H192P/A282P and two single variants, I365L or G506A. Epistasis was determined for several measures of protein stability, including the following: the free-energy barrier to kinetic inactivation, ∆∆G ^‡; thermal transition midpoint temperatures, T _m; and aggregation onset temperatures, T _agg Nonadditive epistasis was observed between neighboring mutations as expected, but also for distant mutations located in the surface and core regions of different domains. Surprisingly, the epistatic behaviors for each measure of stability were often different for any given pairwise recombination, highlighting that kinetic and thermodynamic stabilities do not always depend on the same structural features. Molecular-dynamics simulations and a pairwise cross-correlation analysis revealed that mutations influence the dynamics of their local environment, but also in some cases the dynamics of regions distant in the structure. This effect was found to mediate epistatic interactions between distant mutations and could therefore be exploited in future protein-engineering strategies.

Bioinformatic analysis of the fold type I PLP-dependent enzymes reveals determinants of reaction specificity in l-threonine aldolase from Aeromonas jandaei

Understanding the role of specific amino acid residues in the molecular mechanism of a protein's function is one of the most challenging problems in modern biology. A systematic bioinformatic analysis of protein families and superfamilies can help in the study of structure–function relationships and in the design of improved variants of enzymes/proteins, but represents a methodological challenge. The pyridoxal‐5′‐phosphate (PLP)‐dependent enzymes are catalytically diverse and include the aspartate aminotransferase superfamily which implements a common structural framework known as type fold I. In this work, the recently developed bioinformatic online methods Mustguseal and Zebra were used to collect and study a large representative set of the aspartate aminotransferase superfamily with high structural, but low sequence similarity to l‐threonine aldolase from Aeromonas jandaei (LTAaj), in order to identify conserved positions that provide general properties in the superfamily, and to reveal family‐specific positions (FSPs) responsible for functional diversity. The roles of the identified residues in the catalytic mechanism and reaction specificity of LTAaj were then studied by experimental site‐directed mutagenesis and molecular modelling. It was shown that FSPs determine reaction specificity by coordinating the PLP cofactor in the enzyme's active centre, thus influencing its activation and the tautomeric equilibrium of the intermediates, which can be used as hotspots to modulate the protein's functional properties. Mutagenesis at the selected FSPs in LTAaj led to a reduction in a native catalytic activity and increased the rate of promiscuous reactions. The results provide insight into the structural basis of catalytic promiscuity of the PLP‐dependent enzymes and demonstrate the potential of bioinformatic analysis in studying structure–function relationship in protein superfamilies.

Synthetic Biology to Improve the Production of Lipases and Esterases (Review)

Synthetic biology is an emergent field of research whose aim is to make biology an engineering discipline, thus permitting to design, control, and standardize biological processes. Synthetic biology is therefore expected to boost the development of biotechnological processes such as protein production and enzyme engineering, which can be significantly relevant for lipases and esterases.

The visualCMAT: A web-server to select and interpret correlated mutations/co-evolving residues in protein families

The visualCMAT web-server was designed to assist experimental research in the fields of protein/enzyme biochemistry, protein engineering, and drug discovery by providing an intuitive and easy-to-use interface to the analysis of correlated mutations/co-evolving residues. Sequence and structural information describing homologous proteins are used to predict correlated substitutions by the Mutual information-based CMAT approach, classify them into spatially close co-evolving pairs, which either form a direct physical contact or interact with the same ligand (e.g. a substrate or a crystallographic water molecule), and long-range correlations, annotate and rank binding sites on the protein surface by the presence of statistically significant co-evolving positions. The results of the visualCMAT are organized for a convenient visual analysis and can be downloaded to a local computer as a content-rich all-in-one PyMol session file with multiple layers of annotation corresponding to bioinformatic, statistical and structural analyses of the predicted co-evolution, or further studied online using the built-in interactive analysis tools. The online interactivity is implemented in HTML5 and therefore neither plugins nor Java are required. The visualCMAT web-server is integrated with the Mustguseal web-server capable of constructing large structure-guided sequence alignments of protein families and superfamilies using all available information about their structures and sequences in public databases. The visualCMAT web-server can be used to understand the relationship between structure and function in proteins, implemented at selecting hotspots and compensatory mutations for rational design and directed evolution experiments to produce novel enzymes with improved properties, and employed at studying the mechanism of selective ligand's binding and allosteric communication between topologically independent sites in protein structures. The web-server is freely available at https://biokinet.belozersky.msu.ru/visualcmat and there are no login requirements.

Protein Stability: Enhancement and Measurement

This article defines protein stability, emphasizes its importance and surveys the field of protein stabilization, with summary reference to a selection of 2009-2015 publications. One can enhance stability by, in particular, protein engineering strategies and by chemical modification (including conjugation) in solution. General protocols are set out on how to measure a given protein's (1) kinetic thermal stability, and (2) oxidative stability, and (3) how to undertake chemical modification of a protein in solution.

Creating an Efficient Methanol-Stable Biocatalyst by Protein and Immobilization Engineering Steps towards Efficient Biosynthesis of Biodiesel

Two ternary sol-gel matrices, an octyltriethoxysilane-based aliphatic matrix and a phenyltriethoxysilane (PTEOS)-based aromatic matrix, were used to immobilize a methanol-stable variant of lipase from Geobacillus stearothermophilus T6 for the synthesis of biodiesel from waste oil. Superior thermal stability of the mutant versus the wildtype in methanol was confirmed by intrinsic protein fluorescence measurements. The influence of skim milk and soluble E. coli lysate proteins as bulking and stabilizing agents in conjunction with sol-gel entrapment were investigated. E. coli lysate proteins were better stabilizing agents of the purified lipase mutant than skim milk, as evidenced by reverse engineering of the aromatic-based system. This was also shown for commercial Candida antarctica lipase B (CaLB) and Thermomyces lanuginosus lipase (TLL). Uniform, dense, and nonaggregated particles imaged by scanning electron microscopy and a small particle size of 13 μm pertaining to the system comprising PTEOS and E. coli lysate proteins correlated well with high esterification activity. Combining protein and immobilization engineering resulted in a durable biocatalyst with efficient recycling ability and high biodiesel conversion rates.

Sequence homolog-based molecular engineering for shifting the enzymatic pH optimum

Cell-free synthetic biology system organizes multiple enzymes (parts) from different sources to implement unnatural catalytic functions. Highly adaption between the catalytic parts is crucial for building up efficient artificial biosynthetic systems. Protein engineering is a powerful technology to tailor various enzymatic properties including catalytic efficiency, substrate specificity, temperature adaptation and even achieve new catalytic functions. However, altering enzymatic pH optimum still remains a challenging task. In this study, we proposed a novel sequence homolog-based protein engineering strategy for shifting the enzymatic pH optimum based on statistical analyses of sequence-function relationship data of enzyme family. By two statistical procedures, artificial neural networks (ANNs) and least absolute shrinkage and selection operator (Lasso), five amino acids in GH11 xylanase family were identified to be related to the evolution of enzymatic pH optimum. Site-directed mutagenesis of a thermophilic xylanase from Caldicellulosiruptor bescii revealed that four out of five mutations could alter the enzymatic pH optima toward acidic condition without compromising the catalytic activity and thermostability. Combination of the positive mutants resulted in the best mutant M31 that decreased its pH optimum for 1.5 units and showed increased catalytic activity at pH < 5.0 compared to the wild-type enzyme. Structure analysis revealed that all the mutations are distant from the active center, which may be difficult to be identified by conventional rational design strategy. Interestingly, the four mutation sites are clustered at a certain region of the enzyme, suggesting a potential "hot zone" for regulating the pH optima of xylanases. This study provides an efficient method of modulating enzymatic pH optima based on statistical sequence analyses, which can facilitate the design and optimization of suitable catalytic parts for the construction of complicated cell-free synthetic biology systems.

Parallel workflow manager for non-parallel bioinformatic applications to solve large-scale biological problems on a supercomputer

Rapid expansion of online resources providing access to genomic, structural, and functional information associated with biological macromolecules opens an opportunity to gain a deeper understanding of the mechanisms of biological processes due to systematic analysis of large datasets. This, however, requires novel strategies to optimally utilize computer processing power. Some methods in bioinformatics and molecular modeling require extensive computational resources. Other algorithms have fast implementations which take at most several hours to analyze a common input on a modern desktop station, however, due to multiple invocations for a large number of subtasks the full task requires a significant computing power. Therefore, an efficient computational solution to large-scale biological problems requires both a wise parallel implementation of resource-hungry methods as well as a smart workflow to manage multiple invocations of relatively fast algorithms. In this work, a new computer software mpiWrapper has been developed to accommodate non-parallel implementations of scientific algorithms within the parallel supercomputing environment. The Message Passing Interface has been implemented to exchange information between nodes. Two specialized threads — one for task management and communication, and another for subtask execution — are invoked on each processing unit to avoid deadlock while using blocking calls to MPI. The mpiWrapper can be used to launch all conventional Linux applications without the need to modify their original source codes and supports resubmission of subtasks on node failure. We show that this approach can be used to process huge amounts of biological data efficiently by running non-parallel programs in parallel mode on a supercomputer. The C++ source code and documentation are available from http://biokinet.belozersky.msu.ru/mpiWrapper .

A single residue substitution accounts for the significant difference in thermostability between two isoforms of human cytosolic creatine kinase

Creatine kinase (CK) helps maintain homeostasis of intracellular ATP level by catalyzing the reversible phosphotransfer between ATP and phosphocreatine. In humans, there are two cytosolic CK isoforms, the muscle-type (M) and the brain-type (B), which frequently function as homodimers (hMMCK and hBBCK). Interestingly, these isoenzymes exhibit significantly different thermostabilities, despite high similarity in amino acid sequences and tertiary structures. In order to investigate the mechanism of this phenomenon, in this work, we first used domain swapping and site-directed mutagenesis to search for the key residues responsible for the isoenzyme-specific thermostability. Strikingly, the difference in thermostability was found to principally arise from one single residue substitution at position 36 (Pro in hBBCK vs. Leu in hMMCK). We then engaged the molecular dynamics simulations to study the molecular mechanism. The calculations imply that the P36L substitution introduces additional local interactions around residue 36 and thus further stabilizes the dimer interface through a complex interaction network, which rationalizes the observation that hMMCK is more resistant to thermal inactivation than hBBCK. We finally confirmed this molecular explanation through thermal inactivation assays on Asp36 mutants that were proposed to devastate the local interactions and thus the dimer associations in both isoenzymes.

Sequence and structure-based comparative analysis to assess, identify and improve the thermostability of penicillin G acylases

Penicillin acylases are enzymes employed by the pharmaceutical industry for the manufacture of semi-synthetic penicillins. There is a continuous demand for thermostable and alkalophilic enzymes in such applications. We have carried out a computational analysis of known penicillin G acylases (PGAs) in terms of their thermostable nature using various protein-stabilizing factors. While the presence of disulfide bridges was considered initially to screen putative thermostable PGAs from the database, various other factors such as high arginine to lysine ratio, less content of thermolabile amino acids, presence of proline in β-turns, more number of ion-pair and other non-bonded interactions were also considered for comparison. A modified consensus approach designed could further identify stabilizing residue positions by site-specific comparison between mesostable and thermostable PGAs. A most likely thermostable enzyme identified from the analysis was PGA from Paracoccus denitrificans (PdPGA). This was cloned, expressed and tested for its thermostable nature using biochemical and biophysical experiments. The consensus site-specific sequence-based approach predicted PdPGA to be more thermostable than Escherichia coli PGA, but not as thermostable as the PGA from Achromobacter xylosoxidans. Experimental data showed that PdPGA was comparatively less thermostable than Achromobacter xylosoxidans PGA, although thermostability factors favored a much higher stability. Despite being mesostable, PdPGA being active and stable at alkaline pH is an advantage. Finally, several residue positions could be identified in PdPGA, which upon mutation selectively could improve the thermostability of the enzyme.