Quantum-Mechanical/Molecular-Mechanical Studies of CYP11A1-Catalyzed Biosynthesis of Pregnenolone from Cholesterol Reveal a C-C Bond Cleavage Reaction That Occurs by a Compound I-Mediated Electron Transfer

Electron transfer engineering of artificially designed cell factory for complete biosynthesis of steroids

Biosynthesis of steroids by artificially designed cell factories often involves numerous nicotinamide adenine dinucleotide phosphate (NADPH)-dependent enzymes that mediate electron transfer reactions. However, the unclear mechanisms of electron transfer from regeneration to the final delivery to the NADPH-dependent active centers limit systematically engineering electron transfer to improve steroids production. Here, we elucidate the electron transfer mechanisms of NADPH-dependent enzymes for systematically engineer electron transfer of Saccharomyces cerevisiae, including step-by-step engineering the electron transfer residues of 7-Dehydrocholesterol reductase (DHCR7) and P450 sterol side chain cleaving enzyme (P450scc), electron transfer components for directing carbon flux, and NADPH regeneration pathways, for high-level production of the cholesterol (1.78 g/L) and pregnenolone (0.83 g/L). The electron transfer engineering (ETE) process makes the electron transfer chains shorter and more stable which significantly accelerates deprotonation and proton coupled electron transfer process. This study underscores the significance of ETE strategies in steroids biosynthesis and expands synthetic biology approaches.

Computational Modeling of the Enzymatic Achmatowicz Rearrangement by Heme-Dependent Chloroperoxidase: Reaction Mechanism, Enantiopreference, Regioselectivity, and Substrate Specificity

The chloroperoxidase from Caldariomyces fumago (CfCPO) catalyzes the oxidative ring expansion of α-heterofunctionalized furans via the Achmatowicz rearrangement, providing an elegant tool to convert furan rings into complex-prefunctionalized scaffolds. However, the mechanism of this transformation remains unclear. Herein, the CfCPO-catalyzed reaction of rac-1-(2-furyl)ethanol (1a) is studied by quantum chemical calculations and molecular dynamics simulations. The calculations reveal that the conversion follows the general mechanism of the Achmatowicz reaction. Notably, the binding of 1a to the enzyme's active site influences the Compound I (Cpd I) formation, and the (R)-1a enantiomer binding results in a lower barrier compared to (S)-1a, explaining the observed (R)-enantiopreference toward a racemic substrate. Additionally, due to the weaker steric hindrance between the porphyrin ring and substrate, the nucleophilic attack of Cpd I on the furan core of 1a is preferred at the less-substituted C4=C5 bond, providing a rationale for the experimentally observed regioselectivity. Finally, the bottleneck residues in the substrate delivery channel and also the active site surroundings are proposed to be responsible for the substrate specificity of CfCPO. This study lays a theoretical foundation for the rational design of new CPOs that catalyze the Achmatowicz rearrangement with a broader substrate spectrum or specific stereopreference.

Mechanistic Perspective on C-N and C-S Bond Construction Catalyzed by Cytochrome P450 Enzymes

Cytochrome P450 enzymes catalyze a large number of oxidative transformations that are responsible for natural product synthesis. Recent studies have revealed their unique ability to catalyze the formation of C-N and C-S bonds, broadening their biosynthetic applications. However, the enzymatic mechanisms behind these reactions are still unclear. This review focuses on theoretical insights into the mechanisms of P450-catalyzed C-N and C-S bond formation. The key roles of the conformational dynamics of substrate radicals, influenced by the enzyme environment, in modulating selectivity and reactivity are highlighted. Understanding these reaction mechanisms offers valuable guidance for P450 enzyme engineering and the design of biosynthetic applications.

Biosynthesis of eriodictyol in citrus waster by endowing P450BM3 activity of naringenin hydroxylation

The flavonoid naringenin is abundantly present in pomelo peels, and the unprocessed naringenin in wastes is not friendly for the environment once discarded directly. Fortunately, the hydroxylated product of eriodictyol from naringenin exhibits remarkable antioxidant and anticancer properties. The P450s was suggested promising for the bioconversion of the flavonoids, but less naturally existed P450s show hydroxylation activity to C3' of the naringenin. By well analyzing the catalytic mechanism and the conformations of the naringenin in P450, we proposed that the intermediate Cmpd I ((porphyrin)Fe = O) is more reasonable as key conformation for the hydrolyzation, and the distance between C3'/C5' of naringenin to the O atom of CmpdI determines the hydroxylating activity for the naringenin. Thus, the "flying kite model" that gradually drags the C-H bond of the substrate to the O atom of CmpdI was put forward for rational design. With ab initio design, we successfully endowed the self-sufficient P450-BM3 hydroxylic activity to naringenin and obtained mutant M5-5, with kcat, Km, and kcat/Km values of 230.45 min-1, 310.48 µM, and 0.742 min-1 µM-1, respectively. Furthermore, the mutant M4186 was screened with kcat/Km of 4.28-fold highly improved than the reported M13. The M4186 also exhibited 62.57% yield of eriodictyol, more suitable for the industrial application. This study provided a theoretical guide for the rational design of P450s to the nonnative compounds. KEY POINTS: •The compound I is proposed as the starting point for the rational design of the P450BM3 •"Flying kite model" is proposed based on the distance between O of Cmpd I and C3'/C5' of naringenin •Mutant M15-5 with 1.6-fold of activity than M13 was obtained by ab initio modification.

Progesterone Metabolism in Digitalis and Other Plants-60 Years of Research and Recent Results

5β-Cardenolides are pharmaceutically important metabolites from the specialized metabolism of Digitalis lanata. They were used over decades to treat cardiac insufficiency and supraventricular tachycardia. Since the 1960s, plant scientists have known that progesterone is an essential precursor of cardenolide formation. Therefore, biosynthesis of plant progesterone was mainly analyzed in species of the cardenolide-containing genus Digitalis during the following decades. Today, Digitalis enzymes catalyzing the main steps of progesterone biosynthesis are known. Most of them are found in a broad range of organisms. This review will summarize the findings of 60 years of research on plant progesterone metabolism with particular focus on the recent results in Digitalis lanata and other plants.

Updates on Mechanisms of Cytochrome P450 Catalysis of Complex Steroid Oxidations

Cytochrome P450 (P450) enzymes dominate steroid metabolism. In general, the simple C-hydroxylation reactions are mechanistically straightforward and are generally agreed to involve a perferryl oxygen species (formally FeO3+). Several of the steroid transformations are more complex and involve C-C bond scission. We initiated mechanistic studies with several of these (i.e., 11A1, 17A1, 19A1, and 51A1) and have now established that the dominant modes of catalysis for P450s 19A1 and 51A1 involve a ferric peroxide anion (i.e., Fe3+O2¯) instead of a perferryl ion complex (FeO3+), as demonstrated with 18O incorporation studies. P450 17A1 is less clear. The indicated P450 reactions all involve sequential oxidations, and we have explored the processivity of these multi-step reactions. P450 19A1 is distributive, i.e., intermediate products dissociate and reassociate, but P450s 11A1 and 51A1 are highly processive. P450 17A1 shows intermediate processivity, as expected from the release of 17-hydroxysteroids for the biosynthesis of key molecules, and P450 19A1 is very distributive. P450 11B2 catalyzes a processive multi-step oxidation process with the complexity of a chemical closure of an intermediate to a locked lactol form.

Dry molten globule conformational state of CYP11A1 (SCC) regulates the first step of steroidogenesis in the mitochondrial matrix

Multiple metabolic events occur in mitochondria. Mitochondrial protein translocation from the cytoplasm across compartments depends on the amino acid sequence within the precursor. At the mitochondria associated-ER membrane, misfolding of a mitochondrial targeted protein prior to import ablates metabolism. CYP11A1, cytochrome P450 cholesterol side chain cleavage enzyme (SCC), is imported from the cytoplasm to mitochondrial matrix catalyzing cholesterol to pregnenolone, an essential step for metabolic processes and mammalian survival. Multiple steps regulate the availability of an actively folded SCC; however, the mechanism is unknown. We identified that a dry molten globule state of SCC exists in the matrix by capturing intermediate protein folding steps dictated by its C-terminus. The intermediate dry molten globule state in the mitochondrial matrix of living cells is stable with a limited network of interaction and is inactive. The dry molten globule is activated with hydrogen ions availability, triggering cleavage of cholesterol sidechain, and initiating steroidogenesis.

Selective α-Hydroxyketone Formation and Subsequent C-C Bond Cleavage by Cytochrome P450 Monooxygenase Enzymes

The heme enzymes of the cytochrome P450 superfamily (CYPs) catalyze oxidation reactions with a high level of selectivity. Here, the CYP199A4 enzyme from the bacterium Rhodopseudomonas palustris HaA2 is used to catalyze the hydroxylation of carbonyl-containing compounds to generate α-hydroxyketones. Both 4-propionyl- and 4-(2-oxopropyl)-benzoic acids were regioselectively hydroxylated by this enzyme to generate α-hydroxyketone metabolites, 4-(2-hydroxypropanoyl)benzoic acid and 4-(1-hydroxy-2-oxopropyl)benzoic acid, respectively, with high stereoselectivity. Co-crystallization of CYP199A4 with each substrate allowed high-resolution X-ray crystal structures of the enzyme bound with both to be determined. These provide a rationale for biochemical observations related to substrate binding and activity. As these versatile enzymes have a demonstrated ability to support carbon-carbon (C-C) bond cleavage (lyase) reactions on α-hydroxyketones, we assessed if this activity would be catalyzed by wild-type (WT) CYP199A4. Molecular dynamics (MD) simulations predicted the regioselective hydroxylation of each substrate but indicated that the WT enzyme would not be a good catalyst for lyase activity, in agreement with the experimental observations. The MD simulations also suggested the F182L mutant of CYP199A4 would permit closer approach of the substrate to the ferric-peroxo intermediate, enabling the formation of the lyase transition state. Indeed, this variant was observed to catalyze the cleavage reaction. Furthermore, the F182A variant of CYP199A4 was used to catalyze both the hydroxylation and C-C bond cleavage reactions with both 4-propionyl- and 4-(2-oxopropyl)-benzoic acids using hydrogen peroxide as the oxidant. This dual CYP activity is analogous to that supported by the mammalian CYP17A1 enzyme in steroid biosynthesis.

The multistep oxidation of cholesterol to pregnenolone by human cytochrome P450 11A1 is highly processive

Cytochrome P450 (P450, CYP) 11A1 is the classical cholesterol side chain cleavage enzyme (P450scc) that removes six carbons of the side chain, the first and rate-limiting step in the synthesis of all mammalian steroids. The reaction is a 3-step, 6-electron oxidation that proceeds via formation of 22R-hydroxy (OH) and 20R,22R-(OH)2 cholesterol, yielding pregnenolone. We expressed human P450 11A1 in bacteria, purified the enzyme in the absence of nonionic detergents, and assayed pregnenolone formation by HPLC-mass spectrometry of the dansyl hydrazone. The reaction was inhibited by the nonionic detergent Tween 20, and several lipids did not enhance enzymatic activity. The 22R-OH and 20R,22R-(OH)2 cholesterol intermediates were bound to P450 11A1 relatively tightly, as judged by steady-state optical titrations and koff rates. The electron donor adrenodoxin had little effect on binding; the substrate cholesterol showed a ∼5-fold stimulatory effect on the binding of adrenodoxin to P450 11A1. Presteady-state single-turnover kinetic analysis was consistent with a highly processive reaction with rates of intermediate oxidation steps far exceeding dissociation rates for products and substrates. The presteady-state kinetic analysis revealed a second di-OH cholesterol product, separable by HPLC, in addition to 20R,22R-(OH)2 cholesterol, which we characterized as a rotamer that was also converted to pregnenolone at a similar rate. The first oxidation step (at C-22) is the slowest, limiting the overall rate of cleavage. d3-Cholesterol showed no kinetic deuterium isotope effect on C-22, indicating that C-H bond cleavage is not rate-limiting in the first hydroxylation step.

Computational Characterization of the Reactivity of Compound I in Unspecific Peroxygenases

Unspecific peroxygenases (UPOs) are emerging as promising biocatalysts for selective oxyfunctionalization of unactivated C-H bonds. However, their potential in large-scale synthesis is currently constrained by suboptimal chemical selectivity. Improving the selectivity of UPOs requires a deep understanding of the molecular basis of their catalysis. Recent molecular simulations have sought to unravel UPO's selectivity and inform their design principles. However, most of these studies focused on substrate-binding poses. Few researchers have investigated how the reactivity of CpdI, the principal oxidizing intermediate in the catalytic cycle, influences selectivity in a realistic protein environment. Moreover, the influence of protein electrostatics on the reaction kinetics of CpdI has also been largely overlooked. To bridge this gap, we used multiscale simulations to interpret the regio- and enantioselective hydroxylation of the n-heptane substrate catalyzed by Agrocybe aegerita UPO (AaeUPO). We comprehensively characterized the energetics and kinetics of the hydrogen atom-transfer (HAT) step, initiated by CpdI, and the subsequent oxygen rebound step forming the product. Notably, our approach involved both free energy and potential energy evaluations in a quantum mechanics/molecular mechanics (QM/MM) setting, mitigating the dependence of results on the choice of initial conditions. These calculations illuminate the thermodynamics and kinetics of the HAT and oxygen rebound steps. Our findings highlight that both the conformational selection and the distinct chemical reactivity of different substrate hydrogen atoms together dictate the regio- and enantio-selectivity. Building on our previous study of CpdI's formation in AaeUPO, our results indicate that the HAT step is the rate-limiting step in the overall catalytic cycle. The subsequent oxygen rebound step is swift and retains the selectivity determined by the HAT step. We also pinpointed several polar and charged amino acid residues whose electrostatic potentials considerably influence the reaction barrier of the HAT step. Notably, the Glu196 residue is pivotal for both the CpdI's formation and participation in the HAT step. Our research offers in-depth insights into the catalytic cycle of AaeUPO, which will be instrumental in the rational design of UPOs with enhanced properties.

Molecular cloning and characterization of STC1 gene and its functional analyses in yak (Bos grunniens) cumulus granulosa cells

Cumulus granulosa cells (CGCs), an important type of ovarian somatic cells, carries out various functions related to oogenesis, follicular development, and steroidogenesis. Studying the biological mechanisms involved in the development and function of CGCs makes a great contribution to understanding the reproductive regulation in female animals. Stanniocalcin-1 (STC1) is an important Ca2+-regulated glycoprotein hormone that exhibits high expression levels in ovaries. In this study, we cloned the coding sequence of the yak STC1, predicted the structure of STC1 protein, detected the expression and localization of STC1 in yak ovaries, and analyzed the functions of STC1 in yak CGCs. The CDS (coding sequence) region of yak STC1 gene was found to be 744 bp and encoded 247 amino acids. Homology comparison revealed that STC1 protein was highly conserved among mammals. The STC1 mRNA displayed dynamic expression profiles in different stages of yak ovaries, and the highest expression was found in the follicular phase. Regarding localization, STC1 protein was widely distributed in various kinds of yak ovarian cells, including oocytes, mural granulosa cells, CGCs, and thecal cells. Repressing the expression of STC1 resulted in defective proliferation and survival of yak CGCs. In addition, knockdown the expression of STC1 repressed the secretion of progesterone and promoted the secretion of estrogen. Overexpression of STC1 partially rescued the proliferation of CGCs and resulted in opposite effects on the secretion of progesterone and estrogen. Several apoptosis and steroidogenesis-related genes, including BAX, BCL2, HSD3B1, HSD17B1, CYP11A1 and CYP17A1 showed altered expressions after repressing or increasing the expression of STC1 in yak CGCs. To the best of our knowledge, this study is the first to focus on the role of STC1 in yak CGCs, and the outcomes offer fresh insights into the mechanism governing yak reproduction.

Enantioselective metabolism of novel chiral insecticide Paichongding by human cytochrome P450 3A4: A computational insight

As a novel chiral neonicotinoid insecticide, Paichongding (IPP) has been widely applied in agriculture due to its excellent insecticidal activity. However, the enantioselective metabolism of IPP stereoisomers (5R7R-IPP, 5S7S-IPP, 5R7S-IPP, and 5S7R-IPP) mediated by enzymes in non-target organisms, especially the cytochrome P450s (CYPs), remains unknown. To address this knowledge gap, we developed an integrated computational framework to elucidate the binding interactions and enantioselective metabolism of IPP stereoisomers in human CYP3A4. The results reveal that 5R7R-IPP shows much stronger binding affinity to CYP3A4 than 5S7S-IPP, while enantiomers 5R7S-IPP and 5S7R-IPP have no essential difference in their binding potential, owing to their specific interactions with key CYP3A4 residues. Although enantiomers 5R7R-IPP and 5S7S-IPP feature distinct binding modes resulting from the chiral differences, their transformation activities are slightly different, with C5 and C13 being the primary metabolic sites, respectively. In contrast, CYP3A4 preferably metabolizes 5R7S-IPP over 5S7R-IPP. The metabolism of epimers 5R7R-IPP and 5R7S-IPP share C5-hydroxylation routes due to the conserved 5R-conformaitons, but differ with the transformation routes at C11/C13 and C3 sites. The 7R-chirality of 5S7R-IPP significantly reduces the metabolic potency compared to 5S7S-IPP. CYP3A4-catalyzed hydroxylation and desaturation of IPP stereoisomers generate various chiral metabolites, with C5- and C13-hydroxyIPPs further transforming into depropylated products. Furthermore, the toxicity assessment reveals that IPP, along with the majority of its hydroxylated, desaturated, and depropylated metabolites, can potentially induce adverse effects on human health, specifically hepatotoxicity, respiratory toxicity, and carcinogenicity. This study provides valuable insights into the enantioselective fate of chiral IPP metabolism by CYP3A4, and the identified metabolites can serve as potential biomarkers for monitoring IPP exposure and associated health risk in human body.

QM/MM Modeling Aided Enzyme Engineering in Natural Products Biosynthesis

Natural products and their derivatives are widely used across various industries, particularly pharmaceuticals. Modern engineered biosynthesis provides an alternative way of producing and meeting the growing need for diverse natural products. Natural enzymes, on the other hand, often exhibit unsatisfactory catalytic characteristics and necessitate further enzyme engineering modifications. QM/MM, as a powerful and extensively used computational tool in the field of enzyme catalysis, has been increasingly applied in rational enzyme engineering over the past decade. In this review, we summarize recent advances in QM/MM computational investigation on enzyme catalysis and enzyme engineering for natural product biosynthesis. The challenges and perspectives for future QM/MM applications aided enzyme engineering in natural product biosynthesis will also be discussed.

Multiscale QM/MM modelling of catalytic systems with ChemShell

Hybrid quantum mechanical/molecular mechanical (QM/MM) methods are a powerful computational tool for the investigation of all forms of catalysis, as they allow for an accurate description of reactions occurring at catalytic sites in the context of a complicated electrostatic environment. The scriptable computational chemistry environment ChemShell is a leading software package for QM/MM calculations, providing a flexible, high performance framework for modelling both biomolecular and materials catalysis. We present an overview of recent applications of ChemShell to problems in catalysis and review new functionality introduced into the redeveloped Python-based version of ChemShell to support catalytic modelling. These include a fully guided workflow for biomolecular QM/MM modelling, starting from an experimental structure, a periodic QM/MM embedding scheme to support modelling of metallic materials, and a comprehensive set of tutorials for biomolecular and materials modelling.

Direct Electrophilic Attack of Compound I on the Indole Ring in the Peroxygenase Mechanism of Dehaloperoxidase DHP B in Degrading Haloindole: Electron Transfer Promotes the Reaction

The H2O2-dependent degradation of haloindole catalyzed by the dehaloperoxidase (DHP) from Amphitrite ornate has been reported to employ the peroxygenase mechanism, and the two oxidized products 5-halo-2-oxindole and 5-halo-3-oxindole have a similar amount. According to a previous experimental study, compound I (Cpd I) was suggested to be responsible for triggering the reaction, and the reaction may undergo three possible intermediates; however, the reaction details are still unclear. To clarify the reaction mechanism of DHP, the computational model was constructed on the basis of the high-resolution crystal structure, and a series of the quantum mechanical/molecular mechanical calculations were performed. Based on our calculation results, it is confirmed that the reaction starts from the direct electrophilic attack of Cpd I on the indole ring of the substrate, and the resulted intermediate contains both a carbocation and an oxygen anion, whereas the common hydrogen abstraction by Cpd I was calculated to correspond to a relatively higher barrier. In addition, a net electron transfer from the substrate to the iron center is observed during the attack of Cpd I on the indole ring; therefore, the carbocation/oxygen anion intermediate can easily undergo an intramolecular hydride transfer to form the product 5-halo-2-oxindole or isomerize to the epoxide intermediate which finally generates another product 5-halo-3-oxindole. It is the zwitterionic characteristic of the intermediate that makes the intermolecular hydride transfer quite easy, and it is the high electron affinity of the iron center that promotes the single-electron oxidation of the reaction intermediate. Our calculations well explain the formation of two oxidized products 5-halo-2-oxindole and 5-halo-3-oxindole.

Comparative biochemical characterization of mammalian-derived CYP11A1s with cholesterol side-chain cleavage activities

Steroid drugs, the second largest class of pharmaceuticals after antibiotics, have shown significant anti-inflammatory, anti-allergic, and endocrine-regulating effects. A group of cytochrome P450 enzymes, namely, CYP11A1 isoenzymes from different organisms are capable of converting cholesterol into pregnenolone, which is a pivotal reaction in both steroid metabolism and (bio)synthetic network of steroid products. However, the low activity of CYP11A1s greatly restricts the industrial application of these cholesterol side-chain cleavage enzymes. Herein, we investigate ten CYP11A1 enzymes of different origins and in vitro characterize two CYP11A1s with a relatively higher expression level from Capra hircus and Sus scrofa, together with the CYP11A1s from Homo sapiens and Bos taurus as references. Towards five selected sterol substrates with different side chain structures, S. scrofa CYP11A1 displays relatively higher activities. Through redox partners combination screening, we reveal the optimal redox partner pair of S. scrofa adrenodoxin and C. hircus adrenodoxin reductase. Moreover, the semi-rational mutagenesis for the active sites and substrate entrance channels of human and bovine CYP11A1s is performed based on comparative analysis of their crystal structures. The mutant mBtCYP11A1-Q377A derived from mature B. taurus CYP11A1 shows a 1.46 times higher activity than the wild type enzyme. These results not only demonstrate the tunability of the highly conserved CYP11A1 isoenzymes, but also lay a foundation for the following engineering efforts on these industrially relevant P450 enzymes.

Engineering C-C Bond Cleavage Activity into a P450 Monooxygenase Enzyme

The cytochrome P450 (CYP) superfamily of heme monooxygenases has demonstrated ability to facilitate hydroxylation, desaturation, sulfoxidation, epoxidation, heteroatom dealkylation, and carbon-carbon bond formation and cleavage (lyase) reactions. Seeking to study the carbon-carbon cleavage reaction of α-hydroxy ketones in mechanistic detail using a microbial P450, we synthesized α-hydroxy ketone probes based on the physiological substrate for a well-characterized benzoic acid metabolizing P450, CYP199A4. After observing low activity with wild-type CYP199A4, subsequent assays with an F182L mutant demonstrated enzyme-dependent C-C bond cleavage toward one of the α-hydroxy ketones. This C-C cleavage reaction was subject to an inverse kinetic solvent isotope effect analogous to that observed in the lyase activity of the human P450 CYP17A1, suggesting the involvement of a species earlier than Compound I in the catalytic cycle. Co-crystallization of F182L-CYP199A4 with this α-hydroxy ketone showed that the substrate bound in the active site with a preference for the (S)-enantiomer in a position which could mimic the topology of the lyase reaction in CYP17A1. Molecular dynamics simulations with an oxy-ferrous model of CYP199A4 revealed a displacement of the substrate to allow for oxygen binding and the formation of the lyase transition state proposed for CYP17A1. This demonstration that a correctly positioned α-hydroxy ketone substrate can realize lyase activity with an unusual inverse solvent isotope effect in an engineered microbial system opens the door for further detailed biophysical and structural characterization of CYP catalytic intermediates.

How the Conformational Movement of the Substrate Drives the Regioselective C-N Bond Formation in P450 TleB: Insights from Molecular Dynamics Simulations and Quantum Mechanical/Molecular Mechanical Calculations

P450 TleB catalyzes the oxidative cyclization of the dipeptide N-methylvalyl-tryptophanol into indolactam V through selective intramolecular C-H bond amination at the indole C4 position. Understanding its catalytic mechanism is instrumental for the engineering or design of P450-catalyzed C-H amination reactions. Using multiscale computational methods, we show that the reaction proceeds through a diradical pathway, involving a hydrogen atom transfer (HAT) from N1-H to Cpd I, a conformational transformation of the substrate radical species, and a second HAT from N13-H to Cpd II. Intriguingly, the conformational transformation is found to be the key to enabling efficient and selective C-N coupling between N13 and C4 in the subsequent diradical coupling reaction. The underlined conformational transformation is triggered by the first HAT, which proceeds with an energy-demanding indole ring flip and is followed by the facile approach of the N13-H group to Cpd II. Detailed analysis shows that the internal electric field (IEF) from the protein environment plays key roles in the transformation process, which not only provides the driving force but also stabilizes the flipped conformation of the indole radical. Our simulations provide a clear picture of how the P450 enzyme can smartly modulate the selective C-N coupling reaction. The present findings are in line with all available experimental data, highlighting the crucial role of substrate dynamics in controlling this highly valuable reaction.

Computational Study of the Peroxygenase Mechanism Catalyzed by Hemoglobin Dehaloperoxidase Involved in the Degradation of Chlorophenols

The biochemical evidence showed that hemoglobin dehaloperoxidase (DHP B) from Amphitrite Ornata is a multifunctional hemoprotein that catalyzes both dehalogenation and hydroxylation of halophenols via the peroxidase and peroxygenase mechanism, respectively, which sets the basis for the degradation of halophenols. In the peroxygenase mechanism, the reaction was previously suggested to be triggered either by the hydrogen atom abstraction by the Fe═O center or by the proton abstraction by His55. To illuminate the peroxygenase mechanism of DHP B at the atomistic level, on the basis of the high-resolution crystal structure, computational models were constructed, and a series of quantum mechanical/molecular mechanical calculations have been performed. According to the calculation results, the pathway (Path a) initiated by the H-abstraction by the Fe═O center is feasible. In another pathway (Path b), His55 can abstract the proton from the hydroxyl group of the substrate (4-Cl-o-cresol) to initiate the reaction; however, its feasibility depends on the prior electron transfer from the substrate to the porphyrin group. The rate-limiting step of Path a is the OH-rebound, which corresponds to an energy barrier of 14.7 kcal/mol at the quartet state. His55 acts as an acid-base catalyst and directly involves in the catalysis. Our mutant study indicates that His55 can be replaced by other titratable residues. These findings may provide useful information for further understanding of the catalytic reaction of DHP B and for the design of enzymes in the degradation of pollutants, in particular, halophenols.

Mechanism of Sugar Ring Contraction and Closure Catalyzed by UDP-d-apiose/UDP-d-xylose Synthase (UAXS)

Uridine diphosphate (UDP)-apiose/UDP-xylose synthase (UAXS) is a member of the short-chain dehydrogenase/reductase superfamily (SDR), which catalyzes the ring contraction and closure of UDP-d-glucuronic acid (UDP-GlcA), affording UDP-apiose and UDP-xylose. UAXS is a special enzyme that integrates ring-opening, decarboxylation, rearrangement, and ring closure/contraction in a single active site. Recently, the ternary complex structure of UAXS was crystallized from Arabidopsis thaliana. In this work, to gain insights into the detailed formation mechanism of UDP-apiose and UDP-xylose, an enzyme-substrate reactant model has been constructed and quantum mechanical/molecular mechanical (QM/MM) calculations have been performed. Our calculation results reveal that the reaction starts from the C4-OH oxidation, which is accompanied by the conformational transformation of the sugar ring from chair type to boat type. The sugar ring-opening is prior to decarboxylation, and the deprotonation of the C2-OH group is the prerequisite for sugar ring-opening. Moreover, the keto-enol tautomerization of the decarboxylated intermediate is a necessary step for ring closure/contraction. Based on our calculation results, more UDP-apiose product was expected, which is in line with the experimental observation. Three titratable residues, Tyr185, Cys100, and Cys140, steer the reaction by proton transfer from or to UDP-GlcA, and Arg182, Glu141, and D337 constitute a proton conduit for sugar C2-OH deprotonation. Although Thr139 and Tyr105 are not directly involved in the enzymatic reaction, they are responsible for promoting the catalysis by forming hydrogen-bonding interactions with GlcA. Our calculations may provide useful information for understanding the catalysis of the SDR family.

QM/MM and MM MD simulations on decontamination of the V-type nerve agent VX by phosphotriesterase: Toward a comprehensive understanding of steroselectivity and activity

Due to deadly toxicity and high environmental stability of the nerve agent VX, an efficient decontamination approach is desperately needed in tackling its severe threat to human secu-rity. The enzymatic...

Unexpected Reactions of α,β-Unsaturated Fatty Acids Provide Insight into the Mechanisms of CYP152 Peroxygenases

CYP152 peroxygenases catalyze decarboxylation and hydroxylation of fatty acids using H₂ O₂ as cofactor. To understand the molecular basis for the chemo- and regioselectivity of these unique P450 enzymes, we analyze the activities of three CYP152 peroxygenases (OleT_JE , P450_SPα , P450_BSβ ) towards cis- and trans-dodecenoic acids as substrate probes. The unexpected 6S-hydroxylation of the trans-isomer and 4R-hydroxylation of the cis-isomer by OleT_JE , and molecular docking results suggest that the unprecedented selectivity is due to OleT_JE 's preference of C2-C3 cis-configuration. In addition to the common epoxide products, undecanal is the unexpected major product of P450_SPα and P450_BSβ regardless of the cis/trans-configuration of substrates. The combined H₂ ¹⁸ O₂ tracing experiments, MD simulations, and QM/MM calculations unravel an unusual mechanism for Compound I-mediated aldehyde formation in which the active site water derived from H₂ O₂ activation is involved in the generation of a four-membered ring lactone intermediate. These findings provide new insights into the unusual mechanisms of CYP152 peroxygenases.

Theoretical Investigation on the Elusive Reaction Mechanism of Spirooxindole Formation Mediated by Cytochrome P450s: A Nascent Feasible Charge-Shift C-O Bond Makes a Difference

Spirooxindoles are pivotal biofunctional groups widely distributed in natural products and clinic drugs. However, construction of such subtle chiral skeletons is a long-standing challenge to both organic and bioengineering scientists. The knowledge of enzymatic spirooxindole formation in nature may inspire rational design of new catalysts. To this end, we presented a theoretical investigation on the elusive mechanism of the spiro-ring formation at the 3-position of oxindole mediated by cytochrome P450 enzymes (P450). Our calculated results demonstrated that the electrophilic attack of CpdI, the active species of P450, to the substrate, shows regioselectivity, i.e., the attack at the C9 position forms a tetrahedral intermediate involving an unusual feasible charge-shift C9^δ+-O^δ- bond, while the attack at the C1 position forms an epoxide intermediate. The predominant route is the first route with the charge-shift bonding intermediate due to holding a relatively lower barrier by >5 kcal mol^-1 than the epoxide route, which fits the experimental observations. Such a delocalized charge-shift bond facilitates the formation of a spiro-ring mainly through elongation of the C1-C9 bond to eliminate the aromatization of the tricyclic beta-carboline. Our theoretical results shed profound mechanistic insights for the first time into the elusive spirooxindole formation mediated by P450s.

Multiscale Simulations on the Catalytic Plasticity of CYP76AH1

The catalytic promiscuity and fidelity of cytochrome P450 enzymes are widespread in the skeletal modification of terpenoid natural products and have attracted much attention. CYP76AH1 is involved in key modification reactions in the biosynthetic pathway of tanshinone, a well-known medicinal norditerpenoid. In this work, classical molecular dynamic simulations, metadynamics, and DFT calculations were performed to investigate the protein conformational dynamics, ligand binding poses, and catalytic reaction mechanism in wide-type and mutant CYP76AH1. Our results not only reveal a plausible enzymatic mechanism for mutant CYP76AH1 leading to various products but also provide valuable guidance for rational protein engineering of the CYP76 family.

Steroidogenic cytochrome P450 17A1 structure and function

Cytochrome P450 17A1 (CYP17A1) is a critical steroidogenic enzyme, essential for producing glucocorticoids and sex hormones. This review discusses the complex activity of CYP17A1, looking at its role in both the classical and backdoor steroidogenic pathways and the complex chemistry it carries out to perform both a hydroxylation reaction and a carbon-carbon cleavage, or lyase reaction. Functional and structural investigations have informed our knowledge of these two reactions. This review focuses on a few specific aspects of this discussion: the identities of reaction intermediates, the coordination of hydroxylation and lyase reactions, the effects of cytochrome b5, and conformational selection. These discussions improve understanding of CYP17A1 in a physiological setting, where CYP17A1 is implicated in a variety of steroidogenic diseases. This information can be used to improve ways in which CYP17A1 can be effectively modulated to treat diseases such as prostate and breast cancer, Cushing's syndrome, and glioblastoma.

Computing Metal-Binding Proteins for Therapeutic Benefit

Over one third of biomolecules rely on metal ions to exert their cellular functions. Metal ions can play a structural role by stabilizing the structure of biomolecules, a functional role by promoting a wide variety of biochemical reactions, and a regulatory role by acting as messengers upon binding to proteins regulating cellular metal-homeostasis. These diverse roles in biology ascribe critical implications to metal-binding proteins in the onset of many diseases. Hence, it is of utmost importance to exhaustively unlock the different mechanistic facets of metal-binding proteins and to harness this knowledge to rationally devise novel therapeutic strategies to prevent or cure pathological states associated with metal-dependent cellular dysfunctions. In this compendium, we illustrate how the use of a computational arsenal based on docking, classical, and quantum-classical molecular dynamics simulations can contribute to extricate the minutiae of the catalytic, transport, and inhibition mechanisms of metal-binding proteins at the atomic level. This knowledge represents a fertile ground and an essential prerequisite for selectively targeting metal-binding proteins with small-molecule inhibitors aiming to (i) abrogate deregulated metal-dependent (mis)functions or (ii) leverage metal-dyshomeostasis to selectively trigger harmful cells death.

Conformational Motion of Ferredoxin Enables Efficient Electron Transfer to Heme in the Full-Length P450TT

Cytochrome P450 monooxygenases (P450s) are versatile biocatalysts used in natural products biosynthesis, xenobiotic metabolisms, and biotechnologies. In P450s, the electrons required for O₂ activation are supplied by NAD(P)H through stepwise electron transfers (ETs) mediated by redox partners. While much is known about the machinery of the catalytic cycle of P450s, the mechanisms of long-range ET are largely unknown. Very recently, the first crystal structure of full-length P450_TT was solved. This enables us to decipher the interdomain ET mechanism between the [2Fe-2S]-containing ferredoxin and the heme, by use of molecular dynamics simulations. In contrast to the "distal" conformation characterized in the crystal structure where the [2Fe-2S] cluster is ∼28 Å away from heme-Fe, our simulations demonstrated a "proximal" conformation of [2Fe-2S] that is ∼17 Å [and 13.7 Å edge-to-edge] away from heme-Fe, which may enable the interdomain ET. Key residues involved in ET pathways and interdomain complexation were identified, some of which have already been verified by recent mutation studies. The conformational transit of ferredoxin between "distal" and "proximal" was found to be controlled mostly by the long-range electrostatic interactions between the ferredoxin domain and the other two domains. Furthermore, our simulations show that the full-length P450_TT utilizes a flexible ET pathway that resembles either P450_Scc or P450_cam. Thus, this study provides a uniform picture of the ET process between reductase domains and heme domain.

Frontiers of metal-coordinating drug design

Introduction: The occurrence of metal ions in biomolecules is required to exert vital cellular functions. Metal-containing biomolecules can be modulated by small-molecule inhibitors targeting their metal-moiety. As well, the discovery of cisplatin ushered the rational discovery of metal-containing-drugs. The use of both drug types exploiting metal-ligand interactions is well established to treat distinct pathologies. Therefore, characterizing and leveraging metal-coordinating drugs is a pivotal, yet challenging, part of medicinal chemistry.Area covered: Atomic-level simulations are increasingly employed to overcome the challenges met by traditional drug-discovery approaches and to complement wet-lab experiments in elucidating the mechanisms of drugs' action. Multiscale simulations, allow deciphering the mechanism of metal-binding inhibitors and metallo-containing-drugs, enabling a reliable description of metal-complexes in their biological environment. In this compendium, the authors review selected applications exploiting the metal-ligand interactions by focusing on understanding the mechanism and design of (i) inhibitors targeting iron and zinc-enzymes, and (ii) ruthenium and gold-based anticancer agents targeting the nucleosome and aquaporin protein, respectively.Expert opinion: The showcased applications exemplify the current role and the potential of atomic-level simulations and reveal how their synergic use with experiments can contribute to uncover fundamental mechanistic facets and exploit metal-ligand interactions in medicinal chemistry.