skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: CO Bonding in Hexa- and Pentacoordinate Carboxy-Neuroglobin - A QM/MM and Local Vibrational Mode Study
We present a comprehensive investigation on the different role of CO in carboxy- neuroglobin i) as ligand of the heme group in the active site forming a bond with the heme iron and ii) dissociated from the heme group but still trapped inside the active site, focusing on two specific orientations, one with CO perpendicular to the plane defined by the distal histidine of the enzyme (form A) and one with CO located parallel to that plane (form B). Our study includes wild type carboxy-neuroglobin and nine known protein mutations. Considering that the distal histidine interacting with the heme group can adapt two different tautomeric forms and the two possible orientations of the dissociated CO, a total of 36 protein systems were analyzed in this study. Fully optimized geometries and vibrational frequencies were calculated at the QM/MM level, followed by the local mode analysis, to decode CO bond properties. The intrinsic bond strengths derived from the local mode analysis, complemented with NBO and QTAIM data, reveal that the strength of the CO bond, in the hexacoordinate (where CO is a ligand of the heme group) and pentacoordinate (where CO is dissociated from the heme group) scenarios, is dominated by through bond and through space charge transfer between CO and Fe, fine-tuned by electrostatic and dispersion interactions with the side chain amino acids in the distal heme pocket. Suggestions are made as to advise on how protein modifications can influence the molecular properties of the coordinated or dissociated CO, which could serve the fine-tuning of existing and the design of new neuroglobin models with specific FeC and CO bond strengths.  more » « less
Award ID(s):
2102461
PAR ID:
10342535
Author(s) / Creator(s):
Date Published:
Journal Name:
Journal of computational chemistry
ISSN:
0192-8651
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract Metal–ligand bonding and noncovalent interactions (NCIs), such as hydrogen bonding orπ–πinteractions, play a crucial role in determining the structure, function, and selectivity of both biological and artificial metalloproteins. In this study, we employed a hybrid quantum mechanics/molecular mechanics (QM/MM) approach to investigate the ligation of water or cyanide in a mutated myoglobin system, in which the native heme scaffold was replaced with M-salophen or M-salen Schiff base complexes (M = Cr, Mn, Fe). Using our local vibrational mode analysis, particularly local vibrational mode force constants as intrinsic bond strength parameters, complemented with electron density and natural orbital analyses we explored the role of metal–ligand bonding and NCIs in different environments within the myoglobin pocket. Our analysis revealed that metal–ligand bonding, for both water and cyanide ligands, is strongest in the delta form of distal histidine and favors salophen prosthetic groups, as indicated by an overall increase in metal–ligand bond strength. Hydrogen bonding between the distal histidine and ligand also exhibited greater strength in the delta form; however, this effect was more pronounced with salen prosthetic groups. Additionally, the NCIs within the active pocket of the protein were found to be variable, highlighting the adaptability of local force constants. In summary, our data underscore the potential of computational methodologies in guiding the rational design of artificial metalloproteins for tailored applications, with local vibrational mode analysis serving as a powerful tool for bond strength assessment. 
    more » « less
  2. Abstract For a series of cytochrome b5 proteins isolated from various species, including bacteria, animals, and humans, we analyzed the intrinsic strength of their distal/proximal FeN bonds and the intrinsic stiffness of their axial NFeN bond angles. To assess intrinsic bond strength and bond angle stiffness, we employed local vibrational stretching force constants ka(FeN) and bending force constants ka(NFeN) derived from the local mode theory developed by our group; the ferric and ferrous oxidation states of the heme Fe were considered. All calculations were conducted with the QM/MM methodology. We found that the reduction of the heme Fe from the ferric to the ferrous state makes the FeN axial bonds weaker, longer, less covalent, and less polar. Additionally, the axial NFeN bond angle becomes stiffer and less flexible. Local mode force constants turned out to be far more sensitive to the protein environment than geometries; evaluating force constant trends across diverse protein groups and monitoring changes in the axial heme‐framework revealed redox‐induced changes to the primary coordination sphere of the protein. These results indicate that local mode force constants can serve as useful feature data for training machine learning models that predict cytochrome b5 redox potentials, which currently rely more on geometric data and qualitative descriptors of the protein environment. The insights gained through our investigation also offer valuable guidance for strategically fine‐tuning artificial cytochrome b5 proteins and designing new, versatile variants. 
    more » « less
  3. In our continued investigations of microbial globins, we solved the structure of a truncated hemoglobin from Shewanella benthica, an obligate psychropiezophilic bacterium. The distal side of the heme active site is lined mostly with hydrophobic residues, with the exception of a tyrosine, Tyr34 (CD1) and a histidine, His24 (B13). We found that purified SbHbN, when crystallized in the ferric form with polyethylene glycol as precipitant, turned into a green color over weeks. The electron density obtained from the green crystals accommodated a trans heme d, a chlorin-type derivative featuring a γ-spirolactone and a vicinal hydroxyl group on a pyrroline ring. In solution, exposure of the protein to one equivalent of hydrogen peroxide resulted in a similar green color change, but caused by the formation of multiple products. These were oxidation species released on protein denaturation, likely including heme d, and a species with heme covalently attached to the polypeptide. The Tyr34Phe replacement prevented the formation of both heme d and the covalent linkage. The ready modification of heme b by SbHbN expands the range of chemistries supported by the globin fold and offers a route to a novel heme cofactor. 
    more » « less
  4. Abstract We investigated the intrinsic strength of distal and proximal FeN bonds for both ferric and ferrous oxidation states of bishistidyl hemoproteins from bacteria, animals, human, and plants, including two cytoglobins, ten hemoglobins, two myoglobins, six neuroglobins, and six phytoglobins. As a qualified measure of bond strength, we used local vibrational force constants k(FeN) based on local mode theory developed in our group. All calculations were performed with a hybrid QM/MM ansatz. Starting geometries were taken from available x‐ray structures. k(FeN) values were correlated with FeN bond lengths and covalent bond character. We also investigated the stiffness of the axial NFeN bond angle. Our results highlight that protein effects are sensitively reflected in k(FeN), allowing one to compare trends in diverse protein groups. Moreover, k(NFeN) is a perfect tool to monitor changes in the axial heme framework caused by different protein environments as well as different Fe oxidation states. 
    more » « less
  5. Abstract Ergothioneine (ESH) and ovothiol A (OSHA) are two natural thiol‐histidine derivatives. ESH has been implicated as a longevity vitamin and OSHA inhibits the proliferation of hepatocarcinoma. The key biosynthetic step of ESH and OSHA in the aerobic pathways is the O2‐dependent C−S bond formation catalyzed by non‐heme iron enzymes (e.g., OvoA in ovothiol biosynthesis), but due to the lack of identification of key reactive intermediate the mechanism of this novel reaction is unresolved. In this study, we report the identification and characterization of a kinetically competentS=1 iron(IV) intermediate supported by a four‐histidine ligand environment (three from the protein residues and one from the substrate) in enabling C−S bond formation in OvoA fromMethyloversatilis thermotoleran, which represents the first experimentally observed intermediate spin iron(IV) species in non‐heme iron enzymes. Results reported in this study thus set the stage to further dissect the mechanism of enzymatic oxidative C−S bond formation in the OSHA biosynthesis pathway. They also afford new opportunities to study the structure‐function relationship of high‐valent iron intermediates supported by a histidine rich ligand environment. 
    more » « less