Structural Rigidity and Protein Thermostability in Variants of Lipase A from Bacillus subtilis

Computer-Aided Reconstruction and Application of Bacillus halodurans S7 Xylanase with Heat and Alkali Resistance

β-1,4-Endoxylanase is the most critical hydrolase for xylan degradation during lignocellulosic biomass utilization. However, its poor stability and activity in hot and alkaline environments hinder its widespread application. In this study, BhS7Xyl from Bacillus halodurans S7 was improved using a computer-aided design through isothermal compressibility (βT) perturbation engineering and by combining three thermostability prediction algorithms (ICPE-TPA). The best variant with remarkable improvement in specific activity, heat resistance (70 °C), and alkaline resistance (both pH 9.0 and 70 °C), R69F/E137M/E145L, exhibited a 4.9-fold increase by wild-type in specific activity (1368.6 U/mg), a 39.4-fold increase in temperature half-life (458.1 min), and a 57.6-fold increase in pH half-life (383.1 min). Furthermore, R69F/E137M/E145L was applied to the hydrolysis of agricultural waste (corncob and hardwood pulp) to efficiently obtain a higher yield of high-value xylooligosaccharides. Overall, the ICPE-TPA strategy has the potential to improve the functional performance of enzymes under extreme conditions for the high-value utilization of lignocellulosic biomass.

Lipase A from Bacillus subtilis: Substrate Binding, Conformational Dynamics, and Signatures of a Lid

Protein-ligand binding studies are crucial for understanding the molecular basis of biological processes and for further advancing industrial biocatalysis and drug discovery. Using computational modeling and molecular dynamics simulations, we investigated the binding of a butyrate ester substrate to the lipase A (LipA) enzyme of Bacillus subtilis. Besides obtaining a close agreement of the binding free energy with the experimental value, the study reveals a remarkable reorganization of the catalytic triad upon substrate binding, leading to increased essential hydrogen bond populations. The investigation shows the distortion of the oxyanion hole in both the substrate-bound and unbound states of LipA and highlights the strengthening of the same in the tetrahedral intermediate complex. Principal component analysis of the unbound ensemble reveals the dominant motion in LipA to be the movement of Loop-1 (Tyr129-Arg142) between two states that cover and uncover the active site, mirroring that of a lid prevalent in several lipases. This lid-like motion of Loop-1 is also supported by its tendency to spontaneously open up at an oil-water interface. Overall, this study provides valuable insights into the impact of substrate binding on the structure, flexibility, and conformational dynamics of the LipA enzyme.

A thermostable organic solvent-tolerant lipase from Brevibacillus sp.: production and integrated downstream processing using an alcohol-salt-based aqueous two-phase system

Lipases are used for the synthesis of different compounds in the chemical, pharmaceutical, and food industries. Most of the reactions are carried out in non-aqueous media and often at elevated temperature, requiring the use of organic solvent-tolerant thermostable lipases. However, most known lipases are not stable in the presence of organic solvents and at elevated temperature. In this study, an organic solvent-tolerant thermostable lipase was obtained from Brevibacillus sp. SHI-160, a moderate thermophile isolated from a hot spring in the East African Rift Valley. The enzyme was optimally active at 65°C and retained over 90% of its activity after 1 h of incubation at 70°C. High lipase activity was measured in the pH range of 6.5 to 9.0 with an optimum pH of 8.5. The enzyme was stable in the presence of both polar and non-polar organic solvents. The stability of the enzyme in the presence of polar organic solvents allowed the development of an efficient downstream processing using an alcohol-salt-based aqueous two-phase system (ATPS). Thus, in the presence of 2% salt, over 98% of the enzyme partitioned to the alcohol phase. The ATPS-recovered enzyme was directly immobilized on a solid support through adsorption and successfully used to catalyze a transesterification reaction between paranitrophenyl palmitate and short-chain alcohols in non-aqueous media. This shows the potential of lipase SHI-160 to catalyze reactions in non-aqueous media for the synthesis of valuable compounds. The integrated approach developed for enzyme production and cheap and efficient downstream processing using ATPS could allow a significant reduction in enzyme production costs. The results also show the potential of extreme environments in the East African Rift Valley as sources of valuable microbial genetic resources for the isolation of novel lipases and other industrially important enzymes.

Enzyme adaptation to habitat thermal legacy shapes the thermal plasticity of marine microbiomes

Microbial communities respond to temperature with physiological adaptation and compositional turnover. Whether thermal selection of enzymes explains marine microbiome plasticity in response to temperature remains unresolved. By quantifying the thermal behaviour of seven functionally-independent enzyme classes (esterase, extradiol dioxygenase, phosphatase, beta-galactosidase, nuclease, transaminase, and aldo-keto reductase) in native proteomes of marine sediment microbiomes from the Irish Sea to the southern Red Sea, we record a significant effect of the mean annual temperature (MAT) on enzyme response in all cases. Activity and stability profiles of 228 esterases and 5 extradiol dioxygenases from sediment and seawater across 70 locations worldwide validate this thermal pattern. Modelling the esterase phase transition temperature as a measure of structural flexibility confirms the observed relationship with MAT. Furthermore, when considering temperature variability in sites with non-significantly different MATs, the broadest range of enzyme thermal behaviour and the highest growth plasticity of the enriched heterotrophic bacteria occur in samples with the widest annual thermal variability. These results indicate that temperature-driven enzyme selection shapes microbiome thermal plasticity and that thermal variability finely tunes such processes and should be considered alongside MAT in forecasting microbial community thermal response.

Altering the Regioselectivity of T1 Lipase from Geobacillus zalihae toward sn-3 Acylglycerol Using a Rational Design Approach

The regioselectivity characteristic of lipases facilitate a wide range of novel molecule unit constructions and fat modifications. Lipases can be categorized as sn-1,3, sn-2, and random regiospecific. Geobacillus zalihae T1 lipase catalyzes the hydrolysis of the sn-1,3 acylglycerol chain. The T1 lipase structural analysis shows that the oxyanion hole F16 and its lid domain undergo structural rearrangement upon activation. Site-directed mutagenesis was performed by substituting the lid domain residues (F180G and F181S) and the oxyanion hole residue (F16W) in order to study their effects on the structural changes and regioselectivity. The novel lipase mutant 3M switches the regioselectivity from sn-1,3 to only sn-3. The mutant 3M shifts the optimum pH to 10, alters selectivity toward p-nitrophenyl ester selectivity to C14-C18, and maintains a similar catalytic efficiency of 518.4 × 10−6 (s−1/mM). The secondary structure of 3M lipase comprises 15.8% and 26.3% of the α-helix and β-sheet, respectively, with a predicted melting temperature (Tm) value of 67.8 °C. The in silico analysis was conducted to reveal the structural changes caused by the F180G/F181S/F16W mutations in blocking the binding of the sn-1 acylglycerol chain and orientating the substrate to bond to the sn-3 acylglycerol, which resulted in switching the T1 lipase regioselectivity.

Structural features, temperature adaptation and industrial applications of microbial lipases from psychrophilic, mesophilic and thermophilic origins

Microbial lipases are very prominent biocatalysts because of their ability to catalyze a wide variety of reactions in aqueous and non-aqueous media. Here microbial lipases from different origins (psychrophiles, mesophiles, and thermophiles) have been reviewed. This review emphasizes an update of structural diversity in temperature adaptation and industrial applications, of psychrophilic, mesophilic, and thermophilic lipases. The microbial origins of lipases are logically dynamic, proficient, and also have an extensive range of industrial uses with the manufacturing of altered molecules. It is therefore of interest to understand the molecular mechanisms of adaptation to temperature in occurring lipases. However, lipases from extremophiles (psychrophiles, and thermophiles) are widely used to design biotransformation reactions with higher yields, fewer byproducts, or useful side products and have been predicted to catalyze those reactions also, which otherwise are not possible with the mesophilic lipases. Lipases as a multipurpose biological catalyst have given a favorable vision in meeting the needs of several industries such as biodiesel, foods, and drinks, leather, textile, detergents, pharmaceuticals, and medicals.

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.

Met/Val129 polymorphism of the full-length human prion protein dictates distinct pathways of amyloid formation

Methionine/valine polymorphism at position 129 of the human prion protein, huPrP, is tightly associated with the pathogenic phenotype, disease progress, and age of onset of neurodegenerative diseases such as Creutzfeldt–Jakob disease or Fatal Familial Insomnia. This raises the question of whether and how the amino acid type at position 129 influences the structural properties of huPrP, affecting its folding, stability, and amyloid formation behavior. Here, our detailed biophysical characterization of the 129M and 129V variants of recombinant full-length huPrP(23–230) by amyloid formation kinetics, CD spectroscopy, molecular dynamics simulations, and sedimentation velocity analysis reveals differences in their aggregation propensity and oligomer content, leading to deviating pathways for the conversion into amyloid at acidic pH. We determined that the 129M variant exhibits less secondary structure content before amyloid formation and higher resistance to thermal denaturation compared to the 129V variant, whereas the amyloid conformation of both variants shows similar thermal stability. Additionally, our molecular dynamics simulations and rigidity analyses at the atomistic level identify intramolecular interactions responsible for the enhanced monomer stability of the 129M variant, involving more frequent minimum distances between E196 and R156, forming a salt bridge. Removal of the N-terminal half of the 129M full-length variant diminishes its differences compared to the 129V full-length variant and highlights the relevance of the flexible N terminus in huPrP. Taken together, our findings provide insight into structural properties of huPrP and the effects of the amino acid identity at position 129 on amyloid formation behavior.

Uncoupling of Voltage- and Ligand-Induced Activation in HCN2 Channels by Glycine Inserts

Hyperpolarization-activated cyclic nucleotide-modulated (HCN) channels are tetramers that generate electrical rhythmicity in special brain neurons and cardiomyocytes. The channels are activated by membrane hyperpolarization. The binding of cAMP to the four available cyclic nucleotide-binding domains (CNBD) enhances channel activation. We analyzed in the present study the mechanism of how the effect of cAMP binding is transmitted to the pore domain. Our strategy was to uncouple the C-linker (CL) from the channel core by inserting one to five glycine residues between the S6 gate and the A′-helix (constructs 1G to 5G). We quantified in full-length HCN2 channels the resulting functional effects of the inserted glycines by current activation as well as the structural dynamics and statics using molecular dynamics simulations and Constraint Network Analysis. We show functionally that already in 1G the cAMP effect on activation is lost and that with the exception of 3G and 5G the concentration-activation relationships are shifted to depolarized voltages with respect to HCN2. The strongest effect was found for 4G. Accordingly, the activation kinetics were accelerated by all constructs, again with the strongest effect in 4G. The simulations reveal that the average residue mobility of the CL and CNBD domains is increased in all constructs and that the junction between the S6 and A′-helix is turned into a flexible hinge, resulting in a destabilized gate in all constructs. Moreover, for 3G and 4G, there is a stronger downward displacement of the CL-CNBD than in HCN2 and the other constructs, resulting in an increased kink angle between S6 and A′-helix, which in turn loosens contacts between the S4-helix and the CL. This is suggested to promote a downward movement of the S4-helix, similar to the effect of hyperpolarization. In addition, exclusively in 4G, the selectivity filter in the upper pore region and parts of the S4-helix are destabilized. The results provide new insights into the intricate activation of HCN2 channels.

Structure and Function of Redox-Sensitive Superfolder Green Fluorescent Protein Variant

Aims: Genetically encoded green fluorescent protein (GFP)-based redox biosensors are widely used to monitor specific and dynamic redox processes in living cells. Over the last few years, various biosensors for a variety of applications were engineered and enhanced to match the organism and cellular environments, which should be investigated. In this context, the unicellular intraerythrocytic parasite Plasmodium, the causative agent of malaria, represents a challenge, as the small size of the organism results in weak fluorescence signals that complicate precise measurements, especially for cell compartment-specific observations. To address this, we have functionally and structurally characterized an enhanced redox biosensor superfolder roGFP2 (sfroGFP2). Results: SfroGFP2 retains roGFP2-like behavior, yet with improved fluorescence intensity (FI) in cellulo. SfroGFP2-based redox biosensors are pH insensitive in a physiological pH range and show midpoint potentials comparable with roGFP2-based redox biosensors. Using crystallography and rigidity theory, we identified the superfolding mutations as being responsible for improved structural stability of the biosensor in a redox-sensitive environment, thus explaining the improved FI in cellulo. Innovation: This work provides insight into the structure and function of GFP-based redox biosensors. It describes an improved redox biosensor (sfroGFP2) suitable for measuring oxidizing effects within small cells where applicability of other redox sensor variants is limited. Conclusion: Improved structural stability of sfroGFP2 gives rise to increased FI in cellulo. Fusion to hGrx1 (human glutaredoxin-1) provides the hitherto most suitable biosensor for measuring oxidizing effects in Plasmodium. This sensor is of major interest for studying glutathione redox changes in small cells, as well as subcellular compartments in general. Antioxid. Redox Signal. 37, 1-18.

Viral protein instability enhances host-range evolvability

Viruses are highly evolvable, but what traits endow this property? The high mutation rates of viruses certainly play a role, but factors that act above the genetic code, like protein thermostability, are also expected to contribute. We studied how the thermostability of a model virus, bacteriophage λ, affects its ability to evolve to use a new receptor, a key evolutionary transition that can cause host-range evolution. Using directed evolution and synthetic biology techniques we generated a library of host-recognition protein variants with altered stabilities and then tested their capacity to evolve to use a new receptor. Variants fell within three stability classes: stable, unstable, and catastrophically unstable. The most evolvable were the two unstable variants, whereas seven of eight stable variants were significantly less evolvable, and the two catastrophically unstable variants could not grow. The slowly evolving stable variants were delayed because they required an additional destabilizing mutation. These results are particularly noteworthy because they contradict a widely supported contention that thermostabilizing mutations enhance evolvability of proteins by increasing mutational robustness. Our work suggests that the relationship between thermostability and evolvability is more complex than previously thought, provides evidence for a new molecular model of host-range expansion evolution, and identifies instability as a potential predictor of viral host-range evolution.

Understanding how viruses evolve to infect new hosts is critical for predicting host shifts as well as tuning host-range in phage therapy applications. Yet a mechanistic understanding of the molecular steps required to shift hosts has not been achieved. For this study we examined the evolutionary potential of different strains of a model virus, bacteriophage λ, to gain the ability to use a new receptor, a key step in host shifts. We discovered that λ variants with destabilized host-recognition proteins were more likely to evolve the necessary mutations to use the new receptor than stabilized variants. However, destabilization was only beneficial to a certain point and variants with overly unstable proteins lost all function. These results led us to propose a new molecular model for receptor use evolution in λ; 1) destabilizing mutations evolve that provide protein structural flexibility that allows new protein conformations to form that are able to interact with the new receptors, and 2) mutations evolve that alter the binding surface chemical properties to assist interactions with the new receptor. Our work with a model virus-host system, points to the potential use of viral stability as a phenotypic indicator of the capacity for virus host-range evolution.

Respiratory and C4-photosynthetic NAD-malic enzyme coexist in bundle sheath cell mitochondria and evolved via association of differentially adapted subunits

In plant mitochondria, nicotinamide adenine dinucleotide-malic enzyme (NAD-ME) has a housekeeping function in malate respiration. In different plant lineages, NAD-ME was independently co-opted in C4 photosynthesis. In the C4 Cleome species, Gynandropsis gynandra and Cleome angustifolia, all NAD-ME genes (NAD-MEα, NAD-MEβ1, and NAD-MEβ2) were affected by C4 evolution and are expressed at higher levels than their orthologs in the C3 species Tarenaya hassleriana. In T. hassleriana, the NAD-ME housekeeping function is performed by two heteromers, NAD-MEα/β1 and NAD-MEα/β2, with similar biochemical properties. In both C4 species, this role is restricted to NAD-MEα/β2. In the C4 species, NAD-MEα/β1 is exclusively present in the leaves, where it accounts for most of the enzymatic activity. Gynandropsis gynandra NAD-MEα/β1 (GgNAD-MEα/β1) exhibits high catalytic efficiency and is differentially activated by the C4 intermediate aspartate, confirming its role as the C4-decarboxylase. During C4 evolution, NAD-MEβ1 lost its catalytic activity; its contribution to the enzymatic activity results from a stabilizing effect on the associated α-subunit and the acquisition of regulatory properties. We conclude that in bundle sheath cell mitochondria of C4 species, the functions of NAD-ME as C4 photosynthetic decarboxylase and as a housekeeping enzyme coexist and are performed by isoforms that combine the same α-subunit with differentially adapted β-subunits.

Critical assessment of structure-based approaches to improve protein resistance in aqueous ionic liquids by enzyme-wide saturation mutagenesis

Ionic liquids (IL) and aqueous ionic liquids (aIL) are attractive (co-)solvents for green industrial processes involving biocatalysts, but often reduce enzyme activity. Experimental and computational methods are applied to predict favorable substitution sites and, most often, subsequent site-directed surface charge modifications are introduced to enhance enzyme resistance towards aIL. However, almost no studies evaluate the prediction precision with random mutagenesis or the application of simple data-driven filtering processes. Here, we systematically and rigorously evaluated the performance of 22 previously described structure-based approaches to increase enzyme resistance to aIL based on an experimental complete site-saturation mutagenesis library of Bacillus subtilis Lipase A (BsLipA) screened against four aIL. We show that, surprisingly, most of the approaches yield low gain-in-precision (GiP) values, particularly for predicting relevant positions: 14 approaches perform worse than random mutagenesis. Encouragingly, exploiting experimental information on the thermostability of BsLipA or structural weak spots of BsLipA predicted by rigidity theory yields GiP = 3.03 and 2.39 for relevant variants and GiP = 1.61 and 1.41 for relevant positions. Combining five simple-to-compute physicochemical and evolutionary properties substantially increases the precision of predicting relevant variants and positions, yielding GiP = 3.35 and 1.29. Finally, combining these properties with predictions of structural weak spots identified by rigidity theory additionally improves GiP for relevant variants up to 4-fold to ∼10 and sustains or increases GiP for relevant positions, resulting in a prediction precision of ∼90% compared to ∼9% in random mutagenesis. This combination should be applicable to other enzyme systems for guiding protein engineering approaches towards improved aIL resistance.

Characterization of a novel mesophilic CTP-dependent riboflavin kinase and rational engineering to create its thermostable homologs

Flavins play a central role in metabolism as molecules that catalyze a wide range of redox reactions in living organisms. Several variations in flavin biosynthesis exist among the domains of life, and their analysis has revealed many new structural and mechanistic insights till date. The cytidine triphosphate (CTP)-dependent riboflavin kinase in archaea is one such example - unlike most kinases that use adenosine triphosphate, archaeal riboflavin kinases utilize CTP to phosphorylate riboflavin and produce flavin mononucleotide. In this study, we present the characterization of a new mesophilic archaeal CTP-utilizing riboflavin kinase homolog from Methanococcus maripaludis (MmpRibK), which is linked closely in sequence to the previously characterized thermophilic Methanocaldococcus jannaschii homolog. We reconstitute the activity of MmpRibK, determine its kinetic parameters and molecular factors that contribute to its unique properties, and finally establish the residues that improve its thermostability using computation and a series of experiments.  Our work advances the molecular understanding of flavin biosynthesis in archaea by the characterization of the first mesophilic CTP-dependent riboflavin kinase. Finally, it validates the role of salt bridges and rigidifying amino acid residues in imparting thermostability to this unique structural fold that characterizes archaeal riboflavin kinase enzymes, with implications in enzyme engineering and biotechnological applications.

Aqueous ionic liquids redistribute local enzyme stability via long-range perturbation pathways

Ionic liquids (IL) and aqueous ionic liquids (aIL) are attractive (co-)solvents for biocatalysis due to their unique properties. On the other hand, the incubation of enzymes in IL or aIL often reduces enzyme activity. Recent studies proposed various aIL-induced effects to explain the reduction, classified as direct effects, e.g., local dehydration or competitive inhibition, and indirect effects, e.g., structural perturbations or disturbed catalytic site integrity. However, the molecular origin of indirect effects has largely remained elusive. Here we show by multi-μs long molecular dynamics simulations, free energy computations, and rigidity analyses that aIL favorably interact with specific residues of Bacillus subtilis Lipase A (BsLipA) and modify the local structural stability of this model enzyme by inducing long-range perturbations of noncovalent interactions. The perturbations percolate over neighboring residues and eventually affect the catalytic site and the buried protein core. Validation against a complete experimental site saturation mutagenesis library of BsLipA (3620 variants) reveals that the residues of the perturbation pathways are distinguished sequence positions where substitutions highly likely yield significantly improved residual activity. Our results demonstrate that identifying these perturbation pathways and specific IL ion-residue interactions there effectively predicts focused variant libraries with improved aIL tolerance.

Less Unfavorable Salt Bridges on the Enzyme Surface Result in More Organic Cosolvent Resistance

Biocatalysis for the synthesis of fine chemicals is highly attractive but usually requires organic (co-)solvents (OSs). However, native enzymes often have low activity and resistance in OSs and at elevated temperatures. Herein, we report a smart salt bridge design strategy for simultaneously improving OS resistance and thermostability of the model enzyme, Bacillus subtilits Lipase A (BSLA). We combined comprehensive experimental studies of 3450 BSLA variants and molecular dynamics simulations of 36 systems. Iterative recombination of four beneficial substitutions yielded superior resistant variants with up to 7.6-fold (D64K/D144K) improved resistance toward three OSs while exhibiting significant thermostability (thermal resistance up to 137-fold, and half-life up to 3.3-fold). Molecular dynamics simulations revealed that locally refined flexibility and strengthened hydration jointly govern the highly increased resistance in OSs and at 50-100 °C. The salt bridge redesign provides protein engineers with a powerful and likely general approach to design OSs- and/or thermal-resistant lipases and other α/β-hydrolases.

Promiscuous Esterases Counterintuitively Are Less Flexible than Specific Ones

Understanding mechanisms of promiscuity is increasingly important from a fundamental and application point of view. As to enzyme structural dynamics, more promiscuous enzymes generally have been recognized to also be more flexible. However, examples for the opposite received much less attention. Here, we exploit comprehensive experimental information on the substrate promiscuity of 147 esterases tested against 96 esters together with computationally efficient rigidity analyses to understand the molecular origin of the observed promiscuity range. Unexpectedly, our data reveal that promiscuous esterases are significantly less flexible than specific ones, are significantly more thermostable, and have a significantly increased specific activity. These results may be reconciled with a model according to which structural flexibility in the case of specific esterases serves for conformational proofreading. Our results signify that an esterase sequence space can be screened by rigidity analyses for promiscuous esterases as starting points for further exploration in biotechnology and synthetic chemistry.

An atomistic view of solvent-free protein liquids: the case of Lipase A

Solvent-free enzymes hold the promise of being able to deliver higher activity at elevated temperatures by virtue of them being not limited by the boiling point of the solvent. They have been realized in the liquid phase through a polymer surfactant coating on the protein surface. However, a clear understanding of intermolecular interactions, structure, dynamics, and the behaviour of the minuscule amount of water present in the solvent-free protein liquid is essential to enhance the activity of these biofluids. Using atomistic molecular dynamics simulations, we demonstrate that the scaled spatial correlations between proteins in the hybrid liquid phase of Lipase A enzymes are comparable to the inter-particle correlations in a noble gas fluid. The hydrophilic region of the surfactants forms a coronal layer around each enzyme which percolates throughout the liquid, while the hydrophobic parts are present as disjointed clusters. Inter-surfactant interactions, determined to be attractive and in the range of -200 to -300 kcal mol^-1, stabilize the liquid state. While the protein retains its native state conformational dynamics in the solvent-free form, the fluxionality of its side chains is much reduced; at 333 K, the latter is found to be equivalent to that of the enzyme in an aqueous solution at 249 K. Despite the sluggishness of the solvent-free enzyme, some water molecules exhibit high mobility and transit between enzymes primarily via the interspersed hydrophilic regions. These microscopic insights offer ideas to improve substrate diffusion in the liquid to enable the enhancement of catalytic activity.

Machine Learning Enables Selection of Epistatic Enzyme Mutants for Stability Against Unfolding and Detrimental Aggregation

Machine learning (ML) has pervaded most areas of protein engineering, including stability and stereoselectivity. Using limonene epoxide hydrolase as the model enzyme and innov'SAR as the ML platform, comprising a digital signal process, we achieved high protein robustness that can resist unfolding with concomitant detrimental aggregation. Fourier transform (FT) allows us to take into account the order of the protein sequence and the nonlinear interactions between positions, and thus to grasp epistatic phenomena. The innov'SAR approach is interpolative, extrapolative and makes outside-the-box, predictions not found in other state-of-the-art ML or deep learning approaches. Equally significant is the finding that our approach to ML in the present context, flanked by advanced molecular dynamics simulations, uncovers the connection between epistatic mutational interactions and protein robustness.

Computational design of noncanonical amino acid-based thioether staples at N/C-terminal domains of multi-modular pullulanase for thermostabilization in enzyme catalysis

Enzyme thermostabilization is considered a critical and often obligatory step in biosynthesis, because thermostability is a significant property of enzymes that can be used to evaluate their feasibility for industrial applications. However, conventional strategies for thermostabilizing enzymes generally introduce non-covalent interactions and/or natural covalent bonds caused by natural amino acid substitutions, and the trade-off between the activity and stability of enzymes remains a challenge. Here, we developed a computationally guided strategy for constructing thioether staples by incorporating noncanonical amino acid (ncAA) into the more flexible N/C-terminal domains of the multi-modular pullulanase from Bacillus thermoleovorans (BtPul) to enhance its thermostability. First, potential thioether staples located in the N/C-terminal domains of BtPul were predicted using RosettaMatch. Next, eight variants involving stable thioether staples were precisely predicted using FoldX and Rosetta ddg_monomer. Six positive variants were obtained, of which T73(O2beY)-171C had a 157% longer half-life at 70 °C and an increase of 7.0 °C in T m, when compared with the wild-type (WT). T73(O2beY)-171C/T126F/A72R exhibited an even more improved thermostability, with a 211% increase in half-life at 70 °C and a 44% enhancement in enzyme activity compared with the WT, which was attributed to further optimization of the local interaction network. This work introduces and validates an efficient strategy for enhancing the thermostability and activity of multi-modular enzymes.

Can constraint network analysis guide the identification phase of KnowVolution? A case study on improved thermostability of an endo-β-glucanase

Cellulases are industrially important enzymes, e.g., in the production of bioethanol, in pulp and paper industry, feedstock, and textile. Thermostability is often a prerequisite for high process stability and improving thermostability without affecting specific activities at lower temperatures is challenging and often time-consuming. Protein engineering strategies that combine experimental and computational are emerging in order to reduce experimental screening efforts and speed up enzyme engineering campaigns. Constraint Network Analysis (CNA) is a promising computational method that identifies beneficial positions in enzymes to improve thermostability. In this study, we compare CNA and directed evolution in the identification of beneficial positions in order to evaluate the potential of CNA in protein engineering campaigns (e.g., in the identification phase of KnowVolution). We engineered the industrially relevant endoglucanase EGLII from Penicillium verruculosum towards increased thermostability. From the CNA approach, six variants were obtained with an up to 2-fold improvement in thermostability. The overall experimental burden was reduced to 40% utilizing the CNA method in comparison to directed evolution. On a variant level, the success rate was similar for both strategies, with 0.27% and 0.18% improved variants in the epPCR and CNA-guided library, respectively. In essence, CNA is an effective method for identification of positions that improve thermostability.

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.

Improving the catalytic efficiency of thermostable Geobacillus stearothermophilus xylanase XT6 by single-amino acid substitution

Directed evolution using error-prone polymerase chain reaction was employed in the current study to enhance the catalytic efficiency of a thermostable Geobacillus stearothermophilus xylanase XT6 parent. High-throughput screening identified two variants with enhanced activity. Sequencing analysis revealed the presence of a single-amino acid substitution (P209L or V161L) in each variant. The maximum activity of mutant V161L and P209L was at 85°C and 70°C, respectively. Both mutants exhibited maximum activity at pH 7. The thermal and alkaline tolerance of mutant V161L only were markedly improved. The two mutants were more resistant to ethanol inhibition than the parent. Substrate specificity of the two mutants was shifted from beechwood xylan to birchwood xylan. The potential of the two mutants to hydrolyze rice straw and sugarcane bagasse increased. Both turnover number (kcat) and catalytic efficiency (kcat/kM) increased 12.2- and 5.7-folds for variant P209L and 13- and 6.5-folds for variant V161L, respectively, towards birchwood xylan. Based on the previously published crystal structure of extracellular G. stearothermophilus xylanase XT6, V161L and P209L mutation locate on βα-loops. Conformational changes of the respective loops could potentiate the loop swinging, product release and consequently result in enhancement of the catalytic performance.

iStable 2.0: Predicting protein thermal stability changes by integrating various characteristic modules

Protein mutations can lead to structural changes that affect protein function and result in disease occurrence. In protein engineering, drug design or and optimization industries, mutations are often used to improve protein stability or to change protein properties while maintaining stability. To provide possible candidates for novel protein design, several computational tools for predicting protein stability changes have been developed. Although many prediction tools are available, each tool employs different algorithms and features. This can produce conflicting prediction results that make it difficult for users to decide upon the correct protein design. Therefore, this study proposes an integrated prediction tool, iStable 2.0, which integrates 11 sequence-based and structure-based prediction tools by machine learning and adds protein sequence information as features. Three coding modules are designed for the system, an Online Server Module, a Stand-alone Module and a Sequence Coding Module, to improve the prediction performance of the previous version of the system. The final integrated structure-based classification model has a higher Matthews correlation coefficient than that of the single prediction tool (0.708 vs 0.547, respectively), and the Pearson correlation coefficient of the regression model likewise improves from 0.669 to 0.714. The sequence-based model not only successfully integrates off-the-shelf predictors but also improves the Matthews correlation coefficient of the best single prediction tool by at least 0.161, which is better than the individual structure-based prediction tools. In addition, both the Sequence Coding Module and the Stand-alone Module maintain performance with only a 5% decrease of the Matthews correlation coefficient when the integrated online tools are unavailable. iStable 2.0 is available at http://ncblab.nchu.edu.tw/iStable2.

Systematically Scrutinizing the Impact of Substitution Sites on Thermostability and Detergent Tolerance for Bacillus subtilis Lipase A

Improving an enzyme's (thermo-)stability or tolerance against solvents and detergents is highly relevant in protein engineering and biotechnology. Recent developments have tended toward data-driven approaches, where available knowledge about the protein is used to identify substitution sites with high potential to yield protein variants with improved stability, and subsequently, substitutions are engineered by site-directed or site-saturation (SSM) mutagenesis. However, the development and validation of algorithms for data-driven approaches have been hampered by the lack of availability of large-scale data measured in a uniform way and being unbiased with respect to substitution types and locations. Here, we extend our knowledge on guidelines for protein engineering following a data-driven approach by scrutinizing the impact of substitution sites on thermostability or/and detergent tolerance for Bacillus subtilis lipase A (BsLipA) at very large scale. We systematically analyze a complete experimental SSM library of BsLipA containing all 3439 possible single variants, which was evaluated as to thermostability and tolerances against four detergents under respectively uniform conditions. Our results provide systematic and unbiased reference data at unprecedented scale for a biotechnologically important protein, identify consistently defined hot spot types for evaluating the performance of data-driven protein-engineering approaches, and show that the rigidity theory and ensemble-based approach Constraint Network Analysis yields hot spot predictions with an up to ninefold gain in precision over random classification.

Basal Histamine H4 Receptor Activation: Agonist Mimicry by the Diphenylalanine Motif

Histamine H4 receptor (H4 R) orthologues are G-protein-coupled receptors (GPCRs) that exhibit species-dependent basal activity. In contrast to the basally inactive mouse H4 R (mH4 R), human H4 R (hH4 R) shows a high degree of basal activity. We have performed long-timescale molecular dynamics simulations and rigidity analyses on wild-type hH4 R, the experimentally characterized hH4 R variants S179M, F169V, F169V+S179M, F168A, and on mH4 R to investigate the molecular nature of the differential basal activity. H4 R variant-dependent differences between essential motifs of GPCR activation and structural stabilities correlate with experimentally determined basal activities and provide a molecular explanation for the differences in basal activation. Strikingly, during the MD simulations, F16945.55 dips into the orthosteric binding pocket only in the case of hH4 R, thus adopting the role of an agonist and contributing to the stabilization of the active state. The results shed new light on the molecular mechanism of basal H4 R activation that are of importance for other GPCRs.

The Membrane-Integrated Steric Chaperone Lif Facilitates Active Site Opening of Pseudomonas aeruginosa Lipase A

Lipases are essential and widely used biocatalysts. Hence, the production of lipases requires a detailed understanding of the molecular mechanism of its folding and secretion. Lipase A from Pseudomonas aeruginosa, PaLipA, constitutes a prominent example that has additional relevance because of its role as a virulence factor in many diseases. PaLipA requires the assistance of a membrane-integrated steric chaperone, the lipase-specific foldase Lif, to achieve its enzymatically active state. However, the molecular mechanism of how Lif activates its cognate lipase has remained elusive. Here, we show by molecular dynamics simulations at the atomistic level and potential of mean force computations that Lif catalyzes the activation process of PaLipA by structurally stabilizing an intermediate PaLipA conformation, particularly a β-sheet in the region of residues 17-30, such that the opening of PaLipA's lid domain is facilitated. This opening allows substrate access to PaLipA's catalytic site. A surprising and so far not fully understood aspect of our study is that the open state of PaLipA is unstable compared to the closed one according to our computational and in vitro biochemical results. We thus speculate that further interactions of PaLipA with the Xcp secretion machinery and/or components of the extracellular matrix contribute to the remaining activity of secreted PaLipA. © 2019 Wiley Periodicals, Inc.

Biallelic mutation of human SLC6A6 encoding the taurine transporter TAUT is linked to early retinal degeneration

We previously reported that inactivation of the transmembrane taurine transporter (TauT or solute carrier 6a6) causes early retinal degeneration in mice. Compatible with taurine's indispensability for cell volume homeostasis, protein stabilization, cytoprotection, antioxidation, and immuno- and neuromodulation, mice develop multisystemic dysfunctions (hearing loss; liver fibrosis; and behavioral, heart, and skeletal muscle abnormalities) later on. Here, by genetic, cell biologic, in vivo ¹H-magnetic resonance spectroscopy and molecular dynamics simulation studies, we conducted in-depth characterization of a novel disorder: human TAUT deficiency. Loss of TAUT function due to a homozygous missense mutation caused panretinal degeneration in 2 brothers. TAUT_p.A78E still localized in the plasma membrane but is predicted to impact structural stabilization. ³H-taurine uptake by peripheral blood mononuclear cells was reduced by 95%, and taurine levels were severely reduced in plasma, skeletal muscle, and brain. Extraocular dysfunctions were not yet detected, but significantly increased urinary excretion of 8-oxo-7,8-dihydroguanosine indicated generally enhanced (yet clinically unapparent) oxidative stress and RNA oxidation, warranting continuous broad surveillance.-Preising, M. N., Görg, B., Friedburg, C., Qvartskhava, N., Budde, B. S., Bonus, M., Toliat, M. R., Pfleger, C., Altmüller, J., Herebian, D., Beyer, M., Zöllner, H. J., Wittsack, H.-J., Schaper, J., Klee, D., Zechner, U., Nürnberg, P., Schipper, J., Schnitzler, A., Gohlke, H., Lorenz, B., Häussinger, D., Bolz, H. J. Biallelic mutation of human SLC6A6 encoding the taurine transporter TAUT is linked to early retinal degeneration.

Thermostability of Lipase A and Dynamic Communication Based on Residue Interaction Network

OBJECTIVE

Dynamic communication caused by mutation affects protein stability. The main objective of this study is to explore how mutations affect communication and to provide further insight into the relationship between heat resistance and signal propagation of Bacillus subtilis lipase (Lip A).

METHODS

The relationship between dynamic communication and Lip A thermostability is studied by long-time MD simulation and residue interaction network. The Dijkstra algorithm is used to get the shortest path of each residue pair. Subsequently, time-series frequent paths and spatio-temporal frequent paths are mined through an Apriori-like algorithm.

RESULTS

Time-series frequent paths show that the communication between residue pairs, both in wild-type lipase (WTL) and mutant 6B, becomes chaotic with an increase in temperature; however, more residues in 6B can maintain stable communication at high temperature, which may be associated with the structural rigidity. Furthermore, spatio-temporal frequent paths reflect the interactions among secondary structures. For WTL at 300K, β7, αC, αB, the longest loop, αA and αF contact frequently. The 310-helix between β3 and αA is penetrated by spatio-temporal frequent paths. At 400K, only αC can be frequently transmitted. For 6B, when at 300K, αA and αF are in more tight contact by spatio-temporal frequent paths though I157M and N166Y. Moreover, the rigidity of the active site His156 and the C-terminal of Lip A are increased, as reflected by the spatio-temporal frequent paths. At 400K, αA and αF, 310-helix between β3 and αA, the longest loop, and the loop where the active site Asp133 is located can still maintain stable communication.

CONCLUSION

From the perspective of residue dynamic communication, it is obviously found that mutations cause changes in interactions between secondary structures and enhance the rigidity of the structure, contributing to the thermal stability and functional activity of 6B.

RankProt: A multi criteria-ranking platform to attain protein thermostabilizing mutations and its in vitro applications - Attribute based prediction method on the principles of Analytical Hierarchical Process

Attaining recombinant thermostable proteins is still a challenge for protein engineering. The complexity is the length of time and enormous efforts required to achieve the desired results. Present work proposes a novel and economic strategy of attaining protein thermostability by predicting site-specific mutations at the shortest possible time. The success of the approach can be attributed to Analytical Hierarchical Process and the outcome was a rationalized thermostable mutation(s) prediction tool- RankProt. Briefly the method involved ranking of 17 biophysical protein features as class predictors, derived from 127 pairs of thermostable and mesostable proteins. Among the 17 predictors, ionic interactions and main-chain to main-chain hydrogen bonds were the highest ranked features with eigen value of 0.091. The success of the tool was judged by multi-fold in silico validation tests and it achieved the prediction accuracy of 91% with AUC 0.927. Further, in vitro validation was carried out by predicting thermostabilizing mutations for mesostable Bacillus subtilis lipase and performing the predicted mutations by multi-site directed mutagenesis. The rationalized method was successful to render the lipase thermostable with optimum temperature stability and T_m increase by 20°C and 7°C respectively. Conclusively it can be said that it was the minimum number of mutations in comparison to the number of mutations incorporated to render Bacillus subtilis lipase thermostable, by directed evolution techniques. The present work shows that protein stabilizing mutations can be rationally designed by balancing the biophysical pleiotropy of proteins, in accordance to the selection pressure.

Ensemble- and Rigidity Theory-Based Perturbation Approach To Analyze Dynamic Allostery

Allostery describes the functional coupling between sites in biomolecules. Recently, the role of changes in protein dynamics for allosteric communication has been highlighted. A quantitative and predictive description of allostery is fundamental for understanding biological processes. Here, we integrate an ensemble-based perturbation approach with the analysis of biomolecular rigidity and flexibility to construct a model of dynamic allostery. Our model, by definition, excludes the possibility of conformational changes, evaluates static, not dynamic, properties of molecular systems, and describes allosteric effects due to ligand binding in terms of a novel free-energy measure. We validated our model on three distinct biomolecular systems: eglin c, protein tyrosine phosphatase 1B, and the lymphocyte function-associated antigen 1 domain. In all cases, it successfully identified key residues for signal transmission in very good agreement with the experiment. It correctly and quantitatively discriminated between positively or negatively cooperative effects for one of the systems. Our model should be a promising tool for the rational discovery of novel allosteric drugs.

The mechanism of interactions between tea polyphenols and porcine pancreatic alpha-amylase: Analysis by inhibition kinetics, fluorescence quenching, differential scanning calorimetry and isothermal titration calorimetry

SCOPE

This study aims to use a combination of biochemical and biophysical methods to derive greater mechanistic understanding of the interactions between tea polyphenols and porcine pancreatic α-amylase (PPA).

METHODS AND RESULTS

The interaction mechanism was studied through fluorescence quenching (FQ), differential scanning calorimetry (DSC) and isothermal titration calorimetry (ITC) and compared with inhibition kinetics. The results showed that a higher quenching effect of polyphenols corresponded to a stronger inhibitory activity against PPA. The red-shift of maximum emission wavelength of PPA bound with some polyphenols indicated a potential structural unfolding of PPA. This was also suggested by the decreased thermostability of PPA with these polyphenols in DSC thermograms. Through thermodynamic binding analysis of ITC and inhibition kinetics, the equilibrium of competitive inhibition was shown to result from the binding of particularly galloylated polyphenols with specific sites on PPA. There were positive linear correlations between the reciprocal of competitive inhibition constant (1/Kic ), quenching constant (KFQ ) and binding constant (Kitc ).

CONCLUSION

The combination of inhibition kinetics, FQ, DSC and ITC can reasonably characterize the interactions between tea polyphenols and PPA. The galloyl moiety is an important group in catechins and theaflavins in terms of binding with and inhibiting the activity of PPA.

Contribution of single amino acid and codon substitutions to the production and secretion of a lipase by Bacillus subtilis

Background

Bacillus subtilis produces and secretes proteins in amounts of up to 20 g/l under optimal conditions. However, protein production can be challenging if transcription and cotranslational secretion are negatively affected, or the target protein is degraded by extracellular proteases. This study aims at elucidating the influence of a target protein on its own production by a systematic mutational analysis of the homologous B. subtilis model protein lipase A (LipA). We have covered the full natural diversity of single amino acid substitutions at 155 positions of LipA by site saturation mutagenesis excluding only highly conserved residues and qualitatively and quantitatively screened about 30,000 clones for extracellular LipA production. Identified variants with beneficial effects on production were sequenced and analyzed regarding B. subtilis growth behavior, extracellular lipase activity and amount as well as changes in lipase transcript levels.

Results

In total, 26 LipA variants were identified showing an up to twofold increase in either amount or activity of extracellular lipase. These variants harbor single amino acid or codon substitutions that did not substantially affect B. subtilis growth. Subsequent exemplary combination of beneficial single amino acid substitutions revealed an additive effect solely at the level of extracellular lipase amount; however, lipase amount and activity could not be increased simultaneously.

Conclusions

Single amino acid and codon substitutions can affect LipA secretion and production by B. subtilis. Several codon-related effects were observed that either enhance lipA transcription or promote a more efficient folding of LipA. Single amino acid substitutions could improve LipA production by increasing its secretion or stability in the culture supernatant. Our findings indicate that optimization of the expression system is not sufficient for efficient protein production in B. subtilis. The sequence of the target protein should also be considered as an optimization target for successful protein production. Our results further suggest that variants with improved properties might be identified much faster and easier if mutagenesis is prioritized towards elements that contribute to enzymatic activity or structural integrity.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-017-0772-z) contains supplementary material, which is available to authorized users.

On the potential alternate binding change mechanism in a dimeric structure of Pyruvate Phosphate Dikinase

The pyruvate phosphate dikinase (PPDK) reaction mechanism is characterized by a distinct spatial separation of reaction centers and large conformational changes involving an opening-closing motion of the nucleotide-binding domain (NBD) and a swiveling motion of the central domain (CD). However, why PPDK is active only in a dimeric form and to what extent an alternate binding change mechanism could underlie this fact has remained elusive. We performed unbiased molecular dynamics simulations, configurational free energy computations, and rigidity analysis to address this question. Our results support the hypothesis that PPDK dimerization influences the opening-closing motion of the NBDs, and that this influence is mediated via the CDs of both chains. Such an influence would be a prerequisite for an alternate binding change mechanism to occur. To the best of our knowledge, this is the first time that a possible explanation has been suggested as to why only dimeric PPDK is active.

Complex network analysis of thermostable mutants of Bacillus subtilis Lipase A

Three-dimensional structures of proteins that regulate their functions can be modelled using complex network based approaches for understanding the structure-function relationship. The six mutants of the protein Lipase A from Bacillus subtilis, harbouring 2 to 12 mutations, retain their function at higher temperatures with negligible variation in their overall three-dimensional crystallographic structures. This enhanced thermostability of the mutants questions the structure-function paradigm. In this paper, a coarse-grained complex network approach is used to elucidate the structural basis of enhanced thermostability in the mutant proteins, by uncovering small but significant local changes distributed throughout the structure, rendering stability to the mutants at higher temperatures. Community structure analysis of the six mutant protein networks uncovers the specific reorganisations among the nodes/residues that occur, in absence of overall structural variations, which induce enhanced rigidity underlying the increased thermostability. This study offers a novel and significant application of complex network analysis that proposes to be useful in the understanding and designing of thermostable proteins.

Multifactorial level of extremostability of proteins: can they be exploited for protein engineering?

Research on extremostable proteins has seen immense growth in the past decade owing to their industrial importance. Basic research of attributes related to extreme-stability requires further exploration. Modern mechanistic approaches to engineer such proteins in vitro will have more impact in industrial biotechnology economy. Developing a priori knowledge about the mechanism behind extreme-stability will nurture better understanding of pathways leading to protein molecular evolution and folding. This review is a vivid compilation about all classes of extremostable proteins and the attributes that lead to myriad of adaptations divulged after an extensive study of 6495 articles belonging to extremostable proteins. Along with detailing on the rationale behind extreme-stability of proteins, emphasis has been put on modern approaches that have been utilized to render proteins extremostable by protein engineering. It was understood that each protein shows different approaches to extreme-stability governed by minute differences in their biophysical properties and the milieu in which they exist. Any general rule has not yet been drawn regarding adaptive mechanisms in extreme environments. This review was further instrumental to understand the drawback of the available 14 stabilizing mutation prediction algorithms. Thus, this review lays the foundation to further explore the biophysical pleiotropy of extreme-stable proteins to deduce a global prediction model for predicting the effect of mutations on protein stability.

Application of Rigidity Theory to the Thermostabilization of Lipase A from Bacillus subtilis

Protein thermostability is a crucial factor for biotechnological enzyme applications. Protein engineering studies aimed at improving thermostability have successfully applied both directed evolution and rational design. However, for rational approaches, the major challenge remains the prediction of mutation sites and optimal amino acid substitutions. Recently, we showed that such mutation sites can be identified as structural weak spots by rigidity theory-based thermal unfolding simulations of proteins. Here, we describe and validate a unique, ensemble-based, yet highly efficient strategy to predict optimal amino acid substitutions at structural weak spots for improving a protein’s thermostability. For this, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state has been found as the cause for thermostability. When applied prospectively to lipase A from Bacillus subtilis, we achieved both a high success rate (25% over all experimentally tested mutations, which raises to 60% if small-to-large residue mutations and mutations in the active site are excluded) in predicting significantly thermostabilized lipase variants and a remarkably large increase in those variants’ thermostability (up to 6.6°C) based on single amino acid mutations. When considering negative controls in addition and evaluating the performance of our approach as a binary classifier, the accuracy is 63% and increases to 83% if small-to-large residue mutations and mutations in the active site are excluded. The gain in precision (predictive value for increased thermostability) over random classification is 1.6-fold (2.4-fold). Furthermore, an increase in thermostability predicted by our approach significantly points to increased experimental thermostability (p < 0.05). These results suggest that our strategy is a valuable complement to existing methods for rational protein design aimed at improving thermostability.

Protein thermostability is a crucial factor for biotechnological enzyme applications. However, performance studies of computational approaches for predicting effects of mutations on protein (thermo)stability have suggested that there is still room for improvement. We describe and validate a novel and unique strategy to predict optimal amino acid substitutions at structural weak spots. At variance with other rational approaches, we exploit the fact that in the majority of cases an increased structural rigidity of the folded state is the underlying cause for thermostability. When applied prospectively on lipase LipA from Bacillus subtilis, a high success rate in predicting thermostabilized lipase variants and a remarkably large increase in their thermostability is achieved. This demonstrates the value of the novel strategy, which extends the existing portfolio of methods for rational protein design.

Trading off stability against activity in extremophilic aldolases

Understanding enzyme stability and activity in extremophilic organisms is of great biotechnological interest, but many questions are still unsolved. Using 2-deoxy-D-ribose-5-phosphate aldolase (DERA) as model enzyme, we have evaluated structural and functional characteristics of different orthologs from psychrophilic, mesophilic and hyperthermophilic organisms. We present the first crystal structures of psychrophilic DERAs, revealing a dimeric organization resembling their mesophilic but not their thermophilic counterparts. Conversion into monomeric proteins showed that the native dimer interface contributes to stability only in the hyperthermophilic enzymes. Nevertheless, introduction of a disulfide bridge in the interface of a psychrophilic DERA did confer increased thermostability, suggesting a strategy for rational design of more durable enzyme variants. Constraint network analysis revealed particularly sparse interactions between the substrate pocket and its surrounding α-helices in psychrophilic DERAs, which indicates that a more flexible active center underlies their high turnover numbers.

Protein thermostability engineering

Using structure and sequence based analysis we can engineer proteins to increase their thermal stability.