Challenging Density Functional Theory Calculations with Hemes and Porphyrins

Mechanistic Insights into Photochemical CO2 Reduction to CH4 by a Molecular Iron-Porphyrin Catalyst

Iron tetraphenylporphyrin complex modified with four trimethylammonium groups (Fe-p-TMA) is found to be capable of catalyzing the eight-electron eight-proton reduction of CO₂ to CH₄ photochemically in acetonitrile. In the present work, density functional theory (DFT) calculations have been performed to investigate the reaction mechanism and to rationalize the product selectivity. Our results revealed that the initial catalyst Fe-p-TMA ([Cl-Fe(III)-LR₄]⁴⁺, where L = tetraphenylporphyrin ligand with a total charge of -2, and R₄ = four trimethylammonium groups with a total charge of +4) undergoes three reduction steps, accompanied by the dissociation of the chloride ion to form [Fe(II)-L^••2-R₄]²⁺. [Fe(II)-L^••2-R₄]²⁺, bearing a Fe(II) center ferromagnetically coupled with a tetraphenylporphyrin diradical, performs a nucleophilic attack on CO₂ to produce the ¹η-CO₂ adduct [CO₂^•--Fe(II)-L^•-R₄]²⁺. Two intermolecular proton transfer steps then take place at the CO₂ moiety of [CO₂^•--Fe(II)-L^•-R₄]²⁺, resulting in the cleavage of the C-O bond and the formation of the critical intermediate [Fe(II)-CO]⁴⁺ after releasing a water molecule. Subsequently, [Fe(II)-CO]⁴⁺ accepts three electrons and one proton to generate [CHO-Fe(II)-L^•-R₄]²⁺, which finally undergoes a successive four-electron-five-proton reduction to produce methane without forming formaldehyde, methanol, or formate. Notably, the redox non-innocent tetraphenylporphyrin ligand was found to play an important role in CO₂ reduction since it could accept and transfer electron(s) during catalysis, thus keeping the ferrous ion at a relatively high oxidation state. Hydrogen evolution reaction via the formation of Fe-hydride ([Fe(II)-H]³⁺) turns out to endure a higher total barrier than the CO₂ reduction reaction, therefore providing a reasonable explanation for the origin of the product selectivity.

Amyloid fibril reduction through covalently modified lysine in HEWL and insulin

Proteins possess a variety of nucleophiles, which can carry out different reactions in the functioning cells. Proteins endogenously and synthetically can be modified through their nucleophilic sites. The roles of these chemical modifications have not been completely revealed. These modifications can alter the protein folding process. Protein folding directly affects the function of proteins. If an error in protein folding occurs, it may cause protein malfunction leading to several neurodegenerative disorders such as Alzheimer's and Parkinson's. In this study, Hen Egg White Lysozyme (HEWL) and bovine insulin, as model proteins for studying the amyloid formation, were covalently attached with 5(6)-thiophenolfluorescein. The amyloid formation of the covalently labeled lysozyme and insulin were compared with the native proteins. Interestingly, the results indicated that the covalent attachment of fluorescein slowed down the amyloid formation of HEWL and insulin significantly. The amyloid formation was examined using Thioflavin T (ThT) fluorescence assay, circular dichroism, FTIR, and gel electrophoresis. Tandem mass spectrometry was employed to identify the sites of covalent modifications in HEWL. It turned out that two surface lysine residues (K97 and K 116) in HEWL were modified. Computational studies, including docking and molecular simulations, revealed that 5(6)-thiophenolfluorescein makes several non-covalent interactions with HEWL residues, including Lys 97, leading to the reduction of the β-sheet in the protein. Additionally, AFM analysis confirmed the amyloid fibril reduction of lysine-modified bovine insulin and HEWL. Altogether, our results expand mechanistic insights into preventing amyloid formation by providing an approach for reducing amyloid formation by modifying specific lysine residues in the proteins.

Metal-organic frameworks as O2-selective adsorbents for air separations

Oxygen is a critical gas in numerous industries and is produced globally on a gigatonne scale, primarily through energy-intensive cryogenic distillation of air. The realization of large-scale adsorption-based air separations could enable a significant reduction in associated worldwide energy consumption and would constitute an important component of broader efforts to combat climate change. Certain small-scale air separations are carried out using N2-selective adsorbents, although the low capacities, poor selectivities, and high regeneration energies associated with these materials limit the extent of their usage. In contrast, the realization of O2-selective adsorbents may facilitate more widespread adoption of adsorptive air separations, which could enable the decentralization of O2 production and utilization and advance new uses for O2. Here, we present a detailed evaluation of the potential of metal-organic frameworks (MOFs) to serve as O2-selective adsorbents for air separations. Drawing insights from biological and molecular systems that selectively bind O2, we survey the field of O2-selective MOFs, highlighting progress and identifying promising areas for future exploration. As a guide for further research, the importance of moving beyond the traditional evaluation of O2 adsorption enthalpy, ΔH, is emphasized, and the free energy of O2 adsorption, ΔG, is discussed as the key metric for understanding and predicting MOF performance under practical conditions. Based on a proof-of-concept assessment of O2 binding carried out for eight different MOFs using experimentally derived capacities and thermodynamic parameters, we identify two existing materials and one proposed framework with nearly optimal ΔG values for operation under user-defined conditions. While enhancements are still needed in other material properties, the insights from the assessments herein serve as a guide for future materials design and evaluation. Computational approaches based on density functional theory with periodic boundary conditions are also discussed as complementary to experimental efforts, and new predictions enable identification of additional promising MOF systems for investigation.

Can one use the electronic absorption spectra of metalloporphyrins to benchmark electronic structure methods? A case study on the cobalt porphyrin

In this article, we describe calculations on the absorption spectrum of cobalt(ii) porphyrin, using density functional (DFT) and multireference n-electron valence perturbation (NEVPT) theories. With these calculations, we describe the lowest-energy states of doublet and quartet spin multiplicities, the excited states that originate the Q and B bands of porphyrins, some higher-energy π-π* excitations and charge-transfer states, HOMO-LUMO gaps, and ionisation potentials. Results undoubtedly show that the position of B band is essentially independent on the DFT functional, while the Q band is better described by pure functionals, and these bands do not depend on the initial state of the transition (whether doublet or quartet) as well. However, other excitation energies, orbital energies, and ionisation potentials strongly depend on the functional, in some cases varying more than 2 eV. Based on these results we conclude that one should not use the UV-Vis spectrum of metalloporphyrins to benchmark density functionals, mainly those properties related to coordination with the metallic ion. Furthermore, the results show that functionals that yield correct spectra may be based on an incorrect ground state description. Moreover, we reinforce that one must be skeptical about the reference chosen to benchmark electronic structure calculations, such as DFT functionals and active spaces for multireference calculations.

Bioengineering of Cytochrome P450 OleTJE: How Does Substrate Positioning Affect the Product Distributions?

The cytochromes P450 are versatile enzymes found in all forms of life. Most P450s use dioxygen on a heme center to activate substrates, but one class of P450s utilizes hydrogen peroxide instead. Within the class of P450 peroxygenases, the P450 OleT_JE isozyme binds fatty acid substrates and converts them into a range of products through the α-hydroxylation, β-hydroxylation and decarboxylation of the substrate. The latter produces hydrocarbon products and hence can be used as biofuels. The origin of these product distributions is unclear, and, as such, we decided to investigate substrate positioning in the active site and find out what the effect is on the chemoselectivity of the reaction. In this work we present a detailed computational study on the wild-type and engineered structures of P450 OleT_JE using a combination of density functional theory and quantum mechanics/molecular mechanics methods. We initially explore the wild-type structure with a variety of methods and models and show that various substrate activation transition states are close in energy and hence small perturbations as through the protein may affect product distributions. We then engineered the protein by generating an in silico model of the double mutant Asn242Arg/Arg245Asn that moves the position of an active site Arg residue in the substrate-binding pocket that is known to form a salt-bridge with the substrate. The substrate activation by the iron(IV)-oxo heme cation radical species (Compound I) was again studied using quantum mechanics/molecular mechanics (QM/MM) methods. Dramatic differences in reactivity patterns, barrier heights and structure are seen, which shows the importance of correct substrate positioning in the protein and the effect of the second-coordination sphere on the selectivity and activity of enzymes.

Structure, Properties, and Reactivity of Porphyrins on Surfaces and Nanostructures with Periodic DFT Calculations

Porphyrins are fascinating molecules with applications spanning various scientific fields. In this review we present the use of periodic density functional theory (PDFT) calculations to study the structure, electronic properties, and reactivity of porphyrins on ordered two dimensional surfaces and in the formation of nanostructures. The focus of the review is to describe the application of PDFT calculations for bridging the gaps in experimental studies on porphyrin nanostructures and self-assembly on 2D surfaces. A survey of different DFT functionals used to study the porphyrin-based system as well as their advantages and disadvantages in studying these systems is presented.

Unsymmetrical nonplanar 'push-pull' β-octasubstituted porphyrins: facile synthesis, structural, photophysical and electrochemical redox properties

Mixed substitution at the β-position of porphyrins influences their photophysical and electrochemical redox properties. Two new series of asymmetrically mixed β-octasubstituted porphyrins viz. MTPP(Ph)2Br5X (X = NO2 or Br and M = 2H, Co(ii), Ni(ii), Cu(ii), and Zn(ii)) have been synthesized and characterized by various spectroscopic techniques. The single crystal X-ray structure of H2TPP(NO2)(Ph)2Br5 showed a nonplanar saddle shape conformation of the macrocyclic core. Furthermore, the fully optimized geometries confirmed the saddle shape conformation of H2TPP(Ph)2Br5X (X = NO2 or Br). Electronic spectra revealed a significant bathochromic shift by appending both electron donor and acceptor substituents at the β-position of the meso-tetraphenylporphyrin skeleton, which reflects the following order H2TPP < H2TPP(NO2) < H2TPP(NO2)(Ph)2 < H2TPP(Ph)2Br6 < H2TPP(NO2)(Ph)2Br5. H2TPP(Ph)2Br5X (X = NO2 or Br) exhibited a significant bathochromic shift (Δλmax = 53-61 nm) in the Soret band and (Δλmax = 90-95 nm) in the longest wavelength Qx(0,0) band as compared to H2TPP. Nonplanar conformations and electron withdrawing β-substituents induce higher protonation and deprotonation constants for H2TPP(NO2)(Ph)2Br5 and H2TPP(Ph)2Br6 as compared to precursor porphyrins viz. H2TPP, H2TPP(NO2) and H2TPP(NO2)(Ph)2. The electronic spectral properties and redox potentials of MTPP(Ph)2Br5X (X = NO2 or Br and M = 2H, Co, Ni, Cu and Zn) are affected by β-substituents at the periphery of the porphyrin core. Redox tunability was achieved by appending push-pull substituents at the β-position of the MTPP (M = 2H, CoII, NiII, CuII, and ZnII) skeleton of the macrocycle. CuTPP(Ph)2Br6 and CuTPP(NO2)(Ph)2Br5 exhibited a dramatically reduced HOMO-LUMO gap with a difference of 0.55 V and 0.62 V, respectively as compared to CuTPP due to the push-pull effect of β-substituents and nonplanarity of the porphyrin core.

DFT study on the effect of proximal residues on the Mycobacterium tuberculosis catalase-peroxidase (katG) heme compound I intermediate and its bonding interaction with isoniazid

Isoniazid (INH) is converted into isonicotinyl radical through its interaction with the catalase-peroxidase (katG) enzyme present in the cells of Mycobacterium tuberculosis (M. tb.), the bacteria that causes the tuberculosis disease. This process is important because resistance of M. tb. cells to INH treatment has been associated with the failure of completion of this process. However, this process is poorly understood and there are a variety of conflicting theories about the details of the mechanism of this process. One theory suggests that INH binds to katG and transfers a single electron to the heme while the heme is in its two electron oxidized state, compound I [Fe(iv)Por˙+] (CpdI). In this study, DFT calculations at the UB3LYP/6-31g(d)-LANL2DZ level of theory are used to study the M. tb. katG CpdI molecule. Different models of the M. tb. CpdI molecule were prepared and the calculations revealed the impact of Trp321 on the electronic properties of the heme. Without Trp321 the heme assumed a triradical state with single electrons on the πxy and πyz orbitals of Fe and another on the a2u orbital of the porphyrin ring that can either be coupled with the first two, to assume a quartet state, or decoupled to form a doublet state. With Trp321, however, a transfer of an electron from the πTrp orbital to a2u porphyrin orbital leads to loss of radical character of the porphyrin and leaves the Trp321 group with a radical character. INH was observed to have strong interaction with CpdI and the absence of Trp321 significantly decreased the binding energy by 2 kcal mol-1 explaining the importance of Trp321 in the binding of INH. The results in this study show the importance of Trp321 in the binding of INH and its effect on its electronic properties, which is consistent with previous experimental findings that mutation of Trp321 results in an increase in drug resistance.

Heuristics from Modeling of Spectral Overlap in Förster Resonance Energy Transfer (FRET)

Among the photophysical parameters that underpin Förster resonance energy transfer (FRET), perhaps the least explored is the spectral overlap term ( J). While by definition J increases linearly with acceptor molar absorption coefficient (ε_(A) in M^-1 cm^-1), is proportional to wavelength (λ⁴), and depends on the degree of overlap of the donor fluorescence and acceptor absorption spectra, the question arose as to the value of J for the case of perfect spectral overlap versus that for representative fluorophores with incomplete spectral overlap. Here, Gaussian distributions of absorption and fluorescent spectra have been modeled that encompass varying degrees of overlap, full-width-at-half-maximum (fwhm), and Stokes shift. For ε_(A) = 10⁵ M^-1 cm^-1 and perfect overlap, the J value (in M^-1 cm^-1 nm⁴) ranges from 1.15 × 10¹⁴ (200 nm) to 7.07 × 10¹⁶ (1000 nm), is almost linear with λ⁴ (average of λ_abs and λ_flu), and is nearly independent of fwhm. For visible-region fluorophores with perfectly overlapped Gaussian spectra, the resulting value of J ( J_G-0) is ∼0.71 ε_(A)λ⁴ (M^-1 cm^-1 nm⁴). The experimental J values for homotransfer, as occurs in light-harvesting antennas, were calculated with spectra from a static database of 60 representative compounds (12 groups, 5 compounds each) and found to range from 4.2 × 10¹⁰ ( o-xylene) to 5.3 × 10¹⁶ M^-1 cm^-1 nm⁴ (a naphthalocyanine). The degree of overlap, defined by the ratio of the experimental J to the model J_G-0 for perfectly overlapped spectra, ranges from ∼0.5% (coumarin 151) to 77% (bacteriochlorophyll a). The results provide insights into how a variety of factors affect the resulting J values. The high degree of spectral overlap for (bacterio)chlorophylls prompts brief conjecture concerning the relevance of energy transfer to the question "why chlorophyll".

Mechanistic Insight on the Activity and Substrate Selectivity of Nonheme Iron Dioxygenases

Nonheme iron dioxygenases catalyze vital reactions for human health particularly related to aging processes. They are involved in the biosynthesis of amino acids, but also the biodegradation of toxic compounds. Typically they react with their substrate(s) through oxygen atom transfer, although often with the assistance of a co-substrate like α-ketoglutarate that is converted to succinate and CO₂ . Many reaction processes catalyzed by the nonheme iron dioxygenases are stereoselective or regiospecific and hence understanding the mechanism and protein involvement in the selectivity is important for the design of biotechnological applications of these enzymes. To this end, I will review recent work of our group on nonheme iron dioxygenases and include background information on their general structure and catalytic cycle. Examples of stereoselective and regiospecific reaction mechanisms we elucidated are for the AlkB repair enzyme, prolyl-4-hydroxylase and the ergothioneine biosynthesis enzyme. Finally, I cover an example where we bioengineered S-p-hydroxymandelate synthase into the R-p-hydroxymandelate synthase.

Heme isomers substantially affect heme's electronic structure and function

Inspection of heme protein structures in the protein data bank reveals four isomers of heme characterized by different relative orientations of the vinyl side chains; remarkably, all these have been reported in multiple protein structures. Density functional theory computations explain this as due to similar energy of the isomers but with a sizable (25 kJ mol^-1) barrier to interconversion arising from restricted rotation around the conjugated bonds. The four isomers, EE, EZ, ZE, and ZZ, were then investigated as 4-coordinate hemes, as 5-coordinate deoxyhemes, in 6-coordinate O₂-adducts of globins and as compound I intermediates typical of heme peroxidases. Substantial differences were observed in electronic properties relevant to heme function: notably, the spin state energy gap of O₂-heme adducts, important for fast reversible binding of O₂, depends on the isomer state, and O₂-binding enthalpies change by up to 16 kJ mol^-1; redox potentials change by up to 0.2 V depending on the isomer, and the doublet-quartet energy splitting of compound I, central to "two-state" reactivity, is affected by up to ∼15 kJ mol^-1. These effects are consistently seen with three distinct density functionals, i.e. the effects are not method-dependent. Thus, the nature of the isomer state is an important but overlooked feature of heme chemistry and function, and previous and future studies of hemes may be reconsidered in this new context.

Ultrafast dynamics in co-sensitized photocatalysts under visible and NIR light irradiation

Activity modulation of co-sensitized light harvesting nanohybrids by tuning the ultrafast carrier dynamics under visible and NIR light irradiation.

Ethnopharmacological Approaches for Therapy of Jaundice: Part I

Jaundice is a very common symptom especially in the developing countries. It is associated with several hepatic diseases which are still major causes of death. There are many different approaches to jaundice treatment and the growing number of ethnomedicinal studies shows the plant pharmacology as very promising direction. Many medicinal plants are used for the treatment of jaundice, however a comprehensive review on this subject has not been published. The use of medicinal plants in drug discovery is highly emphasized (based on their traditional and safe uses in different folk medicine systems from ancient times). Many sophisticated analytical techniques are emerging in the pharmaceutical field to validate and discover new biologically active chemical entities derived from plants. Here, we aim to classify and categorize medicinal plants relevant for the treatment of jaundice according to their origin, geographical location, and usage. Our search included various databases like Pubmed, ScienceDirect, Google Scholar. Keywords and phrases used for these searches included: "jaundice," "hyperbilirubinemia," "serum glutamate," "bilirubin," "Ayurveda." The first part of the review focuses on the variety of medicinal plant used for the treatment of jaundice (a total of 207 medicinal plants). In the second part, possible mechanisms of action of biologically active secondary metabolites of plants from five families for jaundice treatment are discussed.

Ethnopharmacological Approaches for Therapy of Jaundice: Part I

Jaundice is a very common symptom especially in the developing countries. It is associated with several hepatic diseases which are still major causes of death. There are many different approaches to jaundice treatment and the growing number of ethnomedicinal studies shows the plant pharmacology as very promising direction. Many medicinal plants are used for the treatment of jaundice, however a comprehensive review on this subject has not been published. The use of medicinal plants in drug discovery is highly emphasized (based on their traditional and safe uses in different folk medicine systems from ancient times). Many sophisticated analytical techniques are emerging in the pharmaceutical field to validate and discover new biologically active chemical entities derived from plants. Here, we aim to classify and categorize medicinal plants relevant for the treatment of jaundice according to their origin, geographical location, and usage. Our search included various databases like Pubmed, ScienceDirect, Google Scholar. Keywords and phrases used for these searches included: “jaundice,” “hyperbilirubinemia,” “serum glutamate,” “bilirubin,” “Ayurveda.” The first part of the review focuses on the variety of medicinal plant used for the treatment of jaundice (a total of 207 medicinal plants). In the second part, possible mechanisms of action of biologically active secondary metabolites of plants from five families for jaundice treatment are discussed.

Singlet versus Triplet Reactivity in an Mn(V)-Oxo Species: Testing Theoretical Predictions Against Experimental Evidence

Discerning the factors that control the reactivity of high-valent metal-oxo species is critical to both an understanding of metalloenzyme reactivity and related transition metal catalysts. Computational studies have suggested that an excited higher spin state in a number of metal-oxo species can provide a lower energy barrier for oxidation reactions, leading to the conclusion that this unobserved higher spin state complex should be considered as the active oxidant. However, testing these computational predictions by experiment is difficult and has rarely been accomplished. Herein, we describe a detailed computational study on the role of spin state in the reactivity of a high-valent manganese(V)-oxo complex with para-Z-substituted thioanisoles and utilize experimental evidence to distinguish between the theoretical results. The calculations show an unusual change in mechanism occurs for the dominant singlet spin state that correlates with the electron-donating property of the para-Z substituent, while this change is not observed on the triplet spin state. Minimum energy crossing point calculations predict small spin-orbit coupling constants making the spin state change from low spin to high spin unlikely. The trends in reactivity for the para-Z-substituted thioanisole derivatives provide an experimental measure for the spin state reactivity in manganese-oxo corrolazine complexes. Hence, the calculations show that the V-shaped Hammett plot is reproduced by the singlet surface but not by the triplet state trend. The substituent effect is explained with valence bond models, which confirm a change from an electrophilic to a nucleophilic mechanism through a change of substituent.