Spin-Forbidden Ligand Binding to the Ferrous−Heme Group: Ab Initio and DFT Studies

Computational simulations uncover enantioselective metabolism of chiral triazole fungicides by human CYP450 enzymes: A case study of tebuconazole

Tebuconazole (TEB), a prominent chiral triazole fungicide, has been extensively utilized for plant pathogen control globally. Despite experimental evidence of TEB metabolism in mammals, the enantioselectivity in the biotransformation of R- and S-TEB enantiomers by specific CYP450s remains elusive. In this work, integrated in silico simulations were employed to unveil the binding interactions and enantioselective metabolic fate of TEB enantiomers within human CYP1A2, 2B6, 2E1, and 3A4. Molecular dynamics (MD) simulations clearly delineated the binding specificity of R- and S-TEB to the four CYP450s, crucially determining their differences in metabolic activity and enantioselectivity. The primary driving force for robust ligand binding was identified as van der Waals interactions with CYP450s, particularly involving the hydrophobic residues. Mechanistic insights derived from quantum mechanics/molecular mechanics (QM/MM) calculations established C2-methyl hydroxylation as the predominant route of R-/S-TEB metabolism, while C6-hydroxylation and triazol epoxidation were deemed kinetically infeasible pathways. Specifically, the resulting hydroxy-R-TEB metabolite primarily originates from R-TEB biotransformation by 1A2, 2E1 and 3A4, whereas hydroxy-S-TEB is preferentially produced by 2B6. These findings significantly contribute to our comprehension of the binding specificity and enantioselective metabolic fate of chiral TEB by CYP450s, potentially informing further research on human health risk assessment associated with TEB exposure.

Gaseous inhibition of the transsulfuration pathway by cystathionine β-synthase

The transsulfuration pathway plays a key role in mammals for maintaining the balance between cysteine and homocysteine, whose concentrations are critical in several biochemical processes. Human cystathionine β-synthase is a heme-containing, pyridoxal 5'-phosphate (PLP)-dependent enzyme found in this pathway. The heme group does not participate directly in catalysis, but has a regulatory function, whereby CO or NO binding inhibits the PLP-dependent reactions. In this study, we explore the detailed structural changes responsible for inhibition using quantum chemical calculations to validate the experimentally observed bonding patterns associated with heme CO and NO binding and molecular dynamics simulations to explore the medium-range structural changes triggered by gas binding and propagating to the PLP active site, which is more than 20 Å distant from the heme group. Our results support a previously proposed mechanical signaling model, whereby the cysteine decoordination associated with gas ligand binding leads to breaking of a hydrogen bond with an arginine residue on a neighbouring helix. In turn, this leads to a shift in position of the helix, and hence also of the PLP cofactor, ultimately disrupting a key hydrogen bond that stabilizes the PLP in its catalytically active form.

Binding of Carbon Monoxide to Hemoglobin in an Oxygen Environment: Force Field Development for Molecular Dynamics

Carbon monoxide (CO) is a byproduct of the incomplete combustion of carbon-based fuels, such as wood, coal, gasoline, or natural gas. As incomplete combustion in a fire accident or in an engine, massively produced CO leads to a serious life threat because CO competes with oxygen (O2) binding to hemoglobin and makes people suffer from hypoxia. Although there is hyperbaric O2 therapy for patients with CO poisoning, the nanoscale mechanism of CO dissociation in the O2-rich environment is not completely understood. In this study, we construct the classical force field parameters compatible with the CHARMM for simulating the coordination interactions between hemoglobin, CO, and O2, and use the force field to reveal the impact of O2 on the binding strength between hemoglobin and CO. Density functional theory and Car-Parrinello molecular dynamics simulations are used to obtain the bond energy and equilibrium geometry, and we used machine learning enabled via a feedforward neural network model to obtain the classical force field parameters. We used steered molecular dynamics simulations with a force field to characterize the mechanical strength of the hemoglobin-CO bond before rupture under different simulated O2-rich environments. The results show that as O2 approaches the Fe2+ of heme at a distance smaller than ∼2.8 Å, the coordination bond between CO and Fe2+ is reduced to 50% bond strength in terms of the peak force observed in the rupture process. This weakening effect is also shown by the free energy landscape measured by our metadynamics simulation. Our work suggests that the O2-rich environment around the hemoglobin-CO bond effectively weakens the bonding, so that designing of O2 delivery vector to the site is helpful for alleviating CO binding, which may shed light on de novo drug design for CO poisoning.

First-principles study of molecular hydrogen binding to heme in competition with O2, NO and CO

Molecular hydrogen shows antioxidant activity and distinct efficacy towards vascular diseases, but the understanding of this is not yet satisfactory at the atomic level. In this work, we study the binding properties of H2 to the heme group in relation with other diatomic molecules (DMs), including O2, NO and CO, and their displacement reactions, using first-principles calculations. We carry out molecular modeling of the heme group, using iron-porphyrin with the imidazole ligand, i.e., FePIm, and smaller models of Fe(CnHn+2N2)2NH3 with n = 3 and 1, and of molecular complexes of heme-DM and -H. Through analysis of optimized geometries and energetics, it is found that the order of binding strength of DMs or H to the Fe of heme is NO > O2 > CO > H > H2 for FePIm-based systems, while it is H > O2 > NO > CO > H2 for model-based systems. We calculate the activation energies for displacement reactions of H2 and H by other DMs, revealing that the H2 displacements occur spontaneously while the H displacements require a large amount of energy. Finally, our calculations corroborate that the rate constants increase with increasing temperature according to the Arrhenius relation.

Insights on cyclophosphamide metabolism and anticancer mechanism of action: A computational study

The oxazaphosphorine cyclophosphamide (CP) is a DNA-alkylating agent commonly used in cancer chemotherapy. This anticancer agent is administered as a prodrug activated by a liver cytochrome P450-catalyzed 4-hydroxylation reaction that yields the active, cytotoxic metabolite. The primary metabolite, 4-hydroxycyclophosphamide, equilibrates with the ring-open aldophosphamide that undergoes β-elimination to yield the therapeutically active DNA cross-linking phosphoramide mustard and the byproduct acrolein. The present paper presents a DFT investigation of the different metabolic phases and an insight into the mechanism by which CP exerts its cytotoxic action. A detailed computational analysis of the energy profiles describing all the involved transformations and the mechanism of DNA alkylation is given with the aim to contribute to an increase of knowledge that, after more than 60 years of unsuccessful attempts, can lead to the design and development of a new generation of oxazaphosphorines.

New insight into biodegradation mechanism of phenylurea herbicides by cytochrome P450 enzymes: Successive N-demethylation mechanism

Phenylurea herbicides (PUHs) present one of the most important herbicides, which have cause serious effects on ecological environment and humans. Nowadays enzyme strategy shows great advantages in degradation of PUHs. Here density functional theory (DFT), quantitative structure - activity relationship (QSAR) and quantum mechanics/molecular mechanics (QM/MM) approaches are used to investigate the degradation mechanism of PUHs catalyzed by P450 enzymes. Two successive N-demethylation pathways are identified and two hydrogen abstraction (H-abstraction) reaction pathways are identified as the rate-determining step through high-throughput DFT calculations. The Boltzmann-weighted average energy barrier of the second H-abstraction pathway (19.95 kcal/mol) is higher than that of the first H-abstraction pathway (16.80 kcal/mol). Two QSAR models are established to predict the energy barriers of the two H-abstraction pathways based on the quantum chemical descriptors and mordred molecular descriptors. The determination coefficient (R2) values of QSAR models are > 0.9, which reveal that the established QSAR models have great predictive capability. QM/MM calculations indicate that human P450 enzymes are more efficient in degradation of PUHs than crop and weed P450 enzymes. Correlations between energy barriers and key structural/charge parameters are revealed and key parameters that have influence on degradation efficiency of PUHs are identified. This study provides lateral insights into the biodegradation strategy and removal method of PUHs and valuable information for designing or engineering of highly efficient degradation enzymes and genetically modified crops.

Benchmarks for transition metal spin-state energetics: why and how to employ experimental reference data?

Accurate prediction of energy differences between alternative spin states of transition metal complexes is essential in computational (bio)inorganic chemistry-for example, in characterization of spin crossover materials and in the theoretical modeling of open-shell reaction mechanisms-but it remains one of the most compelling problems for quantum chemistry methods. A part of this challenge is to obtain reliable reference data for benchmark studies, as even the highest-level applicable methods are known to give divergent results. This Perspective discusses two possible approaches to method benchmarking for spin-state energetics: using either theoretically computed or experiment-derived reference data. With the focus on the latter approach, an extensive general review is provided for the available experimental data of spin-state energetics and their interpretations in the context of benchmark studies, targeting the possibility of back-correcting the vibrational effects and the influence of solvents or crystalline environments. With a growing amount of experience, these effects can be now not only qualitatively understood, but also quantitatively modeled, providing the way to derive nearly chemically accurate estimates of the electronic spin-state gaps to be used as benchmarks and advancing our understanding of the phenomena related to spin states in condensed phases.

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.

Carbon monoxide separation: past, present and future

Large amounts of carbon monoxide are produced by industrial processes such as biomass gasification and steel manufacturing. The CO present in vent streams is often burnt, this produces a large amount of CO2, e.g., oxidation of CO from metallurgic flue gasses is solely responsible for 2.7% of manmade CO2 emissions. The separation of N2 from CO due to their very similar physical properties is very challenging, meaning that numerous energy-intensive steps are required for CO separation, making the CO separation from many process streams uneconomical in spite of CO being a valuable building block in the production of major chemicals through C1 chemistry and the production of linear hydrocarbons by the Fischer-Tropsch process. The development of suitable processes for the separation of carbon monoxide has both industrial and environmental significance. Especially since CO is a main product of electrocatalytic CO2 reduction, an emerging sustainable technology to enable carbon neutrality. This technology also requires an energy-efficient separation process. Therefore, there is a great need to develop energy efficient CO separation processes adequate for these different process streams. As such the urgency of separating carbon monoxide is gaining greater recognition, with research in the field becoming more and more crucial. This review details the principles on which CO separation is based and provides an overview of currently commercialised CO separation processes and their limitations. Adsorption is identified as a technology with the potential for CO separation with high selectivity and energy efficiency. We review the research efforts, mainly seen in the last decades, in developing new materials for CO separation via ad/bsorption and membrane technology. We have geared our review to both traditional CO sources and emerging CO sources, including CO production from CO2 conversion. To that end, a variety of emerging processes as potential CO2-to-CO technologies are discussed and, specifically, the need for CO capture after electrochemical CO2 reduction is highlighted, which is still underexposed in the available literature. Altogether, we aim to highlight the knowledge gaps that could guide future research to improve CO separation performance for industrial implementation.

Emerging Metabolic Profiles of Sulfonamide Antibiotics by Cytochromes P450: A Computational-Experimental Synergy Study on Emerging Pollutants

Metabolism, especially by CYP450 enzymes, is the main reason for mediating the toxification and detoxification of xenobiotics in humans, while some uncommon metabolic pathways, especially for emerging pollutants, probably causing idiosyncratic toxicity are easily overlooked. The pollution of sulfonamide antibiotics in aqueous system has attracted increasing public attention. Hydroxylation of the central amine group can trigger a series of metabolic processes of sulfonamide antibiotics in humans; however, this work parallelly reported the coupling and fragmenting initiated by amino H-abstraction of sulfamethoxazole (SMX) catalyzed by human CYP450 enzymes. Elucidation of the emerging metabolic profiles was mapped via a multistep synergy between computations and experiments, involving preliminary DFT computations and in vitro and in vivo assays, profiling adverse effects, and rationalizing the fundamental factors via targeted computations. Especially, the confirmed SMX dimer was shown to potentially act as a metabolism disruptor in humans, while spin aromatic delocalization resulting in the low electron donor ability of amino radicals was revealed as the fundamental factor to enable coupling of sulfonamide antibiotics by CYP450 through the nonconventional nonrebound pathway. This work may further strengthen the synergistic use of computations prior to experiments to avoid wasteful experimental screening efforts in environmental chemistry and toxicology.

New insight into molecular mechanism of P450-Catalyzed metabolism of emerging contaminants and its consequence for human health: A case study of preservative methylparaben

Hydroxylated metabolites in the living body are considered as a potential biomarker of exposure to emerging contaminations (ECs) and breast cancer, but their formation mechanism has not received enough attention. Besides, the adverse impacts of metabolites during the metabolic transformation of ECs largely remain unknown. In this study, we employed a density functional calculation combing with in-vitro incubation of human liver microsomes to explore the bio-transformation of preservative methylparaben (MPB) in human bodies. Our results showed that hydroxylated metabolites of MPB (OH-MPB) were observed experimentally, while a formation mechanism was revealed at the molecular level. That is, hydroxylated metabolite was exclusively formed via the hydrogen abstraction from the phenolic hydroxyl group of MPB followed by the OH-rebound pathway, rather than the direct hydroxylation on the benzene ring. The increasing of hydroxyl groups on ECs could improve the metabolisms. This was confirmed in the metabolism of ECs without hydroxyl group and with multiple-hydroxyl groups, respectively. Furthermore, toxicity assessments show that compared to parent MPB, the hydroxylated metabolites have increased negative impacts on the gastrointestinal system and liver. A semiquinone product exhibits potential damage in the cardiovascular system and epoxides are toxic to the blood and gastrointestinal system. The findings deepen our insight into the biotransformation of parabens in human health, especially by providing health warnings about the potential impacts caused by semiquinone and epoxides.

Theoretical Kinetic Isotope Effects in Establishing the Precise Biodegradation Mechanisms of Organic Pollutants

Compound-specific isotope analysis (CSIA) for natural isotope ratios has been recognized as a promising tool to elucidate biodegradation pathways of organic pollutants by microbial enzymes by relating reported kinetic isotope effects (KIEs) to apparent KIEs (AKIEs) derived from bulk isotope fractionations (ε_bulk). However, for many environmental reactions, neither are the reference KIE ranges sufficiently narrow nor are the mechanisms elucidated to the point that rate-determining steps have been identified unequivocally. In this work, besides providing reference KIEs and rationalizing AKIEs, good relationships have been explained by DFT computations for diverse biodegradation pathways with known enzymatic models between the theoretical isotope fractionations (ε_bulk^') from intrinsic KIEs on the rate-determining steps and the observed ε_bulk. (1) To confirm the mechanistic details of previously reported pathway-dependent CSIA, it includes isotope changes in MTBE biodegradation between hydroxylation by CYP450 and S_N2 reaction by cobalamin-dependent methyltransferase, the regioselectivity of toluene biodegradation by CYP450, and the rate-determining step in toluene biodegradation by benzylsuccinate synthase. (2) To yield new fundamental insights into some unclear biodegradation pathways, it consists of the oxidative function of toluene dioxygenase in biodegradation of TCE, the epoxidation mode in biodegradation of TCE by toluene 4-monooxygenase, and the weighted average mechanism in biodegradation of cDCE by CYP450.

Theoretical Study of Hydroxylation of α- and β-Pinene by a Cytochrome P450 Monooxygenase Model

Previous studies on biocatalytic transformations of pinenes by cytochrome P450 (CYP) enzymes reveal the formation of different oxygenated products from a single substrate due to the multistate reactivity of CYP and the many reactive sites in the pinene scaffold. Up until now, the detailed mechanism of these biocatalytic transformations of pinenes have not been reported. Hereby, we report a systematic theoretical study of the plausible hydrogen abstraction and hydroxylation reactions of α- and β-pinenes by CYP using the density functional theory (DFT) method. All DFT calculations in this study were based on B3LYP/LAN computational methodology using the Gaussian09 software. We used the B3LYP functional with corrections for dispersive forces, BSSE, and anharmonicity to study the mechanism and thermodynamic properties of these reactions using a bare model (without CYP) and a pinene-CYP model. According to the potential energy surface and Boltzmann distribution for radical conformers, the major reaction products of CYP-catalyzed hydrogen abstraction from β-pinene are the doublet trans (53.4%) and doublet cis (46.1%) radical conformer at delta site. The formation of doublet cis/trans hydroxylated products released a total Gibbs free energy of about 48 kcal/mol. As for alpha pinene, the most stable radicals were trans-doublet (86.4%) and cis-doublet (13.6%) at epsilon sites, and their hydroxylation products released a total of ~50 kcal/mol Gibbs free energy. Our results highlight the likely C-H abstraction and oxygen rebounding sites accounting for the multi-state of CYP (doublet, quartet, and sextet spin states) and the formation of different conformers due to the presence of cis/trans allylic hydrogen in α-pinene and β-pinene molecules.

Common Reactivity and Properties of Heme Peroxidases: A DFT Study of Their Origin

Electronic structure calculations using the density-functional theory (DFT) have been performed to analyse the effect of water molecules and protonation on the heme group of peroxidases in different redox (ferric, ferrous, compounds I and II) and spin states. Shared geometries, spectroscopic properties at the Soret region, and the thermodynamics of peroxidases are discussed. B3LYP and M06-2X density functionals with different basis sets were employed on a common molecular model of the active site (Fe-centred porphine and proximal imidazole). Computed Gibbs free energies indicate that the corresponding aquo complexes are not thermodynamically stable, supporting the five-coordinate Fe(III) centre in native ferric peroxidases, with a water molecule located at a non-bonding distance. Protonation of the ferryl oxygen of compound II is discussed in terms of thermodynamics, Fe-O bond distances, and redox properties. It is demonstrated that this protonation is necessary to account for the experimental data, and computed Gibbs free energies reveal pK_a values of compound II about 8.5-9.0. Computation indicates that the general oxidative properties of peroxidase intermediates, as well as their reactivity towards water and protons and Soret bands, are mainly controlled by the iron porphyrin and its proximal histidine ligand.

The Mechanism of Biochemical NO-Sensing: Insights from Computational Chemistry

The binding of small gas molecules such as NO and CO plays a major role in the signaling routes of the human body. The sole NO-receptor in humans is soluble guanylyl cyclase (sGC) - a histidine-ligated heme protein, which, upon NO binding, activates a downstream signaling cascade. Impairment of NO-signaling is linked, among others, to cardiovascular and inflammatory diseases. In the present work, we use a combination of theoretical tools such as MD simulations, high-level quantum chemical calculations and hybrid QM/MM methods to address various aspects of NO binding and to elucidate the most likely reaction paths and the potential intermediates of the reaction. As a model system, the H-NOX protein from Shewanella oneidensis (So H-NOX) homologous to the NO-binding domain of sGC is used. The signaling route is predicted to involve NO binding to form a six-coordinate intermediate heme-NO complex, followed by relatively facile His decoordination yielding a five-coordinate adduct with NO on the distal side with possible isomerization to the proximal side through binding of a second NO and release of the first one. MD simulations show that the His sidechain can quite easily rotate outward into solvent, with this motion being accompanied in our simulations by shifts in helix positions that are consistent with this decoordination leading to significant conformational change in the protein.

In Silico simulation of Cytochrome P450-Mediated metabolism of aromatic amines: A case study of N-Hydroxylation

Aromatic amines, the widely used raw materials in industry, cause long-term exposure to human bodies. They can be metabolized by cytochrome P450 enzymes to form active electrophilic compounds, which will potentially react with nucleophilic DNA to exert carcinogenesis. The short lifetime and versatility of the oxidant (a high-valent iron (IV)-oxo species, compound I) of P450 enzymes prompts us to use theoretical methods to investigate the metabolism of aromatic amines. In this work, the density functional theory (DFT) has been employed to simulate the hydroxylation metabolism through H-abstraction and to calculate the activation energy of this reaction for 28 aromatic amines. The results indicate that the steric effects, inductive effects and conjugative effects greatly contribute to the metabolism activity of the chemicals. The further correlation reveals that the dissociation energy of -NH₂ (BDE_N-H) can successfully predict the time-consuming calculated activation energy (R² for aromatic and heteroaromatic amines are 0.93 and 0.86, respectively), so BDE_N-H can be taken as a key parameter to characterize the relative stability of aromatic amines in P450 enzymes and further to quickly assess their potential toxicity. The validation results prove such relationship has good statistical performance (q_cv² for aromatic and heteroaromatic amines are 0.95 and 0.90, respectively) and can be used to other aromatic amines in the application domain, greatly reducing computational cost and providing useful support for experimental research.

Computational Investigation of the Bisphenolic Drug Metabolism by Cytochrome P450: What Factors Favor Intramolecular Phenol Coupling

Intramolecular phenol coupling reactions of alkaloids can lead to active metabolites catalyzed by the mammalian cytochrome P450 enzyme (P450); however, the mechanistic knowledge of such an "unusual" process is lacking. This work performs density functional theory computations to reveal the P450-mediated metabolic pathway leading from R-reticuline to the morphine precursor salutaridine by exploring possible intramolecular phenol coupling mechanisms involving diradical coupling, radical addition, and electron transfer. The computed results show that the outer-sphere electron transfer with a high barrier (>20.0 kcal/mol) is unlikely to happen. However, for inter-sphere intramolecular phenol coupling, it reveals that intramolecular phenol coupling of R-reticuline proceeds via the diradical mechanism consecutively by compound I and protonated compound II of P450 rather than the radical addition mechanism. The existence of a much higher radical rebound barrier than that of H-abstraction in the quartet high-spin state can endow the R-reticuline phenoxy radical with a sufficient lifetime to enable intramolecular phenol coupling, while the H-abstraction/radical rebound mode with a negligible rebound barrier leading to phenol hydroxylation can only happen in the doublet low-spin state. Therefore, the ratio [coupling]/[hydroxylation] can be approximately reflected by the relative yield of the high-spin and low-spin H-abstraction by P450, which thus can provide a theoretical ratio of 16:1 for R-reticuline, which is in accordance with previous experimental results. Especially, the high rebound barrier of the phenoxy radical derived from the weak electron-donating ability of the phenoxy radical is revealed as an intrinsic nature. Therefore, the revealed intramolecular phenol coupling mechanism can be potentially extended to several other bisphenolic drugs to infer groups of unexpected metabolites in organisms.

NAST: Nonadiabatic Statistical Theory Package for Predicting Kinetics of Spin-Dependent Processes

We present a nonadiabatic statistical theory (NAST) package for predicting kinetics of spin-dependent processes, such as intersystem crossings, spin-forbidden unimolecular reactions, and spin crossovers. The NAST package can calculate the probabilities and rates of transitions between the electronic states of different spin multiplicities. Both the microcanonical (energy-dependent) and canonical (temperature-dependent) rate constants can be obtained. Quantum effects, including tunneling, zero-point vibrational energy, and reaction path interference, can be accounted for. In the limit of an adiabatic unimolecular reaction proceeding on a single electronic state, NAST reduces to the traditional transition state theory. Because NAST requires molecular properties at only a few points on potential energy surfaces, it can be applied to large molecular systems, used with accurate high-level electronic structure methods, and employed to study slow nonadiabatic processes. The essential NAST input data include the nuclear Hessian at the reactant minimum, as well as the nuclear Hessians, energy gradients, and spin-orbit coupling at the minimum energy crossing point (MECP) between two states. The additional computational tools included in the NAST package can be used to extract the required input data from the output files of electronic structure packages, calculate the effective Hessian at the MECP, and fit the reaction coordinate for more advanced NAST calculations. We describe the theory, its implementation, and three examples of application to different molecular systems.

Computational insight into biotransformation of halophenols by cytochrome P450: Mechanism and reactivity for epoxidation

Halophenols (XPs) have aroused great interests due to their high toxicity and low biodegradability. Previous experimental studies have shown that XPs can be catalytically transformed into epoxides and haloquinones by cytochrome P450 enzymes (CYPs). However, these metabolites have never been detected directly. Moreover, the effects of the reaction site and the type and number of halogen substituents on the biotransformation reactivity of halophenols still remain unknown. In this work, we performed density functional theory (DFT) calculations to simulate the CYP-mediated biotransformation of 36 XPs with mono-, di-, and tri-halogen (F, Cl, and Br) substitutions to unravel the mechanism and relevant kinetics of XPs epoxidation. The whole epoxidation process consists of initial rate-determining O-addition and subsequent ring-closure steps. The simulation results show that the epoxidation in low-spin (LS) state is kinetically preferred over that in high-spin (HS) state, and the formation of epoxide metabolite is strongly exothermic. For all XPs, the epoxidation reactivity follows the order of ortho/para O-addition > meta O-addition. Moreover, the O-addition with higher energy barriers roughly corresponds to chlorophenols and fluorophenols with more halogen atoms. Compared with dichlorophenols, the additional ortho-Cl substitution on trichlorophenols can slightly increase the energy barriers of meta O-addition. By contrast, the additional inclusion of an ortho-Cl to monochlorophenols enhances the meta O-addition reactivity of dichlorophenols. Overall, the present work clarifies the biotransformation routes of XPs to produce epoxides, and identifies the key factors affecting the epoxidation reactivity, which are beneficial in understanding comprehensively the metabolic fate and toxicity of XPs.

Precision Biotransformation of Emerging Pollutants by Human Cytochrome P450 Using Computational-Experimental Synergy: A Case Study of Tris(1,3-dichloro-2-propyl) Phosphate

Precision biotransformation is an envisioned strategy offering detailed insights into biotransformation pathways in real environmental settings using experimentally guided high-accuracy quantum chemistry. Emerging pollutants, whose metabolites are easily overlooked but may cause idiosyncratic toxicity, are important targets of such a strategy. We demonstrate here that complex metabolic reactions of tris(1,3-dichloro-2-propyl) phosphate (TDCIPP) catalyzed by human CYP450 enzymes can be mapped via a three-step synergy strategy: (i) screening the possible metabolites via high-throughout (moderate-accuracy) computations; (ii) analyzing the proposed metabolites in vitro by human liver microsomes and recombinant human CYP450 enzymes; and (iii) rationalizing the experimental data via precise mechanisms using high-level targeted computations. Through the bilateral dialogues from qualitative to semi-quantitative to quantitative levels, we show how TDCIPP metabolism especially by CYP3A4 generates bis(1,3-dichloro-2-propyl) phosphate (BDCIPP) as an O-dealkylation metabolite and bis(1,3-dichloro-2-propyl) 3-chloro-1-hydroxy-2-propyl phosphate (alcohol_β-dehalogen) as a dehalogenation/reduction metabolite via the initial rate-determining H-abstraction from αC- and βC-positions. The relative yield ratio [dehalogenation/reduction]/[O-dealkylation] is derived from the relative barriers of H-abstraction at the βC- and αC-positions by CYP3A4, estimated as 0.002 to 0.23, viz., an in vitro measured ratio of 0.04. Importantly, alcohol_β-dehalogen formation points to a new mechanism involving successive oxidation and reduction functions of CYP450, with its precursor aldehyde_β-dehalogen being a key intermediate detected by trapping assays and rationalized by computations. We conclude that the proposed three-step synergy strategy may meet the increasing challenge of elucidating biotransformation mechanisms of substantial synthesized organic compounds in the future.

Increased adverse effects during metabolic transformation of short-chain chlorinated paraffins by cytochrome P450: A theoretical insight into 1-chlorodecane

Short-chain chlorinated paraffins (SCCPs), frequently detected in human tissues or organs, can result in threat to human health by disturbing normal metabolism. However, their metabolism mechanisms and fates are largely unclear. Therefore, to better understand the impacts of SCCPs and their metabolites on the human health, the metabolic mechanism and kinetics of SCCPs by cytochrome P450 enzymes (CYPs) were explored using density functional theory employed 1-chlorodecane as a model SCCPs. The results show that 1-chlorodecane could be readily metabolized by CYPs, and the rate constant reaches up 42.3 s^-1 in human body. Dechlorination of 1-chlorodecane is unlikely to occur and hydroxylation is dominated via H-abstraction pathways, especially from the intermediate C atom of 1-chlorodecane. The toxicity assessments suggest that the two metabolites, 10-chloro-decan-5-ol and 1-chlorodecanol could exhibit higher bioaccumulation, carcinogenicity and more serious damage on cardiovascular system after the metabolism of 1-chlorodecane. To our knowledge, this is the first study from the viewpoint of theoretical analysis to explore the metabolism of typical SCCPs in human body. It may provide deep insight into the metabolic transformation mechanism of SCCPs and cause the concerns about the adverse effects of their metabolites in human body.

Effect of the Electron Density of the Heme Fe Atom on the Nature of Fe-O2 Bonding in Oxy Myoglobin

Mössbauer spectroscopy has been used to characterize oxygenated myoglobins (oxy Mbs) reconstituted with native and chemically modified ⁵⁷Fe-enriched heme cofactors with different electron densities of the heme Fe atom (ρ_Fe) and to elucidate the effect of a change in the ρ_Fe on the nature of the bond between heme Fe and oxygen (O₂), i.e., the Fe-O₂ bond, in the protein. Quadrupole splitting (ΔE_Q) was found to decrease with decreasing ρ_Fe, and the observed ρ_Fe-dependent ΔE_Q confirmed an increase in the contribution of the ferric-superoxide (Fe³⁺-O₂^-) form to the resonance hybrid of the Fe-O₂ fragment with decreasing ρ_Fe. These observations explicitly accounted for the lowering of O₂ affinity of the protein due to an increase in the O₂ dissociation rate and a decrease in the autoxidation reaction rate of oxy Mb through decreasing H⁺ affinity of the bound ligand with decreasing ρ_Fe. Therefore, the present study demonstrated the mechanism underlying the electronic control of O₂ affinity and the autoxidation of the protein through the heme electronic structure. Carbon monoxide (CO) adducts of reconstituted Mbs (CO-Mbs) were similarly characterized, and we found that the resonance between the two canonical forms of the Fe-CO fragment was also affected by a change in ρ_Fe. Thus, the nature of the Fe-ligand bond in the protein was found to be affected by the ρ_Fe.

Spin-inversion mechanisms in O2 binding to a model heme compound: A perspective from nonadiabatic wave packet calculations

Spin-inversion dynamics in O₂ binding to a model heme complex, which consisted of Fe(II)-porphyrin and imidazole, were studied using nonadiabatic wave packet dynamics calculations. We considered three active nuclear degrees of freedom in the dynamics, including the motions along the Fe-O distance, Fe-O-O angle, and Fe out-of-plane distance. Spin-free potential energy surfaces for the singlet, triplet, quintet, and septet states were developed using density functional theory calculations, and spin-orbit coupling elements were obtained from CASSCF-level electronic structure calculations. The spin-inversion mainly occurred between the singlet state and one of the triplet states due to large spin-orbit couplings and the contributions of other states were extremely small. The present quantum dynamics calculations suggested that the narrow crossing region model plays a dominant role in the O₂ binding dynamics. In addition, the one-dimensional Landau-Zener model underestimated the nonadiabatic transition probability.

Role of the Two Metals in the Active Sites of Heme Copper Oxidases-A Study of NO Reduction in cbb3 Cytochrome c Oxidase

The superfamily of heme copper oxidases reduces molecular oxygen or nitric oxide, and the active sites comprise a high-spin heme group (a₃ or b₃) and a non-heme metal (Cu_B or Fe_B). The cbb₃ C family of cytochrome c oxidases, with the high-spin heme b₃ and Cu_B in the active site, is a subfamily of the heme copper oxidases that can reduce both molecular oxygen, which is the main substrate, and nitric oxide. The mechanism for NO reduction in cbb₃ oxidase is studied here using hybrid density functional theory and compared to other cytochrome c oxidases (A and B families), with a high-spin heme a₃ and Cu_B in the active site, and to cytochrome c dependent NO reductase, with a high-spin heme b₃ and a non-heme Fe_B in the active site. It is found that the reaction mechanism and the detailed reaction energetics of the cbb₃ oxidases are not similar to those of cytochrome c dependent NO reductase, which has the same type of high-spin heme group but a different non-heme metal. This is in contrast to earlier expectations. Instead, the NO reduction mechanism in cbb₃ oxidases is very similar to that in the other cytochrome c oxidases, with the same non-heme metal, Cu_B, and is independent of the type of high-spin heme group. The conclusion is that the type of non-heme metal (Cu_B or Fe_B) in the active site of the heme copper oxidases is more important for the reaction mechanisms than the type of high-spin heme, at least for the NO reduction reaction. The reason is that the proton-coupled reduction potentials of the active site cofactors determine the energetics for the NO reduction reaction, and they depend to a larger extent on the non-heme metal. Observed differences in NO reduction reactivity among the various cytochrome c oxidases may be explained by differences outside the BNC, affecting the rate of proton transfer, rather than in the BNC itself.

Spin-inversion mechanisms in O2 binding to a model heme complex revisited by density function theory calculations

Spin-inversion mechanisms in O₂ binding to a model heme complex, consisting of Fe(II)-porphyrin and imidazole, were investigated using density-functional theory calculations. First, we applied the recently proposed mixed-spin Hamiltonian method to locate spin-inversion structures between different total spin multiplicities. Nine spin-inversion structures were successfully optimized for the singlet-triplet, singlet-quintet, triplet-quintet, and quintet-septet spin-inversion processes. We found that the singlet-triplet spin-inversion points are located around the potential energy surface region at short Fe-O distances, whereas the singlet-quintet and quintet-septet spin-inversion points are located at longer Fe-O distances. This suggests that both narrow and broad crossing models play roles in O₂ binding to the Fe-porphyrin complex. To further understand spin-inversion mechanisms, we performed on-the-fly Born-Oppenheimer molecular dynamics calculations. The reaction coordinates, which are correlated to the spin-inversion dynamics between different spin multiplicities, are also discussed.

Constructing Spin-Adiabatic States for the Modeling of Spin-Crossing Reactions. I. A Shared-Orbital Implementation

In the modeling of spin-crossing reactions, it has become popular to directly explore the spin-adiabatic surfaces. Specifically, through constructing spin-adiabatic states from a two-state Hamiltonian (with spin-orbit coupling matrix elements) at each geometry, one can readily employ advanced geometry optimization algorithms to acquire a "transition state" structure, where the spin crossing occurs. In this work, we report the implementation of a fully-variational spin-adiabatic approach based on Kohn-Sham density functional theory spin states (sharing the same set of molecular orbitals) and the Breit-Pauli one-electron spin-orbit operator. For three model spin-crossing reactions [predissociation of N₂O, singlet-triplet conversion in CH₂, and CO addition to Fe(CO)₄], the spin-crossing points were obtained. Our results also indicated the Breit-Pauli one-electron spin-orbit coupling can vary significantly along the reaction pathway on the spin-adiabatic energy surface. On the other hand, due to the restriction that low-spin and high-spin states share the same set of molecular orbitals, the acquired spin-adiabatic energy surface shows a cusp (i.e. a first-order discontinuity) at the crossing point, which prevents the use of standard geometry optimization algorithms to pinpoint the crossing point. An extension with this restriction removed is being developed to achieve the smoothness of spin-adiabatic surfaces.

Computational Biotransformation Profile of Emerging Phenolic Pollutants by Cytochromes P450: Phenol-Coupling Mechanism

Phenols are ubiquitous environmental pollutants, whose biotransformation involving phenol coupling catalyzed by cytochromes P450 may produce more lipophilic and toxic metabolites. Density functional theory (DFT) computations were performed to explore the debated phenol-coupling mechanisms, taking triclosan as a model substrate. We find that a diradical pathway facilitated by compound I and protonated compound II of P450 is favored vs alternative radical addition or electron-transfer mechanisms. The identified diradical coupling resembles a "two-state reactivity" from compound I characterized by significantly high rebound barriers of the phenoxy radicals, which can be formulated into three equations for calculating the ratio [coupling]/[hydroxylation]. A higher barrier for rebound than for H-abstraction in high-spin triclosan can facilitate the phenoxy radical dissociation and thus enable phenol coupling, while H-abstraction/radical rebound causing phenol hydroxylation via minor rebound barriers mostly occurs via the low-spin state. Therefore, oxidation of triclosan by P450 fits the first equation with a ratio [coupling]/[hydroxylation] of 1:4, consistent with experimental data indicating different extents of triclosan coupling (6-40%). The high rebound barrier of phenoxy radicals, as a key for the mechanistic identification of phenol coupling vs hydroxylation, originates from their weak electron donor ability due to spin aromatic delocalization. We envision that the revealed mechanism can be extended to the cross-coupling reactions between different phenolic pollutants, and the coupling reactions of several other aromatic pollutants, to infer unknown metabolites.

The mechanism for oxygen reduction in the C family cbb3 cytochrome c oxidases - Implications for the proton pumping stoichiometry

Cytochrome c oxidases (CcOs) couple the exergonic reduction of molecular oxygen to proton pumping across the membrane in which they are embedded, thereby conserving a significant part of the free energy. The A family CcOs are known to pump four protons per oxygen molecule, while there is no consensus regarding the proton pumping stoichiometry for the C family cbb₃ oxidases. Hybrid density functional theory is used here to investigate the catalytic mechanism for oxygen reduction in cbb₃ oxidases. A surprising result is that the barrier for O O bond cleavage at the mixed valence reduction level seems to be too high compared to the overall reaction rate of the enzyme. It is therefore suggested that the O O bond is cleaved only after the first proton coupled reduction step, and that this reduction step most likely is not coupled to proton pumping. Furthermore, since the cbb₃ oxidases have only one proton channel leading to the active site, it is proposed that the activated E_H intermediate, suggested to be responsible for proton pumping in one of the reduction steps in the A family, cannot be involved in the catalytic cycle for cbb₃, which results in the lack of proton pumping also in the E to R reduction step. In summary, the calculations indicate that only two protons are pumped per oxygen molecule in cbb₃ oxidases. However, more experimental information on this divergent enzyme is needed, e.g. whether the flow of electrons resembles that in the other more well-studied CcO families.

Water binding to FeIII hemes studied in a cooled ion trap: characterization of a strong 'weak' ligand

The interaction of a water molecule with ferric heme-iron protoporphyrin ([PP Fe^III]⁺) has been investigated in the gas phase in an ion trap and studied theoretically by density functional theory. It is found that the interaction of water with ferric heme leads to a stable [PP-Fe^III-H₂O]⁺ complex in the intermediate spin state (S = 3/2), in the same state as its unligated [PP-Fe^III]⁺ homologue, without spin crossing during water attachment. Using the Van't Hoff equation, the reaction enthalpy for the formation of a Fe-OH₂ bond has been determined for [PP-Fe^III-H₂O]⁺ and [PP-Fe^III-(H₂O)₂]⁺. The corrected binding energy for a single Fe-H₂O bond is -12.2 ± 0.6 kcal mol^-1, while DFT calculations at the OPBE level yield -11.7 kcal mol^-1. The binding energy of the second ligation yielding a six coordinated Fe^III atom is decreased with a bond energy of -9 ± 0.9 kcal mol^-1, well reproduced by calculations as -7.1 kcal mol^-1. However, calculations reveal features of a weaker bond type, such as a rather long Fe-O bond with 2.28 Å for the [PP-Fe^III-H₂O]⁺ complex and the absence of a spin change by complexation. Thus despite a strong bond with H₂O, the Fe^III atom does not show, through theoretical modelling, a strong acceptor character in its half filled 3d_z² orbital. It is also observed that the binding properties of H₂O to hemes seem strikingly specific to ferric heme and we have shown, experimentally and theoretically, that the affinity of H₂O for protonated heme [H PP-Fe]⁺, an intermediate between Fe^III and Fe^II, is strongly reduced compared to that for ferric heme.

A proposed method of enantioselectivity analysis for residual chiral PCBs in gas chromatography

Polychlorinated biphenyls (PCBs) are harmful and persistent organic pollutants. The influence of chiral PCBs acts mainly on the different enantiomer fraction, bioaccumulation and even degradation in environmental media. Solvents and temperatures existed almost everywhere during the analysis of extraction, purification, concentration and detection, which was often underestimated in previous studies. In our study, the configuration stability of the chiral PCBs was examined from solvent and temperature aspects. Transformation phenomena for the analytic stereoisomer monomers of PCB45, PCB95, and PCB149 affected by temperature were observed. We demonstrated that higher inlet temperatures could increase the sensibility for the low-concentration environmental samples, resulting in isomerization of chiral PCBs. Real rice samples were used to verify our analysis method. Combined with density functional theory, we verified the mechanism of isomer conversion with various numbers and sites of the -Cl substituent. PCBs with tetra-ortho substituents (2, 2', 6, 6') were relatively stable and showed the highest rotational barriers (E_a) at approximate 240 kJ mol^-1. Others with trio-ortho substituents (2, 2', 6/6') showed E_a from 170 to 190 kJ mol^-1, whose enantiomeric fractions would be affected by temperature during the analysis process for environmental detection. The method we developed was a promising means to understand the mechanism of isomerization and to predict stabilities of chiral PCBs.

Active Site Midpoint Potentials in Different Cytochrome c Oxidase Families: A Computational Comparison

Cytochrome c oxidase (C cO) is the terminal enzyme in the respiratory electron transport chain, reducing molecular oxygen to water. The binuclear active site in C cO comprises a high-spin heme associated with a Cu_B complex and a redox active tyrosine. The electron transport in the respiratory chain is driven by increasing midpoint potentials of the involved cofactors, resulting in a release of free energy, which is stored by coupling the electron transfer to proton translocation across a membrane, building up an electrochemical gradient. In this context, the midpoint potentials of the active site cofactors in the C cOs are of special interest, since they determine the driving forces for the individual oxygen reduction steps and thereby affect the efficiency of the proton pumping. It has been difficult to obtain useful information on some of these midpoint potentials from experiments. However, since each of the reduction steps in the catalytic cycle of oxygen reduction to water corresponds to the formation of an O-H bond, they can be calculated with a reasonably high accuracy using quantum chemical methods. From the calculated O-H bond strengths, the proton-coupled midpoint potentials of the active site cofactors can be estimated. Using models representing the different families of C cO's (A, B, and C), the calculations give midpoint potentials that should be relevant during catalytic turnover. The calculations also suggest possible explanations for why some experimentally measured potentials deviate significantly from the calculated ones, i.e., for Cu_B in all oxidase families, and for heme b₃ in the C family.

Metabolic Mechanism of Aryl Phosphorus Flame Retardants by Cytochromes P450: A Combined Experimental and Computational Study on Triphenyl Phosphate

Understanding metabolic mechanisms is critical and remains a difficult task in the risk assessment of emerging pollutants. Triphenyl phosphate (TPHP), a widely used aryl phosphorus flame retardant (aryl-PFR), has been frequently detected in the environment, and its major metabolite was considered as diphenyl phosphate (DPHP). However, knowledge of the mechanism for TPHP leading to DPHP and other metabolites is lacking. Our in vitro study shows that TPHP is metabolized into its diester metabolite DPHP and mono- and dihydroxylated metabolites by cytochromes P450 (CYP) in human liver microsomes, while CYP1A2 and CYP2E1 isoforms are mainly involved in such processes. Molecular docking gives the conformation for TPHP binding with the active species Compound I (an iron IV-oxo heme cation radical) in specific CYP isoforms, showing that the aromatic ring of TPHP is likely to undergo metabolism. Quantum chemical calculations have shown that the dominant reaction channel is the O-addition of Compound I onto the aromatic ring of TPHP, followed by a hydrogen-shuttle mechanism leading to ortho-hydroxy-TPHP as the main monohydroxylated metabolite; the subsequent H-abstraction-OH-rebound reaction acting on ortho-hydroxy-TPHP yields the meta- and ipso-position quinol intermediates, while the former of which can be metabolized into dihydroxy-TPHP by fast protonation, and the latter species needs to go through type-I ipso-substitution and fast protonation to be evolved into DPHP. We envision that the identified mechanisms may give inspiration for studying the metabolism of several other aryl-PFRs by CYP.

Computational Insight into the Activation Mechanism of Carcinogenic N'-Nitrosonornicotine (NNN) Catalyzed by Cytochrome P450

Tobacco-specific N'-nitrosonornicotine (NNN), a genotoxic nitrosamine classified as Group 1 carcinogen, is also present in atmospheric particulate matter and has even been detected as a new disinfection byproduct in wastewaters. NNN generally requires metabolic activation by cytochrome P450 enzymes to exert its genotoxicity, but the respective biotransformation pathways have not been described in detail. In this work, we performed density functional theory (DFT) calculations to unravel possible NNN activation pathways including α-hydroxylation, β-hydroxylation, pyridine N-oxidation, and norcotinine formation. The results reveal an initial rate-determining H_α-atom abstraction step for α-hydroxylation, followed by an unexpected kinetic competition between denitrosation and OH rebound, leading to ( iso-)myosmine as a detoxified product and α-hydroxyNNNs as the precursor of carcinogenic diazohydroxides, respectively. Further detoxification routes are given by β-hydroxylation with relative high reaction barrier and N-oxidation with comparable barrier to the toxifying α-hydroxylation. Moreover, we show for the first time how norcotinine can be generated as a minor NNN metabolite that is formed from iso-myosmine through a unique porphyrin-assisted H atom 1,2-transfer mechanism. These results demonstrate that the carcinogenic potential of NNN is subject to a kinetic competition between activating and deactivating metabolic routes, and identify respective biomarkers to inform about the individual risk associated with NNN exposure.

Mechanisms for enzymatic reduction of nitric oxide to nitrous oxide - A comparison between nitric oxide reductase and cytochrome c oxidase

Cytochrome c oxidases (CcO) reduce O₂ to H₂O in the respiratory chain of mitochondria and many aerobic bacteria. In addition, some species of CcO can also reduce NO to N₂O and water while others cannot. Here, the mechanism for NO-reduction in CcO is investigated using quantum mechanical calculations. Comparison is made to the corresponding reaction in a "true" cytochrome c-dependent NO reductase (cNOR). The calculations show that in cNOR, where the reduction potentials are low, the toxic NO molecules are rapidly reduced, while the higher reduction potentials in CcO lead to a slower or even impossible reaction, consistent with experimental observations. In both enzymes the reaction is initiated by addition of two NO molecules to the reduced active site, forming a hyponitrite intermediate. In cNOR, N₂O can then be formed using only the active-site electrons. In contrast, in CcO, one proton-coupled reduction step most likely has to occur before N₂O can be formed, and furthermore, proton transfer is most likely rate-limiting. This can explain why different CcO species with the same heme a₃-Cu active site differ with respect to NO reduction efficiency, since they have a varying number and/or properties of proton channels. Finally, the calculations also indicate that a conserved active site valine plays a role in reducing the rate of NO reduction in CcO.

Molecular Mechanism of Alternative P450-Catalyzed Metabolism of Environmental Phenolic Endocrine-Disrupting Chemicals

Understanding the bioactivation mechanisms to predict toxic metabolites is critical for risk assessment of phenolic endocrine-disrupting chemicals (EDCs). One mechanism involves ipso-substitution, which may contribute to the total turnover of phenolic EDCs, yet the detailed mechanism and its relationship with other mechanisms are unknown. We used density functional theory to investigate the P450-catalyzed ipso-substitution mechanism of the prominent xenoestrogen bisphenol A. The ipso-substitution proceeds via H-abstraction from bisphenol A by Compound I, followed by essentially barrierless OH-rebound onto the ipso-position forming a quinol, which can spontaneously decompose into the carbocation and hydroquinone. This carbocation can further evolve into the highly estrogenic hydroxylated and dimer-type metabolites. The H-abstraction/OH-rebound reaction mechanism has been verified as a general reaction mode for many other phenolic EDCs, such as bisphenol analogues, alkylphenols and chlorophenols. The identified mechanism enables us to effectively distinguish between type I (eliminating-substituent as anion) and type II (eliminating-substituent as cation) ipso-substitution in various phenolic EDCs. We envision that the identified pathways will be applicable for prediction of metabolites from phenolic EDCs whose fate are affected by this alternative type of P450 reactivity, and accordingly enable the screening of these metabolites for endocrine-disrupting activity.

First Principles Calculation of the Reaction Rates for Ligand Binding to Myoglobin: The Cases of NO and CO

Ligand binding by proteins is among the most fundamental processes in nature. Among these processes the binding of small gas molecules, such as O₂ , CO and NO to heme proteins has traditionally received vivid interest, which was further boosted by their recently recognized significant role in gas sensing in the body. At the heart of the binding of these ligands to the heme group is the spinforbidden reaction between high-spin iron(II) and the ligand yielding a low-spin adduct. We use computational means to address the complete mechanism of CO and NO binding by myoglobin. Considering that it involves several steps occurring on different time scales, molecular dynamics simulations were performed to address the diffusion of the ligand through the enzyme, and DFT calculations in combination with statistical rate calculation to investigate the spin-forbidden reaction. The calculations yielded rate constants in qualitative agreement with experiments and revealed that the bottleneck of NO and CO binding is different; for NO, diffusion was found to be rate-limiting, whereas for CO, the spin-forbidden step is the slowest.

The dioxygen adducts of iron and manganese porphyrins: electronic structure and binding energy

The electronic structures of adducts of O₂ and metal porphyrins were thoroughly investigated by highly accurate DMRG-CASPT2.

Ultrafast CO Kinetics in Heme Proteins: Adiabatic Ligand Binding and Heavy Atom Tunneling

The ultrafast kinetics of CO rebinding to carbon monoxide oxidation activator protein (ChCooA) are measured over a wide temperature range and compared with the kinetics of CO binding in other heme systems such as myoglobin (Mb) and hemoglobin (Hb). The Arrhenius prefactor for CO binding to ChCooA and protoheme (∼10¹¹ s^-1) is similar to what is found for spin-allowed NO binding to heme proteins and is several orders of magnitude larger than the prefactor of Mb and Hb (∼10⁹ s^-1). This indicates that the CO binding reaction is adiabatic, in contrast to the commonly held view that it is nonadiabatic due to spin-forbidden (ΔS = 2) selection rules. Under the adiabatic condition, entropic factors, rather than spin-selection rules, are the source of the reduced Arrhenius prefactors associated with CO binding in Mb and Hb. The kinetic response of ChCooA-CO is nonexponential at all temperatures, including 298 K, and is described quantitatively using a distribution of enthalpic rebinding barriers associated with heterogeneity in the heme doming conformation. Above the solvent glass transition (T_g ∼ 180 K), the rebinding progress slows as temperature increases, and this is ascribed to an evolution of the distribution toward increased heme doming and larger enthalpic barriers. Between T_g and ∼60 K, the nonexponential rebinding slows down as the temperature is lowered and the survival fraction follows the predictions expected for a quenched barrier distribution. Below ∼60 K the rebinding kinetics do not follow these predictions unless quantum mechanical tunneling along the heme doming coordinate is also included as an active channel for CO binding.

The mechanism for oxygen reduction in cytochrome c dependent nitric oxide reductase (cNOR) as obtained from a combination of theoretical and experimental results

Bacterial NO-reductases (NOR) belong to the heme-copper oxidase (HCuO) superfamily, in which most members are O₂-reducing, proton-pumping enzymes. This study is one in a series aiming to elucidate the reaction mechanisms of the HCuOs, including the mechanisms for cellular energy conservation. One approach towards this goal is to compare the mechanisms for the different types of HCuOs, cytochrome c oxidase (CcO) and NOR, reducing the two substrates O₂ and NO. Specifically in this study, we describe the mechanism for oxygen reduction in cytochrome c dependent NOR (cNOR). Hybrid density functional calculations were performed on large cluster models of the cNOR binuclear active site. Our results are used, together with published experimental information, to construct a free energy profile for the entire catalytic cycle. Although the overall reaction is quite exergonic, we show that during the reduction of molecular oxygen in cNOR, two of the reduction steps are endergonic with high barriers for proton uptake, which is in contrast to oxygen reduction in CcO, where all reduction steps are exergonic. This difference between the two enzymes is suggested to be important for their differing capabilities for energy conservation. An additional result from this study is that at least three of the four reduction steps are initiated by proton transfer to the active site, which is in contrast to CcO, where electrons always arrive before the protons to the active site. The roles of the non-heme metal ion and the redox-active tyrosine in the active site are also discussed.

Multireference Perturbation Theory with Cholesky Decomposition for the Density Matrix Renormalization Group

We present a second-order N-electron valence state perturbation theory (NEVPT2) based on a density matrix renormalization group (DMRG) reference wave function that exploits a Cholesky decomposition of the two-electron repulsion integrals (CD-DMRG-NEVPT2). With a parameter-free multireference perturbation theory approach at hand, the latter allows us to efficiently describe static and dynamic correlation in large molecular systems. We demonstrate the applicability of CD-DMRG-NEVPT2 for spin-state energetics of spin-crossover complexes involving calculations with more than 1000 atomic basis functions. We first assess, in a study of a heme model, the accuracy of the strongly and partially contracted variant of CD-DMRG-NEVPT2 before embarking on resolving a controversy about the spin ground state of a cobalt tropocoronand complex.