Theoretical Studies on the Mechanism of Thioesterase-Catalyzed Macrocyclization in Erythromycin Biosynthesis

Chemoenzymatic synthesis of macrocyclic peptides and polyketides via thioesterase-catalyzed macrocyclization

Chemoenzymatic strategies that combine synthetic and enzymatic transformations offer efficient approaches to yield target molecules, which have been increasingly employed in the synthesis of bioactive natural products. In the biosynthesis of macrocyclic nonribosomal peptides, polyketides, and their hybrids, thioesterase (TE) domains play a significant role in late-stage macrocyclization. These domains can accept mimics of native substrates in vitro and exhibit potential for use in total synthesis. This review summarizes the recent advances of TE domains in the chemoenzymatic synthesis for these natural products that aim to address the common issues in classical synthetic approaches and increase synthetic efficiencies, which have the potential to facilitate further pharmaceutical research.

Thioesterases as tools for chemoenzymatic synthesis of macrolactones

Macrocycles are a key functional group that can impart unique properties into molecules. Their synthesis has led to the development of many outstanding chemical methodologies and yet still remains challenging. Thioesterase (TE) domains are frequently responsible for macrocyclization in natural product biosynthesis and provide unique strengths for the enzymatic synthesis of macrocycles. In this feature article, we describe our work to characterize the substrate selectivity of TEs and to use these enzymes as biocatalysts. Our efforts have shown that the linear thioester activated substrates are loaded on TEs with limited substrate selectivity to generate acyl-enzyme intermediates. We show that cyclization of the acyl-enzyme intermediates can be highly selective, with competing hydrolysis of the acyl-enzyme intermediates. The mechanisms controlling TE-mediated macrocyclization versus hydrolysis are a significant unsolved problem in TE biochemistry. The potential of TEs as biocatalysts was demonstrated by using them in the chemoenzymatic total synthesis of macrocyclic depsipeptide natural products. This article highlights the strengths and potential of TEs as biocatalysts as well as their limitations, opening exciting research opportunities including TE engineering to optimize these powerful biocatalysts.

Exploring the substrate stereoselectivity and catalytic mechanism of nonribosomal peptide macrocyclization in surugamides biosynthesis

SurE, the first reported penicillin-binding protein-like thioesterase (PBP-like TE), is known as a new off-loading cyclase, which catalyzes heterochiral coupling in nonribosomal peptides (NRPs). However, the structural rationale for substrate stereoselectivity and enzymatic mechanism remains mysterious. Here, computational models, integrating MD simulations and QM/MM methods, unveiled SurE's substrate recognition and catalytic process. An oxyanion hole stabilized the C-terminal D-residue during recognition. Residue R446 anchored the substrate for macrocyclization. A vital hydrogen-bonding network (Y154, K66, N156), verified by mutation results, was responsible for the recognition of N-terminal L-residue and involvement in catalytic process with a calculated 19.4 kcal/mol energy barrier. Four novel-designed peptide precursors were effectively cyclized into cyclopeptides by SurE based on computational analysis. Our results provide a comprehensive understanding of SurE's catalytic mechanism and guiding design of versatile PBP-like TEs for novel macrocyclic NRPs.

Conformational Dynamics in Proteins: Entangled Slow Fluctuations and Nonequilibrium Reaction Events

Proteins exhibit conformational fluctuations and changes over various time scales, ranging from rapid picosecond-scale local atomic motions to slower microsecond-scale global conformational transformations. In the presence of these intricate fluctuations, chemical reactions occur and functions emerge. These conformational fluctuations of proteins are not merely stochastic random motions but possess distinct spatiotemporal characteristics. Moreover, chemical reactions do not always proceed along a single reaction coordinate in a quasi-equilibrium manner. Therefore, it is essential to understand spatiotemporal conformational fluctuations of proteins and the conformational change processes associated with reactions. In this Perspective, we shed light on the complex dynamics of proteins and their role in enzyme catalysis by presenting recent results regarding dynamic couplings and disorder in the conformational dynamics of proteins and rare but rapid enzymatic reaction events obtained from molecular dynamics simulations.

Semi-rational engineering an aldo-keto reductase for stereocomplementary reduction of α-keto amide compounds

Enantio-pure α-hydroxy amides are valuable intermediates for the synthesis of chiral pharmaceuticals. The asymmetric reduction of α-keto amides to generate chiral α-hydroxy amides is a difficult and challenging task in biocatalysis. In this study, iolS, an aldo-keto reductase from Bacillus subtilis 168 was exhibited as a potential biocatalyst, which could catalyze the reduction of diaryl α-keto amide such as 2-oxo-N, 2-diphenyl-acetamide (ONDPA) with moderate S-selectivity (76.1%, ee) and 60.5% conversion. Through semi-rational engineering, two stereocomplementary variants (I57F/F126L and N21A/F126A) were obtained with ee value of 97.6% (S) and 99.9% (R) toward ONDPA (1a), respectively, delivering chiral α-hydroxy amide with > 98% conversions. Moreover, the excellent S- and R-preference variants displayed improved stereoselectivities toward the other α-keto amide compounds. Molecular dynamic and docking analysis revealed that the two key residues at 21 and 126 were identified as the "switch", which specifically controlled the stereopreference of iolS by regulating the shape of substrate binding pocket as well as the substrate orientation. Our results offer an effective strategy to obtain α-hydroxy amides with high optical purity and provide structural insights into altering the stereoselectivity of AKRs.

Molecular Insights into Bifunctional Ambruticin DH3 for Substrate Specificity and Catalytic Mechanism

Dehydratase (DH), a domain located at polyketide synthase (PKS) modules, commonly catalyzes the dehydration of β-hydroxy to an α,β-unsaturated acyl intermediate. As a unique bifunctional dehydratase, AmbDH3 (the DH domain of module 3 of the ambruticin PKS) is verified to be responsible for both dehydration and the following pyran-forming cyclization. Besides, in vitro studies showed that its catalytic efficiency varies with different chiral substrates. However, the detailed molecular mechanism of AmbDH3 remains unclear. In this work, the structural rationale for the substrate specificity (2R/2S- and 6R/6S-substrates) in AmbDH3 was elucidated and the complete reaction pathways including dehydration and cyclization were presented. Both MD simulations and binding free energy calculations indicated AmbDH3 had a stronger preference for 2R-substrates (2R6R-2, 2R6S-3) than 2S-substrates (2S6R-1), and residue H51 and G61 around the catalytic pocket were emphasized by forming stable hydrogen bonds with 2R-substrates. In addition, AmbDH3's mild tolerance at C6 was explained by comparison of substrate conformation and hydrogen bond network in 6S- and 6R-substrate systems. The QM/MM results supported a consecutive one-base dehydration and cyclization mechanism for 2R6S-3 substrate with the energy barrier of 25.2 kcal mol^-1 and 24.5 kcal mol^-1 , respectively. Our computational results uncover the substrate recognition and catalytic process of the first bifunctional dehydratase-cyclase AmbDH3, which will shed light on the application of multifunctional DH domains in PKSs for diverse natural product analogs and benefit the chemoenzymatic synthesis of stereoselective pyran-containing products.

Theoretical Studies of Mutual Effects between 6-m-r Hemiketalization and 26-m-r Lactonization in Pimaricin Thioesterase

Pimaricin is a small polyene macrolide antibiotic and has been broadly used as an antimycotic and antiprotozoal agent in both humans and foods. As a thioesterase in type-I polyketide synthase, pimTE controls the 26-m-r macrolide main chain release in pimaricin biosynthesis. In this work, we sought to determine whether the 6-m-r hemiketal formation was linked to pimTE-catalyzed 26-m-r lactonization. Compared to non-hemiketal TEs, pimTE is characterized by an aspartic acid residue (D179) accessible to the U-turn motif in the acyl-enzyme intermediate. Both the covalent docking and molecular dynamics simulations demonstrate that the reactive conformations for macrocyclic lactonization are drastically promoted by the 6-m-r hemiketal. Moreover, the small-model quantum mechanistic calculations suggest that protic residues can significantly accelerate the 6-m-r hemiketal cyclization. In addition, the post-hemiketal molecular dynamic simulations demonstrate that hydrogen-bonding networks surrounding the substrate U-turn of the hairpin-shaped conformation changes significantly when the 6-m-r hemiketal is formed. In particular, the R-hemiketal intermediate is not only catalyzed by the D179 residue, but also twists the hairpin structure to the 26-m-r lactonizing pre-reaction state. By contrast, the S-hemiketal formation is unlikely catalyzed by D179, which twists the hairpin in an opposite direction. Our results propose that pimTE could be a bi-functional enzyme, which can synergistically catalyze tandem 6-m-r and 26-m-r formations during the main-chain release of pimaricin biosynthesis.

Thioesterase enzyme families: Functions, structures, and mechanisms

Thioesterases are enzymes that hydrolyze thioester bonds in numerous biochemical pathways, for example in fatty acid synthesis. This work reports known functions, structures, and mechanisms of updated thioesterase enzyme families, which are classified into 35 families based on sequence similarity. Each thioesterase family is based on at least one experimentally characterized enzyme, and most families have enzymes that have been crystallized and their tertiary structure resolved. Classifying thioesterases into families allows to predict tertiary structures and infer catalytic residues and mechanisms of all sequences in a family, which is particularly useful because the majority of known protein sequence have no experimental characterization. Phylogenetic analysis of experimentally characterized thioesterases that have structures with the two main structural folds reveal convergent and divergent evolution. Based on tertiary structure superimposition, catalytic residues are predicted.

Exploring the Molecular Basis of Substrate and Product Selectivities of Nocardicin Bifunctional Thioesterase

D-amino acid introduction in peptides can enrich their biological activities and pharmacological properties as potential drugs. This achievement of stereochemical inversion usually owes to an epimerase or racemase. Interestingly, a unique bifunctional thioesterase (NocTE), which is incorporated in nonribosomal peptide synthetase (NRPS) NocA-NocB assembly line for the biosynthesis of monocyclic β-lactam antibiotic nocardicin A, can control the generation of D-products with high stereochemical purity. However, the molecular basis of NocTE selectivity on substrates and products is still unclear. Herein, we constructed a series of systems with different peptides varying in stereochemistry, length, and composition to investigate the substrate selectivity. The studies on binding affinities and loading conformations elucidated the important roles of peptide length and β-lactam ring in substrate selectivity. Through energy decomposition and interaction analyses, some key residues involved in substrate selectivity were captured. On the other hand, natural product undergoing epimerization was found to be liberated from the active pocket more easily in comparison with its diastereomer (epi-nocardicin G), explaining the superiority of nocardicin G. These results provide detailed molecular insights into the exquisite control of substrate and product scopes for NocTE, and encourage to diversification of substrates and final products for NRPS assembly line. The molecular insights into substrate and product selectivities of unique bifunctional thioesterase NocTE were illustrated via several molecular simulations and free energy calculations, contributing to expanding substrate and product scopes of nonribosomal peptide synthetases.

Insights into specificity and catalytic mechanism of amphotericin B/nystatin thioesterase

Polyene polyketides amphotericin B (AMB) and nystatin (NYS) are important antifungal drugs. Thioesterases (TEs), located at the last module of PKS, control the release of polyketides by cyclization or hydrolysis. Intrigued by the tiny structural difference between AMB and NYS, as well as the high sequence identity between AMB TE and NYS TE, we constructed four systems to study the structural characteristics, catalytic mechanism, and product release of AMB TE and NYS TE with combined MD simulations and quantum mechanics/molecular mechanics calculations. The results indicated that compared with AMB TE, NYS TE shows higher specificity on its natural substrate and R26 as well as D186 were proposed to a key role in substrate recognition. The energy barrier of macrocyclization in AMB-TE-Amb and AMB-TE-Nys systems were calculated to be 14.0 and 22.7 kcal/mol, while in NYS-TE-Nys and NYS-TE-Amb systems, their energy barriers were 17.5 and 25.7 kcal/mol, suggesting the cyclization with their natural substrates were more favorable than that with exchanged substrates. At last, the binding free energy obtained with the MM-PBSA.py program suggested that it was easier for natural products to leave TE enzymes after cyclization. And key residues to the departure of polyketide product from the active site were highlighted. We provided a catalytic overview of AMB TE and NYS TE including substrate recognition, catalytic mechanism and product release. These will improve the comprehension of polyene polyketide TEs and benefit for broadening the substrate flexibility of polyketide TEs.

Theoretical study on substrate recognition and catalytic mechanisms of gephyronic acid dehydratase DH1

The bifunctional dehydratase GphF DH1 catalyzes both the dehydration of β-hydroxy and the double bond isomerization with the energy barrier of 27.0 kcal mol⁻¹ and 17.2 kcal mol⁻¹ respectively.

Intrinsic and Extrinsic Programming of Product Chain Length and Release Mode in Fungal Collaborating Iterative Polyketide Synthases

Combinatorial biosynthesis with fungal polyketide synthases (PKSs) promises to produce unprecedented bioactive "unnatural" natural products (uNPs) for drug discovery. Genome mining of the dothideomycete Rhytidhysteron rufulum uncovered a collaborating highly reducing PKS (hrPKS)-nonreducing PKS (nrPKS) pair. These enzymes produce trace amounts of rare S-type benzenediol macrolactone congeners with a phenylacetate core in a heterologous host. However, subunit shuffling and domain swaps with voucher enzymes demonstrated that all PKS domains are highly productive. This contradiction led us to reveal novel programming layers exerted by the starter unit acyltransferase (SAT) and the thioesterase (TE) domains on the PKS system. First, macrocyclic vs linear product formation is dictated by the intrinsic biosynthetic program of the TE domain. Next, the chain length of the hrPKS product is strongly influenced in trans by the off-loading preferences of the nrPKS SAT domain. Last, TE domains are size-selective filters that facilitate or obstruct product formation from certain priming units. Thus, the intrinsic programs of the SAT and TE domains are both part of the extrinsic program of the hrPKS subunit and modulate the observable metaprogram of the whole PKS system. Reconstruction of SAT and TE phylogenies suggests that these domains travel different evolutionary trajectories, with the resulting divergence creating potential conflicts in the PKS metaprogram. Such conflicts often emerge in chimeric PKSs created by combinatorial biosynthesis, reducing biosynthetic efficiency or even incapacitating the system. Understanding the points of failure for such engineered biocatalysts is pivotal to advance the biosynthetic production of uNPs.

In situ and real-time insight into Rhizopus chinensis lipase under high pressure and temperature: Conformational traits and biobehavioural analysis

An in situ and real-time investigation was performed using an optical cell system and in-silico analysis to reveal the impacts of pressure and temperature on the conformational state and behaviours of Rhizopus chinensis lipase (RCL). The fluorescence intensity (FI) of RCL increased remarkably under high pressure, and part of this increase was recovered after depressurization. This result displayed the partially reversible conformational change of RCL, which may be associated with the local change of Trp224 near the catalytic centre. High temperature caused a significant loss of secondary structure, whereas the α-helical segments including the lid were preserved by high pressure even at temperatures over 60 °C. The parameters of enzymatic reaction monitored by UV showed that the hydrolysis rate was remarkably enhanced by the pressure of 200 MPa. In the pressure range of 0.1-200 MPa, the active volume measured by the in situ system decreased from -2.85 to -6.73 mL/mol with the temperature increasing from 20 °C to 40 °C. The high catalytic capacity of the lipase under high pressure and high temperature was primarily attributed to pressure protection on RCL.

Structure and Function of BorB, the Type II Thioesterase from the Borrelidin Biosynthetic Gene Cluster

α/β hydrolases make up a large and diverse protein superfamily. In natural product biosynthesis, cis-acting thioesterase α/β hydrolases can terminate biosynthetic assembly lines and release products by hydrolyzing or cyclizing the biosynthetic intermediate. Thioesterases can also act in trans, removing aberrant intermediates and restarting stalled biosynthesis. Knockout of this "editing" function leads to reduced product titers. The borrelidin biosynthetic gene cluster from Streptomyces parvulus Tü4055 contains a hitherto uncharacterized stand-alone thioesterase, borB. In this work, we demonstrate that purified BorB cleaves acyl substrates with a preference for propionate, which supports the hypothesis that it is also an editing thioesterase. The crystal structure of BorB shows a wedgelike hydrophobic substrate binding crevice that limits substrate length. To investigate the structure-function relationship, we made chimeric BorB variants using loop regions from characterized homologues with different specificities. BorB chimeras slightly reduced activity, arguing that the modified region is a not major determinant of substrate preference. The structure-function relationships described here contribute to the process of elimination for understanding thioesterase specificity and, ultimately, engineering and applying trans-acting thioesterases in biosynthetic assembly lines.

Comprehensive Insights into the Catalytic Mechanism of Middle East Respiratory Syndrome 3C-Like Protease and Severe Acute Respiratory Syndrome 3C-Like Protease

Coronavirus 3C-like protease (3CLPro) is a highly conserved cysteine protease employing a catalytic dyad for its functions. 3CLPro is essential to the viral life cycle and, therefore, is an attractive target for developing antiviral agents. However, the detailed catalytic mechanism of coronavirus 3CLPro remains largely unknown. We took an integrated approach of employing X-ray crystallography, mutational studies, enzyme kinetics study, and inhibitors to gain insights into the mechanism. Such experimental work is supplemented by computational studies, including the prereaction state analysis, the ab initio calculation of the critical catalytic step, and the molecular dynamic simulation of the wild-type and mutant enzymes. Taken together, such studies allowed us to identify a residue pair (Glu-His) and a conserved His as critical for binding; a conserved GSCGS motif as important for the start of catalysis, a partial negative charge cluster (PNCC) formed by Arg-Tyr-Asp as essential for catalysis, and a conserved water molecule mediating the remote interaction between PNCC and catalytic dyad. The data collected and our insights into the detailed mechanism have allowed us to achieve a good understanding of the difference in catalytic efficiency between 3CLPro from SARS and MERS, conduct mutational studies to improve the catalytic activity by 8-fold, optimize existing inhibitors to improve the potency by 4-fold, and identify a potential allosteric site for inhibitor design. All such results reinforce each other to support the overall catalytic mechanism proposed herein.

Evolutionary coupling saturation mutagenesis: Coevolution-guided identification of distant sites influencing Bacillus naganoensis pullulanase activity

Pullulanases are well-known debranching enzymes hydrolyzing α-1,6-glycosidic linkages. To date, engineering of pullulanase is mainly focused on catalytic pocket or domain tailoring based on structure/sequence information. Saturation mutagenesis-involved directed evolution is, however, limited by the low number of mutational sites compatible with combinatorial libraries of feasible size. Using Bacillus naganoensis pullulanase as a target protein, here we introduce the 'evolutionary coupling saturation mutagenesis' (ECSM) approach: residue pair covariances are calculated to identify residues for saturation mutagenesis, focusing directed evolution on residue pairs playing important roles in natural evolution. Evolutionary coupling (EC) analysis identified seven residue pairs as evolutionary mutational hotspots. Subsequent saturation mutagenesis yielded variants with enhanced catalytic activity. The functional pairs apparently represent distant sites affecting enzyme activity.

Insight into Structural Characteristics of Protein-Substrate Interaction in Pimaricin Thioesterase

As a polyene antibiotic of great pharmaceutical significance, pimaricin has been extensively studied to enhance its productivity and effectiveness. In our previous studies, pre-reaction state (PRS) has been validated as one of the significant conformational categories before macrocyclization, and is critical to mutual recognition and catalytic preparation in thioesterase (TE)-catalyzed systems. In our study, molecular dynamics (MD) simulations were conducted on pimaricin TE-polyketide complex and PRS, as well as pre-organization state (POS), a molecular conformation possessing a pivotal intra-molecular hydrogen bond, were detected. Conformational transition between POS and PRS was observed in one of the simulations, and POS was calculated to be energetically more stable than PRS by 4.58 kcal/mol. The structural characteristics of PRS and POS-based hydrogen-bonding, and hydrophobic interactions were uncovered, and additional simulations were carried out to rationalize the functions of several key residues (Q29, M210, and R186). Binding energies, obtained from MM/PBSA calculations, were further decomposed to residues, in order to reveal their roles in product release. Our study advanced a comprehensive understanding of pimaricin TE-catalyzed macrocyclization from the perspectives of conformational change, protein-polyketide recognition, and product release, and provided potential residues for rational modification of pimaricin TE.

How Chorismatases Regulate Distinct Reaction Channels in a Single Conserved Active Pocket: Mechanistic Analysis with QM/MM (ONIOM) Investigations

The FkbO and Hyg5 subfamilies of chorismatases share the same active-site architectures, but perform distinct reaction mechanisms, that is, FkbO employs a hydrolysis reaction whereas Hyg5 proceeds through an intramolecular mechanism. Despite extensive research efforts, the detailed mechanism of the product selectivity in chorismatases need to be further unmasked. In this study, the effects of the A/G residue group (A244^FkbO /G240^Hyg5 ) and the V/Q residue group (V209^FkbO /Q201^Hyg5 ) on the catalytic mechanisms are investigated by employing molecular dynamics simulations and hybrid quantum mechanical/molecular mechanical (QM/MM) calculations of the two wild-type models (FkbO/CHO and Hyg5/CHO; CHO=chorismate) and four mutants models (A244G-FkbO/CHO and G240A-Hyg5/CHO; V209Q-FkbO/CHO and Q201V-Hyg5/CHO). Our results showed that the A/G residue group mentioned by previous works would cause changes in the binding states of the substrate and the orientation of the catalytic glutamate, but only these changes affect the product selectivity in chorismatases limitedly. Interestingly, the distal V/Q residue group, which determines the internal water self-regulating ability at the active site, has significant impact on the selectivity of the catalytic mechanisms. The V/Q residue group is suggested to be an important factor to control the catalytic activities in chorismatases. The results are consistent with biochemical and structural experiments, providing novel insight into the mechanism of product selectivity in chorismatases.

Understanding the R882H mutation effects of DNA methyltransferase DNMT3A: a combination of molecular dynamics simulations and QM/MM calculations

The AML-related high-frequent R882H mutation of DNA (cytosine-5)-methyltransferase 3A (DNMT3A), a key enzyme forde novoepigenetic methylation in human beings, was characterized by a disturbing conformation ofS-adenosylmethionine (SAM).

Why does tautomycetin thioesterase prefer hydrolysis to macrocyclization? Theoretical study on its catalytic mechanism

In this study, we attempted to uncover the reasons why Tautomycetin thioesterase (TMC TE) prefers hydrolysis rather than macrocyclization, and reveal the molecular basis of TE-catalyzed hydrolysis and macrocyclization.

Structural Insight into Enantioselective Inversion of an Alcohol Dehydrogenase Reveals a "Polar Gate" in Stereorecognition of Diaryl Ketones

Diaryl ketones are important building blocks for synthesizing pharmaceuticals and are generally regarded as "difficult-to-reduce" ketones due to the large steric hindrance of their two bulky aromatic side chains. Alcohol dehydrogenase from Kluyveromyces polyspora ( KpADH) has been identified as a robust biocatalyst due to its high conversion of diaryl ketone substrate (4-chlorophenyl)(pyridine-2-yl)ketone (CPMK) with a moderate R-selectivity of 82% ee. To modulate the stereoselectivity of KpADH, a "polarity scanning" strategy was proposed, in which six key residues inside and at the entrance of the substrate binding pocket were identified. After iterative combinatorial mutagenesis, variants Mu-R2 and Mu-S5 with enhanced (99.2% ee, R) and inverted (97.8% ee, S) stereoselectivity were obtained. The crystal structures of KpADH and two mutants in complex with NADPH were resolved to elucidate the evolution of enantioselective inversion. Based on MD simulation, Mu-R2-CPMK_ProR and Mu-S5-CPMK_ProS were more favorable in the formation of prereaction states. Interestingly, a quadrilateral plane formed by α-carbons of four residues (N136, V161, C237, and G214) was identified at the entrance of the substrate binding pocket of Mu-S5; this plane acts as a "polar gate" for substrates. Due to the discrepancy in charge characteristics between chlorophenyl and pyridine substituents, the pro- S orientation of CPMK is defined when it passes through the "polar gate" in Mu-S5, whereas the similar plane in wild-type is blocked by several aromatic residues. Our result paves the way for engineering stereocomplementary ADH toward bulky diaryl ketones and provides structural insight into the mechanism of stereoselective inversion.

A Single Active Site Mutation in the Pikromycin Thioesterase Generates a More Effective Macrocyclization Catalyst

Macrolactonization of natural product analogs presents a significant challenge to both biosynthetic assembly and synthetic chemistry. In the preceding paper , we identified a thioesterase (TE) domain catalytic bottleneck processing unnatural substrates in the pikromycin (Pik) system, preventing the formation of epimerized macrolactones. Here, we perform molecular dynamics simulations showing the epimerized hexaketide was accommodated within the Pik TE active site; however, intrinsic conformational preferences of the substrate resulted in predominately unproductive conformations, in agreement with the observed hydrolysis. Accordingly, we engineered the stereoselective Pik TE to yield a variant (TE_S148C) with improved reaction kinetics and gain-of-function processing of an unnatural, epimerized hexaketide. Quantum mechanical comparison of model TE_S148C and TE_WT reaction coordinate diagrams revealed a change in mechanism from a stepwise addition-elimination (TE_WT) to a lower energy concerted acyl substitution (TE_S148C), accounting for the gain-of-function and improved reaction kinetics. Finally, we introduced the S148C mutation into a polyketide synthase module (PikAIII-TE) to impart increased substrate flexibility, enabling the production of diastereomeric macrolactones.

Towards Precision Engineering of Canonical Polyketide Synthase Domains: Recent Advances and Future Prospects

Modular polyketide synthases (mPKSs) build functionalized polymeric chains, some of which have become blockbuster therapeutics. Organized into repeating clusters (modules) of independently-folding domains, these assembly-line-like megasynthases can be engineered by introducing non-native components. However, poor introduction points and incompatible domain combinations can cause both unintended products and dramatically reduced activity. This limits the engineering and combinatorial potential of mPKSs, precluding access to further potential therapeutics. Different regions on a given mPKS domain determine how it interacts both with its substrate and with other domains. Within the assembly line, these interactions are crucial to the proper ordering of reactions and efficient polyketide construction. Achieving control over these domain functions, through precision engineering at key regions, would greatly expand our catalogue of accessible polyketide products. Canonical mPKS domains, given that they are among the most well-characterized, are excellent candidates for such fine-tuning. The current minireview summarizes recent advances in the mechanistic understanding and subsequent precision engineering of canonical mPKS domains, focusing largely on developments in the past year.