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

Discovery and biochemical characterization of thermostable glycerol oxidases

Alditol oxidases are promising tools for the biocatalytic oxidation of glycerol to more valuable chemicals. By integrating in silico bioprospecting with cell-free protein synthesis and activity screening, an effective pipeline was developed to rapidly identify enzymes that are active on glycerol. Three thermostable alditol oxidases from Actinobacteria Bacterium, Streptomyces thermoviolaceus, and Thermostaphylospora chromogena active on glycerol were discovered. The characterization of these three flavoenzymes demonstrated their glycerol oxidation activities, preference for alkaline conditions, and excellent thermostabilities with melting temperatures higher than 75 °C. Structural elucidation of the alditol oxidase from Actinobacteria Bacterium highlighted a constellation of side chains that engage the substrate through several hydrogen bonds, a histidine residue covalently bound to the FAD prosthetic group, and a tunnel leading to the active site. Upon computational simulations of substrate binding, a double mutant targeting a residue pair at the tunnel entrance was created and found to display an improved thermal stability and catalytic efficiency for glycerol oxidation. The hereby described alditol oxidases form a valuable panel of oxidative biocatalysts that can perform regioselective oxidation of glycerol and other polyols. KEY POINTS: • Rapid pipeline designed to identify putative oxidases • Biochemical and structural characterization of alditol oxidases • Glycerol oxidation to more valuable derivatives.

An engineering strategy to target activated EGFR with CAR T cells

Chimeric antigen receptor (CAR) T cells have shown remarkable response rates in hematological malignancies. In contrast, CAR T cell treatment of solid tumors is associated with several challenges, in particular the expression of most tumor-associated antigens at lower levels in vital organs, resulting in on-target/off-tumor toxicities. Thus, innovative approaches to improve the tumor specificity of CAR T cells are urgently needed. Based on the observation that many human solid tumors activate epidermal growth factor receptor (EGFR) on their surface through secretion of EGFR ligands, we developed an engineering strategy for CAR-binding domains specifically directed against the ligand-activated conformation of EGFR. We show, in several experimental systems, that the generated binding domains indeed enable CAR T cells to distinguish between active and inactive EGFR. We anticipate that this engineering concept will be an important step forward to improve the tumor specificity of CAR T cells directed against EGFR-positive solid cancers.

Molecular Mechanisms Underlying Medium-Chain Free Fatty Acid-Regulated Activity of the Phospholipase PlaF from Pseudomonas aeruginosa

PlaF is a membrane-bound phospholipase A1 from Pseudomonas aeruginosa that is involved in remodeling membrane glycerophospholipids (GPLs) and modulating virulence-associated signaling and metabolic pathways. Previously, we identified the role of medium-chain free fatty acids (FFAs) in inhibiting PlaF activity and promoting homodimerization, yet the underlying molecular mechanism remained elusive. Here, we used unbiased and biased molecular dynamics simulations and free energy computations to assess how PlaF interacts with FFAs localized in the water milieu surrounding the bilayer or within the bilayer and how these interactions regulate PlaF activity. Medium-chain FFAs localized in the upper bilayer leaflet can stabilize inactive dimeric PlaF, likely through interactions with charged surface residues, as has been experimentally validated. Potential of mean force (PMF) computations indicate that membrane-bound FFAs may facilitate the activation of monomeric PlaF by lowering the activation barrier for changing into a tilted, active configuration. We estimated that the coupled equilibria of PlaF monomerization-dimerization and tilting at the physiological concentration of PlaF lead to the majority of PlaF forming inactive dimers when in a cell membrane loaded with decanoic acid (C10). This is in agreement with a suggested in vivo product feedback loop and gas chromatography-mass spectrometry profiling results, indicating that PlaF catalyzes the release of C10 from P. aeruginosa membranes. Additionally, we found that C10 in the water milieu can access the catalytic site of active monomeric PlaF, contributing to the competitive component of C10-mediated PlaF inhibition. Our study provides mechanistic insights into how medium-chain FFAs may regulate the activity of PlaF, a potential bacterial drug target.

In silico design of multipoint mutants for enhanced performance of Thermomyces lanuginosus lipase for efficient biodiesel production

BACKGROUND

Biodiesel, an emerging sustainable and renewable clean energy, has garnered considerable attention as an alternative to fossil fuels. Although lipases are promising catalysts for biodiesel production, their efficiency in industrial-scale application still requires improvement.

RESULTS

In this study, a novel strategy for multi-site mutagenesis in the binding pocket was developed via FuncLib (for mutant enzyme design) and Rosetta Cartesian_ddg (for free energy calculation) to improve the reaction rate and yield of lipase-catalyzed biodiesel production. Thermomyces lanuginosus lipase (TLL) with high activity and thermostability was obtained using the Pichia pastoris expression system. The specific activities of the mutants M11 and M21 (each with 5 and 4 mutations) were 1.50- and 3.10-fold higher, respectively, than those of the wild-type (wt-TLL). Their corresponding melting temperature profiles increased by 10.53 and 6.01 °C, [Formula: see text] (the temperature at which the activity is reduced to 50% after 15 min incubation) increased from 60.88 to 68.46 °C and 66.30 °C, and the optimum temperatures shifted from 45 to 50 °C. After incubation in 60% methanol for 1 h, the mutants M11 and M21 retained more than 60% activity, and 45% higher activity than that of wt-TLL. Molecular dynamics simulations indicated that the increase in thermostability could be explained by reduced atomic fluctuation, and the improved catalytic properties were attributed to a reduced binding free energy and newly formed hydrophobic interaction. Yields of biodiesel production catalyzed by mutants M11 and M21 for 48 h at an elevated temperature (50 °C) were 94.03% and 98.56%, respectively, markedly higher than that of the wt-TLL (88.56%) at its optimal temperature (45 °C) by transesterification of soybean oil.

CONCLUSIONS

An integrating strategy was first adopted to realize the co-evolution of catalytic efficiency and thermostability of lipase. Two promising mutants M11 and M21 with excellent properties exhibited great potential for practical applications for in biodiesel production.

Surface frustration re-patterning underlies the structural landscape and evolvability of fungal orphan candidate effectors

Pathogens secrete effector proteins to subvert host physiology and cause disease. Effectors are engaged in a molecular arms race with the host resulting in conflicting evolutionary constraints to manipulate host cells without triggering immune responses. The molecular mechanisms allowing effectors to be at the same time robust and evolvable remain largely enigmatic. Here, we show that 62 conserved structure-related families encompass the majority of fungal orphan effector candidates in the Pezizomycotina subphylum. These effectors diversified through changes in patterns of thermodynamic frustration at surface residues. The underlying mutations tended to increase the robustness of the overall effector protein structure while switching potential binding interfaces. This mechanism could explain how conserved effector families maintained biological activity over long evolutionary timespans in different host environments and provides a model for the emergence of sequence-unrelated effector families with conserved structures.

Engineering the kinetic stability of a β-trefoil protein by tuning its topological complexity

Kinetic stability, defined as the rate of protein unfolding, is central to determining the functional lifetime of proteins, both in nature and in wide-ranging medical and biotechnological applications. Further, high kinetic stability is generally correlated with high resistance against chemical and thermal denaturation, as well as proteolytic degradation. Despite its significance, specific mechanisms governing kinetic stability remain largely unknown, and few studies address the rational design of kinetic stability. Here, we describe a method for designing protein kinetic stability that uses protein long-range order, absolute contact order, and simulated free energy barriers of unfolding to quantitatively analyze and predict unfolding kinetics. We analyze two β-trefoil proteins: hisactophilin, a quasi-three-fold symmetric natural protein with moderate stability, and ThreeFoil, a designed three-fold symmetric protein with extremely high kinetic stability. The quantitative analysis identifies marked differences in long-range interactions across the protein hydrophobic cores that partially account for the differences in kinetic stability. Swapping the core interactions of ThreeFoil into hisactophilin increases kinetic stability with close agreement between predicted and experimentally measured unfolding rates. These results demonstrate the predictive power of readily applied measures of protein topology for altering kinetic stability and recommend core engineering as a tractable target for rationally designing kinetic stability that may be widely applicable.

Alteration of Chain-Length Selectivity and Thermostability of Rhizopus oryzae Lipase via Virtual Saturation Mutagenesis Coupled with Disulfide Bond Design

Rhizopus oryzae lipase (ROL) is one of the most important enzymes used in the food, biofuel, and pharmaceutical industries. However, the highly demanding conditions of industrial processes can reduce its stability and activity. To seek a feasible method to improve both the catalytic activity and the thermostability of this lipase, first, the structure of ROL was divided into catalytic and noncatalytic regions by identifying critical amino acids in the crevice-like binding pocket. Second, a mutant screening library aimed at improvement of ROL catalytic performance by virtual saturation mutagenesis of residues in the catalytic region was constructed based on Rosetta's Cartesian_ddg protocol. A double mutant, E265V/S267W (with an E-to-V change at residue 265 and an S-to-W change at residue 267), with markedly improved catalytic activity toward diverse chain-length fatty acid esters was identified. Then, computational design of disulfide bonds was conducted for the noncatalytic amino acids of E265V/S267W, and two potential disulfide bonds, S61C-S115C and E190C-E238C, were identified as candidates. Experimental data validated that the variant E265V/S267W/S61C-S115C/E190C-E238C had superior stability, with an increase of 8.5°C in the melting temperature and a half-life of 31.7 min at 60°C, 4.2-fold longer than that of the wild-type enzyme. Moreover, the variant improved the lipase activity toward five 4-nitrophenyl esters by 1.5 to 3.8 times, exhibiting a potential to modify the catalytic efficiency. IMPORTANCE Rhizopus oryzae lipase (ROL) is very attractive in biotechnology and industry as a safe and environmentally friendly biocatalyst. Functional expression of ROL in Escherichia coli facilitates effective high-throughput screening for positive variants. This work highlights a method to improve both selectivity and thermostability based on a combination of virtual saturation mutagenesis in the substrate pocket and disulfide bond prediction in the noncatalytic region. Using the method, ROL thermostability and activity to diverse 4-nitrophenyl esters could be substantially improved. The strategy of rational introduction of multiple mutations in different functional domains of the enzyme is a great prospect in the modification of biocatalysts.

Wide substrate range for a candidate bioremediation enzyme isolated from Nocardioides sp. strain SG-4 G

Narrow substrate ranges can impact heavily on the range of applications and hence commercial viability of candidate bioremediation enzymes. Here we show that an ester hydrolase from Nocardioides strain SG-4 G has potential as a bioremediation agent against various pollutants that can be detoxified by hydrolytic cleavage of some carboxylester, carbamate, or amide linkages. Previously we showed that a radiation-killed, freeze-dried preparation (ZimA) of this strain can rapidly degrade the benzimidazole fungicide carbendazim due to the activity of a specific ester hydrolase, MheI. Here, we report that ZimA also has substantial hydrolytic activity against phthalate diesters (dimethyl, dibutyl, and dioctyl phthalate), anilide (propanil and monalide), and carbamate ester (chlorpropham) herbicides under laboratory conditions. The reaction products are substantially less toxic, or inactive as herbicides, than the parent compounds. Tests of strain SG-4 G and Escherichia coli expressing MheI found they were also able to hydrolyse dimethyl phthalate, propanil, and chlorpropham, indicating that MheI is principally responsible for the above activities.

Proteome-Wide Structural Computations Provide Insights into Empirical Amino Acid Substitution Matrices

The relative contribution of mutation and selection to the amino acid substitution rates observed in empirical matrices is unclear. Herein, we present a neutral continuous fitness-stability model, inspired by the Arrhenius law (qij=aije-ΔΔGij). The model postulates that the rate of amino acid substitution (i→j) is determined by the product of a pre-exponential factor, which is influenced by the genetic code structure, and an exponential term reflecting the relative fitness of the amino acid substitutions. To assess the validity of our model, we computed changes in stability of 14,094 proteins, for which 137,073,638 in silico mutants were analyzed. These site-specific data were summarized into a 20 square matrix, whose entries, ΔΔGij, were obtained after averaging through all the sites in all the proteins. We found a significant positive correlation between these energy values and the disease-causing potential of each substitution, suggesting that the exponential term accurately summarizes the fitness effect. A remarkable observation was that amino acids that were highly destabilizing when acting as the source, tended to have little effect when acting as the destination, and vice versa (source → destination). The Arrhenius model accurately reproduced the pattern of substitution rates collected in the empirical matrices, suggesting a relevant role for the genetic code structure and a tuning role for purifying selection exerted via protein stability.

Affinity and Stability Analysis of Yeast Displayed Proteins

Yeast surface display is a powerful protein engineering technology that is extensively used to improve various properties of proteins, including affinity, specificity, and stability or even to add novel functions (usually ligand binding). Apart from its robustness and versatility as an engineering tool, yeast display offers a further critical advantage: Once the selection campaign is finished, usually resulting in an oligoclonal pool, these enriched protein variants can be analyzed individually on the surface of yeast without the need for any sub-cloning, soluble expression, and purification. Here, we provide detailed protocols for determining both the affinity and the thermal stability of yeast displayed proteins. In addition, we discuss the advantages, challenges, and potential pitfalls associated with affinity and stability analysis using yeast surface display.