Search

Catalysis with engineered enzymes has provided more efficient routes for the production of active pharmaceutical agents. However, the potential of biocatalysis to assist in early-stage drug discovery campaigns remains largely untapped. In this study, we have developed a biocatalytic strategy for the construction of sp3-rich polycyclic compounds via the intramolecular cyclopropanation of benzothiophenes and related heterocycles. Two carbene transferases with complementary regioisomer selectivity were evolved to catalyse the stereoselective cyclization of benzothiophene substrates bearing diazo ester groups at the C2 or C3 position of the heterocycle. The detailed mechanisms of these reactions were elucidated by a combination of crystallographic and computational analyses. Leveraging these insights, the substrate scope of one of the biocatalysts could be expanded to include previously unreactive substrates, highlighting the value of integrating evolutionary and rational strategies to develop enzymes for new-to-nature transformations. The molecular scaffolds accessed here feature a combination of three-dimensional and stereochemical complexity with rule-of-three properties, which should make them highly valuable for fragment-based drug discovery campaigns.

Trifluoromethyl-substituted cyclopropanes (CF3 -CPAs) constitute an important class of compounds for drug discovery. While several methods have been developed for synthesis of trans-CF3 -CPAs, stereoselective production of corresponding cis-diastereomers remains a formidable challenge. We report a biocatalyst for diastereo- and enantio-selective synthesis of cis-CF3 -CPAs with activity on a variety of alkenes. We found that an engineered protoglobin from Aeropyrnum pernix (ApePgb) can catalyze this unusual reaction at preparative scale with low-to-excellent yield (6-55 %) and enantioselectivity (17-99 % ee), depending on the substrate. Computational studies revealed that the steric environment in the active site of the protoglobin forced iron-carbenoid and substrates to adopt a pro-cis near-attack conformation. This work demonstrates the capability of enzyme catalysts to tackle challenging chemistry problems and provides a powerful means to expand the structural diversity of CF3 -CPAs for drug discovery.

New-to-nature radical biocatalysis has recently emerged as a powerful strategy to tame fleeting open-shell intermediates for stereoselective transformations. In 2021, we introduced a novel metalloredox biocatalysis strategy that leverages the innate redox properties of the heme cofactor of P450 enzymes, furnishing new-to-nature atom-transfer radical cyclases (ATRCases) with excellent activity and stereoselectivity. Herein, we report a combined computational and experimental study to shed light on the mechanism and origins of enantioselectivity for this system. Molecular dynamics and quantum mechanics/molecular mechanics (QM/MM) calculations revealed an unexpected role of the key beneficial mutation I263Q. The glutamine residue serves as an essential hydrogen bond donor that engages with the carbonyl moiety of the substrate to promote bromine atom abstraction and enhance the enantioselectivity of radical cyclization. Therefore, the evolved ATRCase is a bifunctional biocatalyst, wherein the heme cofactor enables atom-transfer radical biocatalysis, while the hydrogen bond donor residue further enhances the activity and enantioselectivity. Unlike many enzymatic stereocontrol rationales based on a rigid substrate binding model, our computations demonstrate a high degree of rotational flexibility of the allyl moiety in an enzyme-substrate complex and succeeding intermediates. Therefore, the enantioselectivity is controlled by the radical cyclization transition states rather than the substrate orientation in ground-state complexes in the preceding steps. During radical cyclization, anchoring effects of the Q263 residue and steric interactions with the heme cofactor concurrently control the π-facial selectivity, allowing for highly enantioselective C-C bond formation. Our computational findings are corroborated by experiments with ATRCase mutants generated from site-directed mutagenesis.

Propargyl amines are versatile synthetic intermediates with numerous applications in the pharmaceutical industry. An attractive strategy for efficient preparation of these compounds is nitrene propargylic C(sp3)-H insertion. However, achieving this reaction with good chemo-, regio-, and enantioselective control has proven to be challenging. Here, we report an enzymatic platform for the enantioselective propargylic amination of alkynes using a hydroxylamine derivative as the nitrene precursor. Cytochrome P450 variant PA-G8 catalyzing this transformation was identified after eight rounds of directed evolution. A variety of 1-aryl-2-alkyl alkynes are accepted by PA-G8, including those bearing heteroaromatic rings. This biocatalytic process is efficient and selective (up to 2610 total turnover number (TTN) and 96% ee) and can be performed on preparative scale.

The mechanism by which molecular oxygen is activated by the organic cofactor pyridoxal phosphate (PLP) for oxidation reactions remains poorly understood. Recent work has identified arginine oxidases that catalyze desaturation or hydroxylation reactions. Here, we investigate a desaturase from the Pseudoalteromonas luteoviolacea indolmycin pathway. Our work, combining X-ray crystallographic, biochemical, spectroscopic, and computational studies, supports a shared mechanism with arginine hydroxylases, involving two rounds of single-electron transfer to oxygen and superoxide rebound at the 4' carbon of the PLP cofactor. The precise positioning of a water molecule in the active site is proposed to control the final reaction outcome. This proposed mechanism provides a unified framework to understand how oxygen can be activated by PLP-dependent enzymes for oxidation of arginine and elucidates a shared mechanistic pathway and intertwined evolutionary history for arginine desaturases and hydroxylases.

Iterative P450 enzymes are powerful biocatalysts for selective late-stage C-H oxidation of complex natural product scaffolds. These enzymes represent useful tools for selectivity and cascade reactions, facilitating direct access to core structure diversification. Recently, we reported the structure of the multifunctional bacterial P450 TamI and elucidated the molecular basis of its substrate binding and strict reaction sequence at distinct carbon atoms of the substrate. Here, we report the design and characterization of a toolbox of TamI biocatalysts, generated by mutations at Leu101, Leu244, and/or Leu295, that alter the native selectivity, step sequence, and number of reactions catalyzed, including the engineering of a variant capable of catalyzing a four-step oxidative cascade without the assistance of the flavoprotein and oxidative partner TamL. The tuned enzymes override inherent substrate reactivity, enabling catalyst-controlled C-H functionalization and alkene epoxidation of the tetramic acid-containing natural product tirandamycin. Five bioactive tirandamycin derivatives (6-10) were generated through TamI-mediated enzymatic synthesis. Quantum mechanics calculations and MD simulations provide important insights into the basis of altered selectivity and underlying biocatalytic mechanisms for enhanced continuous oxidation of the iterative P450 TamI.

A cytochrome c heme protein was recently engineered to catalyze the formation of carbon-silicon bonds via carbene insertion into Si-H bonds, a reaction that was not previously known to be catalyzed by a protein. High chemoselectivity toward C-Si bond formation over competing C-N bond formation was achieved, although this trait was not screened for during directed evolution. Using computational and experimental tools, we now establish that activity and chemoselectivity are modulated by conformational dynamics of a protein loop that covers the substrate access to the iron-carbene active species. Mutagenesis of residues computationally predicted to control the loop conformation altered the protein's chemoselectivity from preferred silylation to preferred amination of a substrate containing both N-H and Si-H functionalities. We demonstrate that information on protein structure and conformational dynamics, combined with knowledge of mechanism, leads to understanding of how non-natural and selective chemical transformations can be introduced into the biological world.

Biocatalysis offers an expanding and powerful strategy to construct and diversify complex molecules by C─H bond functionalization. Due to their high selectivity, enzymes have become an essential tool for C─H bond functionalization and offer complementary reactivity to small-molecule catalysts. Hemoproteins, particularly cytochromes P450, have proven effective for selective oxidation of unactivated C─H bonds. Previously, we reported the in vitro characterization of an oxidative tailoring cascade in which TamI, a multifunctional P450 functions co-dependently with the TamL flavoprotein to catalyze regio- and stereoselective hydroxylations and epoxidation to yield tirandamycin A and tirandamycin B. TamI follows a defined order including 1) C10 hydroxylation, 2) C11/C12 epoxidation, and 3) C18 hydroxylation. Here we present a structural, biochemical, and computational investigation of TamI to understand the molecular basis of its substrate binding, diverse reactivity, and specific reaction sequence. The crystal structure of TamI in complex with tirandamycin C together with molecular dynamics simulations and targeted mutagenesis suggest that hydrophobic interactions with the polyene chain of its natural substrate are critical for molecular recognition. QM calculations and molecular dynamics simulations of TamI with variant substrates provided detailed information on the molecular basis of sequential reactivity, and pattern of regio- and stereo-selectivity in catalyzing the three-step oxidative cascade.

Compared to the biological world's rich chemistry for functionalizing carbon, enzymatic transformations of the heavier homologue silicon are rare. We report that a wild-type cytochrome P450 monooxygenase (P450BM3 from Bacillus megaterium, CYP102A1) has promiscuous activity for oxidation of hydrosilanes to give silanols. Directed evolution was applied to enhance this non-native activity and create a highly efficient catalyst for selective silane oxidation under mild conditions with oxygen as the terminal oxidant. The evolved enzyme leaves C-H bonds present in the silane substrates untouched, and this biotransformation does not lead to disiloxane formation, a common problem in silanol syntheses. Computational studies reveal that catalysis proceeds through hydrogen atom abstraction followed by radical rebound, as observed in the native C-H hydroxylation mechanism of the P450 enzyme. This enzymatic silane oxidation extends nature's impressive catalytic repertoire.

Evolution has converged on cation-π interactions for recognition of quaternary alkyl ammonium groups such as trimethyllysine (Kme3). While computational modelling indicates that Trp provides the strongest cation-π interaction of the native aromatic amino acids, there is limited corroborative data from measurements within proteins. Herein we investigate a Tyr to Trp mutation in the binding pocket of the HP1 chromodomain, a reader protein that recognizes Kme3. Binding studies demonstrate that the Trp-mediated cation-π interaction is about -5 kcal mol-1 stronger, and the Y24W crystal structure shows that the mutation is not perturbing. Quantum mechanical calculations indicate that greater enthalpic binding is predominantly due to increased cation-π interactions. NMR studies indicate that differences in the unbound state of the Y24W mutation lead to enthalpy-entropy compensation. These results provide direct experimental quantification of Trp versus Tyr in a cation-π interaction and afford insight into the conservation of aromatic cage residues in Kme3 reader domains.

Cytochromes P450 (P450s) are nature's catalysts of choice for performing demanding and physiologically vital oxidation reactions. Biochemical characterization of these enzymes over the past decades has provided detailed mechanistic insight and highlighted the diversity of substrates P450s accommodate and the spectrum of oxidative transformations they catalyze. Previously, we discovered that the bacterial P450 MycCI from the mycinamicin biosynthetic pathway in Micromonospora griseorubida possesses an unusually broad substrate scope, whereas the homologous P450 from tylosin-producing Streptomyces fradiae (TylHI) exhibits a high degree of specificity for its native substrate. Here, using biochemical, structural, and computational approaches, we aimed to understand the molecular basis for the disparate reactivity profiles of these two P450s. Turnover and equilibrium binding experiments with substrate analogs revealed that TylHI strictly prefers 16-membered ring macrolides bearing the deoxyamino sugar mycaminose. To help rationalize these results, we solved the X-ray crystal structure of TylHI in complex with its native substrate at 1.99-Å resolution and assayed several site-directed mutants. We also conducted molecular dynamics simulations of TylHI and MycCI and biochemically characterized a third P450 homolog from the chalcomycin biosynthetic pathway in Streptomyces bikiniensis These studies provided a basis for constructing P450 chimeras to gain further insight into the features dictating the differences in reaction profile among these structurally and functionally related enzymes, ultimately unveiling the central roles of key loop regions in influencing substrate binding and turnover. Our work highlights the complex nature of P450/substrate interactions and raises interesting questions regarding the evolution of functional diversity among biosynthetic enzymes.

Microbial conversion of aromatic compounds is an emerging and promising strategy for valorization of the plant biopolymer lignin. A critical and often rate-limiting reaction in aromatic catabolism is O-aryl-demethylation of the abundant aromatic methoxy groups in lignin to form diols, which enables subsequent oxidative aromatic ring-opening. Recently, a cytochrome P450 system, GcoAB, was discovered to demethylate guaiacol (2-methoxyphenol), which can be produced from coniferyl alcohol-derived lignin, to form catechol. However, native GcoAB has minimal ability to demethylate syringol (2,6-dimethoxyphenol), the analogous compound that can be produced from sinapyl alcohol-derived lignin. Despite the abundance of sinapyl alcohol-based lignin in plants, no pathway for syringol catabolism has been reported to date. Here we used structure-guided protein engineering to enable microbial syringol utilization with GcoAB. Specifically, a phenylalanine residue (GcoA-F169) interferes with the binding of syringol in the active site, and on mutation to smaller amino acids, efficient syringol O-demethylation is achieved. Crystallography indicates that syringol adopts a productive binding pose in the variant, which molecular dynamics simulations trace to the elimination of steric clash between the highly flexible side chain of GcoA-F169 and the additional methoxy group of syringol. Finally, we demonstrate in vivo syringol turnover in Pseudomonas putida KT2440 with the GcoA-F169A variant. Taken together, our findings highlight the significant potential and plasticity of cytochrome P450 aromatic O-demethylases in the biological conversion of lignin-derived aromatic compounds.

There are few biocatalytic transformations that produce fluorine-containing molecules prevalent in modern pharmaceuticals. To expand the scope of biocatalysis for organofluorine synthesis, we have developed an enzymatic platform for highly enantioselective carbene B-H bond insertion to yield versatile α-trifluoromethylated (α-CF3) organoborons, an important class of organofluorine molecules that contain stereogenic centers bearing both CF3 and boron groups. In contrast to current "carbene transferase" enzymes that use a limited set of simple diazo compounds as carbene precursors, this system based on Rhodothermus marinus cytochrome c (Rma cyt c) can accept a broad range of trifluorodiazo alkanes and deliver versatile chiral α-CF3 organoborons with total turnovers up to 2870 and enantiomeric ratios up to 98.5:1.5. Computational modeling reveals that this broad diazo scope is enabled by an active-site environment that directs the alkyl substituent on the heme CF3-carbene intermediate toward the solvent-exposed face, thereby allowing the protein to accommodate diazo compounds with diverse structural features.

Cytochrome P450 enzymes have gained significant interest as selective oxidants in late-stage chemical synthesis. Their broad substrate scope enables them to be good candidates for their use in non-natural reactivity. Directed evolution evolves new enzyme biocatalysts that promote alternative reactivity for chemical synthesis. While directed evolution has proven useful in developing biocatalysts for specific purposes, this process is very time and labor intensive, and therefore not easily repurposed. Computational analysis of these P450 enzymes provides great insights into the broad substrate scope, the variety of reactions catalyzed, the binding specificity and the study of novel biosynthetic reaction mechanisms. By discovering new P450s and studying their reactivities, we uncover new insights into how this reactivity can be harnessed. We discuss a standard protocol using both DFT calculations and MD simulations to study a variety of cytochrome P450 enzymes. The approach entails theozyme models to study the mechanism and transition states via DFT calculations and subsequent MD simulations to understand the conformational poses and binding mechanisms within the enzyme. We discuss a few examples done in collaboration with the Tang and Sherman/Montgomery groups toward elucidating enzyme mechanisms and rationally designing new enzyme mutants as tools for selective C-H functionalization methods.

Controlling the regioselectivity of Baeyer-Villiger (BV) reactions remains an ongoing issue in organic chemistry, be it by synthetic catalysts or enzymes of the type Baeyer-Villiger monooxygenases (BVMOs). Herein, we address the challenging problem of switching normal to abnormal BVMO regioselectivity by directed evolution using three linear ketones as substrates, which are not structurally biased toward abnormal reactivity. Upon applying iterative saturation mutagenesis at sites lining the binding pocket of the thermostable BVMO from Thermocrispum municipale DSM 44069 (TmCHMO) and using 4-phenyl-2-butanone as substrate, the regioselectivity was reversed from 99:1 (wild-type enzyme in favor of the normal product undergoing 2-phenylethyl migration) to 2:98 in favor of methyl migration when applying the best mutant. This also stands in stark contrast to the respective reaction using the synthetic reagent m-CPBA, which provides solely the normal product. Reversal of regioselectivity was also achieved in the BV reaction of two other linear ketones. Kinetic parameters and melting temperatures revealed that most of the evolved mutants retained catalytic activity, as well as thermostability. In order to shed light on the origin of switched regioselectivity in reactions of 4-phenyl-2-butanone and phenylacetone, extensive QM/MM and MD simulations were performed. It was found that the mutations introduced by directed evolution induce crucial changes in the conformation of the respective Criegee intermediates and transition states in the binding pocket of the enzyme. In mutants that destabilize the normally preferred migration transition state, a reversal of regioselectivity is observed. This conformational control of regioselectivity overrides electronic control, which normally causes preferential migration of the group that is best able to stabilize positive charge. The results can be expected to aid future protein engineering of BVMOs.

Recently, heme proteins have been discovered and engineered by directed evolution to catalyze chemical transformations that are biochemically unprecedented. Many of these nonnatural enzyme-catalyzed reactions are assumed to proceed through a catalytic iron porphyrin carbene (IPC) intermediate, although this intermediate has never been observed in a protein. Using crystallographic, spectroscopic, and computational methods, we have captured and studied a catalytic IPC intermediate in the active site of an enzyme derived from thermostable Rhodothermus marinus (Rma) cytochrome c High-resolution crystal structures and computational methods reveal how directed evolution created an active site for carbene transfer in an electron transfer protein and how the laboratory-evolved enzyme achieves perfect carbene transfer stereoselectivity by holding the catalytic IPC in a single orientation. We also discovered that the IPC in Rma cytochrome c has a singlet ground electronic state and that the protein environment uses geometrical constraints and noncovalent interactions to influence different IPC electronic states. This information helps us to understand the impressive reactivity and selectivity of carbene transfer enzymes and offers insights that will guide and inspire future engineering efforts.

Microbial aromatic catabolism offers a promising approach to convert lignin, a vast source of renewable carbon, into useful products. Aryl-O-demethylation is an essential biochemical reaction to ultimately catabolize coniferyl and sinapyl lignin-derived aromatic compounds, and is often a key bottleneck for both native and engineered bioconversion pathways. Here, we report the comprehensive characterization of a promiscuous P450 aryl-O-demethylase, consisting of a cytochrome P450 protein from the family CYP255A (GcoA) and a three-domain reductase (GcoB) that together represent a new two-component P450 class. Though originally described as converting guaiacol to catechol, we show that this system efficiently demethylates both guaiacol and an unexpectedly wide variety of lignin-relevant monomers. Structural, biochemical, and computational studies of this novel two-component system elucidate the mechanism of its broad substrate specificity, presenting it as a new tool for a critical step in biological lignin conversion.

Duclauxins are dimeric and heptacyclic fungal polyketides with notable bioactivities. We characterized the cascade of redox transformations in the biosynthetic pathway of duclauxin from Talaromyces stipitatus. The redox reaction sequence is initiated by a cupin family dioxygenase DuxM that performs an oxidative cleavage of the peri-fused tricyclic phenalenone and affords a transient hemiketal-oxaphenalenone intermediate. Additional redox enzymes then morph the oxaphenoalenone into either an anhydride or a dihydrocoumarin-containing monomeric building block that is found in dimeric duxlauxins. Oxidative coupling between the monomers to form the initial C-C bond was shown to be catalyzed by a P450 monooxygenase, although the enzyme responsible for the second C-C bond formation was not found in the pathway. Collectively, the number and variety of redox enzymes used in the duclauxin pathway showcase Nature's strategy to generate structural complexity during natural product biosynthesis.

Hapalindole alkaloids are a structurally diverse class of cyanobacterial natural products defined by their varied polycyclic ring systems and diverse biological activities. These complex metabolites are generated from a common biosynthetic intermediate by the Stig cyclases in three mechanistic steps: a rare Cope rearrangement, 6-exo-trig cyclization, and electrophilic aromatic substitution. Here we report the structure of HpiC1, a Stig cyclase that catalyzes the formation of 12-epi-hapalindole U in vitro. The 1.5-Å structure revealed a dimeric assembly with two calcium ions per monomer and with the active sites located at the distal ends of the protein dimer. Mutational analysis and computational methods uncovered key residues for an acid-catalyzed [3,3]-sigmatropic rearrangement, as well as specific determinants that control the position of terminal electrophilic aromatic substitution, leading to a switch from hapalindole to fischerindole alkaloids.

Cytochrome P450 enzymes have gained significant interest as selective oxidants in late-stage chemical synthesis. Their broad substrate scope enables them to be good candidates for their use in non-natural reactivity. Directed evolution evolves new enzyme biocatalysts that promote alternative reactivity for chemical synthesis. While directed evolution has proven useful in developing biocatalysts for specific purposes, this process is very time and labor intensive, and therefore not easily repurposed. Computational analysis of these P450 enzymes provides great insights into the broad substrate scope, the variety of reactions catalyzed, the binding specificity and the study of novel biosynthetic reaction mechanisms. By discovering new P450s and studying their reactivities, we uncover new insights into how this reactivity can be harnessed. We discuss a standard protocol using both DFT calculations and MD simulations to study a variety of cytochrome P450 enzymes. The approach entails theozyme models to study the mechanism and transition states via DFT calculations and subsequent MD simulations to understand the conformational poses and binding mechanisms within the enzyme. We discuss a few examples done in collaboration with the Tang and Sherman/Montgomery groups toward elucidating enzyme mechanisms and rationally designing new enzyme mutants as tools for selective C-H functionalization methods.

Malbrancheamide is a dichlorinated fungal indole alkaloid isolated from both Malbranchea aurantiaca and Malbranchea graminicola that belongs to a family of natural products containing a characteristic bicyclo[2.2.2]diazaoctane core. The introduction of chlorine atoms on the indole ring of malbrancheamide differentiates it from other members of this family and contributes significantly to its biological activity. In this study, we characterized the two flavin-dependent halogenases involved in the late-stage halogenation of malbrancheamide in two different fungal strains. MalA and MalA' catalyze the iterative dichlorination and monobromination of the free substrate premalbrancheamide as the final steps in the malbrancheamide biosynthetic pathway. Two unnatural bromo-chloro-malbrancheamide analogues were generated through MalA-mediated chemoenzymatic synthesis. Structural analysis and computational studies of MalA' in complex with three substrates revealed that the enzyme represents a new class of zinc-binding flavin-dependent halogenases and provides new insights into a potentially unique reaction mechanism.

Hydroalkoxylation is a powerful and efficient method of forming C-O bonds and cyclic ethers in synthetic chemistry. In studying the biosynthesis of the fungal natural product herqueinone, we identified an enzyme that can perform an intramolecular enantioselective hydroalkoxylation reaction. PhnH catalyzes the addition of a phenol to the terminal olefin of a reverse prenyl group to give a dihydrobenzofuran product. The enzyme accelerates the reaction by 3 × 105-fold compared to the uncatalyzed reaction. PhnH belongs to a superfamily of proteins with a domain of unknown function (DUF3237), of which no member has a previously verified function. The discovery of PhnH demonstrates that enzymes can be used to promote the enantioselective hydroalkoxylation reaction and form cyclic ethers.

Epoxides are highly useful synthons and biosynthons for the construction of complex natural products during total synthesis and biosynthesis, respectively. Among enzyme-catalyzed epoxide transformations, a reaction that is notably missing, in regard to the synthetic toolbox, is cationic rearrangement that takes place under strong acid. This is a challenging transformation for enzyme catalysis, as stabilization of the carbocation intermediate upon epoxide cleavage is required. Here, we discovered two Brønsted acid enzymes that can catalyze two unprecedented epoxide transformations in biology. PenF from the penigequinolone pathway catalyzes a cationic epoxide rearrangement under physiological conditions to generate a quaternary carbon center, while AsqO from the aspoquinolone pathway catalyzes a 3-exo-tet cyclization to forge a cyclopropane-tetrahydrofuran ring system. The discovery of these new epoxide-modifying enzymes further highlights the versatility of epoxides in complexity generation during natural product biosynthesis.