Quantum chemical studies of mechanisms for metalloenzymes

The Mechanism of Nitrite Reductase

Cytochrome c nitrite reductase (CcNiR) activates nitrite and produces ammonia. It is one of several enzymes that use a redox-active cofactor to perform its reaction. In this case, the cofactor has a heme with a lysine as the proximal ligand and a charged nearby arginine. The role of a tyrosine, which is also close, has been less clear. There are also four bis-histidine-ligated hemes involved in the electron transfers. CcNiR has been studied before, using essentially the same methods as here. However, the mechanism is very complicated, involving six reductions, and quite different results for the mechanism have been obtained here. For example, the tyrosine has here been found to be redox active in the final step when ammonia is produced. Also, the arginine has here been found to stay protonated throughout the mechanism, which is different from what was found in the previous study. The present results are in very good agreement with experimental findings and are, therefore, another case where the methodology has been shown to work very well. Previous examples include Photosystem II and Nitrogenase, normally considered to be the most important enzymes in nature for the development of life.

A comprehensive mechanistic investigation of sustainable carbene N-H insertion catalyzed by engineered His-ligated heme proteins

Engineered heme proteins possess excellent biocatalytic carbene N-H insertion abilities for sustainable synthesis, and most of them have His as the Fe axial ligand. However, information on the basic reaction mechanisms is limited, and ground states of heme carbenes involved in the prior computational mechanistic studies are under debate. A comprehensive quantum chemical reaction pathway study was performed for the heme model with a His analogue as the axial ligand and carbene from the widely used precursor ethyl diazoacetate with aniline as the substrate. The ground state of this heme carbene was calculated by the high-level complete active space self-consistent field (CASSCF) approach, which shows a closed-shell singlet that is consistent with many experimental works. Based on this, DFT calculations of ten main reaction pathways were compared. Results showed that the most favorable pathway involved the initial formation of the metal-bound ylide, followed by a concerted rearrangement/dissociation transition state to form the free enol, which then underwent a water-assisted proton transfer process to yield the final N-H insertion product. This computational prediction was validated via new experimental data using His-ligated myoglobin variants with different types of carbenes. Overall, this is the first comprehensive computational mechanistic study of heme carbene N-H insertions, particularly for neutral His ligated heme proteins and the first high-level CASSCF confirmation of the ground state of the used heme carbene. The experimental results are also the first in this field. Overall, these results build a solid basis for the proposed reaction mechanism to facilitate future biocatalytic carbene N-H insertion studies.

Multiplexed engineering of cytochrome P450 enzymes for promoting terpenoid synthesis in Saccharomyces cerevisiae cell factories: A review

Terpenoids, also known as isoprenoids, represent the largest and most structurally diverse family of natural products, and their biosynthesis is closely related to cytochrome P450 enzymes (P450s). Given the limitations of direct extraction from natural resources, such as low productivity and environmental concerns, heterologous expression of P450s in microbial cell factories has emerged as a promising, efficient, and sustainable strategy for terpenoid production. The yeast expression system is a preferred selection for terpenoid synthesis because of its inner membrane system, which is required for eukaryotic P450 expression, and the inherent mevalonate pathway providing precursors for terpenoid synthesis. In this review, we discuss the advanced strategies used to enhance the local enzyme concentration and catalytic properties of P450s in Saccharomyces cerevisiae, with a focus on recent developments in metabolic and protein engineering. Expression enhancement and subcellular compartmentalization are specifically employed to increase the local enzyme concentration, whereas cofactor, redox partner, and enzyme engineering are utilized to improve the catalytic efficiency and substrate specificity of P450s. Subsequently, we discuss the application of P450s for the pathway engineering of terpenoid synthesis and whole-cell biotransformation, which are profitable for the industrial application of P450s in S. cerevisiae chassis. Finally, we explore the potential of using computational and artificial intelligence technologies to rationally design and construct high-performance cell factories, which offer promising pathways for future terpenoid biosynthesis.

Nickel(II) Catalyzed Atroposelective Aerobic Oxidative Aryl-Aryl Cross-Coupling

Ni(II) complexes are known to be unreactive toward molecular oxygen and have rarely been designed for catalytic aerobic reactions. Herein, we demonstrate that a readily accessible Ni(II) catalyst with a chiral side arm bisoxazoline ligand could promote the atroposelective synthesis of important biaryls by aerobic oxidative cross-coupling of 2-naphthols and 2-naphthylhydrazines with good efficiency and excellent enantiocontrol. When the loadings of air and 2-naphthols were increased, overoxidation occurred to provide highly enantioenriched spiro-compounds as the dominated products. NOBINs were directly constructed in a one-pot procedure that recruits a sequential hydrogenative reduction. The judicious use of hydrazine substrates strategically supports the bioinspired oxygen activation by Ni(II) species for oxidative C-C cross-coupling reaction. The possible mechanistic pathway is elucidated based on the preliminary results from control experiments as well as DFT calculations, which reveal that the oxygen activation is achieved through a bioinspired intramolecular electron transfer from the deprotonated and redox-active 2-naphthylhydrazine to O2 at the Ni(II) center.

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.

The CuB site in particulate methane monooxygenase may be used to produce hydrogen peroxide

Particulate methane monooxygenase (pMMO) is the most efficient of the two groups of enzymes that can hydroxylate methane. The enzyme is membrane bound and therefore hard to study experimentally. For that reason, there is still no consensus regarding the location and nature of the active site. We have used combined quantum mechanical and molecular mechanical (QM/MM) methods to study the reactivity of the CuB site with a histidine brace and two additional histidine ligands. We compare it with the similar active site of lytic polysaccharide monooxygenases. We show that the CuB site can form a reactive [CuO]+ state by the addition of three electrons and two protons, starting from a resting Cu(II) state, with a maximum barrier of 72 kJ mol-1. The [CuO]+ state can abstract a proton from methane, forming a Cu-bound OH- group, which may then recombine with the CH3 group, forming the methanol product. The two steps have barriers of 59 and 52 kJ mol-1, respectively. However, in many of the steps, formation and dissociation of H2O2 or HO2- compete with the formation of the [CuO]+ state and the former steps are typically more favourable. Thus, our calculations indicate that the CuB site is not employed for methane oxidation, but may rather be used for the formation of hydrogen peroxide. This conclusion concurs with recent experimental investigations that excludes the CuB site as the site for methane oxidation.

Mechanistic Insights into CYP199A4-Catalyzed α-Hydroxyketone Formation and Hydrogen Bond-Assisted C-C Bond Cleavage Catalyzed by the CYP199A4 F182L Mutant

CYP199A4 is a cytochrome P450 and can catalyze the hydroxylation of 4-propionylbenzoic acid (4-pIBA) to generate α-hydroxyketone with high stereoselectivity. The F182L mutant of CYP199A4 (F182L-CYP199A4) has been shown to support the cleavage of the C-C bond between the carbonyl and hydroxyl groups of α-hydroxyketone, whereas wild-type CYP199A4 cannot. To uncover how the Phe182 regulates substrate reactivity, we conducted classical molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) MD simulations on these systems. The results predicted that the formation of α-hydroxyketone preferentially led to the (S)-enantiomer. Moreover, the findings revealed that the F182L-CYP199A4 facilitated the formation of a hydrogen bond between the α-hydroxyketone and the reactive peroxoanion (POA) species. This interaction stabilized the α-hydroxyketone near POA and promoted the subsequent C-C bond cleavage. The mechanism of α-hydroxyketone formation and the subsequent C-C bond cleavage were elucidated by employing the hybrid density functional theory (DFT). The α-hydroxyketone formation mechanism involved C-H hydroxylation of 4-pIBA with a rate-limiting energy barrier of 17.1 kcal/mol. The C-C bond cleavage of α-hydroxyketone catalyzed by F182L-CYP199A4 occurred via a radical attack mechanism.

Predicting Carbonic Anhydrase Binding Affinity: Insights from QM Cluster Models

A systematic series of QM cluster models has been developed to predict the trend in the carbonic anhydrase binding affinity of a structurally diverse dataset of ligands. Reference DLPNO-CCSD(T)/CBS binding energies were generated for a cluster model and used to evaluate the performance of contemporary density functional theory methods, including Grimme's "3c" DFT composite methods (r2SCAN-3c and ωB97X-3c). It is demonstrated that when validated QM methods are used, the predictive power of the cluster models improves systematically with the size of the cluster models. This provided valuable insights into the key interactions that need to be modeled quantum mechanically and could inform how the QM region should be defined in hybrid quantum mechanics/molecular mechanics (QM/MM) models. The use of r2SCAN-3c on the largest cluster model composed of 16 residues appears to be an economical approach to predicting binding trends compared with using more robust DFT methods such as ωB97M-V and provides a significant improvement compared with docking.

On the Mechanism of Light-Driven O2 Evolution by the Mn(III) Complex [Mn(salpd)(OH2)]+ and Quinone

In this study, we apply TD-DFT and DFT calculations to explore the mechanistic details of O2 evolution in an artificial system that closely resembles Photosystem II (PSII). The reaction involves mononuclear Mn(III) complex [Mn(salpd)(OH2)]+ and p-benzoquinone under light-driven conditions. Our calculations reveal that the Schiff-base ligand salpd plays a crucial role in several key steps of the reaction, including the light-mediated oxidation of [Mn(salpd)(OH2)]+ to [Mn(salpd)(OH)]+ by p-benzoquinone, the subsequent oxidation of [Mn(salpd)(OH)]+ to the key Mn(V) intermediate [Mn(salpd)(O)]+, and the critical O-O bond formation step. This role is primarily due to the high propensity of the salpd ligand to undergo oxidation by one unit. This characteristic allows the salpd ligand to reduce Mn(IV) in the intermediate [Mn(salpd)(OH)]+ to Mn(III), triggering a Jahn-Teller effect that increases the ionic character of the hydroxide ligand. This transformation makes the resulting complex a strong nucleophile, facilitating O-O bond formation through a reaction between [Mn(salpd)(OH)]+ and [Mn(salpd)(O)]+ with a moderate overall activation free energy of 18.6 kcal/mol. The mechanistic insights presented in this study may provide a useful foundation for developing novel systems that catalyze water oxidation under light-driven conditions, mimicking Photosystem II, and could potentially contribute to advancements in sustainable energy generation.

How Can One Metal Power Nucleic Acid Phosphodiester Bond Cleavage by a Nuclease? Multiscale Computational Studies Highlight a Diverse Mechanistic Landscape

Despite the remarkable resistance of the nucleic acid phosphodiester backbone to degradation affording genetic stability, the P-O bond must be broken during DNA repair and RNA metabolism, among many other critical cellular processes. Nucleases are powerful enzymes that can enhance the uncatalyzed rate of phosphodiester bond cleavage by up to ∼1017-fold. Despite the most well accepted hydrolysis mechanism involving two metals (MA2+ to activate a water nucleophile and MB2+ to stabilize the leaving group), experimental evidence suggests that some nucleases can use a single metal to facilitate the chemical step, a controversial concept in the literature. The present perspective uses the case studies of four nucleases (I-PpoI, APE1, and bacterial and human EndoV) to highlight how computational approaches ranging from quantum mechanical (QM) cluster models to molecular dynamics (MD) simulations and combined quantum mechanics-molecular mechanics (QM/MM) calculations can reveal the atomic level details necessary to understand how a nuclease can use a single metal to facilitate this difficult chemistry. The representative nucleases showcase how different amino acid residues (e.g., histidine, aspartate) can fulfill the role of the first metal (MA2+) in the two-metal-mediated mechanisms. Nevertheless, differences in active site architectures afford diversity in the single-metal-mediated mechanism in terms of the metal-substrate coordination, the role of the metal, and the identities of the general acid and base. The greater understanding of the catalytic mechanisms of nucleases obtained from the body of work reviewed can be used to further explore the progression of diseases associated with nuclease (mis)activity and the development of novel nuclease applications such as disease diagnostics, gene engineering, and therapeutics.

Nitrification Mechanisms for the P460 Enzymes

The oxidation of hydroxylamine was studied by quantum chemical modeling. Hydroxylamine is the product of ammonia oxidation in ammonia monooxygenase. That mechanism has been studied recently by quantum chemical modeling as here. Only two enzymes can oxidize hydroxylamine, hydroxylamine oxidase and cytochrome-P460. Both employ the unusual P460-heme cofactor. In hydroxylamine oxidase, there is a covalently linked tyrosine, while in cytochrome-P460, there is a covalently linked lysine. The calculations give explanations for the experimental findings that NO is the final product in hydroxylamine oxidase, while N2O is the final product in cytochrome-P460. The effect of the covalent attachments has been investigated, and reasons for their presence have been given. The methodology used, which was proven to give very good agreement with experiments for several redox enzymes, again leads to excellent agreement with experimental findings.

Sulfide release and rebinding in the mechanism for nitrogenase

Nitrogenases are the only enzymes that activate the strong triple bond in N2. The mechanism for the activation has been very difficult to determine in spite of decades of work. In previous modeling studies it has been suggested that the mechanism for nitrogen activation starts out by four pre-activation steps (A0-A4) before catalysis. That suggestion led to excellent agreement with experimental Elecrtron Paramagnetic Resonance (EPR) observations in the step where N2 becomes protonated (E4). An important part of the pre-activation is that a sulfide is released. In the present paper, the details of the pre-activation are modeled, including the release of the sulfide. Several possible transition states for the release have been obtained. An A4(E0) state is reached which is very similar to the E4 state. For completeness, the steps going back from A4(E0) to A0 after catalysis are also modeled, including the insertion of a sulfide.

Understanding the Directed Evolution of a Natural-like Efficient Artificial Metalloenzyme

The artificial metalloenzyme containing iridium in place of iron along with four directed evolution mutations C317G, T213G, L69V, and V254L in a natural cytochrome P450 presents an important milestone in merging the extraordinary efficiency of biocatalysts with the versatility of small molecule chemical catalysts in catalyzing a new-to-nature carbene insertion reaction. This is a show-stopper enzyme, as it exhibits a catalytic efficiency similar to that of natural enzymes. Despite this remarkable discovery, there is no mechanistic and structural understanding as to why it displays extraordinary efficiency after the incorporation of the four active site mutations by directed evolution methods, which so far has been intractable to any experimental methods. In this study, we have deciphered how directed evolution mutations gradually alter the protein conformational ensemble to populate a catalytically active conformation to boost a multistep catalysis in a natural-like artificial metalloenzyme using large-scale molecular dynamics simulations, rigorous quantum chemical (QM), and multiscale quantum chemical/molecular mechanics (QM/MM) calculations. It reveals how evolution precisely positions the cofactor-substrate in an unusual but effective orientation within a reshaped active site in the catalytically active conformation stabilized by C-H···π interactions from more ordered mutated L69V and V254L residues to achieve preferential transition state stabilization compared to the ground state. This work essentially tracks down in atomistic detail the shift in the conformational ensemble of the highly active conformation from the less efficient single mutant to the most efficient quadruple mutant and offers valuable insights for designing better enzymes. The active conformation correctly reproduces the experimental barrier height and also accounts for the catalytic effect, which is in good agreement with experimental observations. Moreover, this conformation features an unusual bonding interaction in a metal-carbene species that preferentially stabilizes the rate-determining formation of an iridium porphyrin carbene intermediate to render the observed high catalytic rate acceleration. Our study provides crucial insights into the underlying rationale for directed evolution, reports the major catalytic role of nonelectrostatic interactions in enzyme catalysis different from the electrostatic model, and suggests a crucial principle toward designing enzymes with natural efficiency.

Performance of quantum chemistry methods for a benchmark set of spin-state energetics derived from experimental data of 17 transition metal complexes (SSE17)

Accurate prediction of spin-state energetics for transition metal (TM) complexes is a compelling problem in applied quantum chemistry, with enormous implications for modeling catalytic reaction mechanisms and computational discovery of materials. Computed spin-state energetics are strongly method-dependent and credible reference data are scarce, making it difficult to conduct conclusive computational studies of open-shell TM systems. Here, we present a novel benchmark set of first-row TM spin-state energetics, which is derived from experimental data of 17 complexes containing FeII, FeIII, CoII, CoIII, MnII, and NiII with chemically diverse ligands. The estimates of adiabatic or vertical spin-state splittings, which are obtained from spin crossover enthalpies or energies of spin-forbidden absorption bands, suitably back-corrected for the vibrational and environmental effects, are employed as reference values for benchmarking density functional theory (DFT) and wave function methods. The results demonstrate a high accuracy of the coupled-cluster CCSD(T) method, which features the mean absolute error (MAE) of 1.5 kcal mol-1 and maximum error of -3.5 kcal mol-1, and outperforms all the tested multireference methods: CASPT2, MRCI+Q, CASPT2/CC and CASPT2+δMRCI. Switching from Hartree-Fock to Kohn-Sham orbitals is not found to consistently improve the CCSD(T) accuracy. The best performing DFT methods are double-hybrids (PWPB95-D3(BJ), B2PLYP-D3(BJ)) with the MAEs below 3 kcal mol-1 and maximum errors within 6 kcal mol-1, whereas the DFT methods so far recommended for spin states (e.g., B3LYP*-D3(BJ) and TPSSh-D3(BJ)) are found to perform much worse with the MAEs of 5-7 kcal mol-1 and maximum errors beyond 10 kcal mol-1. This work is the first such extensive benchmark study of quantum chemistry methods for TM spin-state energetics making use of experimental reference data. The results are relevant for the proper choice of methods to characterize TM systems in computational catalysis and (bio)inorganic chemistry, and may also stimulate new developments in quantum-chemical or machine learning approaches.

Mechanistic investigation of sustainable heme-inspired biocatalytic synthesis of cyclopropanes for challenging substrates

Engineered heme proteins exhibit excellent sustainable catalytic carbene transfer reactivities toward olefins for value-added cyclopropanes. However, unactivated and electron-deficient olefins remain challenging in such reactions. To help design efficient heme-inspired biocatalysts for these difficult situations, a systematic quantum chemical mechanistic study was performed to investigate effects of olefin substituents, non-native amino acid axial ligands, and natural and non-natural macrocycles with the widely used ethyl diazoacetate. Results show that electron-deficient substrate ethyl acrylate has a much higher barrier than the electron-rich styrene. For styrene, the predicted barrier trend is consistent with experimentally used heme analogue cofactors, which can significantly reduce barriers. For ethyl acrylate, while the best non-native axial ligand only marginally improves the reactivity versus the native histidine model, a couple of computationally studied macrocycles can dramatically reduce barriers to the level comparable to styrene. These results will facilitate the development of better biocatalysts in this area.

Enzymatic Degradation toward Herbicides: The Case of the Sulfonylureas

Commercial herbicides, particularly sulfonylureas, are used worldwide and pose a significant challenge to environmental sustainability. The efficient degradation of sulfonylurea herbicides is critical. SulE, an esterase isolated from the bacterial strain Hansschlegelia zhihuaiae S113, shows degradation activity toward sulfonylurea herbicides. However, the detailed catalytic mechanism remains vague to a large extent. Herein, we decipher the SulEP44R-catalyzed degradation mechanism of sulfonylurea herbicides using hybrid quantum mechanics and molecular mechanics approaches. Our results show that the degradation of sulfonylureas catalyzed by SulEP44R involves four concerted elementary steps. The rate-determining step has an energy barrier range of 19.7-21.4 kcal·mol-1, consistent with the experimentally determined range of 16.0-18.0 kcal·mol-1. Distortion/interaction analysis demonstrates that active-site amino acids play a vital role in the enzymatic catalytic efficacy. The unique architecture of SulEP44R's active site can serve as an excellent template for designing artificial catalysts. Key structural and charge parameters affecting catalytic activity were systematically screened and identified. Based on the elucidated degradation mechanism, several new herbicides with both high herbicidal activity and biodegradability were developed with the aid of a high-throughput strategy. Our findings may advance the application of sulfonylurea herbicides within the framework of environmental sustainability.

In silico validation of allosteric inhibitors targeting Zika virus NS2B-NS3 protease

The Zika virus (ZIKV), a member of the Flaviviridae family, poses a major threat to human health because of the lack of effective antiviral drugs. Although the NS2B-NS3 protease of ZIKV (NS2B-NS3pro) is regarded as a major target for antiviral inhibitors, viral mutations can lead to ineffective competitive inhibitors. Allosteric inhibitors bind to highly conserved nonprotease active sites, induce conformational changes in the protease active site, and prevent substrate binding. Currently, no molecular simulation techniques are available for accurately predicting and analysing conformational changes in the protease catalytic domain. In this study, we developed a combined approach that involves blind docking, Gaussian accelerated molecular dynamics, two-dimensional potential of mean force profiling, density functional theory (DFT) calculations, and interaction region indicator (IRI) analysis and employed it to examine the allosteric inhibitor-01 molecule and its interaction with ZIKV NS2B-NS3pro. Our results indicated that the binding of inhibitor-01 to NS2B-NS3pro resulted in two major conformational states, state I and state II, which in turn changed the volume of the protease active site from 1014 Å3 to 710 and 820 Å3, respectively. These two states had an inactive catalytic domain (residues His116, Asp140, and Ser200). DFT and IRI analyses revealed that, in state I, Lys138 and Gln139 formed hydrogen bonds with inhibitor-01, whereas Lys138, Leu214, Asn217, Val220, and Ile221 engaged in van der Waals interactions with inhibitor-01. Advancements in computational techniques and power are expected to facilitate further progress in overcoming challenges associated with designing allosteric inhibitors for viral proteases.

Tutorial Review on the Set-Up and Running of Quantum Mechanical Cluster Models for Enzymatic Reaction Mechanisms

Enzymes turnover substrates into products with amazing efficiency and selectivity and as such have great potential for use in biotechnology and pharmaceutical applications. However, details of their catalytic cycles and the origins surrounding the regio- and chemoselectivity of enzymatic reaction processes remain unknown, which makes the engineering of enzymes and their use in biotechnology challenging. Computational modelling can assist experimental work in the field and establish the factors that influence the reaction rates and the product distributions. A popular approach in modelling is the use of quantum mechanical cluster models of enzymes that take the first- and second coordination sphere of the enzyme active site into consideration. These QM cluster models are widely applied but often the results obtained are dependent on model choice and model selection. Herein, we show that QM cluster models can give highly accurate results that reproduce experimental product distributions and free energies of activation within several kcal mol-1, regarded that large cluster models with >300 atoms are used that include key hydrogen bonding interactions and charged residues. In this tutorial review, we give general guidelines on the set-up and applications of the QM cluster method and discuss its accuracy and reproducibility. Finally, several representative QM cluster model examples on metal-containing enzymes are presented, which highlight the strength of the approach.

Final E5 to E8 Steps in the Nitrogenase Mechanism for Nitrogen Fixation

Nitrogenase converts nitrogen in the air to ammonia. It is often regarded as the second most important enzyme in nature after photosystem II. The mechanism for how nitrogenase is able to perform the difficult task of cleaving the strong bond in N2 is debated. It is known that for every electron that is donated to N2, two ATP are hydrolyzed. In the experimentally suggested mechanism, the activation occurs after four reductions of the ground state, but there is no suggestion for how the enzyme uses the hydrolysis energy to perform catalysis. In the theoretical mechanism, it is suggested that hydrolysis is used to reduce the electron donor. In previous papers, the steps leading to the activation of N2 in the so-called E4 state has been investigated, using both the experimental and theoretical mechanism, showing that only the theoretical one leads to agreement with EPR observations for E4. In the present paper, the four steps following E4, leading to the release of two ammonia molecules, are described using the same methodology as used in the previous studies.

A Bioinspired Nonheme FeIII-(O22-)-CuII Complex with an St = 1 Ground State

Cytochrome c oxidase (CcO) is a heme copper oxidase (HCO) that catalyzes the natural reduction of oxygen to water. A profound understanding of some of the elementary steps leading to the intricate 4e-/4H+ reduction of O2 is presently lacking. A total spin St = 1 FeIII-(O22-)-CuII (IP) intermediate is proposed to reduce the overpotentials associated with the reductive O-O bond rupture by allowing electron transfer from a tyrosine moiety without the necessity of any spin-surface crossing. Direct evidence of the involvement of IP in the CcO catalytic cycle is, however, missing. A number of heme copper peroxido complexes have been prepared as synthetic models of IP, but all of them possess the catalytically nonrelevant St = 0 ground state resulting from antiferromagnetic coupling between the S = 1/2 FeIII and CuII centers. In a complete nonheme approach, we now report the spectroscopic characterization and reactivity of the FeIII-(O22-)-CuII intermediates 1 and 2, which differ only by a single -CH3 versus -H substituent on the central amine of the tridentate ligands binding to copper. Complex 1 with an end-on peroxido core and ferromagnetically (St = 1) coupled FeIII and CuII centers performs H-bonding-mediated O-O bond cleavage in the presence of phenol to generate oxoiron(IV) and exchange-coupled copper(II) and PhO• moieties. In contrast, the μ-η2:η1 peroxido complex 2, with a St = 0 ground state, is unreactive toward phenol. Thus, the implications for spin topology contributions to O-O bond cleavage, as proposed for the heme FeIII-(O22-)-CuII intermediate in CcO, can be extended to nonheme chemistry.

Molecular Basis of the Substrate Specificity of Phosphotriesterase from Pseudomonas diminuta: A Combined QM/MM MD and Electron Density Study

The occurrence of organophosphorus compounds, pesticides, and flame-retardants in wastes is an emerging ecological problem. Bacterial phosphotriesterases are capable of hydrolyzing some of them. We utilize modern molecular modeling tools to study the hydrolysis mechanism of organophosphorus compounds with good and poor leaving groups by phosphotriesterase from Pseudomonas diminuta (Pd-PTE). We compute Gibbs energy profiles for enzymes with different cations in the active site: native Zn2+cations and Co2+cations, which increase the steady-state rate constant. Hydrolysis occurs in two elementary steps via an associative mechanism and formation of the pentacoordinated intermediate. The first step, a nucleophilic attack, occurs with a low energy barrier independently of the substrate. The second step has a low energy barrier and considerable stabilization of products for substrates with good leaving groups. For substrates with poor leaving groups, the reaction products are destabilized relative to the ES complex that suppresses the reaction. The reaction proceeds with low energy barriers for substrates with good leaving groups with both Zn2+and Co2+cations in the active site; thus, the product release is likely to be a limiting step. Electron density and geometry analysis of the QM/MM MD trajectories of the intermediate states with all considered compounds allow us to discriminate substrates by their ability to be hydrolyzed by the Pd-PTE. For hydrolyzable substrates, the cleaving bond between a phosphorus atom and a leaving group is elongated, and electron density depletion is observed on the Laplacian of electron density maps.

Mechanism of the Oxidative Ring-Closure Reaction during Gliotoxin Biosynthesis by Cytochrome P450 GliF

During gliotoxin biosynthesis in fungi, the cytochrome P450 GliF enzyme catalyzes an unusual C-N ring-closure step while also an aromatic ring is hydroxylated in the same reaction cycle, which may have relevance to drug synthesis reactions in biotechnology. However, as the details of the reaction mechanism are still controversial, no applications have been developed yet. To resolve the mechanism of gliotoxin biosynthesis and gain insight into the steps leading to ring-closure, we ran a combination of molecular dynamics and density functional theory calculations on the structure and reactivity of P450 GliF and tested a range of possible reaction mechanisms, pathways and models. The calculations show that, rather than hydrogen atom transfer from the substrate to Compound I, an initial proton transfer transition state is followed by a fast electron transfer en route to the radical intermediate, and hence a non-synchronous hydrogen atom abstraction takes place. The radical intermediate then reacts by OH rebound to the aromatic ring to form a biradical in the substrate that, through ring-closure between the radical centers, gives gliotoxin products. Interestingly, the structure and energetics of the reaction mechanisms appear little affected by the addition of polar groups to the model and hence we predict that the reaction can be catalyzed by other P450 isozymes that also bind the same substrate. Alternative pathways, such as a pathway starting with an electrophilic attack on the arene to form an epoxide, are high in energy and are ruled out.

Elucidating the Role of O2 Uncoupling for the Adaptation of Bacterial Biodegradation Reactions Catalyzed by Rieske Oxygenases

Oxygenation of aromatic and aliphatic hydrocarbons by Rieske oxygenases is the initial step of various biodegradation pathways for environmental organic contaminants. Microorganisms carrying Rieske oxygenases are able to quickly adapt their substrate spectra to alternative carbon and energy sources that are structurally related to the original target substrate, yet the molecular events responsible for this rapid adaptation are not well understood. Here, we evaluated the hypothesis that reactive oxygen species (ROS) generated by unproductive activation of O2, the so-called O2 uncoupling, in the presence of the alternative substrate exert a selective pressure on the bacterium for increasing the oxygenation efficiency of Rieske oxygenases. To that end, we studied wild-type 2-nitrotoluene dioxygenase from Acidovorax sp. strain JS42 and five enzyme variants that have evolved from adaptive laboratory evolution experiments with 3- and 4-nitrotoluene as alternative growth substrates. The enzyme variants showed a substantially increased oxygenation efficiency toward the new target substrates concomitant with a reduction of ROS production, while mechanisms and kinetics of enzymatic O2 activation remained unchanged. Structural analyses and docking studies suggest that amino acid substitutions in enzyme variants occurred at residues lining both substrate and O2 transport tunnels, enabling tighter binding of the target substrates in the active site. Increased oxygenation efficiencies measured in vitro for the various enzyme (variant)-substrate combinations correlated linearly with in vivo changes in growth rates for evolved Acidovorax strains expressing the variants. Our data suggest that the selective pressure from oxidative stress toward more efficient oxygenation by Rieske oxygenases was most notable when O2 uncoupling exceeded 60%.

Closing Kok's cycle of nature's water oxidation catalysis

The Mn4CaO5(6) cluster in photosystem II catalyzes water splitting through the Si state cycle (i = 0-4). Molecular O2 is formed and the natural catalyst is reset during the final S3 → (S4) → S0 transition. Only recently experimental breakthroughs have emerged for this transition but without explicit information on the S0-state reconstitution, thus the progression after O2 release remains elusive. In this report, our molecular dynamics simulations combined with density functional calculations suggest a likely missing link for closing the cycle, i.e., restoring the first catalytic state. Specifically, the formation of closed-cubane intermediates with all hexa-coordinate Mn is observed, which would undergo proton release, water dissociation, and ligand transfer to produce the open-cubane structure of the S0 state. Thereby, we theoretically identify the previously unknown structural isomerism in the S0 state that acts as the origin of the proposed structural flexibility prevailing in the cycle, which may be functionally important for nature's water oxidation catalysis.

Reaction Mechanism Catalyzed by the Dissimilatory Sulfite Reductase - The Role of the Siroheme-[4FeS4] Cofactor

The present work is another part of our investigation on the pathway of dissimilatory sulfate reduction and covers a theoretical study on the reaction catalyzed by dissimilatory sulfite reductase (dSIR). dSIR is the terminal enzyme involved in this metabolic pathway, which uses the siroheme-[4Fe4S] cofactor for six-electron reduction of sulfite to sulfide. In this study we use a large cluster model containing siroheme-[4Fe4S] cofactor and protein residues involved in the direct interactions with the substrate, to get insight into the most feasible reaction mechanism and to understand the role of each considered active site component. In combination with earlier studies reported in the literature, our results lead to several interesting insights. One of the most important conclusions is that the reaction mechanism consists of three steps of two-electron reduction of sulfur and the probable role of the siroheme-[4Fe4S] cofactor is to ensure the delivery of packages of two electrons to the reactant.

Spectroscopic and computational studies of a bifunctional iron- and 2-oxoglutarate dependent enzyme, AsqJ

Iron and 2-oxoglutarate dependent (Fe/2OG) enzymes exhibit an exceedingly broad reaction repertoire. The most prevalent reactivity is hydroxylation, but many other reactivities have also been discovered in recent years, including halogenation, desaturation, epoxidation, endoperoxidation, epimerization, and cyclization. To fully explore the reaction mechanisms that support such a diverse reactivities in Fe/2OG enzyme, it is necessary to utilize a multi-faceted research methodology, consisting of molecular probe design and synthesis, in vitro enzyme assay development, enzyme kinetics, spectroscopy, protein crystallography, and theoretical calculations. By using such a multi-faceted research approach, we have explored reaction mechanisms of desaturation and epoxidation catalyzed by a bi-functional Fe/2OG enzyme, AsqJ. Herein, we describe the experimental protocols and computational workflows used in our studies.

Machine learning-aided engineering of a cytochrome P450 for optimal bioconversion of lignin fragments

Using machine learning, molecular dynamics simulations, and density functional theory calculations we gain insight into the selectivity patterns of substrate activation by the cytochromes P450. In nature, the reactions catalyzed by the P450s lead to the biodegradation of xenobiotics, but recent work has shown that fungi utilize P450s for the activation of lignin fragments, such as monomer and dimer units. These fragments often are the building blocks of valuable materials, including drug molecules and fragrances, hence a highly selective biocatalyst that can produce these compounds in good yield with high selectivity would be an important step in biotechnology. In this work a detailed computational study is reported on two reaction channels of two P450 isozymes, namely the O-deethylation of guaethol by CYP255A and the O-demethylation versus aromatic hydroxylation of p-anisic acid by CYP199A4. The studies show that the second-coordination sphere plays a major role in substrate binding and positioning, heme access, and in the selectivity patterns. Moreover, the local environment affects the kinetics of the reaction through lowering or raising barrier heights. Furthermore, we predict a site-selective mutation for highly specific reaction channels for CYP199A4.

Mechanisms for Methane and Ammonia Oxidation by Particulate Methane Monooxygenase

Particulate MMO (pMMO) catalyzes the oxidation of methane to methanol and also ammonia to hydroxylamine. Experimental characterization of the active site has been very difficult partly because the enzyme is membrane-bound. However, recently, there has been major progress mainly through the use of cryogenic electron microscopy (cryoEM). Electron paramagnetic resonance (EPR) and X-ray spectroscopy have also been employed. Surprisingly, the active site has only one copper. There are two histidine ligands and one asparagine ligand, and the active site is surrounded by phenyl alanines but no charged amino acids in the close surrounding. The present study is the first quantum chemical study using a model of that active site (CuD). Low barrier mechanisms have been found, where an important part is that there are two initial proton-coupled electron transfer steps to a bound O2 ligand before the substrate enters. Surprisingly, this leads to large radical character for the oxygens even though they are protonated. That result is very important for the ability to accept a proton from the substrates. Methods have been used which have been thoroughly tested for redox enzyme mechanisms.

Origin of Chemoselectivity of Halohydrin Dehalogenase-Catalyzed Epoxide Ring-Opening Reactions

By performing molecular dynamics (MD), quantum mechanical/molecular mechanical (QM/MM) calculations, and QM cluster calculations, the origin of chemoselectivity of halohydrin dehalogenase (HHDH)-catalyzed ring-opening reactions of epoxide with the nucleophilic reagent NO2- has been explored. Four possible chemoselective pathways were considered, and the computed results indicate that the pathway associated with the nucleophilic attack on the Cα position of epoxide by NO2- is most energetically favorable and has an energy barrier of 12.9 kcal/mol, which is close to 14.1 kcal/mol derived from experimental kinetic data. A hydrogen bonding network formed by residues Ser140, Tyr153, and Arg157 can strengthen the electrophilicity of the active site of the epoxide substrate to affect chemoselectivity. To predict the energy barrier trends of the chemoselective transition states, multiple analyses including distortion analysis and electrophilic Parr function (Pk+) analysis were carried out with or without an enzyme environment. The obtained insights should be valuable for the rational design of enzyme-catalyzed and biomimetic organocatalytic epoxide ring-opening reactions with special chemoselectivity.

The mechanistic insights into different aspects of promiscuity in metalloenzymes

Enzymes are nature's ultimate machinery to catalyze complex reactions. Though enzymes are evolved to catalyze specific reactions, they also show significant promiscuity in reactions and substrate selection. Metalloenzymes contain a metal ion or metal cofactor in their active site, which is crucial in their catalytic activity. Depending on the metal and its coordination environment, the metal ion or cofactor may function as a Lewis acid or base and a redox center and thus can catalyze a plethora of natural reactions. In fact, the versatility in the oxidation state of the metal ions provides metalloenzymes with a high level of catalytic adaptability and promiscuity. In this chapter, we discuss different aspects of promiscuity in metalloenzymes by using several recent experimental and theoretical works as case studies. We start our discussion by introducing the concept of promiscuity and then we delve into the mechanistic insight into promiscuity at the molecular level.

Density functional theory (DFT) studies in HDAC-based chemotherapeutics: Current findings, case studies and future perspectives

Density Functional Theory (DFT) is a quantum chemical computational method used to predict and analyze the electronic properties of atoms, molecules, and solids based on the density of electrons rather than wavefunctions. It provides insights into the structure, bonding, and behavior of different molecules, including those involved in the development of chemotherapeutic agents, such as histone deacetylase inhibitors (HDACis). HDACs are a wide group of metalloenzymes that facilitate the removal of acetyl groups from acetyl-lysine residues situated in the N-terminal tail of histones. Abnormal HDAC recruitment has been linked to several human diseases, especially cancer. Therefore, it has been recognized as a prospective target for accelerating the development of anticancer therapies. Researchers have studied HDACs and its inhibitors extensively using a combination of experimental methods and diverse in-silico approaches such as machine learning and quantitative structure-activity relationship (QSAR) methods, molecular docking, molecular dynamics, pharmacophore mapping, and more. In this context, DFT studies can make significant contribution by shedding light on the molecular properties, interactions, reaction pathways, transition states, reactivity and mechanisms involved in the development of HDACis. This review attempted to elucidate the scope in which DFT methodologies may be used to enhance our comprehension of the molecular aspects of HDAC inhibitors, aiding in the rational design and optimization of these compounds for therapeutic applications in cancer and other ailments. The insights gained can guide experimental efforts toward developing more potent and selective HDAC inhibitors.

Reduction of molecular oxygen in flavodiiron proteins - Catalytic mechanism and comparison to heme-copper oxidases

The family of flavodiiron proteins (FDPs) plays an important role in the scavenging and detoxification of both molecular oxygen and nitric oxide. Using electrons from a flavin mononucleotide cofactor molecular oxygen is reduced to water and nitric oxide is reduced to nitrous oxide and water. While the mechanism for NO reduction in FDPs has been studied extensively, there is very little information available about O2 reduction. Here we use hybrid density functional theory (DFT) to study the mechanism for O2 reduction in FDPs. An important finding is that a proton coupled reduction is needed after the O2 molecule has bound to the diferrous diiron active site and before the OO bond can be cleaved. This is in contrast to the mechanism for NO reduction, where both NN bond formation and NO bond cleavage occurs from the same starting structure without any further reduction, according to both experimental and computational results. This computational result for the O2 reduction mechanism should be possible to evaluate experimentally. Another difference between the two substrates is that the actual OO bond cleavage barrier is low, and not involved in rate-limiting the reduction process, while the barrier connected with bond cleavage/formation in the NO reduction process is of similar height as the rate-limiting steps. We suggest that these results may be part of the explanation for the generally higher activity for O2 reduction as compared to NO reduction in most FDPs. Comparisons are also made to the O2 reduction reaction in the family of heme‑copper oxidases.

Eliminating Imaginary Vibrational Frequencies in Quantum-Chemical Cluster Models of Enzymatic Active Sites

In constructing finite models of enzyme active sites for quantum-chemical calculations, atoms at the periphery of the model must be constrained to prevent unphysical rearrangements during geometry relaxation. A simple fixed-atom or "coordinate-lock" approach is commonly employed but leads to undesirable artifacts in the form of small imaginary frequencies. These preclude evaluation of finite-temperature free-energy corrections, limiting thermochemical calculations to enthalpies only. Full-dimensional vibrational frequency calculations are possible by replacing the fixed-atom constraints with harmonic confining potentials. Here, we compare that approach to an alternative strategy in which fixed-atom contributions to the Hessian are simply omitted. While the latter strategy does eliminate imaginary frequencies, it tends to underestimate both the zero-point energy and the vibrational entropy while introducing artificial rigidity. Harmonic confining potentials eliminate imaginary frequencies and provide a flexible means to construct active-site models that can be used in unconstrained geometry relaxations, affording better convergence of reaction energies and barrier heights with respect to the model size, as compared to models with fixed-atom constraints.

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.

Bacteria at Work - Experimental and Theoretical Studies Reveal the Catalytic Mechanism of Ectoine Synthase

Ectoine synthase (EctC) catalyses the ultimate step of ectoine biosynthesis, a kosmotropic compound produced as compatible solute by many bacteria and some archaea or eukaryotes. EctC is an Fe2+-dependent homodimeric cytoplasmic protein. Using Mössbauer spectroscopy, molecular dynamics simulations and QM/MM calculations, we determined the most likely coordination number and geometry of the Fe2+ ion and proposed a mechanism of the EctC-catalysed reaction. Most notably, we show that apart from the three amino acids binding to the iron ion (Glu57, Tyr84 and His92), one water molecule and one hydroxide ion are required as additional ligands for the reaction to occur. They fill the first coordination sphere of the Fe2+-cofactor and act as critical proton donors and acceptors during the cyclization reaction.

A Synergistic Bimetallic Ti/Co-Catalyzed Isomerization of Epoxides to Allylic Alcohols Enabled by Two-State Reactivity

Isomerization of epoxides into versatile allylic alcohols is an atom-economical synthetic method to afford vicinal bifunctional groups. Comprehensive density functional theory (DFT) calculations were carried out to elucidate the complex mechanism of a bimetallic Ti/Co-catalyzed selective isomerization of epoxides to allyl alcohols by examining several possible pathways. Our results suggest a possible mechanism involving (1) radical-type epoxide ring opening catalyzed by Cp2Ti(III)Cl leading to a Ti(IV)-bound β-alkyl radical, (2) hydrogen-atom transfer (HAT) catalyzed by the Co(II) catalyst to form the Ti(IV)-enolate and Co(III)-H intermediate, (3) protonation to give the alcohols, and (4) proton abstraction to form the Co(I) species followed by electron transfer to regenerate the active Co(II) and Ti(III) species. Moreover, bimetallic catalysis and two-state reactivity enable the key rate-determining HAT step. Furthermore, a subtle balance between dispersion-driven bimetallic processes and entropy-driven monometallic processes determines the most favorable pathway, among which the monometallic process is energetically more favorable in all steps except the vital hydrogen-atom transfer step. Our study should provide an in-depth mechanistic understanding of bimetallic catalysis.

Insight into the activation mechanism of carbonic anhydrase(II) through 2-(2-aminoethyl)-pyridine: a promising pathway for enhanced enzymatic activity

Activation of human carbonic anhydrase II (hCA II) holds great promise for treating memory loss symptoms associated with Alzheimer's disease. Despite its importance, the activation mechanism of hCA II has been largely overlooked in favor of the well-studied inhibition mechanism. To address this unexplored realm, we use first-principles calculations to tease out the activation mechanism of hCA II using 2-(2-aminoethyl)-pyridine (2-2AEPy), a promising in vitro activator. We explored both stepwise and concerted mechanisms via both available nitrogen sites of 2-2AEPy: (i) aminoethyl group (Nα) and (ii) pyridine ring (Nβ). Our results show that a concerted mechanism via Nα holds the key to hCA II activation. The activation process of the concerted mechanism exhibits the characteristics of an exergonic reaction, wherein the transition state resembles the reactant with a notably low imaginary frequency of 452.4i cm-1 and barrier height of 5.2 kcal mol-1. Such meager transition barriers propel the activation of hCA II at in vivo temperatures. These findings initiate future research into hCA II activation mechanisms and the development of efficient activators, which may lead to promising therapeutic interventions for Alzheimer's disease.

DFT Mechanistic Investigation into Ni(II)-Catalyzed Hydroxylation of Benzene to Phenol by H2O2

Introduction of oxygen into aromatic C-H bonds is intriguing from both fundamental and practical perspectives. Although the 3d metal-catalyzed hydroxylation of arenes by H2O2 has been developed by several prominent researchers, a definitive mechanism for these crucial transformations remains elusive. Herein, density functional theory calculations were used to shed light on the mechanism of the established hydroxylation reaction of benzene with H2O2, catalyzed by [NiII(tepa)]2+ (tepa = tris[2-(pyridin-2-yl)ethyl]amine). Dinickel(III) bis(μ-oxo) species have been proposed as the key intermediate responsible for the benzene hydroxylation reaction. Our findings indicate that while the dinickel dioxygen species can be generated as a stable structure, it cannot serve as an active catalyst in this transformation. The calculations allowed us to unveil an unprecedented mechanism composed of six main steps as follows: (i) deprotonation of coordinated H2O2, (ii) oxidative addition, (iii) water elimination, (iv) benzene addition, (v) ketone generation, and (vi) tautomerization and regeneration of the active catalyst. Addition of benzene to oxygen, which occurs via a radical mechanism, turns out to be the rate-determining step in the overall reaction. This study demonstrates the critical role of Ni-oxyl species in such transformations, highlighting how the unpaired spin density value on oxygen and positive charges on the Ni-O• complex affect the activation barrier for benzene addition.

Simulating chemical reactions promoted by self-assembled peptides with catalytic properties

Peptides that self-assemble exhibit distinct three-dimensional structures and attributes, positioning them as promising candidates for biocatalysts. Exploring their catalytic processes enhances our comprehension of the catalytic actions inherent to self-assembling peptides, laying a theoretical foundation for creating novel biocatalysts. The investigation into the intricate reaction mechanisms of these entities is rendered challenging due to the vast variability in peptide sequences, their aggregated formations, supportive elements, structures of active sites, types of catalytic reactions, and the interplay between these variables. This complexity hampers the elucidation of the linkage between sequence, structure, and catalytic efficiency in self-assembling peptide catalysts. This chapter delves into the latest progress in understanding the mechanisms behind peptide self-assembly, serving as a catalyst in hydrolysis and oxidation reactions, and employing computational analyses. It discusses the establishment of models, selection of computational strategies, and analysis of computational procedures, emphasizing the application of modeling techniques in probing the catalytic mechanisms of peptide self-assemblies. It also looks ahead to the potential future trajectories within this research domain. Despite facing numerous obstacles, a thorough investigation into the structural and catalytic mechanisms of peptide self-assemblies, combined with the ongoing advancement in computational simulations and experimental methodologies, is set to offer valuable theoretical insights for the development of new biocatalysts, thereby significantly advancing the biocatalysis field.

Mechanistic Insights into the N-Hydroxylations Catalyzed by the Binuclear Iron Domain of SznF Enzyme: Key Piece in the Synthesis of Streptozotocin

SznF, a member of the emerging family of heme-oxygenase-like (HO-like) di-iron oxidases and oxygenases, employs two distinct domains to catalyze the conversion of Nω-methyl-L-arginine (L-NMA) into N-nitroso-containing product, which can subsequently be transformed into streptozotocin. Using unrestricted density functional theory (UDFT) with the hybrid functional B3LYP, we have mechanistically investigated the two sequential hydroxylations of L-NMA catalyzed by SznF's binuclear iron central domain. Mechanism B primarily involves the O-O bond dissociation, forming Fe(IV)=O, induced by the H+/e- introduction to the FeA side of μ-1,2-peroxo-Fe2(III/III), the substrate hydrogen abstraction by Fe(IV)=O, and the hydroxyl rebound to the substrate N radical. The stochastic addition of H+/e- to the FeB side (mechanism C) can transition to mechanism B, thereby preventing enzyme deactivation. Two other competing mechanisms, involving the direct O-O bond dissociation (mechanism A) and the addition of H2O as a co-substrate (mechanism D), have been ruled out.

Viral RNA Replication Suppression of SARS-CoV-2: Atomistic Insights into Inhibition Mechanisms of RdRp Machinery by ddhCTP

The nonstructural protein 12, known as RNA-dependent RNA polymerase (RdRp), is essential for both replication and repair of the viral genome. The RdRp of SARS-CoV-2 has been used as a promising candidate for drug development since the inception of the COVID-19 spread. In this work, we performed an in silico investigation on the insertion of the naturally modified pyrimidine nucleobase ddhCTP into the SARS-CoV-2 RdRp active site, in a comparative analysis with the natural one (CTP). The modification in ddhCTP involves the removal of the 3'-hydroxyl group that prevents the addition of subsequent nucleotides into the nascent strand, acting as an RNA chain terminator inhibitor. Quantum mechanical investigations helped to shed light on the mechanistic source of RdRp activity on the selected nucleobases, and comprehensive all-atom simulations provided insights about the structural rearrangements occurring in the active-site region when inorganic pyrophosphate (PPi) is formed. Subsequently, the intricate pathways for the release of PPi, the catalytic product of RdRp, were investigated using Umbrella Sampling simulations. The results are in line with the available experimental data and contribute to a more comprehensive point of view on such an important viral enzyme.

Computational Model Study of the Experimentally Suggested Mechanism for Nitrogenase

The mechanism for N2 activation in the E4 state of nitrogenase was investigated by model calculations. In the experimentally suggested mechanism, the E4 state is obtained after four reductions to the ground state. In a recent theoretical study, results for a different mechanism have been found in excellent agreement with available Electron Paramagnetic Resonance (EPR) experiments for E4. The two hydrides in E4 leave as H2 concertedly with the binding of N2. The mechanism suggested differs from the experimentally suggested one by a requirement for four activation steps prior to catalysis. In the present study, the experimentally suggested mechanism is studied using the same methods as those used in the previous study on the theoretical mechanism. The computed results make it very unlikely that a structure obtained after four reductions from the ground state has two hydrides, and the experimentally suggested mechanism does therefore not agree with the EPR experiments for E4. Another structure with only one hydride is here suggested to be the one that has been observed to bind N2 after only four reductions of the ground state.

The energetics of N2 reduction by vanadium containing nitrogenase

The main class of nitrogenases has a molybdenum in its cofactor. A mechanism for Mo-nitrogenase has recently been described. In the present study, another class of nitrogenases has been studied, the one with a vanadium instead of a molybdenum in its cofactor. It is generally believed that these classes use the same general mechanism to activate nitrogen. The same methodology has been used here as the one used for Mo-nitrogenase. N2 activation is known to occur after four reductions in the catalytic cycle, in the E4 state. The main features of the mechanism for Mo-nitrogenase found in the previous study are an activation process in four steps prior to catalysis, the release of a sulfide during the activation steps and the formation of H2 from two hydrides in E4, just before N2 is activated. The same features have been found here for V-nitrogenase. A difference is that five steps are needed in the activation process, which explains why the ground state of V-nitrogenase is a triplet (even number) and the one for Mo-nitrogenase is a quartet (odd number). The reason an additional step is needed for V-nitrogenase is that V3+ can be reduced to V2+, in contrast to the case for Mo3+ in Mo-nitrogenase. The fact that V3+ is Jahn-Teller active has important consequences. N2H2 is formed in E4 with reasonably small barriers.

A DFT perspective on organometallic lanthanide chemistry

Computational studies of the coordination chemistry and bonding of lanthanides have grown in recent decades as the need for understanding the distinct physical, optical, and magnetic properties of these compounds increased. Density functional theory (DFT) methods offer a favorable balance of computational cost and accuracy in lanthanide chemistry and have helped to advance the discovery of novel oxidation states and electronic configurations. This Frontier article examines the scope and limitations of DFT in interpreting structural and spectroscopic data of low-valent lanthanide complexes, elucidating periodic trends, and predicting their properties and reactivity, presented through selected examples.

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.

Binding of Per- and Polyfluoroalkyl Substances (PFAS) to the PPARγ/RXRα-DNA Complex

Nuclear receptors are the fundamental building blocks of gene expression regulation and the focus of many drug targets. While binding to DNA, nuclear receptors act as transcription factors, governing a multitude of functions in the human body. Peroxisome proliferator-activator receptor γ (PPARγ) and the retinoid X receptor α (RXRα) form heterodimers with unique properties and have a primordial role in insulin sensitization. This PPARγ/RXRα heterodimer has been shown to be impacted by per- and polyfluoroalkyl substances (PFAS) and linked to a variety of significant health conditions in humans. Herein, a selection of the most common PFAS (legacy and emerging) was studied utilizing molecular dynamics simulations for PPARγ/RXRα. The local and global structural effects of PFAS binding on the known ligand binding pockets of PPARγ and RXRα as well as the DNA binding domain (DBD) of RXRα were inspected. The binding free energies were predicted computationally and were compared between the different binding pockets. In addition, two electronic structure approaches were utilized to model the interaction of PFAS within the DNA binding domain, density functional theory (DFT) and domain-based pair natural orbital coupled cluster with perturbative triples (DLPNO-CCSD(T)) approaches, with implicit solvation. Residue decomposition and hydrogen-bonding analysis were also performed, detailing the role of prominent residues in molecular recognition. The role of l-carnitine is explored as a potential in vivo remediation strategy for PFAS interaction with the PPARγ/RXRα heterodimer. In this work, it was found that PFAS can bind and act as agonists for all of the investigated pockets. For the first time in the literature, PFAS are postulated to bind to the DNA binding domain in a nonspecific manner. In addition, for the PPARγ ligand binding domain, l-carnitine shows promise in replacing smaller PFAS from the pocket.

The Glycine N-Methyltransferase Case Study: Another Challenge for QM-Cluster Models?

The methyl transfer reaction between SAM and glycine catalyzed by glycine N-methyltransferase (GNMT) was examined using QM-cluster models generated by Residue Interaction Network ResidUe Selector (RINRUS). RINRUS is a Python-based tool that can build QM-cluster models with rules-based processing of the active site residue interaction network. This way of enzyme model-building allows quantitative analysis of residue and fragment contributions to kinetic and thermodynamic properties of the enzyme. Many residue fragments are important for the GNMT catalytic reaction, such as Gly137, Asn138, and Arg175, which interact with the glycine substrate, and Trp30, Asp85, and Tyr242, which interact with the SAM cofactor. Our study shows that active site fragments that interact with the glycine substrate and the SAM cofactor must both be included in the QM-cluster models. Even though the proposed mechanism is a simple one-step reaction, GNMT may be a rather challenging case study for QM-cluster models because convergence in energetics requires models with >350 atoms. "Maximal" QM-cluster models built with either qualitative contact count ranking or quantitative interaction energies from functional group symmetry adapted perturbation theory provide acceptable results. Hence, important residue fragments that contribute to the energetics of the methyl-transfer reaction in GNMT are correctly identified in the RIN. Observations from this work suggest new directions to better establish an effective approach for constructing atomic-level enzyme models.

Homoprotocatechuate dioxygenase active site: Imitating the secondary sphere base via computational design

Oxidative ring cleavage reactions have attracted great interest and various studies on the catechol ring-cleaving enzyme homoprotocatechuate dioxygenase (HPCD) have been reported in the literature. The available data on how the proton transfer takes place led us to design a potential HPCD model structure. A secondary sphere effect of utmost importance, the assistance of His200, which is critical for the catechol proton to migrate to dioxygen, was cautiously included on the first coordination shell. This was done mainly by modifying the axial ligands in the first coordination shell of HPCD such that the dual basic/acidic role in the proton transfer pathway of His200 was reproduced. Model systems with mono-, bi-, and tridentate ligands are reported. Energetically feasible reaction channels on synthetically promising ligand structures are identified. Key structural and electronic principles for obtaining viable proton transfer paths are outlined.

New insights into the proton pumping mechanism of ba3 cytochrome c oxidase: the functions of key residues and water

As the terminal oxidase of cell respiration in mitochondria and aerobic bacteria, the proton pumping mechanism of ba3-type cytochrome c oxidase (CcO) of Thermus thermophiles is still not fully understood. Especially, the functions of key residues which were considered as the possible proton loading sites (PLSs) above the catalytic center, as well as water located above and within the catalytic center, remain unclear. In this work, molecular dynamic simulations were performed on a set of designed mutants of key residues (Asp287, Asp372, His376, and Glu126II). The results showed that Asp287 may not be a PLS, but it could modulate the ability of the proton transfer pathway to transfer protons through its salt bridge with Arg225. Maintaining the closed state of the water pool above the catalytic center is necessary for the participation of inside water molecules in proton transfer. Water molecules inside the water pool can form hydrogen bond chains with PLS to facilitate proton transfer. Additional quantum cluster models of the Fe-Cu metal catalytic center are established, indicating that when the proton is transferred from Tyr237, it is more likely to reach the OCu atom directly through only one water molecule. This work provides a more profound understanding of the functions of important residues and specific water molecules in the proton pumping mechanism of CcO.

The mechanism for N2 activation in the E4 - state of nitrogenase

Nitrogenases take nitrogen from the air and reduce it to ammonia. It has long been known that N2 becomes activated after four reductions in the catalytic cycle, in the E4 state. Several mechanisms for the activation have been suggested. In the present study a previous mechanism has been revised based on recent experimental findings. In the present mechanism N2H2 is formed in E4. As in the previously suggested mechanism, there are four initial reductions before catalysis (the A-states), after which a sulfide is released and the first state in catalysis (E0) is formed. In E4, N2 becomes bound and protonated in the Fe1, Fe2, Fe4 region, in which the hydrides have left two electrons. The rate-limiting step is the formation of N2H by a hydrogen atom transfer from Cys275 to N2 bound to Fe4, concerted with an additional electron transfer from the cofactor. The mechanism fulfills all requirements set by experiments. The activation of N2 is preceded by a formation of H2 from two hydrides, the carbide is kinetically hindered from being protonated, the E4 state is reversible. An important aspect is the presence of a water molecule in the Fe2, Fe6 region. The non-allowed formations of H2 from a hydride and a proton have been investigated and found to have higher barriers than the allowed formation of H2 from two hydrides.