Computational design of carboxylase for the synthesis of 4-hydroxyisophthalic acid from p-hydroxybenzoic acid by fixing CO2

Carbon dioxide (CO2) emissions constitute the primary contribution to global climate change. Synthetic CO2 fixation represents an exceptionally appealing and sustainable method for carbon neutralization. Unlike the limitations of chemical catalysis, biological CO2 fixation displays high selectivity and the ability to operate under mild conditions. The superfamily of amidohydrolases has demonstrated the ability to synthesize a range of aromatic monocarboxylic acids. However, there is a scarcity of reported carboxylases capable of synthesizing aromatic dicarboxylic acids. Among these, 4-hydroxyisophthalic acid holds significant potential for applications across various fields, yet no enzyme has been reported for its synthesis. In this study, we developed for the first time that exhibits starting activity in fixing CO2 to synthesize 4-hydroxyisophthalic acid. Furthermore, we have devised a computational strategy that effectively enhances the catalytic activity of this enzyme. A focused library comprising only 13 variants was generated. Experimental validation confirmed a threefold improvement in the carboxylation activity of the optimal variant (L47M). The computational enzyme design strategy proposed in this paper demonstrates broad applicability in developing carboxylases for synthesizing other aromatic dicarboxylic acids. This lays the groundwork for leveraging biocatalysis in industrial synthesis for CO2 fixation.

Decoding CRISPR-Cas PAM recognition with UniDesign

The critical first step in Clustered Regularly Interspaced Short Palindromic Repeats (CRISPR)-associated (CRISPR-Cas) protein-mediated gene editing is recognizing a preferred protospacer adjacent motif (PAM) on target DNAs by the protein's PAM-interacting amino acids (PIAAs). Thus, accurate computational modeling of PAM recognition is useful in assisting CRISPR-Cas engineering to relax or tighten PAM requirements for subsequent applications. Here, we describe a universal computational protein design framework (UniDesign) for designing protein-nucleic acid interactions. As a proof of concept, we applied UniDesign to decode the PAM-PIAA interactions for eight Cas9 and two Cas12a proteins. We show that, given native PIAAs, the UniDesign-predicted PAMs are largely identical to the natural PAMs of all Cas proteins. In turn, given natural PAMs, the computationally redesigned PIAA residues largely recapitulated the native PIAAs (74% and 86% in terms of identity and similarity, respectively). These results demonstrate that UniDesign faithfully captures the mutual preference between natural PAMs and native PIAAs, suggesting it is a useful tool for engineering CRISPR-Cas and other nucleic acid-interacting proteins. UniDesign is open-sourced at https://github.com/tommyhuangthu/UniDesign.

Between Protein Fold and Nucleophile Identity: Multiscale Modeling of the TEV Protease Enzyme-Substrate Complex

The cysteine protease from the tobacco etch virus (TEVp) is a well-known and widely utilized enzyme. TEVp's chymotrypsin-like fold is generally associated with serine catalytic triads that differ in terms of a reaction mechanism from the most well-studied papain-like cysteine proteases. The question of what dominates the TEVp mechanism, nucleophile identity, or structural composition has never been previously addressed. Here, we use enhanced sampling multiscale modeling to uncover that TEVp combines the features of two worlds in such a way that potentially hampers its activity. We show that TEVp cysteine is strictly in the anionic form in a free enzyme similar to papain. Peptide binding shifts the equilibrium toward the nucleophile's protonated form, characteristic of chymotrypsin-like proteases, although the cysteinyl anion form is still present and interconversion is rapid. This way cysteine protonation generates enzyme states that are a diversion from the most effective course of action, with only 13.2% of Michaelis complex sub-states able to initiate the reaction. As a result, we propose an updated view on the reaction mechanism catalyzed by TEVp. We also demonstrate that AlphaFold is able to construct protease-substrate complexes with high accuracy. We propose that our findings open a way for its industrious use in enzymological tasks. Unique features of TEVp discovered in this work open a discussion on the evolutionary history and trade-offs of optimizing serine triad-associated folds to cysteine as a nucleophile.

Enhancing cellulosic ethanol production through coevolution of multiple enzymatic characteristics of β-glucosidase from Penicillium oxalicum 16

In previous studies, we isolated a novel β-glucosidase from Penicillium oxalicum 16 (16BGL), which is useful for producing cellulosic ethanol. However, 16BGL has a relatively low enzyme activity and product tolerance; besides, a huge gap exists between the optimum temperature of 16BGL (70 °C) and the fermentation temperature for producing cellulosic ethanol (40 °C). Here, we present a directed evolution-based study, which combines one-round error-prone PCR with three rounds of high-throughput screening, for coevolving multiple enzymatic characteristics of 16BGL. We identified an improved variant Y-1-B1 with a triple mutant (G414S/D421V/T441S). Y-1-B1 had an optimum temperature of 50 °C, much closer to the fermentation temperature. The catalytic efficiency of Y-1-B1 for hydrolyzing pNPG was 1355 mM^-1 s^-1 at 50 °C and pH 5, significantly higher than that of 16BGL (807 mM^-1 s^-1). Y-1-B1 also achieved a slightly reduced strength of product inhibition of 1.1 at a glucose concentration of 20 mM, compared with the ratio of 1.3 for 16BGL. A maximum titer of 6.9 g/L for ethanol production was achieved in the reaction with Y-1-B1, which was 22% higher than that achieved with 16BGL. Structure modeling revealed that the mutations are distant from the active-site pocket. Therefore, we performed molecular dynamics (MD) simulations to understand why these mutations can improve catalytic efficiency. MD simulation revealed that the nucleophilic residue Asp261 had a much closer contact with the glucosidic center of pNPG in the simulation with Y-1-B1 than that with 16BGL, suggesting that the mutant is more favorable for catalysis. KEY POINTS: • Multiple enzymatic properties of Penicillium oxalicum 16 BGL were coevolved. • A catalytically efficient triple mutant G414S/D421V/T441S was reported. • Molecular dynamics simulation supported the enhanced catalytic activity.

Computer-aided engineering of adipyl-CoA synthetase for enhancing adipic acid synthesis

OBJECTIVE

To enhance adipic acid production, a computer-aided approach was employed to engineer the adipyl-CoA synthetase from Thermobifida fusca by combining sequence analysis, protein structure modeling, in silico site-directed mutagenesis, and molecular dynamics simulation.

RESULTS

Two single mutants of T. fusca adipyl-CoA synthetase, E210βN and E210βQ, achieved a specific enzyme activity of 1.95 and 1.84 U/mg, respectively, which compared favorably with the 1.48 U/mg for the wild-type. The laboratory-level fermentation experiments showed that E210βN and E210βQ achieved a maximum adipic acid titer of 0.32 and 0.3 g/L. In contrast, the wild-type enzyme yielded a titer of 0.15 g/L under the same conditions. Molecular dynamics (MD) simulations revealed that the mutants (E210βN and E210βQ) could accelerate the dephosphorylation process in catalysis and enhance enzyme activity.

CONCLUSIONS

The combined computational-experimental approach provides an effective strategy for enhancing enzymatic characteristics, and the mutants may find a useful application for producing adipic acid.

Co-evolution of β-glucosidase activity and product tolerance for increasing cellulosic ethanol yield

β-Glucosidase (BGL) plays a key role in cellulose hydrolysis. However, it is still a great challenge to enhance product tolerance and enzyme activity of BGL simultaneously. Here, we utilized one round error-prone PCR to engineer the Penicillium oxalicum 16 BGL (16BGL) for improving the cellulosic ethanol yield. We identified a new variant (L-6C), a triple mutant (M280T/V484L/D589E), with enhanced catalytic efficiency ([Formula: see text]) for hydrolyzing pNPG and reduced strength of inhibition ([Formula: see text]) by glucose. To be specific, L-6C achieved a [Formula: see text] of 0.35 at a glucose concentration of 20 mM, which was 3.63 times lower than that attained by 16BGL. The catalytic efficiency for L-6C to hydrolyze pNPG was determined to be 983.68 mM^-1 s^-1, which was 22% higher than that for 16BGL. However, experiments showed that L-6C had reduced binding affinity (2.88 mM) to pNGP compared with 16BGL (1.69 mM). L-6C produced 6.15 g/L ethanol whose yield increased by about 10% than 16BGL. We performed molecular docking and molecular dynamics (MD) simulation, and binding free energy calculation using the Molecular Mechanics/Poisson Boltzmann surface area (MM/PBSA) method. MD simulation together with the MM/PBSA calculation suggested that L-6C had reduced binding free energy to pNPG, which was consistent with the experimental data.

Strategy for the Biosynthesis of Short Oligopeptides: Green and Sustainable Chemistry

Short oligopeptides are some of the most promising and functionally important amide bond-containing components, with widespread applications. Biosynthesis of these oligopeptides may potentially become the ultimate strategy because it has better cost efficiency and environmental-friendliness than conventional solid phase peptide synthesis and chemo-enzymatic synthesis. To successfully apply this strategy for the biosynthesis of structurally diverse amide bond-containing components, the identification and selection of specific biocatalysts is extremely important. Given that perspective, this review focuses on the current knowledge about the typical enzymes that might be potentially used for the synthesis of short oligopeptides. Moreover, novel enzymatic methods of producing desired peptides via metabolic engineering are highlighted. It is believed that this review will be helpful for technological innovation in the production of desired peptides.

Modelling of substrate access and substrate binding to cephalosporin acylases

Semisynthetic cephalosporins are widely used antibiotics currently produced by different chemical steps under harsh conditions, which results in a considerable amount of toxic waste. Biocatalytic synthesis by the cephalosporin acylase from Pseudomonas sp. strain N176 is a promising alternative. Despite intensive engineering of the enzyme, the catalytic activity is still too low for a commercially viable process. To identify the bottlenecks which limit the success of protein engineering efforts, a series of MD simulations was performed to study for two acylase variants (WT, M6) the access of the substrate cephalosporin C from the bulk to the active site and the stability of the enzyme-substrate complex. In both variants, cephalosporin C was binding to a non-productive substrate binding site (E86α, S369β, S460β) at the entrance to the binding pocket, preventing substrate access. A second non-productive binding site (G372β, W376β, L457β) was identified within the binding pocket, which competes with the active site for substrate binding. Noteworthy, substrate binding to the protein surface followed a Langmuir model resulting in binding constants K = 7.4 and 9.2 mM for WT and M6, respectively, which were similar to the experimentally determined Michaelis constants KM = 11.0 and 8.1 mM, respectively.