Does compound I vary significantly between isoforms of cytochrome P450?

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.

An In Crystallo Reaction with an Engineered Cytochrome P450 Peroxygenase

The cytochrome P450 monooxygenases (CYPs) are a class of heme-thiolate enzymes that insert oxygen into unactivated C-H bonds. These enzymes can be converted into peroxygenases via protein engineering, which enables their activity to occur using hydrogen peroxide (H2 O2 ) without the requirement for additional nicotinamide co-factors or partner proteins. Here, we demonstrate that soaking crystals of an engineered P450 peroxygenase with H2 O2 enables the enzymatic reaction to occur within the crystal. Crystals of the designed P450 peroxygenase, the T252E mutant of CYP199A4, in complex with 4-methoxybenzoic acid were soaked with different concentrations of H2 O2 for varying times to initiate the in crystallo O-demethylation reaction. Crystal structures of T252E-CYP199A4 showed a distinct loss of electron density that was consistent with the O-demethylated metabolite, 4-hydroxybenzoic acid. A new X-ray crystal structure of this enzyme with the 4-hydroxybenzoic acid product was obtained to enable comparison alongside the existing substrate-bound structure. The visualisation of enzymatic catalysis in action is challenging in structural biology and the ability to initiate the reactions of P450 enzymes, in crystallo by simply soaking crystals with H2 O2 will enable new structural biology methods and techniques to be applied to study their mechanism of action.

Atomic-level design of metalloenzyme-like active pockets in metal-organic frameworks for bioinspired catalysis

Natural metalloenzymes with astonishing reaction activity and specificity underpin essential life transformations. Nevertheless, enzymes only operate under mild conditions to keep sophisticated structures active, limiting their potential applications. Artificial metalloenzymes that recapitulate the catalytic activity of enzymes can not only circumvent the enzymatic fragility but also bring versatile functions into practice. Among them, metal-organic frameworks (MOFs) featuring diverse and site-isolated metal sites and supramolecular structures have emerged as promising candidates for metalloenzymes to move toward unparalleled properties and behaviour of enzymes. In this review, we systematically summarize the significant advances in MOF-based metalloenzyme mimics with a special emphasis on active pocket engineering at the atomic level, including primary catalytic sites and secondary coordination spheres. Then, the deep understanding of catalytic mechanisms and their advanced applications are discussed. Finally, a perspective on this emerging frontier research is provided to advance bioinspired catalysis.

Insights into Cytochrome P450 Enzyme Catalyzed Defluorination of Aromatic Fluorides

Density functional calculations establish a novel mechanism of aromatic defluorination by P450 Compound I. This is achieved via either an initial epoxide intermediate or through a 1,2-fluorine shift in an electrophilic intermediate, which highlights that the P450s can defluorinate fluoroarenes. However, in the absence of a proton donor a strong Fe-F bond can be obtained as shown from the calculations.

Melatonin Activation by Human Cytochrome P450 Enzymes: A Comparison between Different Isozymes

Cytochrome P450 enzymes in the human body play a pivotal role in both the biosynthesis and the degradation of the hormone melatonin. Melatonin plays a key role in circadian rhythms in the body, but its concentration is also linked to mood fluctuations as well as emotional well-being. In the present study, we present a computational analysis of the binding and activation of melatonin by various P450 isozymes that are known to yield different products and product distributions. In particular, the P450 isozymes 1A1, 1A2, and 1B1 generally react with melatonin to provide dominant aromatic hydroxylation at the C6-position, whereas the P450 2C19 isozyme mostly provides O-demethylation products. To gain insight into the origin of these product distributions of the P450 isozymes, we performed a comprehensive computational study of P450 2C19 isozymes and compared our work with previous studies on alternative isozymes. The work covers molecular mechanics, molecular dynamics and quantum mechanics approaches. Our work highlights major differences in the size and shape of the substrate binding pocket amongst the different P450 isozymes. Consequently, substrate binding and positioning in the active site varies substantially within the P450 isozymes. Thus, in P450 2C19, the substrate is oriented with its methoxy group pointing towards the heme, and therefore reacts favorably through hydrogen atom abstraction, leading to the production of O-demethylation products. On the other hand, the substrate-binding pockets in P450 1A1, 1A2, and 1B1 are tighter, direct the methoxy group away from the heme, and consequently activate an alternative site and lead to aromatic hydroxylation instead.

The Oxidation of Oxygen and Sulfur-Containing Heterocycles by Cytochrome P450 Enzymes

The cytochrome P450 (CYP) superfamily of monooxygenase enzymes play important roles in the metabolism of molecules which contain heterocyclic, aromatic functional groups. Here we study how oxygen- and sulfur-containing heterocyclic groups interact with and are oxidized using the bacterial enzyme CYP199A4. This enzyme oxidized both 4-(thiophen-2-yl)benzoic acid and 4-(thiophen-3-yl)benzoic acid almost exclusively via sulfoxidation. The thiophene oxides produced were activated towards Diels-Alder dimerization after sulfoxidation, forming dimeric metabolites. Despite X-ray crystal structures demonstrating that the aromatic carbon atoms of the thiophene ring were located closer to the heme than the sulfur, sulfoxidation was still favoured with 4-(thiophen-3-yl)benzoic acid. These results highlight a preference of this cytochrome P450 enzyme for sulfoxidation over aromatic hydroxylation. Calculations predict a strong preference for homodimerization of the enantiomers of the thiophene oxides and the formation of a single major product, in broad agreement with the experimental data. 4-(Furan-2-yl)benzoic acid was oxidized to 4-(4'-hydroxybutanoyl)benzoic acid using a whole-cell system. This reaction proceeded via a γ-keto-α,β-unsaturated aldehyde species which could be trapped in vitro using semicarbazide to generate a pyridazine species. The combination of the enzyme structures, the biochemical data and theoretical calculations provides detailed insight into the formation of the metabolites formed from these heterocyclic compounds.

Unraveling the Major Differences between the Trinuclear Cyclopentadienylmetal Carbonyl Chemistry of Cobalt and That of Nickel-A Theoretical Study

The geometries and energetics of the trinuclear cyclopentadienylmetal carbonyls Cp₃M₃(CO)_n (Cp = η⁵-C₅H₅); M = Co, Ni; n = 3, 2, 1, 0) have been investigated by density functional theory. The cobalt and nickel systems are found to be rather different owing to the different electronic configurations of the metal atoms. For cobalt, the small calculated energy separation of 5.0 kcal/mol between the two lowest-energy singlet Cp₃Co₃(μ₃-CO)(μ-CO)₂ and Cp₃Co₃(μ-CO)₃ tricarbonyl structures accounts for the experimental results of both isomers as stable species that can be isolated and structurally characterized by X-ray crystallography. The corresponding Cp₃Ni₃(CO)₃ species in the nickel system are predicted not to be viable owing to exothermic CO dissociation to give the experimentally observed very stable Cp₃Ni₃(μ-CO)₂, which is found to be the lowest-energy isomer by a substantial margin of ∼25 kcal/mol. In all of the low-energy Cp₃M₃(CO)_n (n = 2, 1) structures, including that of the experimentally known triplet spin state Cp₃Co₃(μ₃-CO)₂, all of the carbonyl groups are face-bridging or face-semi-bridging μ₃-CO groups bonded to all three metal atoms of the M₃ triangle. In the lowest-energy carbonyl-free Cp₃M₃ (M = Co, Ni) structures, agostic C-H-M interactions are found using hydrogens of the Cp rings. In addition, the lowest-energy Cp₃Ni₃ is the only structure among all of the low-energy Cp₃M₃(CO)_n (M = Co, Ni; n = 3, 2, 1, 0) structures in which each Cp ring is a bridging rather than terminal ligand.

Melatonin Activation by Cytochrome P450 Isozymes: How Does CYP1A2 Compare to CYP1A1?

Cytochrome P450 enzymes are versatile enzymes found in most biosystems that catalyze mono-oxygenation reactions as a means of biosynthesis and biodegradation steps. In the liver, they metabolize xenobiotics, but there are a range of isozymes with differences in three-dimensional structure and protein chain. Consequently, the various P450 isozymes react with substrates differently and give varying product distributions. To understand how melatonin is activated by the P450s in the liver, we did a thorough molecular dynamics and quantum mechanics study on cytochrome P450 1A2 activation of melatonin forming 6-hydroxymelatonin and N-acetylserotonin products through aromatic hydroxylation and O-demethylation pathways, respectively. We started from crystal structure coordinates and docked substrate into the model, and obtained ten strong binding conformations with the substrate in the active site. Subsequently, for each of the ten substrate orientations, long (up to 1 μs) molecular dynamics simulations were run. We then analyzed the orientations of the substrate with respect to the heme for all snapshots. Interestingly, the shortest distance does not correspond to the group that is expected to be activated. However, the substrate positioning gives insight into the protein residues it interacts with. Thereafter, quantum chemical cluster models were created and the substrate hydroxylation pathways calculated with density functional theory. These relative barrier heights confirm the experimental product distributions and highlight why certain products are obtained. We make a detailed comparison with previous results on CYP1A1 and identify their reactivity differences with melatonin.

How does multiple substrate binding lead to substrate inhibition of CYP2D6 metabolizing dextromethorphan? A theoretical study

CYP2D6 is one of the most important metalloenzymes involved in the biodegradation of many drug molecules in the human body. It has been found that multiple substrate binding can lead to substrate inhibition of CYP2D6 metabolizing dextromethorphan (DM), but the corresponding theoretical mechanism is rarely reported. Therefore, we chose DM as the probe and performed molecular dynamics simulations and quantum mechanical calculations on CYP2D6-DM systems to investigate the mechanism of how the multiple substrate binding leads to the substrate inhibition of CYP2D6 metabolizing substrates. According to our results, three gate residues (Arg221, Val374, and Phe483) for the catalytic pocket are determined. We also found that the multiple substrate binding can lead to substrate inhibition by reducing the stability of CYP2D6 binding DM and increasing the reactive activation energy of the rate-determining step. Our findings would help to understand the substrate inhibition of CYP2D6 metabolizing the DM and enrich the knowledge of the drug-drug interactions for the cytochrome P450 superfamily.

Investigating the Active Oxidants Involved in Cytochrome P450 Catalyzed Sulfoxidation Reactions

Cytochrome P450 (CYP) heme-thiolate monooxygenases catalyze the hydroxylation of the C-H bonds of organic molecules. This reaction is initiated by a ferryl-oxo heme radical cation (Cpd I). These enzymes can also catalyze sulfoxidation reactions and the ferric-hydroperoxy complex (Cpd 0) and the Fe(III)-H₂ O₂ complex have been proposed as alternative oxidants for this transformation. To investigate this, the oxidation of 4-alkylthiobenzoic acids and 4-methoxybenzoic acid by the CYP199A4 enzyme from Rhodopseudomonas palustris HaA2 was compared using both monooxygenase and peroxygenase pathways. By examining mutants at the mechanistically important, conserved acid alcohol-pair (D251N, T252A and T252E) the relative amounts of the reactive intermediates that would form in these reactions were disturbed. Substrate binding and X-ray crystal structures helped to understand changes in the activity and enabled an attempt to evaluate whether multiple oxidants can participate in these reactions. In peroxygenase reactions the T252E mutant had higher activity towards sulfoxidation than O-demethylation but in the monooxygenase reactions with the WT enzyme the activity of both reactions was similar. The peroxygenase activity of the T252A mutant was greater for sulfoxidation reactions than the WT enzyme, which is the reverse of the activity changes observed for O-demethylation. The monooxygenase activity and coupling efficiency of sulfoxidation and oxidative demethylation were reduced by similar degrees with the T252A mutant. These observations infer that while Cpd I is required for O-dealkylation, another oxidant may contribute to sulfoxidation. Based on the activity of the CYP199A4 mutants it is proposed that this is the Fe(III)-H₂ O₂ complex which would be more abundant in the peroxide-driven reactions.

Spin state dependent peroxidase activity of heme bound amyloid β peptides relevant to Alzheimer's disease

The colocalization of heme rich deposits in the senile plaque of Aβ in the cerebral cortex of the Alzheimer's disease (AD) brain along with altered heme homeostasis and heme deficiency symptoms in AD patients has invoked the association of heme in AD pathology. Heme bound Aβ complexes, depending on the concentration of the complex or peptide to heme ratio, exhibit an equilibrium between a high-spin mono-His bound peroxidase-type active site and a low-spin bis-His bound cytochrome b type active site. The high-spin heme-Aβ complex shows higher peroxidase activity than free heme, where compound I is the reactive oxidant. It is also capable of oxidizing neurotransmitters like serotonin in the presence of peroxide, owing to the formation of compound I. The low-spin bis-His heme-Aβ complex on the other hand shows enhanced peroxidase activity relative to high-spin heme-Aβ. It reacts with H₂O₂ to produce two stable intermediates, compound 0 and compound I, which are characterized by absorption, EPR and resonance Raman spectroscopy. The stability of compound I of low-spin heme-Aβ is accountable for its enhanced peroxidase activity and oxidation of the neurotransmitter serotonin. The effect of the second sphere Tyr10 residue of Aβ on the formation and stability of the intermediates of low-spin heme-Aβ has also been investigated. The higher stability of compound I for low-spin heme-Aβ is likely due to H-bonding interactions involving Tyr10 in the distal pocket.

Mechanism of Melatonin Metabolism by CYP1A1: What Determines the Bifurcation Pathways of Hydroxylation versus Deformylation?

Melatonin, a widely applied cosmetic active ingredient, has a variety of uses as a skin protector through antioxidant and anti-inflammatory functions as well as giving the body UV-induced defenses and immune system support. In the body, melatonin is synthesized from a tryptophan amino acid in a cascade of reactions, but as melatonin is toxic at high concentrations, it is metabolized in the human skin by the cytochrome P450 enzymes. The P450s are diverse heme-based mono-oxygenases that catalyze oxygen atom-transfer processes that trigger metabolism and detoxification reactions in the body. In the catalytic cycle of the P450s, a short-lived high-valent iron(IV)-oxo heme cation radical is formed that has been proposed to be the active oxidant. How and why it activates melatonin in the human body and what the origin of the product distributions is, are unknown. This encouraged us to do a detailed computational study on a typical human P450 isozyme, namely CYP1A1. We initially did a series of molecular dynamics simulations with substrate docked into several orientations. These simulations reveal a number of stable substrate-bound positions in the active site, which may lead to differences in substrate activation channels. Using tunneling analysis on the full protein structures, we show that two of the four binding conformations lead to open substrate-binding pockets. As a result, in these open pockets, the substrate is not tightly bound and can escape back into the solution. In the closed conformations, in contrast, the substrate is mainly oriented with the methoxy group pointing toward the heme, although under a different angle. We then created large quantum cluster models of the enzyme and focused on the chemical reaction mechanisms for melatonin activation, leading to competitive O-demethylation and C⁶-aromatic hydroxylation pathways. The calculations show that active site positioning determines the product distributions, but the bond that is activated is not necessarily closest to the heme in the enzyme-substrate complex. As such, the docking and molecular dynamics positioning of the substrate versus oxidant can give misleading predictions on product distributions. In particular, in quantum mechanics cluster model I, we observe that through a tight hydrogen bonding network, a preferential 6-hydroxylation of melatonin is obtained. However, O-demethylation becomes possible in alternative substrate-binding orientations that have the C⁶-aromatic ring position shielded. Finally, we investigated enzymatic and non-enzymatic O-demethylation processes and show that the hydrogen bonding network in the substrate-binding pocket can assist and perform this step prior to product release from the enzyme.

Reliably assessing the electronic structure of cytochrome P450 on today's classical computers and tomorrow's quantum computers

Chemical simulation is one of the most promising applications for future quantum computers. It is thought that quantum computers may enable accurate simulation for complex molecules that are otherwise impossible to simulate classically; that is, it displays quantum advantage. To better understand quantum advantage in chemical simulation, we explore what quantum and classical resources are required to simulate a series of pharmaceutically relevant molecules. Using classical methods, we show that reliable classical simulation of these molecules requires significant resources and therefore is a promising candidate for quantum simulation. We estimate the quantum resources, both in overall simulation time and the size. The insights from this study pave the way for future quantum simulation of complex molecules.

An accurate assessment of how quantum computers can be used for chemical simulation, especially their potential computational advantages, provides important context on how to deploy these future devices. To perform this assessment reliably, quantum resource estimates must be coupled with classical computations attempting to answer relevant chemical questions and to define the classical algorithms simulation frontier. Herein, we explore the quantum computation and classical computation resources required to assess the electronic structure of cytochrome P450 enzymes (CYPs) and thus define a classical–quantum advantage boundary. This is accomplished by analyzing the convergence of density matrix renormalization group plus n-electron valence state perturbation theory (DMRG+NEVPT2) and coupled-cluster singles doubles with noniterative triples [CCSD(T)] calculations for spin gaps in models of the CYP catalytic cycle that indicate multireference character. The quantum resources required to perform phase estimation using qubitized quantum walks are calculated for the same systems. Compilation into the surface code provides runtime estimates to compare directly to DMRG runtimes and to evaluate potential quantum advantage. Both classical and quantum resource estimates suggest that simulation of CYP models at scales large enough to balance dynamic and multiconfigurational electron correlation has the potential to be a quantum advantage problem and emphasizes the important interplay between classical computations and quantum algorithms development for chemical simulation.

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

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

Tris(Butadiene) Compounds versus Butadiene Oligomerization in Second-Row Transition Metal Chemistry: Effects of Increased Ligand Fields

The geometries, energetics, and preferred spin states of the second-row transition metal tris(butadiene) complexes (C₄H₆)₃M (M = Zr–Pd) and their isomers, including the experimentally known very stable molybdenum derivative (C₄H₆)₃Mo, have been examined by density functional theory. Such low-energy structures are found to have low-spin singlet and doublet spin states in contrast to the corresponding derivatives of the first-row transition metals. The three butadiene ligands in the lowest-energy (C₄H₆)₃M structures of the late second-row transition metals couple to form a C₁₂H₁₈ ligand that binds to the central metal atom as a hexahapto ligand for M = Pd but as an octahapto ligand for M = Rh and Ru. However, the lowest-energy (C₄H₆)₃M structures of the early transition metals have three separate tetrahapto butadiene ligands for M = Zr, Nb, and Mo or two tetrahapto butadiene ligands and one dihapto butadiene ligand for M = Tc. The low energy of the experimentally known singlet (C₄H₆)₃Mo structure contrasts with the very high energy of its experimentally unknown singlet chromium (C₄H₆)₃Cr analog relative to quintet (C₁₂H₁₈)Cr isomers with an open-chain C₁₂H₁₈ ligand.

Deciphering Reaction Determinants of Altered-Activity CYP2D6 Variants by Well-Tempered Metadynamics Simulation and QM/MM Calculations

The xenobiotic metabolizing enzyme CYP2D6 is the P450 cytochrome family member with the highest rate of polymorphism. This causes changes in the enzyme activity and specificity, which can ultimately lead to adverse reactions during drug treatment. To avoid or lower CYP-related toxicity risks, prediction of the most likely positions within a molecule where a metabolic reaction might occur is paramount. In order to obtain accurate predictions, it is crucial to understand all phenomena within the active site of the enzyme that contribute to an efficient substrate recognition and the subsequent catalytic reaction together with their relative weight within the overall thermodynamic context. This study aims to define the weight of the driving forces upon the C-H bond activation within CYP2D6 wild-type and a clinically relevant allelic variant with increased activity (CYP2D6*53) featuring two amino acid mutations in close vicinity of the heme. First, we investigated the steric and electrostatic complementarity of the substrate bufuralol using well-tempered metadynamics simulations with the aim to obtain the free energy profiles for each site of metabolism (SoM) within the different active sites. Second, the stereoelectronic complementarity was determined for each SoM within the two different active-site environments. Relying on the well-tempered metadynamics simulation energy profiles of each SoM, we identified the binding mode that was closest to the preferred transition-state geometry for efficient C-H bond activation. The binding modes were then used as starting structures for the quantum mechanics/molecular mechanics calculations performed to quantify the corresponding activation barriers. Our results show the relevance of the steric component in orienting the SoM in an energetically accessible position toward the heme. However, the corresponding intrinsic reactivity and electronic complementarity within the active site must be accurately evaluated in order to obtain a meaningful reaction prediction, from which the predominant SoM can be determined. The F120I mutation lowered the activation barrier for the major site and one of the minor SoMs. However, it had an impact neither on the CYP2D6 enantioselectivity preference of the oxidation reaction nor on the stereoselectivity from the substrate point of view.

CYP154C5 Regioselectivity in Steroid Hydroxylation Explored by Substrate Modifications and Protein Engineering*

CYP154C5 from Nocardia farcinica is a P450 monooxygenase able to hydroxylate a range of steroids with high regio- and stereoselectivity at the 16α-position. Using protein engineering and substrate modifications based on the crystal structure of CYP154C5, an altered regioselectivity of the enzyme in steroid hydroxylation had been achieved. Thus, conversion of progesterone by mutant CYP154C5 F92A resulted in formation of the corresponding 21-hydroxylated product 11-deoxycorticosterone in addition to 16α-hydroxylation. Using MD simulation, this altered regioselectivity appeared to result from an alternative binding mode of the steroid in the active site of mutant F92A. MD simulation further suggested that the entrance of water to the active site caused higher uncoupling in this mutant. Moreover, exclusive 15α-hydroxylation was observed for wild-type CYP154C5 in the conversion of 5α-androstan-3-one, lacking an oxy-functional group at C17. Overall, our data give valuable insight into the structure-function relationship of this cytochrome P450 monooxygenase for steroid hydroxylation.

Lignin Biodegradation by a Cytochrome P450 Enzyme: A Computational Study into Syringol Activation by GcoA

A recently characterized cytochrome P450 isozyme GcoA activates lignin components through a selective O-demethylation or alternatively an acetal formation reaction. These are important reactions in biotechnology and, because lignin is readily available; it being the main component in plant cell walls. In this work we present a density functional theory study on a large active site model of GcoA to investigate syringol activation by an iron(IV)-oxo heme cation radical oxidant (Compound I) leading to hemiacetal and acetal products. Several substrate-binding positions were tested and full energy landscapes calculated. The study shows that substrate positioning determines the product distributions. Thus, with the phenol group pointing away from the heme, an O-demethylation is predicted, whereas an initial hydrogen-atom abstraction of the weak phenolic O-H group would trigger a pathway leading to ring-closure to form acetal products. Predictions on how to engineer P450 GcoA to get more selective product distributions are given.

Perfluoroolefin complexes versus perfluorometallacycles and perfluorocarbene complexes in cyclopentadienylcobalt chemistry

Fluorocarbons have been shown experimentally by Baker and coworkers to combine with the cyclopentadienylcobalt (CpCo) moiety to form fluoroolefin and fluorocarbene complexes as well as fluorinated cobaltacyclic rings. In this connection density functional theory (DFT) studies on the cyclopentadienylcobalt fluorocarbon complexes CpCo(L)(C_nF_2n) (L = CO, PMe₃; n = 3 and 4) indicate structures with perfluoroolefin ligands to be the lowest energy structures followed by perfluorometallacycle structures and finally by structures with perfluorocarbene ligands. Thus, for the CpCo(L)(C₃F₆) (L = CO, PMe₃) complexes, the perfluoropropene structure has the lowest energy, followed by the perfluorocobaltacyclobutane structure and the perfluoroisopropylidene structure less stable by 8 to 11 kcal mol^-1, and the highest energy perfluoropropylidene structure less stable by more than 12 kcal mol^-1. For the two metal carbene structures Cp(L)Co[double bond, length as m-dash]C(CF₃)₂ and Cp(L)Co[double bond, length as m-dash]CF(C₂F₅), the former is more stable than the latter, even though the latter has Fischer carbene character. For the CpCo(L)(C₄F₈) (L = CO, PMe₃) complexes, the perfluoroolefin complex structures have the lowest energies, followed by the perfluorometallacycle structures at 10 to 20 kcal mol^-1, and the structures with perfluorocarbene ligands at yet higher energies more than 20 kcal mol^-1 above the lowest energy structure. This is consistent with the experimentally observed isomerization of the perfluorinated cobaltacyclobutane complexes CpCo(PPh₂Me)(-CFR-CF₂-CF₂-) (R = F, CF₃) to the perfluoroolefin complexes CpCo(PPh₂Me)(RCF[double bond, length as m-dash]CF₂) in the presence of catalytic quantities of HN(SO₂CF₃)₂. Further refinement of the relative energies by the state-of-the-art DLPNO-CCSD(T) method gives results essentially consistent with the DFT results summarized above.

Preclinical Pharmacokinetics of Lamivudine and Its Interaction with Schisandra chinensis Extract in Rats

Schisandra chinensis (Turcz.) Baill. (S. chinensis) extract and its active ingredient, schizandrin, have been used as a botanical medicine and dietary supplement for the treatment of hepatitis. Lamivudine is an antiretroviral drug and is used to treat hepatitis B viral infection. The aim of this study was to develop an ultrahigh-performance liquid chromatography-tandem mass spectrometry (UHPLC-MS/MS) method for the measurement of lamivudine and to determine the pharmacokinetic behaviors of an aqueous-ethanol extract of S. chinensis in rats. The separation was performed on a phenyl column maintained at 40 °C. The experimental animals were distributed into three groups: (1) lamivudine alone (10 mg/kg, i.v.); (2) lamivudine (10 mg/kg, i.v.) + pretreatment with S. chinensis (3 g/kg, p.o.); and (3) lamivudine (10 mg/kg, i.v.) + pretreatment with S. chinensis (10 g/kg, p.o.). The experimental results indicated that neither treatment with lamivudine alone nor pretreatment with S. chinensis (3 or 10 g/kg) significantly changed the pharmacokinetic parameters. In conclusion, based on the above preclinical experimental model, the combination of lamivudine with the herbal extract of S. chinensis did not exhibit significant pharmacokinetic interactions. These data offer useful information for assessing the preclinical safety of nutritional supplementation with lamivudine.

Investigation of the requirements for efficient and selective cytochrome P450 monooxygenase catalysis across different reactions

The cytochrome P450 metalloenzyme (CYP) CYP199A4 from Rhodopseudomonas palustris HaA2 catalyzes the highly efficient oxidation of para-substituted benzoic acids. Here we determined crystal structures of CYP199A4, and the binding and turnover parameters, with different meta-substituted benzoic acids in order to establish which criteria are important for efficient catalysis. When compared to the para isomers, the meta-substituted benzoic acids were less efficiently oxidized. For example, 3-formylbenzoic acid was oxidized with lower activity than the equivalent para isomer and 3-methoxybenzoic acid did not undergo O-demethylation by CYP199A4. The structural data highlighted that the meta-substituted benzoic acids bound in the enzyme active site in a modified position with incomplete loss of the distal water ligand of the heme moiety. However, for both sets of isomers the meta- or para-substituent pointed towards, and was in close proximity, to the heme iron. The absence of oxidation activity with 3-methoxybenzoic acid was assigned to the observation that the CH bonds of this molecule point away from the heme iron. In contrast, in the para isomer they are in an ideal location for abstraction. These findings were confirmed by using the bulkier 3-ethoxybenzoic acid as a substrate which removed the water ligand and reoriented the meta-substituent so that the methylene hydrogens pointed towards the heme, enabling more efficient oxidation. Overall we show relatively small changes in substrate structure and position in the active site can have a dramatic effect on the activity.

In Silico Predictions of Endocrine Disruptors Properties

Endocrine-disrupting chemicals (EDCs) are a broad class of molecules present in our environment that are suspected to cause adverse effects in the endocrine system by interfering with the synthesis, transport, degradation, or action of endogenous ligands. The characterization of the harmful interaction between environmental compounds and their potential cellular targets and the development of robust in vivo, in vitro, and in silico screening methods are important for assessment of the toxic potential of large numbers of chemicals. In this context, computer-aided technologies that will allow for activity prediction of endocrine disruptors and environmental risk assessments are being developed. These technologies must be able to cope with diverse data and connect chemistry at the atomic level with the biological activity at the cellular, organ, and organism levels. Quantitative structure-activity relationship methods became popular for toxicity issues. They correlate the chemical structure of compounds with biological activity through a number of molecular descriptors (e.g., molecular weight and parameters to account for hydrophobicity, topology, or electronic properties). Chemical structure analysis is a first step; however, modeling intermolecular interactions and cellular behavior will also be essential. The increasing number of three-dimensional crystal structures of EDCs' targets has provided a wealth of structural information that can be used to predict their interactions with EDCs using docking and scoring procedures. In the present review, we have described the various computer-assisted approaches that use ligands and targets properties to predict endocrine disruptor activities.

Proton relay network in P450cam formed upon docking of putidaredoxin

Cytochromes P450 are versatile heme-based enzymes responsible for vital life processes. Of these, P450cam (substrate camphor) has been most studied. Despite this, precise mechanisms of the key O─O cleavage step remain partly elusive to date; effects observed in various enzyme mutants remain partly unexplained. We have carried out extended (to 1000 ns) MM-MD and follow-on quantum mechanics/molecular mechanics computations, both on the well-studied FeOO state and on Cpd(0) (compound 0). Our simulations include (all camphor-bound): (a) WT (wild type), FeOO state. (b) WT, Cpd(0). (c) Pdx (Putidaredoxin, redox partner of P450)-docked-WT, FeOO state. (d) Pdx-docked WT, Cpd(0). (e) Pdx-docked T252A mutant, Cpd(0). Among our key findings: (a) Effect of Pdx docking appears to go far beyond that indicated in prior studies: it leads to specific alterations in secondary structure that create the crucial proton relay network. (b) Specific proton relay networks we identify are: FeOO(H)⋯T252⋯nH ₂ O⋯D251 in WT; FeOO(H)⋯nH ₂ O⋯D251 in T252A mutant; both occur with Pdx docking. (c) Direct interaction of D251 with -FeOOH is, respectively, rare/frequent in WT/T252A mutant. (d) In WT, T252 is in the proton relay network. (e) Positioning of camphor appears significant: when camphor is part of H-bonding network, second protonation appears to be facilitated.

Tris(butadiene) Metal Complexes of the First-Row Transition Metals versus Coupling of Butadiene to Eight- and Twelve-Carbon Hydrocarbon Chains

The role that zerovalent nickel plays in catalyzing the trimerization of butadiene to 1,5,9-cyclododecatriene conveys interest in the properties of the tris(butadiene)metal complexes (C₄H₆)₃M. In this connection the complexes (C₄H₆)₃M (M = Ti-Ni) of the first-row transition metals have been investigated by density functional theory. The intermediate C₁₂H₁₈Ni which has been isolated in the nickel-catalyzed trimerization of butadiene but is too unstable for X-ray crystallography is suggested here to have an open-chain hexahapto η^3,3-C₁₂H₁₈ ligand rather than the octahapto such ligand suggested by some investigators. The lowest energy (C₄H₆)₃M structures of the other first-row transition metals from vanadium to cobalt are found to have related structures with open-chain C₁₂H₁₈ ligands having hapticities ranging from four to eight with hexahapto structures being most common. The nickel and cobalt (C₁₂H₁₈)M derivatives favor low-spin singlet and doublet spin states, respectively, whereas the manganese derivative (C₁₂H₁₈)Mn favors the high-spin sextet state corresponding to the half-filled d⁵ shell of Mn(II). A (C₄H₆)₃Cr structure with three separate tetrahapto butadiene ligands analogous to the very stable (η⁴-C₄H₆)₃M (M = Mo, W) with the favored 18-electron metal configuration is found to be a very high energy structure relative to isomers containing an open-chain C₁₂H₁₈ ligand.

Coupling of fluoroborylene ligands in manganese carbonyl chemistry to give a difluorodiborene ligand

A bridging difluorodiborene μ-B₂F₂ ligand has been observed in two of the three lowest energy Mn₂(BF)₂(CO)₇ structures as well as in the lowest energy Mn₂(BF)₂(CO)₆ structure.

Alternative modes of bonding of C4F8 units in mononuclear and binuclear iron carbonyl complexes

The lowest energy C₄F₈Fe(CO)₄ structure is not the experimentally known ferracyclopentane complex but instead isomeric (perfluorobutene)iron tetracarbonyls. However, activation energies for the fluorine shifts required to form the latter isomers are very high. The lowest energy (C₄F₈)₂Fe₂(CO)_n (n = 7, 6) structures have bridging perfluorocarbene and terminal perfluoroolefin ligands.

Higher spin states in some low-energy bis(tetramethyl-1,2-diaza-3,5-diborolyl) sandwich compounds of the first row transition metals: boraza analogues of the metallocenes

Density functional studies on (Me₄B₂N₂CH)₂M (M = Ti, V, Cr, Mn, Fe, Co, Ni) show low-energy sandwich structures for all seven metals. The lowest-energy such Cr and Mn derivatives have higher spin states than the corresponding metallocenes.

Quantum chemical simulations revealed the toxicokinetic mechanisms of organic phosphorus flame retardants catalyzed by P450 enzymes

The metabolic fate and toxicokinetics of organic phosphorus flame retardants catalyzed by cytochrome P450 enzymes (CYPs) are here investigated by in silico simulations, leveraging an active center model to mimic the CYPs, triphenyl phosphate (TPHP), tris(2-butoxyethyl) phosphate and tris(1,3-dichloro-2-propyl) phosphate as substrates. Our calculations elucidated key main pathways and predicted products, which were corroborated by current in vitro data. Results showed that alkyl OPFRs are eliminated faster than aryl and halogenated alkyl-substituted OPFRs. In addition, we discovered a proton shuttle pathway for aryl hydroxylation of TPHP and P = O bond-assisted H-transfer mechanisms (rather than nonenzymatic hydrolysis) that lead to O-dealkylation/dearylation of phosphotriesters.

Butadiene as a ligand in open sandwich compounds

Density functional theory shows the lowest energy bis(butadiene)metal structures (C₄H₆)₂M (M = Ti to Ni) to have a staggered orientation of the two butadiene ligands corresponding to a tetrahedral coordination of the central metal atom.

Binuclear Cyclopentadienylmetal Methylene Sulfur Dioxide Complexes of Rhodium and Iridium Related to a Photochromic Metal Dithionite Complex

The photochromic dithionite complex Cp*₂Rh₂(μ-CH₂)₂(μ-O₂SSO₂) (Cp* = η⁵-Me₅C₅) is of interest because it undergoes an unusual fully reversible unimolecular photochemical rearrangement to the isodithionite complex Cp*₂Rh₂(μ-CH₂)₂(μ-O₂SOSO). In order to obtain more insight into these systems, a comprehensive density functional theory study has been carried out on isomeric Cp₂M₂(CH₂)₂(SO₂)₂ (M = Rh, Ir) derivatives. The experimentally observed rhodium complexes with coupled sulfur dioxide (SO₂) units to give dithionite or isodithionite ligands are surprisingly high-energy kinetic isomers in our analysis, reflecting the need for dithionite rather than SO₂ for their synthesis. Many isomeric structures containing two separate SO₂ ligands are found to lie at lower energies than these dithionite and isodithionite complexes. In the lowest-energy Cp₂M₂(CH₂)₂(SO₂)₂ isomers, the two methylene groups couple to form an ethylene ligand that can be either terminal or bridging. In slightly higher energy structures, a formal hydrogen shift is predicted to occur within the ethylene ligand to give a methylcarbene CH₃CH ligand. Isomers with a bridging methylcarbene ligand are energetically preferred over isomers with a terminal methylcarbene ligand. Generation of the lower-energy Cp₂Rh₂(CH₂)₂(SO₂)₂ isomers containing separate SO₂ ligands should be achievable through reactions of SO₂ with more highly reduced cyclopentadienylrhodium methylene complexes such as Cp*₂Rh₂(μ-CH₂)₂.