Structural and catalytic effects of proline substitution and surface loop deletion in the extended active site of human carbonic anhydrase II

Impact of single-residue mutations on protein thermal stability: The case of threonine 83 of BC2L-CN lectin

The thermal stability of trimeric lectin BC2L-CN was investigated and found to be considerably altered when mutating residue 83, originally a threonine, located at the fucose-binding loop. Mutants were analyzed using differential scanning calorimetry and isothermal microcalorimetry. Although most mutations decreased the affinity of the protein for oligosaccharide H type 1, six mutations increased the melting temperature (Tm) by >5 °C; one mutation, T83P, increased the Tm value by 18.2 °C(T83P, Tm = 96.3 °C). In molecular dynamic simulations, the investigated thermostable mutants, T83P, T83A, and T83S, had decreased fluctuations in the loop containing residue 83. In the T83S mutation, the side-chain hydroxyl group of serine formed a hydrogen bond with a nearby residue, suggesting that the restricted movement of the side-chain resulted in fewer fluctuations and enhanced thermal stability. Residue 83 is located at the interface and near the upstream end of the equivalent loop in a different protomer; therefore, fluctuations by this residue likely propagate throughout the loop. Our study of the dramatic change in thermal stability by a single amino acid mutation provides useful insights into the rational design of protein structures, especially the structures of oligomeric proteins.

Bioprocess development to produce a hyperthermostable S-methyl-5'-thioadenosine phosphorylase in Escherichia coli

Nucleoside phosphorylases are important biocatalysts for the chemo-enzymatic synthesis of nucleosides and their analogs which are, among others, used for the treatment of viral infections or cancer. S-methyl-5'-thioadenosine phosphorylases (MTAP) are a group of nucleoside phosphorylases and the thermostable MTAP of Aeropyrum pernix (ApMTAP) was described to accept a wide range of modified nucleosides as substrates. Therefore, it is an interesting biocatalyst for the synthesis of nucleoside analogs for industrial and therapeutic applications. To date, thermostable nucleoside phosphorylases were produced in shake flask cultivations using complex media. The drawback of this approach is low volumetric protein yields which hamper the wide-spread application of the thermostable nucleoside phosphorylases in large scale. High cell density (HCD) cultivations allow the production of recombinant proteins with high volumetric yields, as final optical densities >100 can be achieved. Therefore, in this study, we developed a suitable protocol for HCD cultivations of ApMTAP. Initially, optimum expression conditions were determined in 24-well plates using a fed-batch medium. Subsequently, HCD cultivations were performed using E. coli BL21-Gold cells, by employing a glucose-limited fed-batch strategy. Comparing different growth rates in stirred-tank bioreactors, cultivations revealed that growth at maximum growth rates until induction resulted in the highest yields of ApMTAP. On a 500-mL scale, final cell dry weights of 87.1-90.1 g L-1 were observed together with an overproduction of ApMTAP in a 1.9%-3.8% ratio of total protein. Compared to initially applied shake flask cultivations with terrific broth (TB) medium the volumetric yield increased by a factor of 136. After the purification of ApMTAP via heat treatment and affinity chromatography, a purity of more than 90% was determined. Activity testing revealed specific activities in the range of 0.21 ± 0.11 (low growth rate) to 3.99 ± 1.02 U mg-1 (growth at maximum growth rate). Hence, growth at maximum growth rate led to both an increased expression of the target protein and an increased specific enzyme activity. This study paves the way towards the application of thermostable nucleoside phosphorylases in industrial applications due to an improved heterologous expression in Escherichia coli.

Improvement of the stability and catalytic efficiency of heparan sulfate N-sulfotransferase for preparing N-sulfated heparosan

The chemo-enzymatic and enzymatic synthesis of heparan sulfate and heparin are considered as an attractive alternative to the extraction of heparin from animal tissues. Sulfation of the hydroxyl group at position 2 of the deacetylated glucosamine is a prerequisite for subsequent enzymatic modifications. In this study, multiple strategies, including truncation mutagenesis based on B-factor values, site-directed mutagenesis guided by multiple sequence alignment, and structural analysis were performed to improve the stability and activity of human N-sulfotransferase. Eventually, a combined variant Mut02 (MBP-hNST-NΔ599-602/S637P/S741P/E839P/L842P/K779N/R782V) was successfully constructed, whose half-life at 37°C and catalytic activity were increased by 105-fold and 1.35-fold, respectively. After efficient overexpression using the Escherichia coli expression system, the variant Mut02 was applied to N-sulfation of the chemically deacetylated heparosan. The N-sulfation content reached around 82.87% which was nearly 1.88-fold higher than that of the wild-type. The variant Mut02 with high stability and catalytic efficiency has great potential for heparin biomanufacturing.

Targeting the BspC-vimentin interaction to develop anti-virulence therapies during Group B streptococcal meningitis

Bacterial infections are a major cause of morbidity and mortality worldwide and the rise of antibiotic resistance necessitates development of alternative treatments. Pathogen adhesins that bind to host cells initiate disease pathogenesis and represent potential therapeutic targets. We have shown previously that the BspC adhesin in Group B Streptococcus (GBS), the leading cause of bacterial neonatal meningitis, interacts with host vimentin to promote attachment to brain endothelium and disease development. Here we determined that the BspC variable (V-) domain contains the vimentin binding site and promotes GBS adherence to brain endothelium. Site directed mutagenesis identified a binding pocket necessary for GBS host cell interaction and development of meningitis. Using a virtual structure-based drug screen we identified compounds that targeted the V-domain binding pocket, which blocked GBS adherence and entry into the brain in vivo. These data indicate the utility of targeting the pathogen-host interface to develop anti-virulence therapeutics.

Genome-Wide Identification of Aldehyde Oxidase Genes in Moths and Butterflies Suggests New Insights Into Their Function as Odorant-Degrading Enzymes

Aldehyde oxidases (AOXs) are common detoxifying enzymes in several organisms. In insects, AOXs act in xenobiotic metabolism and as odorant-degrading enzymes (ODEs). These last appear as crucial enzymes in the life cycle of insects, helping to reset their olfactory system, particularly in lepidopterans, which fulfill important ecological roles (e.g., pollination or destructive life cycles). A comprehensive understanding of their olfactory system has provided opportunities to study key chemosensory proteins. However, no significant advance has been made around lepidopteran AOXs research, and even less around butterflies, a recently evolved lineage. In this study we have identified novel AOX gene families in moths and butterflies in order to understand their role as ODEs. Eighteen genomes from both moths and butterflies were used for phylogenetics, molecular evolution and sequence analyses. We identified 164 AOXs, from which 91 are new. Their phylogeny showed two main clades that are potentially related to odorant-degrading function, where both moths and butterflies have AOXs. A first ODE-related clade seems to have a non-ditrysian origin, likely related to plant volatiles. A second ODE-related clade could be more pheromone-biased. Molecular evolution analysis suggests a slight purifying selection process, though a number of sites appeared under positive selection. ODE-related AOXs have changed a phenylalanine residue by proline in the active site. Finally, this study could serve as a reference for further evolutionary and functional studies around Lepidopteran AOXs.

Tolerance to Glutaraldehyde in Escherichia coli Mediated by Overexpression of the Aldehyde Reductase YqhD by YqhC

Glutaraldehyde is a widely used biocide on the market for about 50 years. Despite its broad application, several reports on the emergence of bacterial resistance, and occasional outbreaks caused by poorly disinfection, there is a gap of knowledge on the bacterial adaptation, tolerance, and resistance mechanisms to glutaraldehyde. Here, we analyze the effects of the independent selection of mutations in the transcriptional regulator yqhC for biological replicates of Escherichia coli cells subjected to adaptive laboratory evolution (ALE) in the presence of glutaraldehyde. The evolved strains showed improved survival in the biocide (11-26% increase in fitness) as a result of mutations in the activator yqhC, which led to the overexpression of the yqhD aldehyde reductase gene by 8 to over 30-fold (3.1-5.2 log2FC range). The protective effect was exclusive to yqhD as other aldehyde reductase genes of E. coli, such as yahK, ybbO, yghA, and ahr did not offer protection against the biocide. We describe a novel mechanism of tolerance to glutaraldehyde based on the activation of the aldehyde reductase YqhD by YqhC and bring attention to the potential for the selection of such tolerance mechanism outside the laboratory, given the existence of YqhD homologs in various pathogenic and opportunistic bacterial species.

Engineering stable carbonic anhydrases for CO2 capture: a critical review

Targeted inhibition of misregulated protein-protein interactions (PPIs) has been a promising area of investigation in drug discovery and development for human diseases. However, many constraints remain, including shallow binding surfaces and dynamic conformation changes upon interaction. A particularly challenging aspect is the undesirable off-target effects caused by inherent structural similarity among the protein families. To tackle this problem, phage display has been used to engineer PPIs for high-specificity binders with improved binding affinity and greatly reduced undesirable interactions with closely related proteins. Although general steps of phage display are standardized, library design is highly variable depending on experimental contexts. Here in this review, we examined recent advances in the structure-based combinatorial library design and the advantages and limitations of different approaches. The strategies described here can be explored for other protein-protein interactions and aid in designing new libraries or improving on previous libraries.

Bioengineering oligomerization and monomerization of enzymes: learning from natural evolution to matching the demands for industrial applications

It is generally accepted that oligomeric enzymes evolve from their monomeric ancestors, and the evolution process generates superior structural benefits for functional advantages. Furthermore, adjusting the transition between different oligomeric states is an important mechanism for natural enzymes to regulate their catalytic functions for adapting environmental fluctuations in nature, which inspires researchers to mimic such a strategy to develop artificially oligomerized enzymes through protein engineering for improved performance under specific conditions. On the other hand, transforming oligomeric enzymes into their monomers is needed in fundamental research for deciphering catalytic mechanisms as well as exploring their catalytic capacities for better industrial applications. In this article, strategies for developing artificially oligomerized and monomerized enzymes are reviewed and highlighted by their applications. Furthermore, advances in the computational prediction of oligomeric structures are introduced, which would accelerate the systematic design of oligomeric and monomeric enzymes. Finally, the current challenges and future directions in this field are discussed.

Catalytic and structural effects of flexible loop deletion in organophosphorus hydrolase enzyme: A thermostability improvement mechanism

Thermostability improvement of enzymes used industrially or commercially would develop their capacity and commercial potential due to increased enzymatic competence and cost-effectiveness. Several stabilizing factors have been suggested to be the base of thermal stability, like proline replacements, disulfide bonds, surface loop truncation and ionic pair networks creation. This research evaluated the mechanism of increasing the rigidity of organophosphorus hydrolase enzyme by flexible loop truncation. Bioinformatics analysis revealed that the mutated protein retains its stability after loop truncation (five amino acids deleted). The thermostability of the wild-type (OPH-wt) and mutated (OPH-D5) enzymes were investigated by half-life, Delta Gi, and fluorescence and far-UV CD analysis. Results demonstrated an increase half-life and Delta Gi in OPH-D5 compared to OPH-wt. These results were confirmed by extrinsic fluorescence and circular dichroism (CD) spectrometry experiments, therefore, as rigidity increased in OPHD5 after loop truncation, half-life and Delta Gi also increased. Based on these findings, a strong case is presented for thermostability improvement of OPH enzyme by flexible loop truncation after bioinformatics analysis.

Finding the generalized molecular principles of protein thermal stability

Are there any generalized molecular principles of thermal adaptation? Here, integrating the concepts of structural bioinformatics, sequence analysis, and classical knot theory, we develop a robust computational framework that seeks for mechanisms of thermal adaptation by comparing orthologous mesophilic-thermophilic and mesophilic-hyperthermophilic proteins of remarkable structural and topological similarities, and still leads us to context-independent results. A comprehensive analysis of 4741 high-resolution, non-redundant X-ray crystallographic structures collected from 11 hyperthermophilic, 32 thermophilic and 53 mesophilic prokaryotes unravels at least five "nearly universal" signatures of thermal adaptation, irrespective of the enormous sequence, structure, and functional diversity of the proteins compared. A careful investigation further extracts a set of amino acid changes that can potentially enhance protein thermal stability, and remarkably, these mutations are overrepresented in protein crystallization experiments, in disorder-to-order transitions and in engineered thermostable variants of existing mesophilic proteins. These results could be helpful to find a precise, global picture of thermal adaptation.

Enhancing subtilisin thermostability through a modified normalized B-factor analysis and loop-grafting strategy

Rational design-guided improvement of protein thermostability typically requires identification of residues or regions contributing to instability and introduction of mutations into these residues or regions. One popular method, B-FIT, utilizes B-factors to identify unstable residues or regions and combines them with other strategies, such as directed evolution. Here, we performed structure-based engineering to improve the thermostability of the subtilisin E-S7 (SES7) peptidase. The B-value of each residue was redefined in a normalized B-factor calculation, which was implemented with a refined bioinformatics analysis strategy to identify the critical area (loop 158-162) related to flexibility and to screen for suitable thermostable motif sequences in the Protein Data Bank that can act as transplant loops. In total, we analyzed 445 structures and identified 29 thermostable motifs as candidates. Using these motifs as a starting point, we performed iterative homologous modeling to obtain a desirable chimera loop and introduced five different mutations into this loop to construct thermostable SES7 proteins. Differential scanning fluorimetry revealed increases of 7.3 °C in the melting temperature of an SES7 variant designated M5 compared with the WT. The X-ray crystallographic structure of this variant was resolved at 1.96 Å resolution. The crystal structure disclosed that M5 forms more hydrogen bonds than the WT protein, consistent with design and molecular dynamics simulation results. In summary, the modified B-FIT strategy reported here has yielded a subtilisin variant with improved thermostability and promising industrial applications, supporting the notion that this modified method is a powerful tool for protein engineering.

Thermodynamic, kinetic, and structural parameterization of human carbonic anhydrase interactions toward enhanced inhibitor design

The aim of rational drug design is to develop small molecules using a quantitative approach to optimize affinity. This should enhance the development of chemical compounds that would specifically, selectively, reversibly, and with high affinity interact with a target protein. It is not yet possible to develop such compounds using computational (i.e., in silico) approach and instead the lead molecules are discovered in high-throughput screening searches of large compound libraries. The main reason why in silico methods are not capable to deliver is our poor understanding of the compound structure-thermodynamics and structure-kinetics correlations. There is a need for databases of intrinsic binding parameters (e.g., the change upon binding in standard Gibbs energy (ΔGint), enthalpy (ΔHint), entropy (ΔSint), volume (ΔVintr), heat capacity (ΔCp,int), association rate (ka,int), and dissociation rate (kd,int)) between a series of closely related proteins and a chemically diverse, but pharmacophoric group-guided library of compounds together with the co-crystal structures that could help explain the structure-energetics correlations and rationally design novel compounds. Assembly of these data will facilitate attempts to provide correlations and train data for modeling of compound binding. Here, we report large datasets of the intrinsic thermodynamic and kinetic data including over 400 primary sulfonamide compound binding to a family of 12 catalytically active human carbonic anhydrases (CA). Thermodynamic parameters have been determined by the fluorescent thermal shift assay, isothermal titration calorimetry, and by the stopped-flow assay of the inhibition of enzymatic activity. Kinetic measurements were performed using surface plasmon resonance. Intrinsic thermodynamic and kinetic parameters of binding were determined by dissecting the binding-linked protonation reactions of the protein and sulfonamide. The compound structure-thermodynamics and kinetics correlations reported here helped to discover compounds that exhibited picomolar affinities, hour-long residence times, and million-fold selectivities over non-target CA isoforms. Drug-lead compounds are suggested for anticancer target CA IX and CA XII, antiglaucoma CA IV, antiobesity CA VA and CA VB, and other isoforms. Together with 85 X-ray crystallographic structures of 60 compounds bound to six CA isoforms, the database should be of help to continue developing the principles of rational target-based drug design.

Computation-Aided Rational Deletion of C-Terminal Region Improved the Stability, Activity, and Expression Level of GH2 β-Glucuronidase

In this study, computation-aided design on the basis of structural analysis was employed to rationally identify a highly dynamic C-terminal region that regulates the stability, expression level, and activity of a GH2 fungal glucuronidase from Aspergillus oryzae Li-3 (PGUS). Then, four mutants with a precisely truncated C-terminal region in different lengths were constructed; among them, mutant D591-604 with a 3.8-fold increase in half-life at 65 °C and a 6.8 kJ/mol increase in Gibbs free energy showed obviously improved kinetic and thermodynamic stability in comparison to PGUS. Mutants D590-604 and D591-604 both showed approximately 2.4-fold increases in the catalytic efficiency k_cat/ K_m and 1.8-fold increases in the expression level. Additionally, the expression level of PGUS was doubled through a C-terminal region swap with bacterial GUS from E. coli (EGUS). Finally, the robust PGUS mutants D590-604 and D591-604 were applied in the preparation of glycyrrhetinic acid with 4.0- and 4.4-fold increases in concentration through glycyrrhizin hydrolysis by a fed-batch process.

Enhanced trypsin thermostability in Pichia pastoris through truncating the flexible region

Background

High thermostability is required for trypsin to have wider industrial applications. Target mutagenesis at flexible regions has been proved to be an efficient protein engineering method to enhance the protein thermostability.

Results

The flexible regions in porcine trypsin were predicted using the methods including molecular dynamic simulation, FlexPred, and FoldUnfold. The amino acids 78–90 was predicted to be the highly flexible region simultaneously by the three methods and hence selected to be the mutation target. We constructed five variants (D3, D5, D7, D9, and D11) by truncating the region. And the variant D9 showed higher thermostability, with a 5 °C increase in T_opt, 5.8 °C rise in \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$T_{50}^{10}$$\end{document}T5010, and a 4.5 °C rise in T_m, compared to the wild-type. Moreover, the half-life value of the variant D9 was also found to be dramatically improved by 46 min. Circular dichroism and intrinsic fluorescence indicated that the structures had no significant change between the variant D9 and the wild-type. The surface hydrophobicity of D9 was measured to be lower than that of wild-type, indicating the increased hydrophobic interaction, which could have contributed to the improved thermostability of D9.

Conclusions

These results showed that the thermostability of variant D9 was increased. The variant D9 could be expected to be a promising tool enzyme for its wider industrial applications. The method of truncating the flexible region used in our study has the potential to be used for enhancing the thermostability of other proteins.

Electronic supplementary material

The online version of this article (10.1186/s12934-018-1012-x) contains supplementary material, which is available to authorized users.

Physical and molecular bases of protein thermal stability and cold adaptation

The molecular bases of thermal and cold stability and adaptation, which allow proteins to remain folded and functional in the temperature ranges in which their host organisms live and grow, are still only partially elucidated. Indeed, both experimental and computational studies fail to yield a fully precise and global physical picture, essentially because all effects are context-dependent and thus quite intricate to unravel. We present a snapshot of the current state of knowledge of this highly complex and challenging issue, whose resolution would enable large-scale rational protein design.

Engineering de novo disulfide bond in bacterial α-type carbonic anhydrase for thermostable carbon sequestration

Exploiting carbonic anhydrase (CA), an enzyme that rapidly catalyzes carbon dioxide hydration, is an attractive biomimetic route for carbon sequestration due to its environmental compatibility and potential economic viability. However, the industrial applications of CA are strongly hampered by the unstable nature of enzymes. In this work, we introduced in silico designed, de novo disulfide bond in a bacterial α-type CA to enhance thermostability. Three variants were selected and expressed in Escherichia coli with an additional disulfide bridge. One of the variants showed great enhancement in terms of both kinetic and thermodynamic stabilities. This improvement could be attributed to the loss of conformational entropy of the unfolded state, showing increased rigidity. The variant showed an upward-shifted optimal temperature and appeared to be thermoactivated, which compensated for the lowered activity at 25 °C. Collectively, the variant constructed by the rapid and effective de novo disulfide engineering can be used as an efficient biocatalyst for carbon sequestration under high temperature conditions.

Activity and anion inhibition studies of the α-carbonic anhydrase from Thiomicrospira crunogena XCL-2 Gammaproteobacterium

Thiomicrospira crunogena XCL-2 expresses an α-carbonic anhydrase (TcruCA). Sequence alignments reveal that TcruCA displays a high sequence identity (>30%) relative to other α-CAs. This includes three conserved histidines that coordinate the active site zinc, a histidine proton shuttling residue, and opposing hydrophilic and hydrophobic sides that line the active site. The catalytic efficiency of TcruCA is considered moderate relative to other α-CAs (k(cat)/K(M)=1.1×10(7) M(-1) s(-1)), being a factor of ten less efficient than the most active α-CAs. TcruCA is also inhibited by anions with Cl(-), Br(-), and I(-), all showing Ki values in the millimolar range (53-361 mM). Hydrogen sulfide (HS(-)) revealed the highest affinity for TcruCA with a Ki of 1.1 μM. It is predicted that inhibition of TcruCA by HS(-) (an anion commonly found in the environment where Thiomicrospira crunogena is located) is a way for Thiomicrospira crunogena to regulate its carbon-concentrating mechanism (CCM) and thus the organism's metabolic functions. Results from this study provide preliminary insights into the role of TcruCA in the general metabolism of Thiomicrospira crunogena.

Structural and biophysical characterization of the α-carbonic anhydrase from the gammaproteobacterium Thiomicrospira crunogena XCL-2: insights into engineering thermostable enzymes for CO2 sequestration

Biocatalytic CO2 sequestration to reduce greenhouse-gas emissions from industrial processes is an active area of research. Carbonic anhydrases (CAs) are attractive enzymes for this process. However, the most active CAs display limited thermal and pH stability, making them less than ideal. As a result, there is an ongoing effort to engineer and/or find a thermostable CA to fulfill these needs. Here, the kinetic and thermal characterization is presented of an α-CA recently discovered in the mesophilic hydrothermal vent-isolate extremophile Thiomicrospira crunogena XCL-2 (TcruCA), which has a significantly higher thermostability compared with human CA II (melting temperature of 71.9°C versus 59.5°C, respectively) but with a tenfold decrease in the catalytic efficiency. The X-ray crystallographic structure of the dimeric TcruCA shows that it has a highly conserved yet compact structure compared with other α-CAs. In addition, TcruCA contains an intramolecular disulfide bond that stabilizes the enzyme. These features are thought to contribute significantly to the thermostability and pH stability of the enzyme and may be exploited to engineer α-CAs for applications in industrial CO2 sequestration.

Thermostable Carbonic Anhydrases in Biotechnological Applications

Carbonic anhydrases are ubiquitous metallo-enzymes which catalyze the reversible hydration of carbon dioxide in bicarbonate ions and protons. Recent years have seen an increasing interest in the utilization of these enzymes in CO₂ capture and storage processes. However, since this use is greatly limited by the harsh conditions required in these processes, the employment of thermostable enzymes, both those isolated by thermophilic organisms and those obtained by protein engineering techniques, represents an interesting possibility. In this review we will provide an extensive description of the thermostable carbonic anhydrases so far reported and the main processes in which these enzymes have found an application.