Protein effects in non-heme iron enzyme catalysis: insights from multiscale models

Implications of non-native metal substitution in carbonic anhydrase - engineered enzymes and models

The enzyme carbonic anhydrase has been intensely studied over decades as a means to understand the role of zinc in hydrating CO2. The naturally occurring enzyme has also been immobilized on distinct heterogeneous platforms, which results in a different hybrid class of catalysts that are useful for the adsorption and hydration of CO2. However, the reusability and robustness of such natural and immobilized systems are substantially affected when tested under industrial conditions, such as high temperature and high flow rate. This led to the generation of model systems in the form of metal-coordination complexes, metal-organic frameworks, metallo-peptide self-assembled supramolecules and nanomaterials that mimic the primary, and, to some extent, secondary coordination sphere of the active site of the natural carbonic anhydrase enzymes. Furthermore, the effects of zinc-substitution by other relevant transition metals in both the naturally occurring enzymes and model systems has been reported. It has been observed that some other transition metal ions in the active site of carbonic anhydrase and its models can also accomplish similar activity, established by various reaction probes and ideas. Herein, we present a comprehensive highlight about substituting zinc in the active site of the modified enzymes and its biomimetic model systems with non-native metal ions and review how they affect the structural orientation and reactivity towards CO2 hydration. In addition, the utility of artificially engineered carbonic anhydrases towards a number of non-natural reactions is also discussed.

A Hemicyanine-Based Highly Sensitive and Selective Near-Infrared Fluorescent Probe for Fe3+ in Aqueous Media

Ferric ion is widely distributed in human cells and an important component of hemoglobin, which can promote the transportation of blood in the human body. Its ability to bind with oxygen is essential for participating in oxidation reactions and enzymatic reactions. Deficiency of iron (III) and excessive accumulation of iron in the blood can cause many health problems in the body, such as anemia, loss of appetite, fibrosis, reduced work routines, and decreased immunity. As a fluorescent dye, cyanine has the advantages of high fluorescence quantum yield, good optical and chemical stability, excitation and emission wavelengths in the near-infrared region, high molar extinction coefficient, and small influence on pH variation. It is a fluorescent dye that modern scholars choose more. A new near-IR hemicyanine derivative bearing 1,3-dithiane moiety was developed as an efficient fluorescence probe 1 for the detection of Fe3+ in aqueous medium (THF/H2O, 1:1, v:v). This near-IR fluorescence probe 1 exhibits high sensitivity and selectivity toward Fe3+ sensing with detection limit of 0.5 μM under lager scope of pH value (pH = 2-12).

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

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

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

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

Fabrication of P and N Co-Doped Carbon Dots for Fe3+ Detection in Serum and Lysosomal Tracking in Living Cells

Doping with heteroatoms allows the retention of the general characteristics of carbon dots while allowing their physicochemical and photochemical properties to be effectively modulated. In this work, we report the preparation of ultrastable P and N co-doped carbon dots (PNCDs) that can be used for the highly selective detection of Fe³⁺ and the tracking of lysosomes in living cells. Fluorescent PNCDs were facilely prepared via a hydrothermal treatment of ethylenediamine and phytic acid, and they exhibited a high quantum yield of 22.0%. The strong coordination interaction between the phosphorus groups of PNCDs and Fe³⁺ rendered them efficient probes for use in selective Fe³⁺ detection, with a detection limit of 0.39 μM, and we demonstrated their practicability by accurately detecting the Fe³⁺ contents in bio-samples. At the same time, PNCDs exhibited high lysosomal location specificity in different cell lines due to surface lipophilic amino groups, and real-time tracking of the lysosome morphology in HeLa cells was achieved. The present work suggests that the fabrication of heteroatom-doped CDs might be an effective strategy to provide promising tools for cytology, such as organelle tracking.

Second coordination sphere effects on the mechanistic pathways for dioxygen activation by a ferritin: involvement of a Tyr radical and the identification of a cation binding site

Ferritins are ubiquitous diiron enzymes involved in iron(II) detoxification and oxidative stress responses and can act as metabolic iron stores. The overall reaction mechanisms of ferritin enzymes are still unclear, particularly concerning the role of the conserved, near catalytic center Tyr residue. Thus, we carried out a computational study of a ferritin using a large cluster model of well over 300 atoms including its first- and second-coordination sphere. The calculations reveal important insight into the structure and reactivity of ferritins. Specifically, the active site Tyr residue delivers a proton and electron in the catalytic cycle prior to iron(II) oxidation. In addition, the calculations highlight a likely cation binding site at Asp65, which through long-range electrostatic interactions, influences the electronic configuration and charge distributions of the metal center. The results are consistent with experimental observations but reveal novel detail of early mechanistic steps that lead to an unusual mixed-valent iron(III)-iron(II) center.

Conformational Dynamics Underlies Different Functions of Human KDM7 Histone Demethylases

The human KDM7 subfamily histone H3 Nϵ-methyl lysine demethylases PHF8 (KDM7B) and KIAA1718 (KDM7A) have different substrate selectivities and are linked to genetic diseases and cancer. We describe experimentally based computational studies revealing that flexibility of the region linking the PHD finger and JmjC domains in PHF8 and KIAA1718 regulates interdomain interactions, the nature of correlated motions, and ultimately H3 binding and demethylation site selectivity. F279S an X-linked mental retardation mutation in PHF8 is involved in correlated motions with the iron ligands and second sphere residues. The calculations reveal key roles of a flexible protein environment in productive formation of enzyme-substrate complexes and suggest targeting the flexible KDM7 linker region is of interest from a medicinal chemistry perspective.

Revisiting the catalytic mechanism of Mo-Cu carbon monoxide dehydrogenase using QM/MM and DFT calculations

Previous density functional theory (DFT) studies have shown that the release of the produced carbon dioxide (CO2) from an active-site cluster is a thermodynamically or kinetically difficult step in the enzymatic carbon monoxide (CO) oxidation catalyzed by Mo-Cu carbon monoxide dehydrogenase (Mo-Cu CODH). To better understand the effect of the protein environment on this difficult CO2 release step as well as other reaction steps, we applied hybrid quantum mechanics and molecular mechanics (QM/MM) calculations to the Mo-Cu CODH enzyme. The results show that in the first step, the equatorial Mo[double bond, length as m-dash]O group in the active-site cluster attacks the nearby CO molecule bound to the Cu site. Afterward, a stable thiocarbonate intermediate is formed in which the CO2 molecule is embedded and the copper-S(μ-sulfido) bond is broken. A free CO2 molecule, i.e., the final product, is then released from the active-site cluster, not directly from the thiocarbonate intermediate but via a previously formed intermediate that also contains CO2 but retains the Cu-S(μ-sulfido) bond. In contrast to the previous DFT results, the calculated barrier for this process was low in our QM/MM calculations. An additional QM/MM analysis of the barrier height showed that the effect of the protein environment on this barrier lowering is not very large. We found that the reason for the low barrier obtained by QM/MM is that the barrier for CO2 release is already not high at the DFT level. These results allow us to conclude that the CO oxidation reaction passes through the formation of a thiocarbonate intermediate, and that the subsequent CO2 release is kinetically not difficult. Nevertheless, the protein environment has an important role to play in making the latter process thermodynamically favored. No low-barrier pathway for the product release could be obtained for the reaction of n-butylisocyanide, which is consistent with the experimental fact that n-butylisocyanide inhibits Mo-Cu CODH.

Computational approaches for deciphering the equilibrium and kinetic properties of iron transport proteins

With the advances in three-dimensional structure determination techniques, high quality structures of the iron transport proteins transferrin and the bacterial ferric binding protein (FbpA) have been deposited in the past decade. These are proteins of relatively large size, and developments in hardware and software have only recently made it possible to study their dynamics using standard computational resources. We review computational techniques towards understanding the equilibrium and kinetic properties of iron transport proteins under different environmental conditions. At the level of detail that requires quantum chemical treatments, the octahedral geometry around iron has been scrutinized and it has been established that the iron coordinating tyrosines are in an unusual deprotonated state. At the atomistic level, both the N-lobe and the full bilobal structure of transferrin have been studied under varying conditions of pH, ionic strength and binding of other metal ions by molecular dynamics (MD) simulations. These studies have allowed questions to be answered, among others, on the function of second shell residues in iron release, the role of synergistic anions in preparing the active site for iron binding, and the differences between the kinetics of the N- and the C-lobe. MD simulations on FbpA have led to the detailed observation of the binding kinetics of phosphate to the apo form, and to the conformational preferences of the holo form under conditions mimicking the environmental niches provided by the periplasmic space. To study the dynamics of these proteins with their receptors, one must resort to coarse-grained methodologies, since these systems are prohibitively large for atomistic simulations. A study of the complex of human transferrin (hTf) with its pathogenic receptor by such methods has revealed a potential mechanistic explanation for the defense mechanism that arises in evolutionary warfare. Meanwhile, the motions in the transferrin receptor bound hTf have been shown to disfavor apo hTf dissociation, explaining why the two proteins remain in complex during the recycling process from the endosome to the cell surface. Open problems and possible technological applications related to metal ion binding-release in iron transport proteins that may be handled by hybrid use of quantum mechanical, MD and coarse-grained approaches are discussed.