Saturation Mutagenesis Research Articles

Engineering a human P2X2 receptor with altered ligand selectivity in yeast

P2X receptors are a family of ligand gated ion channels found in a range of eukaryotic species including humans but are not naturally present in the yeast Saccharomyces cerevisiae. We demonstrate the first recombinant expression and functional gating of the P2X2 receptor in baker’s yeast. We leverage the yeast host for facile genetic screens of mutant P2X2 by performing site saturation mutagenesis at residues of interest, including SNPs implicated in deafness and at residues involved in native binding. Deep mutational analysis and rounds of genetic engineering yield mutant P2X2 F303Y A304W, which has altered ligand selectivity towards the ATP analog AMP-PNP. The F303Y A304 variant shows over 100-fold increased intracellular calcium amplitudes with AMP-PNP compared to the WT receptor and has a much lower desensitization rate. Since AMP-PNP does not naturally activate P2X receptors, the F303Y A304 P2X2 may be a starting point for downstream applications in chemogenetic cellular control. Interestingly, the A304W mutation selectively destabilizes the desensitized state, which may provide a mechanistic basis for receptor opening with suboptimal agonists. The yeast system represents an inexpensive, scalable platform for ion channel characterization and engineering by circumventing the more expensive and time-consuming methodologies involving mammalian hosts.

Abstract 377: Saturation mutagenesis of KEAP1 to identify pathogenic mutations in lung cancer patients and gain insight into KEAP1 molecular function

Abstract Introduction: KEAP1 and NFE2L2 (NFR2) mutations occur in roughly 20% of non-small cell lung cancer (NSCLC) patients and are associated with resistance to a broad range of therapies including radiotherapy, chemotherapy, targeted therapies, and immunotherapy. In prior work, we demonstrated that not all KEAP1 mutations observed in NSCLC patients are pathogenic. These benign mutations have a wild-type phenotype that does not impart therapy resistance. However, only roughly 100 of ~13,000 (<1%) possible KEAP1 single amino acid substitution mutations have been characterized. Here, we performed saturation mutagenesis of KEAP1 to comprehensively evaluate all possible point mutation phenotypes in vitro. Methods: Saturation mutagenesis was accomplished using the High-throughput Mutagenesis by Integrated TilEs (MITE) method. In brief, the KEAP1 sequence was divided into 21 tiles each containing a 90 bp variable region flanked by a 30 bp constant region. At each codon position, mutations representing the 20 possible natural amino acids and a stop codon were included for a total of 13,076 KEAP1 variants. The tiles were assembled such that each full length KEAP1 construct contained a single desired mutation. The KEAP1 MITE lentivirus library was transduced into H1299-KEAP1NULL cells that were then treated with vehicle or H2O2 for 24 days. The KEAP1 ORF was amplified from gDNA and mutations were quantified by next generation sequencing. Mutation enrichment analysis was used to derive functional annotations. Results: Oxidative stress resulted in a bimodal distribution of KEAP1 mutations. Cells with KEAP1 nonsense mutations were enriched while cells with silent mutations were depleted. Using a threshold of 95% specificity for silent mutations to distinguish benign from pathogenic mutations, 94% of nonsense mutations and 88% of previously characterized pathogenic mutations were identified as pathogenic. While the majority of KEAP1 mutations were benign (64%), analysis of mutations found in TCGA NSCLC patients showed an enrichment of pathogenic mutations (72%). Overlaying mutation calls onto the KEAP1 structure showed enrichment of pathogenic mutations in the BTB homodimerization interface, hydrophobic core of the oxidative stress sensing IVR domain, and NRF2 binding pocket of the Kelch domain. Conclusions: Saturation mutagenesis of KEAP1 generated a comprehensive list of pathogenic KEAP1 single amino acid substitutions. Assignments strongly agreed with biologic function and prior literature. The assessment of KEAP1 mutation phenotypes is important as benign mutations occur in a significant subset of NSCLC patients and are not expected to impart therapy resistance. Our results will allow identification of benign and pathogenic KEAP1 mutations in clinical studies, which will enable testing of personalized therapeutic strategies for patients with KEAP1 mutations. Citation Format: Noah Kastelowitz, Shashank Shrishrimal, Soyeong Jun, Anni M. Zhang, Rui Wang, Ilayda Ilerten, Maximilian Diehn. Saturation mutagenesis of KEAP1 to identify pathogenic mutations in lung cancer patients and gain insight into KEAP1 molecular function [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 377.

Studies on Insertion/Deletion Residues for Functional Analysis and Improved Amination Activity in Meso-DAPDH from Corynebacterium glutamicum

Meso-diaminopimelate dehydrogenase (meso-DAPDH) from Corynebacterium glutamicum ATCC13032 (CgDAPDH) is a type I meso-DAPDH that shows obvious preference toward meso-diaminopimelate (meso-DAP) and exhibits almost no amination activity toward 2-keto acids. There are seven distinct conserved insertions and deletions (indels) between type I and type II meso-DAPDH. The current functional analysis of indels is not comprehensive in meso-DAPDH. Continuing from our previous work on these indels, we first examined the functions of the other indels shown as insertion residues in type I CgDAPDH. Alanine mutations in M216, T240, K289, and Q290 lost at least 40% of their activity, highlighting the importance of these four sites in CgDAPDH. Molecular dynamic analysis indicated that the four non-active sites altered the dynamic network of interactions within the protein. Subsequently, these four sites together with the previously identified indel-related residues R180, L176, and H193 were targeted by site-saturation mutagenesis to improve the amination ability of CgDAPDH toward pyruvic acid. The most significant improvement was observed with the mutant CgL176R, which showed a six-fold increase toward pyruvic acid in kcat/Km compared to wild-type CgDAPDH. Overall, our study provides new hotspots and ideas for the subsequent protein engineering of CgDAPDH, which may also be applied to other meso-DAPDHs.

Abstract 1689: CRISPR functional genetic screening for immune enhancers of inhibitor receptors in human primary CD8 T cells

Abstract The restrained efficacy of T cell therapies in cancer treatment is associated with T cell exhaustion (TEX), an impaired effector function state that is marked by distinct epigenetic landscape. These unique epigenetic features promoted the progression of TEX, hindered reinvigoration of TEX state and controlled the abnormal expression of inhibitor receptors (IRs), which could serve as target sites in cancer treatment. Thus, deciphering the regulatory elements in primary human T cells will deeply reveal the mechanism controlling TEX and provide novel therapeutic strategies to promote clinical response. However, how enhancers can be effectively applied in anti-cancer treatment is still poorly understood because high throughput identifying the functional enhancers of inhibitor receptors at scale is technically limited in human T cells. Here, we first developed a pooled screening platform that employs a compact CRISPRi system and enabled efficient one-step delivery of pooled screening libraries to primary human T cells. Leveraging this platform, we screened large genomic regions (~300kb) at saturation around key TEX genes, two IRs—PDCD1, HAVCR2 and one transcription factor (TF)—TBX21 to interrogate functional enhancers in human primary CD8 T cells. We identified two enhancers for PDCD1, three enhancers for HAVCR2 and three enhancers for TBX21, all of which are located over 12kb away from each gene’s transcription start sites. Then, to parse the essential sequences that fulfill enhancer function, we employed Cas9-based in-situ saturation mutagenesis screening to the identified enhancers of PDCD1 and HAVCR2. We found disrupt essential sequences impaired the entire enhancer function, quantitatively phenocopied the effect of targeting the entire enhancer with CRISPRi. Transcription factor motifs analyses of these sequences revealed potential regulation of ZFX for controlling PDCD1 expression and STAT for HAVCR2 expression. Finally, we explored the functional consequences of applying enhancers to anti-cancer treatment. Indeed, targeting the identified enhancers of PDCD1 and HAVCR2 in CAR T cells lead to enhanced cancer cell killing efficiency in vitro. Together, our study provides a platform and paradigm to interrogate functional enhancers in human T cells, illuminated the great potential of manipulating enhancers in T cell therapies. Citation Format: Ping Wang, Josephine R. Giles, Hua Huang, Junwei Shi, E.John Wherry. CRISPR functional genetic screening for immune enhancers of inhibitor receptors in human primary CD8 T cells [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1689.

Abstract 1220: Deep mutational scanning of EGFR reveals potential domain-specific TKI sensitivities in lung cancer and glioblastoma

Abstract The epidermal growth factor receptor, EGFR, is frequently activated in lung cancer and glioblastoma by genomic alterations including missense mutations. The different mutation spectra in these diseases are reflected in divergent responses to EGFR inhibition: significant patient benefit in lung cancer, but limited in glioblastoma. Here, we report a comprehensive mutational analysis of EGFR function. We performed saturation mutagenesis of EGFR and assessed function of ~22,500 variants in a human EGFR-dependent lung cancer cell line. This approach revealed enrichment of erlotinib-insensitive variants of known and unknown significance in the dimerization, transmembrane, and kinase domains. Multiple EGFR extracellular domain variants, lacking approved targeted therapies, were sensitive to dacomitinib. Strikingly, at least two glioblastoma patients with G598V dimerization domain mutations showed responses to dacomitinib treatment together with subsequent off-target resistance in one case. In summary, this comprehensive screen reveals novel functional EGFR variants and suggests broader clinical investigation of EGFR inhibition for cancers harboring extracellular domain mutations. Citation Format: Tikvah Hayes, Elisa Aquilanti, Erica Kim, Nicole Persky, Lisa Brenan, Amy Goodale, Ted Sharpe, Doug Alan, Robert Shue, Lior Golomb, Brianna Silverman, Yvonne Li, Lindsay Westlake, Federica Piccioni, Andrew Cherniack, Xiaoping Yang Yang, David Root, J. Kevin Hicks, Andrew Chi, Jorg Dietrich, Daniel Cahill, Cory Johannessen, Tracy Batchelor, Matthew Meyerson. Deep mutational scanning of EGFR reveals potential domain-specific TKI sensitivities in lung cancer and glioblastoma [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2024; Part 1 (Regular Abstracts); 2024 Apr 5-10; San Diego, CA. Philadelphia (PA): AACR; Cancer Res 2024;84(6_Suppl):Abstract nr 1220.

Engineering the activity and thermostability of a carboxylic acid reductase in the conversion of vanillic acid to vanillin

Vanillin is a valuable natural product that can be used as a fragrance and additive. Recent research in the biosynthesis of vanillin has brought attention to a key enzyme, carboxylic acid reductase (CAR), which catalyzes the reduction of vanillic acid to vanillin. Nevertheless, the biosynthesis of vanillin is hampered by the low activity and stability of CAR. As such, a rational design campaign was conducted on a well-documented carboxylic acid reductase from Segniliparus rugosus (SrCAR), using vanillic acid as the model substrate. After combined active site saturation and iterative site-specific mutagenesis, the best quadruple mutant N292H/K524S/A627L/E1121W (M3) was successfully obtained. In comparison to the wildtype SrCAR, M3 demonstrated a 4.2-fold increase in catalytic efficiency (kcat/Km), and its half-life (t1/2) was enhanced by 3.8 times up to 385.08 minutes at 40 °C. In silico docking and molecular dynamics simulation provided insights into the improved activity and stability. In the subsequent preparative-scale reaction with 100 mM (16.8 g L−1) vanillic acid, the whole cell catalysis utilizing M3 produced 10.15 g·L−1 of vanillin and 1.11 g·L−1 of vanillyl alcohol, respectively. This work demonstrates a dual improvement in the activity and thermal stability of SrCAR, thereby potentially facilitating the application of carboxylic acid reductase in the biosynthesis of vanillin.

Plant enhancers exhibit both cooperative and additive interactions among their functional elements.

Enhancers are cis-regulatory elements that shape gene expression in response to numerous developmental and environmental cues. In animals, several models have been proposed to explain how enhancers integrate the activity of multiple transcription factors. However, it remains largely unclear how plant enhancers integrate transcription factor activity. Here, we use Plant STARR-seq to characterize 3 light-responsive plant enhancers-AB80, Cab-1, and rbcS-E9-derived from genes associated with photosynthesis. Saturation mutagenesis revealed mutations, many of which clustered in short regions, that strongly reduced enhancer activity in the light, in the dark, or in both conditions. When tested in the light, these mutation-sensitive regions did not function on their own; rather, cooperative interactions with other such regions were required for full activity. Epistatic interactions occurred between mutations in adjacent mutation-sensitive regions, and the spacing and order of mutation-sensitive regions in synthetic enhancers affected enhancer activity. In contrast, when tested in the dark, mutation-sensitive regions acted independently and additively in conferring enhancer activity. Taken together, this work demonstrates that plant enhancers show evidence for both cooperative and additive interactions among their functional elements. This knowledge can be harnessed to design strong, condition-specific synthetic enhancers.

Rationally Designed Novel Antimicrobial Peptides Targeting Chitin Synthase for Combating Soybean Phytophthora Blight.

Soybean phytophthora blight is a severe menace to global agriculture, causing annual losses surpassing USD 1 billion. Present crop loss mitigation strategies primarily rely on chemical pesticides and disease-resistant breeding, frequently surpassed by the pathogens' quick adaptive evolution. In this urgent scenario, our research delves into innovative antimicrobial peptides characterized by low drug resistance and environmental friendliness. Inhibiting chitin synthase gene activity in Phytophthora sojae impairs vital functions such as growth and sporulation, presenting an effective method to reduce its pathogenic impact. In our study, we screened 16 previously tested peptides to evaluate their antimicrobial effects against Phytophthora using structure-guided drug design, which involves molecular docking, saturation mutagenesis, molecular dynamics, and toxicity prediction. The in silico analysis identified AMP_04 with potential inhibitory activity against Phytophthora sojae's chitin synthase. Through three rounds of saturation mutagenesis, we pin-pointed the most effective triple mutant, TP (D10K, G11I, S14L). Molecular dynamic simulations revealed TP's stability in the chitin synthase-TP complex and its transmembrane mechanism, employing an all-atom force field. Our findings demonstrate the efficacy of TP in occupying the substrate-binding pocket and translocation catalytic channel. Effective inhibition of the chitin synthase enzyme can be achieved. Specifically, the triple mutant demonstrates enhanced antimicrobial potency and decreased toxicity relative to the wild-type AMP_04, utilizing a mechanism akin to the barrel-stave model during membrane translocation. Collectively, our study provides a new strategy that could be used as a potent antimicrobial agent in combatting soybean blight, contributing to sustainable agricultural practices.

Study of dual binding specificity of aptamer to ochratoxin A and norfloxacin and the development of fluorescent aptasensor in milk detection

Target specificity, one of aptamer characteristics that determine recognition efficiency of biosensors, is generally considered to be an intrinsic property of aptamer. However, a high-affinity aptamer may have additional target binding specificity, little is known about the specificity of aptamer binding to multiple targets, which may result in false-positive results that hinder the accuracy of detection. Herein, an aptamer OBA3 with dual target ochratoxin A (OTA) and norfloxacin (NOR) was used as an example to explore the binding specificity mechanism and developed rapid fluorescent aptasensing methods. The nucleotide 15th T of aptamer OBA3 was demonstrated to be critical for specificity and affinity binding of target OTA via site-saturation mutagenesis. Substituting the 15th T base for C base could directly improve recognition specificity of aptamer for NOR and remove the binding affinity for OTA. The combination of π-π stacking interactions, salt bridges and hydrogen bonds between loop pocket of aptamer and quinolone skeleton, piperazinyl group may contributes to the fluoroquinolone antibiotics (NOR and difloxacin)-aptamer recognition interaction. Based on this understanding, a dual-aptamer fluorescent biosensor was fabricated for simultaneous detection of OTA and NOR, which has a linear detection range of 50–6000 nM with a detection limit of 31 nM for OTA and NOR. Combined with T15C biosensor for eliminating interference of OTA, the assay was applied to milk samples with satisfactory recovery (94.06–100.93%), which can achieve detection of OTA and NOR individually within 40 min.

Integrated engineering at distal site and active center regulates stereoselectivity and activity of carbonyl reductase towards N-Boc-pyrrolidone

Green and efficient biosynthesis of two enantiomeric products of N-Boc-pyrrolidinol (NPBL) as key pharmaceutical blocks for treatment of cancer and HIV has received much attention. The precursor, N-Boc-pyrrolidone (NPBO), is considered as a difficult-to-reduce ketone by carbonyl reductases, because it has a bulky Boc group and sterically similar substitutions on either side of the carbonyl group resulting in low stereoselectivity and conversion. Moreover, activity enhancement concomitant with reversal of stereoselectivity is challenging. Carbonyl reductase (CpCR) from Candida parapsilosis shows 90 %ee for (S)-NBPL, but low catalytic efficiency. In this study, we considered both activity and selectivity, and proposed an integrated engineering strategy. The distal site was introduced from an activity perspective and the active pocket sites were introduced by virtual saturation mutagenesis from a stereoselectivity perspective. Mutants L34A/W116A (95 %ee (S)) and L34A/W116T/F285C/W286S (94 %ee (R)) were obtained with 26.8-fold and 2.9-fold higher catalytic efficiencies, respectively. Molecular dynamics simulations revealed the mechanism of stereoselective flip-flop and activity enhancement; mutation altered the substrate-binding mode and changed the shape and size of the cavities, thus contributing to the change in the ratio of active conformations in the pre-reactive state. This study, which differs from the traditional method of exchanging large and small pockets, provides a viable approach for the rational design of carbonyl reductases with high stereoselectivity for target substrates.

Engineering Site 228 of Streptomyces coelicolor Laccase for Optimizing Catalytic Activity.

Malachite green (MG) poses a formidable threat to ecosystems and human health. Laccase emerges as a promising candidate for MG degradation, prompting an investigation into the catalytic activity modulation of a small laccase (SLAC) from Streptomyces coelicolor, with a focus on amino acid position 228. Through saturation mutagenesis, five mutants with a 50% increase in the specific activity were generated. Characterization revealed notable properties, Km of E228F was 8.8% of the wild type (WT), and E288T exhibited a 133% kcat compared to WT. Structural analyses indicated improved hydrophobicity and electrostatic potential on the mutants' surfaces, with the stable E228F-ABTS complex exhibiting reduced flexibility, possibly contributing to the observed decrease in turnover rate. Mutants demonstrated enhanced MG decolorization, particularly E228G. Site 228 acts as a crucial functional control switch, suggesting its potential role in SLAC engineering. This study provides insights into laccase modulation and offers promising avenues for enzymatic bioremediation applications.

Engineering the Propeptide of Microbial Transglutaminase Zymogen: Enabling Substrate-Dependent Activation for Bioconjugation Applications.

Microbial transglutaminase (MTG) from Streptomyces mobaraensis is a powerful biocatalytic glue for site-specific cross-linking of a range of biomolecules and synthetic molecules that have an MTG-reactive moiety. The preparation of active recombinant MTG requires post-translational proteolytic digestion of a propeptide that functions as an intramolecular chaperone to assist the correct folding of the MTG zymogen (MTGz) in the biosynthesis. Herein, we report engineered active zymogen of MTG (EzMTG) that is expressed in soluble form in the host Escherichia coli cytosol and exhibits cross-linking activity without limited proteolysis of the propeptide. We found that the saturation mutagenesis of residues K10 or Y12 in the propeptide domain generated several active MTGz mutants. In particular, the K10D/Y12G mutant exhibited catalytic activity comparable to that of mature MTG. However, the expression level was low, possibly because of decreased chaperone activity and/or the promiscuous substrate specificity of MTG, which is potentially harmful to the host cells. The K10R/Y12A mutant exhibited specific substrate-dependent reactivity toward peptidyl substrates. Quantitative analysis of the binding affinity of the mutated propeptides to the active site of MTG suggested an inverse relationship between the binding affinity and the catalytic activity of EzMTG. Our proof-of-concept study provides insights into the design of a new biocatalyst using the MTGz as a scaffold and a potential route to high-throughput screening of EzMTG mutants for bioconjugation applications.

Increase in the solubility of uvsY using a site saturation mutagenesis library for application in a lyophilized reagent for recombinase polymerase amplification.

Recombinase uvsY from bacteriophage T4, along with uvsX, is a key enzyme for recombinase polymerase amplification (RPA), which is used to amplify a target DNA sequence at a constant temperature. uvsY, though essential, poses solubility challenges, complicating the lyophilization of RPA reagents. This study aimed to enhance uvsY solubility. Our hypothesis centered on the C-terminal region of uvsY influencing solubility. To test this, we generated a site-saturation mutagenesis library for amino acid residues Lys91-Glu134 of the N-terminal (His)6-tagged uvsY. Screening 480 clones identified A116H as the variant with superior solubility. Lyophilized RPA reagents featuring the uvsY variant A116H demonstrated enhanced performance compared to those with wild-type uvsY. The uvsY variant A116H emerges as an appealing choice for RPA applications, offering improved solubility and heightened lyophilization feasibility.

Engineering A-type Dye-Decolorizing Peroxidases by Modification of a Conserved Glutamate Residue.

Dye-decolorizing peroxidases (DyPs) are recently identified microbial enzymes that have been used in several Biotechnology applications from wastewater treatment to lignin valorization. However, their properties and mechanism of action still have many open questions. Their heme-containing active site is buried by three conserved flexible loops with a putative role in modulating substrate access and enzyme catalysis. Here, we investigated the role of a conserved glutamate residue in stabilizing interactions in loop 2 of A-type DyPs. First, we did site saturation mutagenesis of this residue, replacing it with all possible amino acids in bacterial DyPs from Bacillus subtilis (BsDyP) and from Kitasatospora aureofaciens (KaDyP1), the latter being characterized here for the first time. We screened the resulting libraries of variants for activity towards ABTS and identified variants with increased catalytic efficiency. The selected variants were purified and characterized for activity and stability. We furthermore used Molecular Dynamics simulations to rationalize the increased catalytic efficiency and found that the main reason is the electron channeling becoming easier from surface-exposed tryptophans. Based on our findings, we also propose that this glutamate could work as a pH switch in the wild-type enzyme, preventing intracellular damage.

Mutability landscape guided engineering of a promiscuous microbial glycosyltransferase for regioselective synthesis of salidroside and icariside D2

Microbial glycosyltransferases efficiently synthesize glucosides and have garnered increasing interest. However, limited regioselectivity has impeded their broad application, particularly in the pharmaceutical industry. In this study, the UDP-glycosyltransferase YjiC from Bacillus licheniformis (BlYjiC) was engineered to achieve the bidirectional regioselective glycosylation of tyrosol and its derivatives. Initially, site-directed saturation mutagenesis was performed on two newly identified substrate-binding cavities in the acceptor pocket of BlYjiC to provide a comprehensive blueprint of the interplay between mutations and function (mutability landscape). Iterative saturation mutagenesis was performed, guided by the mutability landscape. Two highly regioselective mutants M6 (M112L/I325Y/L70R/Q136E/I67E/M77R) and M2’ (M112D/I62L) were generated, exhibiting >99 % regioselectivity toward the alcoholic and phenolic hydroxyl of tyrosol, respectively, compared with the wild-type (product mixture: 51:49 %). Both mutants exhibited excellent regioselectivity toward several dihydroxy phenolic substrates, offering valuable biocatalysts for the regioselective synthesis of glucosides. Their application was confirmed in a short synthesis of salidroside (3.6 g/L) and icariside D2 (2.4 g/L), which exhibited near-perfect regioselectivity. This study provides valuable insights into future protein engineering of similar enzymes and opens new avenues for their practical applications.

Roles of the N-terminal motif in improving the activity and soluble expression of phenylalanine ammonia lyases in Escherichia coli

Phenylalanine ammonia-lyase (PAL) has various applications in fine chemical manufacturing and the pharmaceutical industry. In particular, PAL derived from Anabaena variabilis (AvPAL) is used as a therapeutic agent to the treat phenylketonuria in clinical settings. In this study, we aligned the amino acid sequences of AvPAL and PAL derived from Nostoc punctiforme (NpPAL) to obtain several mutants with enhanced activity, expression yield, and thermal stability via amino acid substitution and saturation mutagenesis at the N-terminal position. Enzyme kinetic experiments revealed that the kcat values of NpPAL-N2K, NpPAL-I3T, and NpPAL-T4L mutants were increased to 3.2-, 2.8-, and 3.3-fold that of the wild-type, respectively. Saturation mutagenesis of the fourth amino acid in AvPAL revealed that the kcat values of AvPAL-L4N, AvPAL-L4P, AvPAL-L4Q and AvPAL-L4S increased to 4.0-, 3.7-, 3.6-, and 3.2-fold, respectively. Additionally, the soluble protein yield of AvPAL-L4K increased to approximately 14 mg/L, which is approximately 3.5-fold that of AvPAL. Molecular dynamics studies further revealed that maintaining the attacking state of the reaction and N-terminal structure increased the rate of catalytic reaction and improved the solubility of proteins. These findings provide new insights for the rational design of PAL in the future.

Thermostability Enhancement of GH 62 α-l-Arabinofuranosidase by Directed Evolution and Rational Design.

GH 62 arabinofuranosidases are known for their excellent specificity for arabinoxylan of agroindustrial residues and their synergism with endoxylanases and other hemicellulases. However, the low thermostability of some GH enzymes hampers potential industrial applications. Protein engineering research highly desires mutations that can enhance thermostability. Therefore, we employed directed evolution using one round of error-prone PCR and site-saturation mutagenesis for thermostability enhancement of GH 62 arabinofuranosidase from Aspergillus fumigatus. Single mutants with enhanced thermostability showed significant ΔΔG changes (<-2.5 kcal/mol) and improvements in perplexity scores from evolutionary scale modeling inverse folding. The best mutant, G205K, increased the melting temperature by 5 °C and the energy of denaturation by 41.3%. We discussed the functional mechanisms for improved stability. Analyzing the adjustments in α-helices, β-sheets, and loops resulting from point mutations, we have obtained significant knowledge regarding the potential impacts on protein stability, folding, and overall structural integrity.

Potts Hamiltonian Models and Molecular Dynamics Free Energy Simulations for Predicting the Impact of Mutations on Protein Kinase Stability.

Single-point mutations in kinase proteins can affect their stability and fitness, and computational analysis of these effects can provide insights into the relationships among protein sequence, structure, and function for this enzyme family. To assess the impact of mutations on protein stability, we used a sequence-based Potts Hamiltonian model trained on a kinase family multiple-sequence alignment (MSA) to calculate the statistical energy (fitness) effects of mutations and compared these against relative folding free energies (ΔΔGs) calculated from all-atom molecular dynamics free energy perturbation (FEP) simulations in explicit solvent. The fitness effects of mutations in the Potts model (ΔEs) showed good agreement with experimental thermostability data (Pearson r = 0.68), similar to the correlation we observed with ΔΔGs predicted from structure-based relative FEP simulations. Recognizing the possible advantages of using Potts models to rapidly estimate protein stability effects of kinase mutations seen in cancer genomics data, we used the Potts statistical energy model to estimate the stability effects of 65 conservative and nonconservative mutations across three distinct kinases (Wee1, Abl1, and Cdc7) with somatic mutations reported in the Genomic Data Commons (GDC) database. The ΔEs of these mutations calculated from the Potts model are consistent with the corresponding ΔΔGs from FEP simulations (Pearson ratio of 0.72). The agreement between these methods suggests that the Potts model may be used as a sequence-based tool for high-throughput screening of mutational effects as part of a computational pipeline for predicting the stability effects of mutations. We also demonstrate how the scalability of the fitness-based Potts model calculations permits analyses that are not easily accessed using FEP simulations. To this end, we employed site-saturation mutagenesis in the Potts model in order to investigate the relative stability effects of mutations seen in different cancer evolutionary scenarios. We used this approach to analyze the effects of drug pressure in Abl kinase by contrasting the relative fitness penalties of somatic mutations seen in miscellaneous cancer types with those calculated for mutations associated with cancer drug resistance. We observed that, in contrast to somatic mutations of Abl seen in various tumors that appear to have evolved neutrally, cancer mutations that evolved under drug pressure in Abl-targeted therapies tend to preserve enzyme stability.

Controlling Monoterpene Isomerization by Guiding Challenging Carbocation Rearrangement Reactions in Engineered Squalene-Hopene Cyclases.

The interconversion of monoterpenes is facilitated by a complex network of carbocation rearrangement pathways. Controlling these isomerization pathways is challenging when using common Brønsted and Lewis acid catalysts, which often produce product mixtures that are difficult to separate. In contrast, natural monoterpene cyclases exhibit high control over the carbocation rearrangement reactions but are reliant on phosphorylated substrates. In this study, we present engineered squalene-hopene cyclases from Alicyclobacillus acidocaldarius (AacSHC) that catalyze the challenging isomerization of monoterpenes with unprecedented precision. Starting from a promiscuous isomerization of (+)-β-pinene, we first demonstrate noticeable shifts in the product distribution solely by introducing single point mutations. Furthermore, we showcase the tuneable cation steering by enhancing (+)-borneol selectivity from 1 % to >90 % (>99 % de) aided by iterative saturation mutagenesis. Our combined experimental and computational data suggest that the reorganization of key aromatic residues leads to the restructuring of the water network that facilitates the selective termination of the secondary isobornyl cation. This work expands our mechanistic understanding of carbocation rearrangements and sets the stage for target-oriented skeletal reorganization of broadly abundant terpenes.

Controlling Monoterpene Isomerization by Guiding Challenging Carbocation Rearrangement Reactions in Engineered Squalene‐Hopene Cyclases

AbstractThe interconversion of monoterpenes is facilitated by a complex network of carbocation rearrangement pathways. Controlling these isomerization pathways is challenging when using common Brønsted and Lewis acid catalysts, which often produce product mixtures that are difficult to separate. In contrast, natural monoterpene cyclases exhibit high control over the carbocation rearrangement reactions but are reliant on phosphorylated substrates. In this study, we present engineered squalene‐hopene cyclases from Alicyclobacillus acidocaldarius (AacSHC) that catalyze the challenging isomerization of monoterpenes with unprecedented precision. Starting from a promiscuous isomerization of (+)‐β‐pinene, we first demonstrate noticeable shifts in the product distribution solely by introducing single point mutations. Furthermore, we showcase the tuneable cation steering by enhancing (+)‐borneol selectivity from 1 % to &gt;90 % (&gt;99 % de) aided by iterative saturation mutagenesis. Our combined experimental and computational data suggest that the reorganization of key aromatic residues leads to the restructuring of the water network that facilitates the selective termination of the secondary isobornyl cation. This work expands our mechanistic understanding of carbocation rearrangements and sets the stage for target‐oriented skeletal reorganization of broadly abundant terpenes.