Density functional theory studies on the conversion of hydroxyheme to iron-verdoheme in the presence of dioxygen

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.

How to correct relative voxel scale factors for calculations of vector-difference Fourier maps in cryo-EM

The atomic coordinates derived from cryo-electron microscopy (cryo-EM) maps can be inaccurate when the voxel scaling factors are not properly calibrated. Here, we describe a method for correcting relative voxel scaling factors between pairs of cryo-EM maps for the same or similar structures that are expanded or contracted relative to each other. We find that the correction of scaling factors reduces the amplitude differences of Fourier-inverted structure factors from voxel-rescaled maps by up to 20-30%, as shown by two cryo-EM maps of the SARS-CoV-2 spike protein measured at pH 4.0 and pH 8.0. This allows for the calculation of the difference map after properly scaling, revealing differences between the two structures for individual amino acid residues. Unexpectedly, the analysis uncovers two previously overlooked differences of amino acid residues in structures and their local structural changes. Furthermore, we demonstrate the method as applied to two cryo-EM maps of monomeric apo-photosystem II from the cyanobacteria Synechocystis sp. PCC 6803 and Thermosynechococcus elongatus. The resulting difference maps reveal many changes in the peripheral transmembrane PsbX subunit between the two species.

Ionic bond in hydrogen transferring of the ferrous and/or ferric human/mouse verdoheme oxygenase

Formation of five coordinated ferric (ferrous) verdoheme oxygenase complexes have been investigated at ωB97X-D/6-31G(d) level of theory. This process was carried out by adsorption of imidazole and human/mouse verdoheme oxygenase (VO) compounds. Global reactivity indexes show electrophile and nucleophile roles of the VO complexes and Imidazole, respectively. This result confirms their interaction, molecular electrostatic potential (MEP) maps, and low HOMO_FRVMO-LUMO_Imidazole gap. These interactions can cause in adsorption and five coordinated of the VO complexes. More negative value (-64.3 kJ mol^-1) of adsorption energy (E_ads) in the FRVMO complex shows better adsorption strength and stable configuration. Significant point of this interaction is hydrogen transfer from imidazole to the nearest oxygen of the VO complexes; this issue is approved using quantum theory of atom in molecule (QTAIM) and natural bond orbital (NBO) analysis. QTAIM calculations confirm ionic bonding between the transferred hydrogen and the oxygen atom of the VO. The 312.2-kcal mol^-1 s order stabilization energies in this complex are confirmation for strong donation and better formation of five coordinated complex in electron view point.

Theoretical kinetics and thermodynamics study: Peripheral substituent effects on the hydrolysis of verdoheme

The heme oxygenase (HO) enzyme is a free heme protein that binds to heme in the body. Heme acts as both a cofactor and a substrate in this enzyme. The catabolism of heme into biliverdin, monoxide carbon, and free-iron, catalyzed by heme oxygenase via three consecutive oxygenation steps, in which the heme group functions as the prosthetic group as well as the substrate. Investigations of the reactions of the peripheral substituent on the heme ring with 5-oxaporphyrin iron complexes (verdohemes) have been assumed to provide models and largely unknown for the primary step in the hydrolysis of verdohemes. In this work, a theoretical kinetics and thermodynamics study of the degradation reactions of verdohemes was performed, and calculations show that the [Formula: see text] in the hydrolysis of verdohemes with non-peripheral substituents is more negative than hydrolysis of verdohemes with peripheral substituents. In other words, the hydrolysis of verdohemes with non-peripheral substituents is more energy-efficient than verdohemes with a peripheral substituents. Equilibrium constant calculations show that hydrolysis of verdohemes with non-peripheral substituents is much faster than that of verdohemes with peripheral substituents, which is due to a more convenient nucleophilic attack on the cationic ring than the anionic ring. To acquire a good molecular understanding, peripheral substituent effects on the hydrolysis of verdoheme’s inhibitory role was studied using the DFT method.

ONIOM investigations of the heme degradation mechanism by MhuD: the critical function of heme ruffling

A unique ruffling conformation of hydroxyheme in MhuD inhibits its “on-site” monooxygenation but induces “remote-site” dioxygenation.

On the Difference Between Additive and Subtractive QM/MM Calculations

The combined quantum mechanical (QM) and molecular mechanical (MM) approach (QM/MM) is a popular method to study reactions in biochemical macromolecules. Even if the general procedure of using QM for a small, but interesting part of the system and MM for the rest is common to all approaches, the details of the implementations vary extensively, especially the treatment of the interface between the two systems. For example, QM/MM can use either additive or subtractive schemes, of which the former is often said to be preferable, although the two schemes are often mixed up with mechanical and electrostatic embedding. In this article, we clarify the similarities and differences of the two approaches. We show that inherently, the two approaches should be identical and in practice require the same sets of parameters. However, the subtractive scheme provides an opportunity to correct errors introduced by the truncation of the QM system, i.e., the link atoms, but such corrections require additional MM parameters for the QM system. We describe and test three types of link-atom correction, viz. for van der Waals, electrostatic, and bonded interactions. The calculations show that electrostatic and bonded link-atom corrections often give rise to problems in the geometries and energies. The van der Waals link-atom corrections are quite small and give results similar to a pure additive QM/MM scheme. Therefore, both approaches can be recommended.

A novel mechanism of heme degradation to biliverdin studied by QM/MM and QM calculations

Heme degradation by heme oxygenase enzymes is important for maintaining iron homeostasis and prevention of oxidative stress.

QM/MM Study of the Conversion of Oxophlorin into Verdoheme by Heme Oxygenase

Heme oxygenase is an enzyme that degrades heme, thereby recycling iron in most organisms, including humans. Pervious density functional theory (DFT) calculations have suggested that iron(III) hydroxyheme, an intermediate generated in the first step of heme degradation by heme oxygenase, is converted to iron(III) superoxo oxophlorin in the presence of dioxygen. In this article, we have studied the detailed mechanism of conversion of iron(III) superoxo oxophlorin to verdoheme by using combined quantum mechanics and molecular mechanics (QM/MM) calculations. The calculations employed the B3LYP method and the def2-QZVP basis set, considering dispersion effects with the DFT-D3 approach, obtaining accurate energies with large QM regions of almost 1000 atoms. The reaction was found to be exothermic by -35 kcal/mol, with a rate-determining barrier of 19 kcal/mol in the doublet state. The protein environment and especially water in the enzyme pocket significantly affects the reaction by decreasing the reaction activation energies and changing the structures by providing strategic hydrogen bonds.