Saturation Mutagenesis Research Articles

Revealing the critical role of Leucine145 of α-glucosidase AglA for enhancing α-arbutin production

The α-glucosidase AglA from Xanthomonas campestris suffers from low catalytic activity for α-arbutin production, which leads to unsatisfied productivity. To address this issue, herein, the critical amino acid residue Leucine145 of α-glucosidase AglA, which locates on the loop adjacent to maltose binding domain, is identified with potential capability of affecting the enzyme activity with the aid of computational analysis. Site-directed and saturation mutagenesis on this key residue is performed, and the best mutant L145V with 7.2-fold improvement in catalytic efficiency (kcat/Km) is achieved. Following maltose hydrolysis test and α-arbutin production by whole cell biotransformation further demonstrate that L145V has the most excellent catalytic performance. The production of α-arbutin produced by Escherichia coli BL21(DE3) cells harboring mutant L145V reaches 63.1 mM within 3 h, which is 4.7-fold higher than that of wild type AglA. Finally, molecular dynamics simulations studies provide molecular insights into the possible maltose binding poses, and shed light on the possible reasons behind the greatly enhanced activity of mutant L145V. These results reveal critical roles of residue L145 of α-glucosidase AglA in regulating catalytic activity, which will be very useful for guiding the further engineering of other glucosidases.

Improving the activity and expression level of a phthalate-degrading enzyme by a combination of mutagenesis strategies and strong promoter replacement.

Phthalic acid esters (PAEs) are widely used plasticizers found in consumer products, which enter the environment and pose severe threats to human health. Here, a new PAE-degrading enzyme EstJ6 was modified by combining mutagenesis strategies and a strong promoter replacement to improve its catalytic activity and expression level. Four mutants with enhanced activity were obtained by random mutation, among which EstJ6M1.1 exhibited the highest catalytic activity with an increase in catalytic activity by 2.9-fold toward dibutyl phthalate (DBP) than that of the wild-type (WT) enzyme. With these mutants as a template, a variant EstJ6M2 with 3.1-fold higher catalytic activity and 4.61 times higher catalytic efficiency (Kcat/Km) was identified by staggered extension PCR. Targeting four mutation sites of EstJ6M2, a variant EstJ6M3.1 was gained by site-directed saturation mutagenesis and displayed 4.3-fold higher activity and 5.97 times higher Kcat/Km than WT. The expression level of three mutants EstJ6M1.1, EstJ6M2, and EstJ6M3.1, as well as the WT, increased nearly threefold after a strong promoter replacement. These results provide a proof-theoretical basis and practicable pipeline for applying PAE-degrading enzymes.

Engineering of Hydroxymandelate Oxidase and Cascade Reactions for High-Yielding Conversion of Racemic Mandelic Acids to Phenylglyoxylic Acids and (R)- and (S)-Phenylglycines

An alcohol oxidase (AOx) for the (S)-enantioselective oxidation of 4-hydroxymandelic acid (4-HMA) to 4-hydroxyphenylglyoxylic acid (4-HPGA) with enhanced activity was developed by directed evolution of hydroxymandelate oxidase (HMO) through three rounds of iterative saturation mutagenesis. The engineered HMO mutant A80G-T159S-T162Q (HMOTM) catalyzed the oxidation of (S)-4-HMA to 4-HPGA with a 23-fold enhancement in catalytic efficiency (kcat/KM). Substrate docking simulation on HMOTM suggested that A80G reoriented FMN, while T159S and T162Q formed hydrogen bonds with the carboxylic group of the substrate, thus facilitating substrate binding and catalysis. (S)-Enantioselective HMOTM was used together with mandelic acid racemase (MR) and catalase (KatE) to achieve high-yielding oxidation of rac-4-HMA to 4-HPGA with either purified enzymes (up to 93% yield and 426 mM) or Escherichia coli (HMC) cells expressing the three enzymes (up to 93% yield and 140 mM). The HMOTM-MR-KatE cascades were applied for the oxidation of seven other substituted MAs, producing the corresponding phenylglyoxylic acids with 90–99% conversions using purified enzymes or whole cells. Efficient conversion of racemic α-hydroxy acids to (S)- or (R)-α-amino acids was achieved by combining HMOTM-MR-KatE with (S)-enantioselective transaminase (EcαTA) or (R)-enantioselective transaminase (DpgAT), together with glutamate dehydrogenase (GluDH). Coupling of E. coli (HMC) with E. coli (E-G) expressing EcαTA and GluDH for one-pot biotransformation of five rac-MAs produced the corresponding (S)-phenylglycines (PGs) in 91–99% ee with 73–98% conversions. Using E. coli (HMC) and E. coli (D-G) expressing DpgAT and GluDH enabled the production of five (R)-PGs with 93–99% ee and 73–98% conversions from the corresponding rac-MAs. The engineered AOx, AOx-MR-KatE cascades, and AOx-MR-Kat-GluDH-(S)- or (R)-TA cascades provide useful tools for highly active and enantioselective oxidation of racemic hydroxyacids and high-yielding conversion of racemic hydroxyacids to ketoacids and enantiopure (S)- or (R)-aminoacids, respectively, which are challenging and useful chemical reactions.

Versatile Strategy for the Construction of a Transcription Factor-Based Orthogonal Gene Expression Toolbox in Monascus spp.

Gene expression is needed to be conducted in an orthogonal manner and controllable independently from the host's native regulatory system. However, there is a shortage of gene expression regulatory toolboxes that function orthogonally from each other and toward the host. Herein, we developed a strategy based on the mutant library to generate orthogonal gene expression toolboxes. A transcription factor, MaR, located in the Monascus azaphilone biosynthetic gene cluster, was taken as a typical example. Nine DNA-binding residues of MaR were identified by molecular simulation and site-directed mutagenesis. We created five MaR multi-site saturation mutagenesis libraries consisting of 10743 MaR variants on the basis of five cognate promoters. A functional analysis revealed that all five tested promoters were orthogonally regulated by five different MaR variants, respectively. Furthermore, fine gene expression tunability and high signal sensitivity of this toolbox are demonstrated by introducing chemically inducible expression modules, designing synthetic promoter elements, and creating protein-protein interaction between MaRs. This study paves the way for a bottom-up approach to build orthogonal gene expression toolboxes.

Analysis of critical protein–protein interactions of SARS-CoV-2 capping and proofreading molecular machineries towards designing dual target inhibitory peptides

In recent years, the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), as the cause of the coronavirus disease (COVID-19) global pandemic, and its variants, especially those with higher transmissibility and substantial immune evasion, have highlighted the imperative for developing novel therapeutics as sustainable solutions other than vaccination to combat coronaviruses (CoVs). Beside receptor recognition and virus entry, members of the SARS-CoV-2 replication/transcription complex are promising targets for designing antivirals. Here, the interacting residues that mediate protein–protein interactions (PPIs) of nsp10 with nsp16 and nsp14 were comprehensively analyzed, and the key residues’ interaction maps, interaction energies, structural networks, and dynamics were investigated. Nsp10 stimulates both nsp14’s exoribonuclease (ExoN) and nsp16’s 2′O-methyltransferase (2′O-MTase). Nsp14 ExoN is an RNA proofreading enzyme that supports replication fidelity. Nsp16 2′O-MTase is responsible for the completion of RNA capping to ensure efficient replication and translation and escape from the host cell’s innate immune system. The results of the PPIs analysis proposed crucial information with implications for designing SARS-CoV-2 antiviral drugs. Based on the predicted shared protein–protein interfaces of the nsp16-nsp10 and nsp14-nsp10 interactions, a set of dual-target peptide inhibitors was designed. The designed peptides were evaluated by molecular docking, peptide–protein interaction analysis, and free energy calculations, and then further optimized by in silico saturation mutagenesis. Based on the predicted evolutionary conservation of the interacted target residues among CoVs, the designed peptides have the potential to be developed as dual target pan-coronavirus inhibitors.

Rationally Controlling Selective Steroid Hydroxylation via Scaffold Sampling of a P450 Family

Many steroids are important pharmaceutically active compounds, while cytochrome P450 monooxygenases (CYPs) are attractive enzymes for applications in steroidal drug synthesis. However, the catalytic efficiency of existing P450s is not routinely high enough, as well as the molecular basis for selectivity control is unclear, which severely restrict their real applications. Here, a 16β steroid-hydroxylase CYP109B4 from Bacillus sonorensis is identified with excellent selectivity and activity. The crystallization and structural analysis of CYP109B4 reveal potential three “hotspot” residues (V84, V292, and S387) responsible for selectivity control. Then, guided by the sequence–function relationships revealed from the mutability landscape construction on the three residues, focused rational iterative site-specific mutagenesis (FRISM) and limited iterative saturation mutagenesis were performed, which provide variant B4-M7 (L240V/S387F/V84L/V292S/I291T/M290F/F294I) with completely switched regioselectivity from 16β to 15β. The subsequent computational analysis uncovers insights into the substrate binding modes in CYP109B4 and its variants, which further confirms the critical role of the “hotspot” residues for selectivity control. Finally, the generality of conserved-“hotspots”-mediated selectivity control is demonstrated by performing scaffold sampling between a panel of CYP109B members. Overall, in addition to the present chemical results, our study provides guidance in rationally designing more excellent P450 biocatalysts for potential practical (industrial) applications.

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.

Structure of an Alkaline Pectate Lyase and Rational Engineering with Improved Thermo-Alkaline Stability for Efficient Ramie Degumming.

Alkaline pectate lyases have biotechnological applications in plant fiber processing, such as ramie degumming. Previously, we characterized an alkaline pectate lyase from Bacillus clausii S10, named BacPelA, which showed potential for enzymatic ramie degumming because of its high cleavage activity toward methylated pectins in alkaline conditions. However, BacPelA displayed poor thermo-alkaline stability. Here, we report the 1.78 Å resolution crystal structure of BacPelA in apo form. The enzyme has the characteristic right-handed β-helix fold of members of the polysaccharide lyase 1 family and shows overall structural similarity to them, but it displays some differences in the details of the secondary structure and Ca2+-binding site. On the basis of the structure, 10 sites located in flexible regions and showing high B-factor and positive ΔTm values were selected for mutation, aiming to improve the thermo-alkaline stability of the enzyme. Following site-directed saturation mutagenesis and screening, mutants A238C, R150G, and R216H showed an increase in the T5015 value at pH 10.0 of 3.0 °C, 6.5 °C, and 7.0 °C, respectively, compared with the wild-type enzyme, interestingly accompanied by a 24.5%, 46.6%, and 61.9% increase in activity. The combined mutant R150G/R216H/A238C showed an 8.5 °C increase in the T5015 value at pH 10.0, and an 86.1% increase in the specific activity at 60 °C, with approximately doubled catalytic efficiency, compared with the wild-type enzyme. Moreover, this mutant retained 86.2% activity after incubation in ramie degumming conditions (4 h, 60 °C, pH 10.0), compared with only 3.4% for wild-type BacPelA. The combined mutant increased the weight loss of ramie fibers in degumming by 30.2% compared with wild-type BacPelA. This work provides a thermo-alkaline stable, highly active pectate lyase with great potential for application in the textile industry, and also illustrates an effective strategy for rational design and improvement of pectate lyases.

Potent De Novo Macrocyclic Peptides That Inhibit O‐GlcNAc Transferase through an Allosteric Mechanism

AbstractGlycosyltransferases are a superfamily of enzymes that are notoriously difficult to inhibit. Here we apply an mRNA display technology integrated with genetic code reprogramming, referred to as the RaPID (random non‐standard peptides integrated discovery) system, to identify macrocyclic peptides with high binding affinities for O‐GlcNAc transferase (OGT). These macrocycles inhibit OGT activity through an allosteric mechanism that is driven by their binding to the tetratricopeptide repeats of OGT. Saturation mutagenesis in a maturation screen using 39 amino acids, including 22 non‐canonical residues, led to an improved unnatural macrocycle that is ≈40 times more potent than the parent compound (Kiapp=1.5 nM). Subsequent derivatization delivered a biotinylated derivative that enabled one‐step affinity purification of OGT from complex samples. The high potency and novel mechanism of action of these OGT ligands should enable new approaches to elucidate the specificity and regulation of OGT.

Potent De Novo Macrocyclic Peptides That Inhibit O-GlcNAc Transferase through an Allosteric Mechanism.

Glycosyltransferases are a superfamily of enzymes that are notoriously difficult to inhibit. Here we apply an mRNA display technology integrated with genetic code reprogramming, referred to as the RaPID (random non-standard peptides integrated discovery) system, to identify macrocyclic peptides with high binding affinities for O-GlcNAc transferase (OGT). These macrocycles inhibit OGT activity through an allosteric mechanism that is driven by their binding to the tetratricopeptide repeats of OGT. Saturation mutagenesis in a maturation screen using 39 amino acids, including 22 non-canonical residues, led to an improved unnatural macrocycle that is ≈40 times more potent than the parent compound (Ki app =1.5 nM). Subsequent derivatization delivered a biotinylated derivative that enabled one-step affinity purification of OGT from complex samples. The high potency and novel mechanism of action of these OGT ligands should enable new approaches to elucidate the specificity and regulation of OGT.

Improving the activity of creatinase from Alcaligenes sp. KS-85 through semi-rational design

Creatinine levels in biological fluids are important indicators for the clinical evaluation of renal function. Creatinase (CRE, EC3.5.3.3) is one of the key enzymes in the enzymatic measurement of creatinine concentration, and it is also the rate-limiting enzyme in the whole enzymatic cascade system. The poor catalytic activity of CRE severely limits its clinical and industrial applications. To address this issue, a semi-rational design is applied to increase the activity of a creatinase from Alcaligenes sp. KS-85 (Al-CRE). By high-throughput screen of saturation mutagenesis libraries on the selected hotspot mutations, multiple variant enzymes with increased activity are obtained. The five-point best variant enzyme (I304L/F395V/K351V/Y63S/Q88A) were further obtained by recombine the improved mutations sites that to showed a 2.18-fold increased specific activity. Additionally, structure analysis is conducted to understand the mechanism of the activity change. This study paves the way for a better practical application of creatinase and may help further understand its catalytic mechanism.

Ter-Seq: A high-throughput method to stabilize transient ternary complexes and measure associated kinetics.

Regulation of biological processes by proteins often involves the formation of transient, multimeric complexes whose characterization is mechanistically important but challenging. The bacterial toxin CcdB binds and poisons DNA Gyrase. The corresponding antitoxin CcdA extracts CcdB from its complex with Gyrase through the formation of a transient ternary complex, thus rejuvenating Gyrase. We describe a high throughput methodology called Ter-Seq to stabilize probable ternary complexes and measure associated kinetics using the CcdA-CcdB-GyrA14 ternary complex as a model system. The method involves screening a yeast surface display (YSD) saturation mutagenesis library of one partner (CcdB) for mutants that show enhanced ternary complex formation. We also isolated CcdB mutants that were either resistant or sensitive to rejuvenation, and used surface plasmon resonance (SPR) with purified proteins to validate the kinetics measured using the surface display. Positions, where CcdB mutations lead to slower rejuvenation rates, are largely involved in CcdA-binding, though there were several notable exceptions suggesting allostery. Mutations at these positions reduce the affinity towards CcdA, thereby slowing down the rejuvenation process. Mutations at GyrA14-interacting positions significantly enhanced rejuvenation rates, either due to reduced affinity or complete loss of CcdB binding to GyrA14. We examined the effect of different parameters (CcdA affinity, GyrA14 affinity, surface accessibilities, evolutionary conservation) on the rate of rejuvenation. Finally, we further validated the Ter-Seq results by monitoring the kinetics of ternary complex formation for individual CcdB mutants in solution by fluorescence resonance energy transfer (FRET) studies.

Directed evolution of a β-N-acetylhexosaminidase from Haloferula sp. for lacto-N-triose II and lacto-N-neotetraose synthesis from chitin

In our previous study, a β-N-acetylhexosaminidase (HaHex74) from Haloferula sp. showing high human milk oligosaccharides (HMOs) synthesis ability was identified and characterized. In this study, HaHex74 was further engineered by directed evolution and site-saturation mutagenesis to improve its transglycosylation activity for HMOs synthesis. A mutant (mHaHex74) with improved transglycosylation activity (HaHex74-Asn401Ile/His394Leu) was obtained and characterized. mHaHex74 exhibited maximal activity at pH 5.5 and 35 °C, respectively, which were distinct from that of HaHex74 (pH 6.5 and 45 °C). Moreover, mHaHex74 showed the highest LNT2 conversion ratio of 28.2% from N,N’-diacetyl chitobiose (GlcNAc2), which is 2.2 folds higher than that of HaHex74. A three-enzyme cascade reaction for the synthesis of LNT2 and LNnT from chitin was performed in a 5–L reactor, and the contents of LNT2 and LNnT reached up to 15.0 g L1 and 4.9 g L1, respectively. Therefore, mHaHex74 maybe a good candidate for enzymatic synthesis of HMOs.

Improving the Cβ stereoselectivity of L-threonine aldolase for the preparation of L-threo-3,4-dihydroxyphenylserine, a powerful anti-Parkinson's disease drug

L-threo-DOPS (L-threo-3,4-dihydroxyphenylserine) is a powerful anti-Parkinson's disease drug approved by FDA in 2014. The industrial production of optical L-threo-DOPS relied on optical resolution after chemical reactions, which did not fit into the atom-economy principle. L-Threonine aldolase (L-TA) could catalyze the synthesis of L-threo-DOPS from glycine and 3,4-dihydroxybenzaldehyde in one step, which is a promising alternative route. However, the substrate 3,4-dihydroxybenzaldehyde bearing an electron-rich catechol group is inactive for aldol condensation reaction, resulting in poor Cβ diastereoselectivity and low yields in enzymatic synthesis. In this work, virtual protein mutagenesis revealed the importance of His128 of L-TA to its binding with L-threo-DOPS intermediate. Further, the site-saturation mutagenesis of His128 provided a high-diastereoselectivity mutant H128N with 92.9 % de in the synthesis of L-threo-DOPS. With the whole cells harboring H128N, L-threo-DOPS was synthesized in a great de value of 95.4 % and a yield of 1.6 g/L in 1 h, which was significantly superior to the best-ever enzyme (55.4 % de). Based on the mutagenesis result and in silico docking, it was suggested that the enlarged volume and active hydrogen donor of the side chain of the amino acid 128 were disadvantageous, while proper hydrogen receptors could enhance substrate binding, leading to improved conversion.

Generation of an anti-idiotypic affibody-based masking domain for conditional activation of EGFR-targeting

Conditional activation of engineered affinity proteins by proteolytic processing is an interesting approach for a wide range of applications. We have generated an anti-idiotypic masking domain with specificity for the binding surface of an EGFR-targeting affibody molecule using an in-house developed staphylococcal display method. The masking domain could specifically abrogate EGFR-binding on cancer cells when fused to the EGFR-targeting affibody molecule via a linker comprising a protease cleavage site. EGFR-binding was restored by proteolytic cleavage of the linker region resulting in release of the masking domain. A saturation mutagenesis study provided detailed information on the interaction between the EGFR-targeting affibody molecule and the masking domain. Introducing an anti-idiotypic masking affibody domain is a viable approach for blocking EGFR-binding and allows for conditional activation by proteolytic processing. The results warrant further studies evaluating the therapeutic and diagnostic applicability both in vitro and in vivo.

Reshaping the Binding Pocket of Cellobiose 2-Epimerase for Improved Substrate Affinity and Isomerization Activity for Enabling Green Synthesis of Lactulose.

Enzymatic isomerization of lactose into lactulose via cellobiose 2-epimerase (CE) could provide an eco-friendly route for the industrial production of lactulose, a valuable food prebiotic. However, poor substrate affinity for lactose and preference for epimerization over isomerization hinder this application. Previous studies on CE improvement have focused on random mutagenesis or active site rational design; little is known about the relationship between substrate binding and enzyme efficacy, which was hence the subject of this study. First, residues 372W and 308W were identified as key for disaccharide recognition in CEs based on crystal structure alignment of the N-acetyl-glucosamine 2-epimerase superfamily and site-directed mutation. This binding domain was then reshaped through site saturation mutagenesis, resulting in seven mutants with enhanced isomerization activity. The optimal mutant CsCE/Q371E had significantly enhanced substrate affinity (Km, 269.65 mM vs Km, 417.5 mM), reduced epimerization activity, and 3.3-fold increased isomerization activity over the original CsCE. Molecular dynamics simulation further revealed that substituting Gln-371 with Glu strengthened the hydrogen-bonding network and altered the active site-substrate interactions, increasing the substrate stability and shifting the catalytic direction. This study uncovered new information about the substrate binding region and its mechanisms and impact on CE catalytic performance, paving the way for potential commercial applications.

Mutability‐Landscape‐Guided Engineering of l‐Threonine Aldolase Revealing the Prelog Rule in Mediating Diastereoselectivity of C−C Bond Formation

Abstractl‐threonine aldolase (LTA) catalyzes C−C bond synthesis with moderate diastereoselectivity. In this study, with LTA from Cellulosilyticum sp (CpLTA) as an object, a mutability landscape was first constructed by performing saturation mutagenesis at substrate access tunnel amino acids. The combinatorial active‐site saturation test/iterative saturation mutation (CAST/ISM) strategy was then used to tune diastereoselectivity. As a result, the diastereoselectivity of mutant H305L/Y8H/V143R was improved from 37.2 %syn to 99.4 %syn. Furthermore, the diastereoselectivity of mutant H305Y/Y8I/W307E was inverted to 97.2 %anti. Based on insight provided by molecular dynamics simulations and coevolution analysis, the Prelog rule was employed to illustrate the diastereoselectivity regulation mechanism of LTA, holding that the asymmetric formation of the C−C bond was caused by electrons attacking the carbonyl carbon atom of the substrate aldehyde from the re or si face. The study would be useful to expand LTA applications and guide engineering of other C−C bond‐forming enzymes.

Mutability-Landscape-Guided Engineering of l-Threonine Aldolase Revealing the Prelog Rule in Mediating Diastereoselectivity of C-C Bond Formation.

l-threonine aldolase (LTA) catalyzes C-C bond synthesis with moderate diastereoselectivity. In this study, with LTA from Cellulosilyticum sp (CpLTA) as an object, a mutability landscape was first constructed by performing saturation mutagenesis at substrate access tunnel amino acids. The combinatorial active-site saturation test/iterative saturation mutation (CAST/ISM) strategy was then used to tune diastereoselectivity. As a result, the diastereoselectivity of mutant H305L/Y8H/V143R was improved from 37.2 %syn to 99.4 %syn . Furthermore, the diastereoselectivity of mutant H305Y/Y8I/W307E was inverted to 97.2 %anti . Based on insight provided by molecular dynamics simulations and coevolution analysis, the Prelog rule was employed to illustrate the diastereoselectivity regulation mechanism of LTA, holding that the asymmetric formation of the C-C bond was caused by electrons attacking the carbonyl carbon atom of the substrate aldehyde from the re or si face. The study would be useful to expand LTA applications and guide engineering of other C-C bond-forming enzymes.

Effect of naturally-occurring mutations on the stability and function of cancer-associated NQO1: Comparison of experiments and computation.

Recent advances in DNA sequencing technologies are revealing a large individual variability of the human genome. Our capacity to establish genotype-phenotype correlations in such large-scale is, however, limited. This task is particularly challenging due to the multifunctional nature of many proteins. Here we describe an extensive analysis of the stability and function of naturally-occurring variants (found in the COSMIC and gnomAD databases) of the cancer-associated human NAD(P)H:quinone oxidoreductase 1 (NQO1). First, we performed in silico saturation mutagenesis studies (>5,000 substitutions) aimed to identify regions in NQO1 important for stability and function. We then experimentally characterized twenty-two naturally-occurring variants in terms of protein levels during bacterial expression, solubility, thermal stability, and coenzyme binding. These studies showed a good overall correlation between experimental analysis and computational predictions; also the magnitude of the effects of the substitutions are similarly distributed in variants from the COSMIC and gnomAD databases. Outliers in these experimental-computational genotype-phenotype correlations remain, and we discuss these on the grounds and limitations of our approaches. Our work represents a further step to characterize the mutational landscape of NQO1 in the human genome and may help to improve high-throughput in silico tools for genotype-phenotype correlations in this multifunctional protein associated with disease.

The XPO7/Npat Axis Is a Potential Therapeutic Target for TP53-Mutated AML

Acute myeloid leukemia (AML) is a devastating disease with a long-term survival rate below 30%. While AML is a genetically heterogeneous disease, TP53 mutations are among the most powerful risk factors in AML, which underscores the critical need to devise a novel therapeutic strategy for TP53-mutated AML. To identify genes/pathways whose loss is vulnerable to TP53 deficiency in AML cells, we performed genome-wide CRISPR-Cas9 screens using Trp53-knockout (KO) and wild-type (WT) mouse AML cells. To generate AML lines with a relatively "clean" genetic background, we established mouse AML lines harboring WT Trp53 with normal karyotype (Yamauchi et al. Cancer Cell, 2018). We then generated Trp53-KO AML lines using a single guide RNA (sgRNA) targeting Trp53. Genome-wide CRISPR-Cas9 dropout screens were performed using these lines to identify genes/pathways whose loss is vulnerable to TP53 deficiency. We identified Xpo7, a putative nuclear/cytoplasmic transporter, as a factor necessary for the survival of Trp53-KO AML cells. sgRNAs targeting Xpo7 were enriched in Trp53-WT AML cells after the 16-day culture, while they were significantly depleted in Trp53-KO cells. Trp53-KO cells were vulnerable to Xpo7 depletion, while Xpo7 functioned as a Trp53-dependent tumor suppressor in Trp53-WT AML cells. As expected, CRISPR/Cas9-mediated Xpo7 depletion and shRNA-mediated Xpo7 knockdown significantly delayed cell cycle progression and suppressed cell growth only in Trp53-KO AML cells, but not in WT cells. We next performed CRISPR-Cas9 saturation mutagenesis scan targeting all Xpo7 exons using Trp53-WT and -KO mouse AML cells. Strikingly, sgRNAs targeting Xpo7 coding regions were significantly enriched only in the Trp53-WT cells, namely in the presence of an MDM2 inhibitor. In stark contrast, the same sgRNAs which target coding exons including the C-terminal importin-β domain were mostly depleted in Trp53-KO AMLs. Since Xpo7 presumably mediates the nuclear import and export of proteins, we hypothesized that Xpo7 promotes retention of Trp53 protein in the nucleus in Trp53-WT AML cells. As expected, Trp53 protein levels in the nucleus and mRNA levels of Trp53-targeting genes were significantly decreased upon Xpo7 depletion. Furthermore, Xpo7 depletion resulted in resistance to Selinexor, a XPO1 inhibitor that reportedly blocks Trp53 protein export from the nucleus. These data suggest that Xpo7 retains WT-Trp53 in the nucleus and functions as a Trp53-dependent tumor suppressor in Trp53-WT AML. Next, to elucidate the mechanisms underlying the toxic effects observed in Trp53-KO AMLs upon Xpo7 depletion, we set out to identify proteins whose cytoplasmic/nuclear localization was potentially regulated by Xpo7 in Trp53-KO AML cells. We first performed IP-mass spectrometry using an anti-Xpo7 antibody to identify Xpo7 interacting proteins. Then, we obtained both cytoplasmic and nuclear protein extracts from both Xpo7-WT and -KO AML cells and identified proteins whose localization was significantly altered upon Xpo7 deletion by mass spectrometry. We found that Npat, an activator of histone transcription in the G1/S transition, interacted with Xpo7 and localized more in the cytoplasm upon Xpo7 depletion in Trp53-KO AMLs to a significant extent. Strikingly, Npat protein levels were robustly upregulated following Trp53 depletion, and Npat depletion significantly delayed the progression of the G1/S transition and proliferation in Trp53-KO AML cells. Finally, to explore the functional significance of the XPO7/NPAT axis in human TP53-mutated AML cells, we assessed the correlation between XPO7 expression levels and AML subtypes and/or genetic background using publicly-available datasets. XPO7 and NPAT mRNA levels were significantly upregulated in TP53-mutated AMLs in the TCGA datasets (Ley et al. NEJM, 2013). Moreover, XPO7 and NPAT mRNAs levels were remarkably high in acute erythroid leukemia (AEL) cases, where TP53 mutations are frequently observed (Tyner et al. Nature, 2018 and Iacobucci et al. Nat Genet. 2019). In fact, shRNA-mediated XPO7 or NPAT knockdown significantly suppressed proliferation of HEL, a human AEL cell line harboring TP53 mutation, in vitro and in vivo. In summary, we identified a synthetic lethal relationship between TP53 and XPO7. Our data suggest that XPO7/NPAT inactivation is a therapeutic vulnerability for TP53-mutated AML, such as AEL.