Molecular Basis for Chemoselectivity Control in Oxidations of Internal Aryl-Alkenes Catalyzed by Laboratory Evolved P450s

P450 enzymes naturally perform selective oxidations of unfunctionalized hydrocarbon substrates, among other reactions. The adaptation of P450 enzymes to particular oxidative reactions involving alkenes is of great interest for the design of new biocatalysts. However, the mechanism that these enzymes utilize to precisely modulate the chemoselectivity and distinguishing between competing alkene double bond epoxidations and allylic C-H hydroxylations is sometimes not clear, which hampers the rational design of specific biocatalysts. In a previous work, P450LA1 was engineered in the laboratory using directed evolution to catalyze the direct oxidation of trans-b-methylstyrene to phenylacetone. The final variant, KS, was able to overcome the intrinsic preference for alkene epoxidation to directly generate a ketone product via the formation of a highly reactive carbocation intermediate. Here, additional library screening along this evolutionary lineage permitted to serendipitously detect a mutation that overcomes epoxidation and carbonyl formation by exhibiting a large selectivity of 94% towards allylic C-H hydroxylation. A multiscalar computational methodology was applied to reveal the molecular basis towards this hydroxylation preference. Enzyme modelling suggests that introduction of bulky substitution dramatically changes accessible conformations of the substrate in the active site, thus modifying the enzymatic selectivity towards terminal hydroxylation and avoiding the competing epoxidation pathway.

Efficient Transferase Engineering for SAM Analog Synthesis from Iodoalkanes

S-Adenosyl-l-methionine (SAM) is an important cosubstrate in various biochemical processes, including selective methyl transfer reactions. Simple methods for the (re)generation of SAM analogs could expand the chemistry accessible with SAM-dependent transferases and go beyond methylation reactions. Here we present an efficient enzyme engineering strategy to synthesize different SAM analogs from "off-the-shelf" iodoalkanes through enzymatic alkylation of S-adenosyl-l-homocysteine (SAH). This was achieved by mutating multiple hydrophobic and structurally dynamic amino acids simultaneously. Combinatorial mutagenesis was guided by the natural amino acid diversity and generated a highly functional mutant library. This approach increased the speed as well as the scale of enzyme engineering by providing a panel of optimized enzymes with orders of magnitude higher activities for multiple substrates in just one round of enzyme engineering. The optimized enzymes exhibit catalytic efficiencies up to 31 M-1 s-1, convert various iodoalkanes, including substrates bearing cyclopropyl or aromatic moieties, and catalyze S-alkylation of SAH with very high stereoselectivities (>99 % de). We further report a high throughput chromatographic screening system for reliable and rapid SAM analog analysis. We believe that the methods and enzymes described herein will further advance the field of selective biocatalytic alkylation chemistry by enabling SAM analog regeneration with "off-the-shelf" reagents.

Biocatalytic strategy for the construction of sp3-rich polycyclic compounds from directed evolution and computational modelling

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.

Enantio- and Diastereodivergent Cyclopropanation of Allenes by Directed Evolution of an Iridium-Containing Cytochrome

Alkylidene cyclopropanes (ACPs) are valuable synthetic intermediates because of their constrained structure and opportunities for further diversification. Although routes to ACPs are known, preparations of ACPs with control of both the configuration of the cyclopropyl (R vs S) group and the geometry of the alkene (E vs Z) are unknown. We describe enzymatic cyclopropanation of allenes with ethyl diazoacetate (EDA) catalyzed by an iridium-containing cytochrome (Ir(Me)-CYP119) that controls both stereochemical elements. Two mutants of Ir(Me)-CYP119 identified by 6-codon (6c, VILAFG) saturation mutagenesis catalyze the formation of (E)-ACPs with -93% to >99% ee and >99:1 E/Z ratio with just three rounds of 96 mutants. By four additional rounds of mutagenesis, an enzyme variant was identified that forms (Z)-ACPs with up to 94% ee and a 28:72 E/Z ratio. Computational studies show that the orientation of the carbene unit dictated by the mutated positions accounts for the stereoselectivity.

Reversing the Enantioselectivity of Enzymatic Carbene N-H Insertion Through Mechanism-Guided Protein Engineering

We report a computationally driven approach to access enantiodivergent enzymatic carbene N-H insertions catalyzed by P411 enzymes. Computational modeling was employed to rationally guide engineering efforts to control the accessible conformations of a key lactone-carbene (LAC) intermediate in the enzyme active site by installing a new H-bond anchoring point. This H-bonding interaction controls the relative orientation of the reactive carbene intermediate, orienting it for an enantioselective N-nucleophilic attack by the amine substrate. By combining MD simulations and site-saturation mutagenesis and screening targeted to only two key residues, we were able to reverse the stereoselectivity of previously engineered S-selective P411 enzymes. The resulting variant, L5_FL-B3, accepts a broad scope of amine substrates for N-H insertion with excellent yields (up to >99 %), high efficiency (up to 12 300 TTN), and good enantiocontrol (up to 7 : 93 er).

Protoglobin-Catalyzed Formation of cis-Trifluoromethyl-Substituted Cyclopropanes by Carbene Transfer

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.

Enzymatic Control over Reactive Intermediates Enables Direct Oxidation of Alkenes to Carbonyls by a P450 Iron-Oxo Species

The aerobic oxidation of alkenes to carbonyls is an important and challenging transformation in synthesis. Recently, a new P450-based enzyme (aMOx) has been evolved in the laboratory to directly oxidize styrenes to their corresponding aldehydes with high activity and selectivity. The enzyme utilizes a heme-based, high-valent iron-oxo species as a catalytic oxidant that normally epoxidizes alkenes, similar to other catalysts. How the evolved aMOx enzyme suppresses the commonly preferred epoxidation and catalyzes direct carbonyl formation is currently not well understood. Here, we combine computational modelling together with mechanistic experiments to study the reaction mechanism and unravel the molecular basis behind the selectivity achieved by aMOx. Our results describe that although both pathways are energetically accessible diverging from a common covalent radical intermediate, intrinsic dynamic effects determine the strong preference for epoxidation. We discovered that aMOx overrides these intrinsic preferences by controlling the accessible conformations of the covalent radical intermediate. This disfavors epoxidation and facilitates the formation of a carbocation intermediate that generates the aldehyde product through a fast 1,2-hydride migration. Electrostatic preorganization of the enzyme active site also contributes to the stabilization of the carbocation intermediate. Computations predicted that the hydride migration is stereoselective due to the enzymatic conformational control over the intermediate species. These predictions were corroborated by experiments using deuterated styrene substrates, which proved that the hydride migration is cis- and enantioselective. Our results demonstrate that directed evolution tailored a highly specific active site that imposes strong steric control over key fleeting biocatalytic intermediates, which is essential for accessing the carbonyl forming pathway and preventing competing epoxidation.

Engineered P450 Atom-Transfer Radical Cyclases are Bifunctional Biocatalysts: Reaction Mechanism and Origin of Enantioselectivity

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.

Changes in Protonation States of In-Pathway Residues can Alter Ligand Binding Pathways Obtained From Spontaneous Binding Molecular Dynamics Simulations

Protein-ligand binding processes often involve changes in protonation states that can be key to recognize and orient the ligand in the binding site. The pathways through which (bio)molecules interplay to attain productively bound complexes are intricate and involve a series of interconnected intermediate and transition states. Molecular dynamics (MD) simulations and enhanced sampling techniques are commonly used to characterize the spontaneous binding of a ligand to its receptor. However, the effect of protonation state changes of in-pathway residues in spontaneous binding MD simulations remained mostly unexplored. Here, we used molecular dynamics simulations to reconstruct the trypsin-benzamidine binding pathway considering different protonation states of His57. This residue is part of the trypsin catalytic triad and is located more than 10 Å away from Asp189, which is responsible for benzamidine binding in the trypsin S1 pocket. Our MD simulations showed that the binding pathways that benzamidine follow to target the S1 binding site are critically dependent on the His57 protonation state. Binding of benzamidine frequently occurs when His57 is protonated in the delta nitrogen while the binding process is significantly less frequent when His57 is positively charged. Constant-pH MD simulations retrieved the equilibrium populations of His57 protonation states at trypsin active pH offering a clearer picture of benzamidine recognition and binding. These results indicate that properly accounting for protonation states of distal residues can be important in spontaneous binding MD simulations.

Directed evolution of nonheme iron enzymes to access abiological radical-relay C(sp3)-H azidation

We report the reprogramming of nonheme iron enzymes to catalyze an abiological C(sp³)‒H azidation reaction through iron-catalyzed radical relay. This biocatalytic transformation uses amidyl radicals as hydrogen atom abstractors and Fe(III)‒N₃ intermediates as radical trapping agents. We established a high-throughput screening platform based on click chemistry for rapid evolution of the catalytic performance of identified enzymes. The final optimized variants deliver a range of azidation products with up to 10,600 total turnovers and 93% enantiomeric excess. Given the prevalence of radical relay reactions in organic synthesis and the diversity of nonheme iron enzymes, we envision that this discovery will stimulate future development of metalloenzyme catalysts for synthetically useful transformations unexplored by natural evolution.

An Enzymatic Platform for Primary Amination of 1-Aryl-2-alkyl Alkynes

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(sp³)-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 Unexplored Importance of Fleeting Chiral Intermediates in Enzyme-Catalyzed Reactions

Decades of extensive research efforts by biochemists, organic chemists, and protein engineers have led to an understanding of the basic mechanisms of essentially all known types of enzymes, but in a formidable number of cases an essential aspect has been overlooked. The occurrence of short-lived chiral intermediates formed by symmetry-breaking of prochiral precursors in enzyme catalyzed reactions has been systematically neglected. We designate these elusive species as fleeting chiral intermediates and analyze such crucial questions as "Do such intermediates occur in homochiral form?" If so, what is the absolute configuration, and why did Nature choose that particular stereoisomeric form, even when the isolable final product may be achiral? Does the absolute configuration of a chiral product depend in any way on the absolute configuration of the fleeting chiral precursor? How does this affect the catalytic proficiency of the enzyme? If these issues continue to be unexplored, then an understanding of the mechanisms of many enzyme types remains incomplete. We have systematized the occurrence of these chiral intermediates according to their structures and enzyme types. This is followed by critical analyses of selected case studies and by final conclusions and perspectives. We hope that the fascinating concept of fleeting chiral intermediates will attract the attention of scientists, thereby opening an exciting new research field.

Engineering P450 TamI as an Iterative Biocatalyst for Selective Late-Stage C-H Functionalization and Epoxidation of Tirandamycin Antibiotics

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.

Accessing Chemo- and Regioselective Benzylic and Aromatic Oxidations by Protein Engineering of an Unspecific Peroxygenase

Unspecific peroxygenases (UPOs) enable oxyfunctionalizations of a broad substrate range with unparalleled activities. Tailoring these enzymes for chemo- and regioselective transformations represents a grand challenge due to the difficulties in their heterologous productions. Herein, we performed protein engineering in Saccharomyces cerevisiae using the MthUPO from Myceliophthora thermophila. More than 5300 transformants were screened. This protein engineering led to a significant reshaping of the active site as elucidated by computational modelling. The reshaping was responsible for the increased oxyfunctionalization activity, with improved k_cat/K_m values of up to 16.5-fold for the model substrate 5-nitro-1,3-benzodioxole. Moreover, variants were identified with high chemo- and regioselectivities in the oxyfunctionalization of aromatic and benzylic carbons, respectively. The benzylic hydroxylation was demonstrated to perform with enantioselectivities of up to 95% ee. The proposed evolutionary protocol and rationalization of the enhanced activities and selectivities acquired by MthUPO variants represent a step forward toward the use and implementation of UPOs in biocatalytic synthetic pathways of industrial interest.

Origin and Control of Chemoselectivity in Cytochrome c Catalyzed Carbene Transfer into Si-H and N-H bonds

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.

Simultaneous screening of multiple substrates with an unspecific peroxygenase enabled modified alkane and alkene oxyfunctionalisations

Protein engineering of an unspecific peroxygenase (UPO) was performed with three substrates and six products in parallel by a high throughput GC-MS setup. Modified chemo- and regioselective variants were identified for aliphatic substrates.

Molecular Basis of Iterative C─H Oxidation by TamI, a Multifunctional P450 monooxygenase from the Tirandamycin Biosynthetic Pathway

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.

Selective Enzymatic Oxidation of Silanes to Silanols

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.

Regioselective Synthesis and Characterization of Tris- and Tetra-Prato Adducts of M3N@C80 (M = Y, Gd)

The tris- and tetra-adducts of M₃N@I_h-C₈₀ metallofullerenes were synthesized and characterized for the first time. The 1,3-dipolar cycloaddition (Prato reaction) of Y₃N@I_h-C₈₀ and Gd₃N@I_h-C₈₀ with an excess of N-ethylglycine and formaldehyde provided tris- and tetra-fulleropyrrolidine adducts in a regioselective manner. Purification by HPLC and analyses of the isolated peaks by NMR, MS, and vis-NIR spectra revealed that the major products were four tris- and one tetra-isomers for both Y₃N@I_h-C₈₀ and Gd₃N@I_h-C₈₀. Considering the large number of possible isomers (e.g., at least 1140 isomers for the tris-adduct), the limited number of isomers obtained indicated that the reactions proceeded with high regioselectivity. NMR analyses of the Y₃N@I_h-C₈₀ adducts found that the tris-adducts were all-[6,6]- or [6,6][6,6][5,6]-isomers and that some showed mutual isomerization or remained intact at room temperature. The tetra-adduct obtained as a major product was all-[6,6] and stable. For the structural elucidation of Gd₃N@I_h-C₈₀ tris- and tetra-adducts, density functional theory (DFT) calculations were performed to estimate the relative stabilities of tris- and tetra-adducts formed upon Prato functionalization of the most pyramidalized regions of the fullerene structure. The most stable structures corresponded to additions on the most pyramidalized (i.e., strained) bonds. Taking together the experimental vis-NIR spectra, NMR assignments, and the computed relative DFT stabilities of the potential tris- and tetra-adducts, the structures of the isolated adducts were elucidated. Electron resonance (ESR) measurements measurements of pristine, bis-, and tris-adducts of Gd₃N@C₈₀ suggested that the rotation of the endohedral metal cluster slowed upon increase of the addition numbers to C₈₀ cage, which is favored for accommodating the Gd atoms of the relatively large Gd₃N cluster inner space at the sp³ addition sites. This is presumably related to the high regioselectivity in the Prato addition reaction driven by the strain release of the Gd₃N@C₈₀ fullerene structure.

In Vivo Selection for Formate Dehydrogenases with High Efficiency and Specificity toward NADP. ACS Catal 2020;10:7512-7525. [PMID: 32733773 PMCID: PMC7384739 DOI: 10.1021/acscatal.0c01487]

The efficient regeneration of cofactors is vital for the establishment of biocatalytic processes. Formate is an ideal electron donor for cofactor regeneration due to its general availability, low reduction potential, and benign byproduct (CO₂). However, formate dehydrogenases (FDHs) are usually specific to NAD⁺, such that NADPH regeneration with formate is challenging. Previous studies reported naturally occurring FDHs or engineered FDHs that accept NADP⁺, but these enzymes show low kinetic efficiencies and specificities. Here, we harness the power of natural selection to engineer FDH variants to simultaneously optimize three properties: kinetic efficiency with NADP⁺, specificity toward NADP⁺, and affinity toward formate. By simultaneously mutating multiple residues of FDH from Pseudomonas sp. 101, which exhibits practically no activity toward NADP⁺, we generate a library of >10⁶ variants. We introduce this library into an E. coli strain that cannot produce NADPH. By selecting for growth with formate as the sole NADPH source, we isolate several enzyme variants that support efficient NADPH regeneration. We find that the kinetically superior enzyme variant, harboring five mutations, has 5-fold higher efficiency and 14-fold higher specificity in comparison to the best enzyme previously engineered, while retaining high affinity toward formate. By using molecular dynamics simulations, we reveal the contribution of each mutation to the superior kinetics of this variant. We further determine how nonadditive epistatic effects improve multiple parameters simultaneously. Our work demonstrates the capacity of in vivo selection to identify highly proficient enzyme variants carrying multiple mutations which would be almost impossible to find using conventional screening methods.

Thermodynamic consequences of Tyr to Trp mutations in the cation-π-mediated binding of trimethyllysine by the HP1 chromodomain

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⁻¹ 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.

In this work, we experimentally validate that tryptophan provides the strongest cation–π binding interaction among aromatic amino acids and also lend insight into the importance of residue identity in trimethyllysine recognition by reader proteins.

Exploring the molecular basis for substrate specificity in homologous macrolide biosynthetic cytochromes P450

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.

Site-Selectivity of Prato Additions to C70: Experimental and Theoretical Studies of a New Thermodynamic Product at the dd-[5,6]-Junction

Three Prato monoadduct isomers were synthesized and structurally characterized by ¹H, ¹³C NMR spectra and single-crystal X-ray diffraction, and one adduct on the dd-[5,6]-bond was found as the first example of a Prato [5,6]-adduct of C₇₀. To investigate the mechanism in the generation of this dd-[5,6]-adduct, computational studies were employed to show that it was thermodynamically obtained by sigmatropic rearrangement from the presumed initial kinetic product de-[6,6]-adduct.

Structures of Gd3N@C80 Prato Bis-Adducts: Crystal Structure, Thermal Isomerization, and Computational Study

The structures of two bis-ethylpyrrolidinoadducts of Gd₃N@I_h-C₈₀, obtained by regioselective 1,3-dipolar cycloadditions, were elucidated by single crystal X-ray, visible-near infrared (vis-NIR) spectra, studies on their thermal isomerization, and theoretical calculations. The structure of the minor-bis-adduct reveals a C₂-symmetric carbon cage with [6,6][6,6]-addition sites and with an endohedral Gd₃N cluster that is completely flattened. This is the first example of a crystal structure of Gd₃N@I_h-C₈₀ derivatives. The structure of the major-bis-adduct was inferred by the vis-NIR spectrum being corresponded to the structure of a previously reported major-bis-adduct of Y₃N@I_h-C₈₀ known to have an asymmetric [6,6][6,6]-structure. Based on experimental results showing that the minor-bis-adduct of Gd₃N@I_h-C₈₀ isomerized to the major-adduct, a possible second addition site was elucidated with support from density functional theory calculations.

A Biocatalytic Platform for Synthesis of Chiral α-Trifluoromethylated Organoborons

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.

Size-selective encapsulation of C60 and C60-derivatives within an adaptable naphthalene-based tetragonal prismatic supramolecular nanocapsule

A novel naphthalene-based 5·(BArF)8 capsule allows for the size-selective inclusion of C60 from fullerene mixtures. Its size selectivity towards C60 has been rationalized by its dynamic adaptability in solution that has been investigated by molecular dynamics. Additionally, 5·(BArF)8 encapsulates C60-derivatives such as C60-PCBM and N-methylpyrrolidine-C60. The latter can be separated from C60 since 5·(BArF)8 displays distinct affinity for them.

Computational Protocol to Understand P450 Mechanisms and Design of Efficient and Selective Biocatalysts

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.

Overriding Traditional Electronic Effects in Biocatalytic Baeyer-Villiger Reactions by Directed Evolution

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.

Catalytic iron-carbene intermediate revealed in a cytochrome c carbene transferase

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.

A promiscuous cytochrome P450 aromatic O-demethylase for lignin bioconversion

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.

Biosynthesis of Heptacyclic Duclauxins Requires Extensive Redox Modifications of the Phenalenone Aromatic Polyketide

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.

On the regioselectivity of the Diels-Alder cycloaddition to C60 in high spin states

Controlling the regioselectivity in the exohedral functionalization of fullerenes and endohedral metallofullerenes is essential to produce specific desired fullerene derivatives. In this work, using density functional theory (DFT) calculations, we show that the regioselectivity of the Diels-Alder (DA) cycloaddition of cyclopentadiene to 2S+1C60 changes from the usual [6,6] addition in the singlet ground state to the [5,6] attack in high spin states of C60. Changes in the aromaticity of the five- and six-membered rings when going from singlet to high spin C60 provide a rationale to understand this regioselectivity change. Experimentally, however, we find that the DA cycloaddition of isoindene to triplet C60 yields the usual [6,6] adduct. Further DFT calculations and computational analysis give an explanation to this unanticipated experimental result by showing the presence of an intersystem crossing close to the formed triplet biradical intermediate.

Structural basis of the Cope rearrangement and cyclization in hapalindole biogenesis

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.

Effect of incarcerated HF on the exohedral chemical reactivity of HF@C60

The first chemical modification on the brand new endohedral HF@C₆₀ is reported. In particular, the isomerization from optically pure (2S,5S)-cis-pyrrolidino[3,4:1,2][60]fullerene 2b to (2S,5R)-trans-pyrrolidino[3,4:1,2][60]fullerene 2b has been studied and compared with empty C₆₀ (2a) and endohedral H₂O@C₆₀ (3). The comparative study shows a kinetic order for the isomerization process of H₂O@C₆₀ > HF@C₆₀ > C₆₀, thus confirming the effect of the incarcerated species on the zwitterionic intermediate stability.

Exploring the origins of selectivity in soluble epoxide hydrolase from Bacillus megaterium

Epoxide hydrolase (EH) enzymes catalyze the hydration of racemic epoxides to yield their corresponding vicinal diols. In this work, the Bacillus megaterium epoxide hydrolase (BmEH)-mediated hydrolysis of racemic styrene oxide (rac-SO) and its para-nitro styrene oxide (rac-p-NSO) derivative are computationally investigated using density functional theory (DFT).

Epoxide hydrolase (EH) enzymes catalyze the hydration of racemic epoxides to yield their corresponding vicinal diols. These enzymes present different enantio- and regioselectivity depending upon either the substrate structure or the substitution pattern of the epoxide ring. In this study, we computationally investigate the Bacillus megaterium epoxide hydrolase (BmEH)-mediated hydrolysis of racemic styrene oxide (rac-SO) and its para-nitro styrene oxide (rac-p-NSO) derivative using density functional theory (DFT) and an active site cluster model consisting of 195 and 197 atoms, respectively. Full reaction mechanisms for epoxide ring opening were evaluated considering the attack at both oxirane carbons and considering two possible orientations of the substrate at the BmEH active site. Our results indicate that for both SO and p-NSO substrates the BmEH enantio- and regioselectivity is opposite to the inherent (R)-BmEH selectivity, the attack at the benzylic position (C1) of the (S)-enantiomer being the most favoured chemical outcome.

Role of Conformational Dynamics in the Evolution of Retro-Aldolase Activity

Enzymes exist as ensembles of conformations that are important for function. Tuning these populations of conformational states through mutation enables evolution toward additional activities. Here we computationally evaluate the population shifts induced by distal and active site mutations in a family of computationally designed and experimentally optimized retro-aldolases. The conformational landscape of these enzymes was significantly altered during evolution, as pre-existing catalytically active conformational substates became major states in the most evolved variants. We further demonstrate that key residues responsible for these substate conversions can be predicted computationally. Significantly, the identified residues coincide with those positions mutated in the laboratory evolution experiments. This study establishes that distal mutations that affect enzyme catalytic activity can be predicted computationally and thus provides the enzyme (re)design field with a rational strategy to determine promising sites for enhancing activity through mutation.

Computational Design of Protein Function

The computational design of enzymes is a tremendous challenge for both chemistry and biochemistry. The ability to design stable and functional biocatalysts that could operate under different conditions to perform chemical reactions without precedent in nature, allowing the large-scale production of chemicals à la carte, would revolutionise both synthetic, pharmacologic and materials chemistry. Despite the great advances achieved, this highly multidisciplinary area of research is still in its infancy. This chapter describes the ‘inside-out’ protocol for computational enzyme design and both the achievements and limitations of the current technology are highlighted. Furthermore, molecular dynamics simulations have proved to be invaluable in the enzyme design process, constituting an important tool for discovering elusive catalytically relevant conformations of the engineered or designed enzyme. As a complement to the ‘inside-out’ design protocol, different examples where hybrid QM/MM approaches have been directly applied to discover beneficial mutations in rational computational enzyme design are described.

Function and Structure of MalA/MalA', Iterative Halogenases for Late-Stage C-H Functionalization of Indole Alkaloids

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.

Rationalizing the relative abundances of trimetallic nitride template-based endohedral metallofullerenes from aromaticity measures

The synthesis of endohedral metallofullerenes (EMFs) from a carbon soot sample of an arc discharge leads to a variety of EMFs that are obtained in different relative abundances. In the present work, we show that these abundances can be predicted from aromaticity calculations. In particular, we use the normalized Additive Local Aromaticity (ALAN) index. Our results show that the most abundant Sc3N-based and Y3N-based EMFs in fullerene soot are the most aromatic. This study reinforces the idea that aromaticity plays a key role in determining the stability of EMFs.

Enzyme-Catalyzed Intramolecular Enantioselective Hydroalkoxylation

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 × 10⁵-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.

Computational tools for the evaluation of laboratory-engineered biocatalysts

Biocatalysis is based on the application of natural catalysts for new purposes, for which enzymes were not designed. Although the first examples of biocatalysis were reported more than a century ago, biocatalysis was revolutionized after the discovery of an in vitro version of Darwinian evolution called Directed Evolution (DE). Despite the recent advances in the field, major challenges remain to be addressed. Currently, the best experimental approach consists of creating multiple mutations simultaneously while limiting the choices using statistical methods. Still, tens of thousands of variants need to be tested experimentally, and little information is available on how these mutations lead to enhanced enzyme proficiency. This review aims to provide a brief description of the available computational techniques to unveil the molecular basis of improved catalysis achieved by DE. An overview of the strengths and weaknesses of current computational strategies is explored with some recent representative examples. The understanding of how this powerful technique is able to obtain highly active variants is important for the future development of more robust computational methods to predict amino-acid changes needed for activity.

Reactivity of Single-Walled Carbon Nanotubes in the Diels-Alder Cycloaddition Reaction: Distortion-Interaction Analysis along the Reaction Pathway

Diels-Alder cycloaddition is one of the most powerful tools for the functionalization of single-walled carbon nanotubes (SWCNTs). Density functional theory at the B3-LYP level of theory has been used to investigate the reactivity of different-diameter SWCNTs (4-9,5) in Diels-Alder reactions with 1,3-butadiene; the reactivity was found to decrease with increasing SWCNT diameter. Distortion/interaction analysis along the whole reaction pathway was found to be a better way to explore the reactivity of this type of reaction. The difference in interaction energy along the reaction pathway is larger than that of the corresponding distortion energy. However, the distortion energy plots for these reactions show the same trend. Therefore, the formation of the transition state can be determined from the interaction energy. A lower interaction energy leads to an earlier transition state, which indicates a lower activation energy. The computational results also indicate that the original distortion of the SWCNTs leads to an increase in the reactivity of the SWCNTs.

On the physical origins of interaction-induced vibrational (hyper)polarizabilities

This paper presents the results of a pioneering exploration of the physical origins of vibrational contributions to the interaction-induced electric properties of molecular complexes. In order to analyze the excess nuclear relaxation (hyper)polarizabilities, a new scheme was proposed which relies on the computationally efficient Bishop-Hasan-Kirtman method for determining the nuclear relaxation contributions to electric properties. The extension presented herein is general and can be used with any interaction-energy partitioning method. As an example, in this study we employed the variational-perturbational interaction-energy decomposition scheme (at the MP2/aug-cc-pVQZ level) and the extended transition state method by employing three exchange-correlation functionals (BLYP, LC-BLYP, and LC-BLYP-dDsC) to study the excess properties of the HCN dimer. It was observed that the first-order electrostatic contribution to the excess nuclear relaxation polarizability cancels with the negative exchange repulsion term out to a large extent, resulting in a positive value of Δα(nr) due to the contributions from the delocalization and the dispersion terms. In the case of the excess nuclear relaxation first hyperpolarizability, the pattern of interaction contributions is very similar to that for Δα(nr), both in terms of their sign as well as relative magnitude. Finally, our results show that the LC-BLYP and LC-BLYP-dDsC functionals, which yield smaller values of the orbital relaxation term than BLYP, are more successful in predicting excess properties.

Reaction Mechanism and Regioselectivity of the Bingel-Hirsch Addition of Dimethyl Bromomalonate to La@C2v -C82

We quantum chemically explore the thermodynamics and kinetics of all 65 possible mechanistic pathways of the Bingel-Hirsch addition of dimethyl bromomalonate to the endohedral metallofullerene La@C2v -C82 that result from the combination of 24 nonequivalent carbon atoms and 35 different bonds present in La@C2v -C82 by using dispersion-corrected DFT calculations. Experimentally, this reaction leads to four singly bonded derivatives and one fulleroid adduct. Of these five products, only the singly bonded derivative on C23 could be experimentally identified unambiguously. Our calculations show that La@C2v -C82 is not particularly regioselective under Bingel-Hirsch conditions. From the obtained results, however, it is possible to make a tentative assignment of the products observed experimentally. We propose that the observed fulleroid adduct results from the attack at bond 19 and that the singly bonded derivatives correspond to the C2, C19, C21, and C23 initial attacks. However, other possibilities cannot be ruled out completely.