An official website of the United States government
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
This content will become publicly available on January 27, 2026
Iron‐Histidine Coordination in Cytochrome b5: A Local Vibrational Mode Study
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
- Award ID(s):
- 2102461
- PAR ID:
- 10576959
- Publisher / Repository:
- Wiley Blackwell (John Wiley & Sons)
- Date Published:
- Journal Name:
- ChemPhysChem
- Volume:
- 26
- Issue:
- 7
- ISSN:
- 1439-4235
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
ABSTRACT Using the QM/MM methodology and a local mode analysis, we investigated a character and a strength of FeS bonds of heme groups in oxidized and reduced forms of Bacterioferritin fromAzotobacter vinelandii. The strength of the FeS bonds was correlated with a bond length, an energy density at a bond critical point, and a charge difference of the F and S atoms. Changing the oxidation state from ferrous to ferric generally makes the FeS bonds weaker, longer, more covalent, and more polar. We also investigated the SFeS bond bending and found that the stronger FeS bond, generally makes the SFeS bond bending stiffer, which could play a key role in the balance between ferric and ferrous oxidation states and related biological activities.more » « less
-
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
-
Abstract Iron oxide is commonly found in natural or industrial glass compositions and can exist as ferrous (Fe2+) or ferric (Fe3+) species, with their ratios depending on glass composition, temperature, pressure and the redox reactions during the glass forming process. The iron redox ratio plays an important role on silicate glass structures and consequently various properties. This work aims to study the effect of iron oxide, and particularly the iron redox ratio, on the structures of borosilicate and boroaluminosilicate glasses using molecular dynamics simulations with newly developed iron potential parameters that are compatible with the borosilicate potentials. The results provide detailed cation coordination states of both iron species and the effect of redox ratio on boron coordination and other structural features. Particularly, competition for charge compensating modifier cations (such as Na+) among the fourfold‐coordinated cations such as B3+, Al3+, and Fe3+is investigated by calculating the cation–cation pair distribution functions and coordination preferential ratios. The results show that the trivalent ferric ions, with a shorter Fe–O bond distance and better defined first coordiation shell with mainly four‐fold coordination, act as a glass former whereas the divalent ferrous ions mainly play the role of glass modifier. The ferrous/ferric ratio (Fe2+/Fe3+) was found to affect the glass chemistry and hence glass properties by regulating the amount of four‐coordinated boron, the fraction of non‐briding oxygen and other features. The results are compared with available experimental data to gain insights of the complex structures and charge compensation schemes of the glass system.more » « less
-
Information resulting from a comprehensive investigation into the intrinsic strengths of hydrated divalent magnesium clusters is useful for elucidating the role of aqueous solvents on the Mg2+ ion, which can be related to those in bulk aqueous solution. However, the intrinsic Mg–O and intermolecular hydrogen bond interactions of hydrated magnesium ion clusters have yet to be quantitatively measured. In this work, we investigated a set of 17 hydrated divalent magnesium clusters by means of local vibrational mode force constants calculated at the ωB97X-D/6-311++G(d,p) level of theory, where the nature of the ion–solvent and solvent–solvent interactions were interpreted from topological electron density analysis and natural population analysis. We found the intrinsic strength of inner shell Mg–O interactions for [Mg(H2O)n]2+ (n = 1–6) clusters to relate to the electron density at the bond critical point in Mg–O bonds. From the application of a secondary hydration shell to [Mg(H2O)n]2+ (n = 5–6) clusters, stronger Mg–O interactions were observed to correspond to larger instances of charge transfer between the lp(O) orbitals of the inner hydration shell and the unfilled valence shell of Mg. As the charge transfer between water molecules of the first and second solvent shell increased, so did the strength of their intermolecular hydrogen bonds (HBs). Cumulative local vibrational mode force constants of explicitly solvated Mg2+, having an outer hydration shell, reveal a CN of 5, rather than a CN of 6, to yield slightly more stable configurations in some instances. However, the cumulative local mode stretching force constants of implicitly solvated Mg2+ show the six-coordinated cluster to be the most stable. These results show that such intrinsic bond strength measures for Mg–O and HBs offer an effective way for determining the coordination number of hydrated magnesium ion clusters.more » « less
