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

Using the Theozyme Model to Study the Dynamical Mechanism of the Post-Transition State Bifurcation Reaction by NgnD Enzyme

Post-transition state bifurcation (PTSB) is a fundamental process in which a single transition state leads to multiple products. This phenomenon is important in both biological and chemical contexts and offers valuable insights into reaction mechanisms and their applications. The theozyme model, which focuses on key residues within enzymes, offers a computationally efficient method for studying these processes while preserving the enzyme's catalytic properties. This approach enhances our understanding of how enzymes stabilize and direct the transition state, thereby influencing product distribution and selectivity. In this study, we investigate the dynamics and regulatory mechanisms of the PTSB reaction catalyzed by the enzyme NgnD. The enzyme NgnD facilitates a cycloaddition reaction that produces both [6 + 4] and [4 + 2] adducts, with a preference for the [6 + 4] adduct. By analyzing the potential energy surface, bond length distribution, and interactions between the theozyme and the ambimodal transition state, we elucidate the role of the enzyme's active site residues in determining product selectivity. We illustrate how these key residues contribute to the formation of different adducts, providing insights from various perspectives. Using theozyme models, we propose how the four most influential active residues collectively might control the direction of adduct formation through their cumulative effects.

Harnessing bio and (Photo)catalysts for microplastics degradation and remediation in soil environment

Soil pollution by microplastics (MPs) is an escalating environmental crisis with far-reaching consequences. However, current research on the degradation and/or remediation of MPs has mainly focused on water-simulated environments, with little attention given to soil MPs. Therefore, the review explores such terrestrial territory, exploring the potential of biodegradation and novel photocatalytic technologies for MPs degradation/remediation in soil. This review comprehensively investigates the potential of biological and photocatalytic approaches for soil MPs degradation and remediation. A temporal analysis of research from 2004 to 2024 highlights the increasing focus on this critical issue. The review explores the biocatalytic roles of diverse enzymes, including cutinase, PETase, MHETase, hydrolase, lipase, laccase, lignin peroxidase, and Mn-peroxidase, in MPs degradation. Strategies for enzyme engineering, such as protein engineering and immobilization, are explored to enhance catalytic efficiency. The potential for developing enzyme consortia for optimized MP degradation is also discussed. Photocatalytic remediation using TiO2, ZnO, clay, hydrogel, and other photocatalysts is examined, emphasizing their mechanisms and effectiveness. Computational modeling is proposed to deepen understanding of soil MPs-catalyst interactions, primarily aiming to develop novel catalysts tailored for soil environments for environmental safety and sustainable restoration. A comparative analysis of biological and photocatalytic approaches evaluates their environmental implications and the potential for synergistic combinations, with emphasis on soil quality protection, restoration and impact on soil ecosystems. Hence, this review accentuates the urgent need for innovative solutions to address MPs pollution in soil and provides a foundational understanding of the current knowledge gaps, as well as paves the way for future research and development.

Role of mutations in a chemoenzymatic enantiodivergent C(sp3)-H insertion: exploring the mechanism and origin of stereoselectivity

New-to-nature enzymes have emerged as powerful catalysts in recent years for streamlining various stereoselective organic transformations. While synthetic strategies employing engineered enzymes have witnessed proliferating success, there is limited clarity on the mechanistic front and more so when considering molecular-level insights into the role of selected mutations, dramatically escalating catalytic competency and selectivity. We have investigated the mechanism and correlation between mutations and exquisite stereoselectivity of a lactone carbene insertion into the C(sp3)-H bond of substituted aniline, catalyzed by two mutants of a cytochrome P450 variant, "P411" (engineered through directed evolution) in which the axial cysteine has been mutated to serine, utilizing various computational tools. The pivotal role of S264 and L/R328 mutations in the active site has been delineated computationally using two cluster models, thus rationalizing the enantiodivergence. This report provides much-needed insights into the origin of enantiodivergence, furnishing a mechanistic framework for understanding the anchoring effects of H-bond donor residues with the lactone ring. This study is expected to have important implications in the rational design of stereodivergent enzymes and toward successful in silico enzyme designing.

Computation-Aided Engineering of Cytochrome P450 for the Production of Pravastatin

CYP105AS1 is a cytochrome P450 from Amycolatopsis orientalis that catalyzes monooxygenation of compactin to 6-epi-pravastatin. For fermentative production of the cholesterol-lowering drug pravastatin, the stereoselectivity of the enzyme needs to be inverted, which has been partially achieved by error-prone PCR mutagenesis and screening. In the current study, we report further optimization of the stereoselectivity by a computationally aided approach. Using the CoupledMoves protocol of Rosetta, a virtual library of mutants was designed to bind compactin in a pro-pravastatin orientation. By examining the frequency of occurrence of beneficial substitutions and rational inspection of their interactions, a small set of eight mutants was predicted to show the desired selectivity and these variants were tested experimentally. The best CYP105AS1 variant gave >99% stereoselective hydroxylation of compactin to pravastatin, with complete elimination of the unwanted 6-epi-pravastatin diastereomer. The enzyme-substrate complexes were also examined by ultrashort molecular dynamics simulations of 50 × 100 ps and 5 × 22 ns, which revealed that the frequency of occurrence of near-attack conformations agreed with the experimentally observed stereoselectivity. These results show that a combination of computational methods and rational inspection could improve CYP105AS1 stereoselectivity beyond what was obtained by directed evolution. Moreover, the work lays out a general in silico framework for specificity engineering of enzymes of known structure.

C-N Bond Forming Radical Rebound Is the Enantioselectivity-Determining Step in P411-Catalyzed Enantioselective C(sp3)-H Amination: A Combined Computational and Experimental Investigation

Engineered metalloenzymes represent promising catalysts for stereoselective C-H functionalization reactions. Recently, P450 enzymes have been evolved to allow for new-to-nature intramolecular C(sp3)-H amination reactions via a nitrene transfer mechanism, giving rise to diamine derivatives with excellent enantiocontrol. To shed light on the origin of enantioselectivity, a combined computational and experimental study was carried out. Hybrid quantum mechanics/molecular mechanics calculations were performed to investigate the activation energies and enantioselectivities of both the hydrogen atom transfer (HAT) and the subsequent C-N bond forming radical rebound steps. Contrary to previously hypothesized enantioinduction mechanisms, our calculations show that the radical rebound step is enantioselectivity-determining, whereas the preceding HAT step is only moderately stereoselective. Furthermore, the selectivity in the initial HAT is ablated by rapid conformational change of the radical intermediate prior to C-N bond formation. This finding is corroborated by our experimental study using a set of enantiomerically pure, monodeuterated substrates. Furthermore, classical and ab initio molecular dynamics simulations were carried out to investigate the conformational flexibility of the carbon-centered radical intermediate. This key radical species undergoes a facile conformational change in the enzyme active site from the pro-(R) to the pro-(S) configuration, whereas the radical rebound is slower due to the spin-state change and ring strain of the cyclization process, thereby allowing stereoablative C-N bond formation. Together, these studies revealed an underappreciated enantioinduction mechanism in biocatalytic C(sp3)-H functionalizations involving radical intermediates, opening up new avenues for the development of other challenging asymmetric C(sp3)-H functionalizations.

Oriented External Electric Fields Regurating the Reaction Mechanism of CH4 Oxidation Catalyzed by Fe(IV)-Oxo-Corrolazine: Insight from Density Functional Calculations

Methane is the simplest alkane and can be used as an alternative energy source for oil and coal, but the greenhouse effect caused by its leakage into the air is not negligible, and its conversion into liquid methanol not only facilitates transportation, but also contributes to carbon neutrality. In order to find an efficient method for converting methane to methanol, CH4 oxidation catalyzed by Fe(IV)-Oxo-corrolazine (Fe(IV)-Oxo-Cz) and its reaction mechanism regulation by oriented external electric fields (OEEFs) are systematically studied by density functional calculations. The calculations show that Fe(IV)-Oxo-Cz can abstract one H atom from CH4 to form the intermediate with OH group connecting on the corrolazine ring, with the energy barrier of 25.44 kcal mol-1. And then the product methanol is formed through the following rebound reaction. Moreover, the energy barrier can be reduced to 20.72 kcal mol-1 through a two-state reaction pathway. Furthermore, the effect of OEEFs on the reaction is investigated. We found that OEEFs can effectively regulate the reaction by adjusting the stability of the reactant and the transition state through the interaction of electric field-molecular dipole moment. When the electric field is negative, the energy barrier of the reaction decreases with the increase of electric intensity. Moreover, the OEEF aligned along the intrinsic Fe‒O reaction axis can effectively regulate the ability of forming the OH on the corrolazine ring by adjusting the charges of O and H atoms. When the electric field intensity is -0.010 a.u., the OH can be directly rebounded to the CH3· before it is connecting on the corrolazine ring, thus forming the product directly from the transition state without passing through the intermediate with only an energy barrier of 17.34 kcal mol-1, which greatly improves the selectivity of the reaction.

Modeling the Ligand Effect on the Structure of CYP 450 Within the Density Functional Theory

An improved understanding of the P450 structure is relevant to the development of biomimetic catalysts and inhibitors for controlled CH-bond activation, an outstanding challenge of synthetic chemistry. Motivated by the experimental findings of an unusually short Fe-S bond of 2.18 Å for the wild-type (WT) OleT P450 decarboxylase relative to a cysteine pocket mutant form (A369P), a computational model that captures the effect of the thiolate axial ligand on the iron-sulfur distance is presented. With the computational efficiency and streamlined analysis in mind, this model combines a cluster representation of the enzyme─40-110 atoms, depending on the heme and ligand truncation level─with a density functional theory (DFT) description of the electronic structure (ES) and is calibrated against the experimental data. The optimized Fe-S distances show a difference of 0.25 Å between the low and high spin states, in agreement with the crystallographic structures of the OleT WT and mutant forms. We speculate that this difference is attributable to the packing of the ligand; the mutant is bulkier due to an alanine-to-proline replacement, meaning that it is excluded from the energetically favored low-spin minimum because of steric constraints. The presence of pure spin-state pairs and the intersection of the low/high spin states for the enzyme model is indicative of the limitations of single-reference ES methods in such systems and emphasizes the significance of using the proper state when modeling the hydrogen atom transfer (HAT) reaction catalyzed by OleT. At the same time, the correct characterization of both the short and long Fe-S bonds within a small DFT-based model of 42 atoms paves the way for quantum dynamics modeling of the HAT step, which initiates the OleT decarboxylation reaction.

Stereodivergent atom-transfer radical cyclization by engineered cytochromes P450

Naturally occurring enzymes can be a source of unnatural reactivity that can be molded by directed evolution to generate efficient biocatalysts with valuable activities. Owing to the lack of exploitable stereocontrol elements in synthetic systems, steering the absolute and relative stereochemistry of free-radical processes is notoriously difficult in asymmetric catalysis. Inspired by the innate redox properties of first-row transition-metal cofactors, we repurposed cytochromes P450 to catalyze stereoselective atom-transfer radical cyclization. A set of metalloenzymes was engineered to impose substantial stereocontrol over the radical addition step and the halogen rebound step in these unnatural processes, allowing enantio- and diastereodivergent radical catalysis. This evolvable metalloenzyme platform represents a promising solution to tame fleeting radical intermediates for asymmetric catalysis.

Computer-aided understanding and engineering of enzymatic selectivity

Enzymes offering chemo-, regio-, and stereoselectivity enable the asymmetric synthesis of high-value chiral molecules. Unfortunately, the drawback that naturally occurring enzymes are often inefficient or have undesired selectivity toward non-native substrates hinders the broadening of biocatalytic applications. To match the demands of specific selectivity in asymmetric synthesis, biochemists have implemented various computer-aided strategies in understanding and engineering enzymatic selectivity, diversifying the available repository of artificial enzymes. Here, given that the entire asymmetric catalytic cycle, involving precise interactions within the active pocket and substrate transport in the enzyme channel, could affect the enzymatic efficiency and selectivity, we presented a comprehensive overview of the computer-aided workflow for enzymatic selectivity. This review includes a mechanistic understanding of enzymatic selectivity based on quantum mechanical calculations, rational design of enzymatic selectivity guided by enzyme-substrate interactions, and enzymatic selectivity regulation via enzyme channel engineering. Finally, we discussed the computational paradigm for designing enzyme selectivity in silico to facilitate the advancement of asymmetric biosynthesis.

Natural Compounds as Pharmaceuticals: The Key Role of Cytochromes P450 Reactivity

The design of drugs from natural products is a re-emerging area due to the need for bioactive compounds. The exploitation of natural products and their derivatives obtained by biocatalysis is in line with the higher attention given today to new sustainable technologies that better preserve the environment (green chemistry). The research field of cytochromes P450 (CYPs) is continuously providing new enzymes and mutants that produce metabolites suitable for late-stage functionalization for new potential drugs. This review provides an overview of the exploitation of CYPs as biocatalysts in drug synthesis. Additionally, recent progress in protein and metabolic engineering is provided to show how these enzymes offer a toolbox that can be combined with other biocatalytic or chemical processes to build new platforms for the green production of new drugs.

Advances in Sustainable Catalysis: A Computational Perspective

The enormous challenge of moving our societies to a more sustainable future offers several exciting opportunities for computational chemists. The first principles approach to "catalysis by design" will enable new and much greener chemical routes to produce vital fuels and fine chemicals. This prospective outlines a wide variety of case studies to underscore how the use of theoretical techniques, from QM/MM to unrestricted DFT and periodic boundary conditions, can be applied to biocatalysis and to both homogeneous and heterogenous catalysts of all sizes and morphologies to provide invaluable insights into the reaction mechanisms they catalyze.

Use of Computational Biochemistry for Elucidating Molecular Mechanisms of Nitric Oxide Synthase

Nitric oxide (NO) is an essential signaling molecule in the regulation of multiple cellular processes. It is endogenously synthesized by NO synthase (NOS) as the product of L-arginine oxidation to L-citrulline, requiring NADPH, molecular oxygen, and a pterin cofactor. Two NOS isoforms are constitutively present in cells, nNOS and eNOS, and a third is inducible (iNOS). Despite their biological relevance, the details of their complex structural features and reactivity mechanisms are still unclear. In this review, we summarized the contribution of computational biochemistry to research on NOS molecular mechanisms. We described in detail its use in studying aspects of structure, dynamics and reactivity. We also focus on the numerous outstanding questions in the field that could benefit from more extensive computational investigations.