Abstract
Carbon monoxide (CO) is a byproduct of the incomplete combustion of carbon-based fuels, such as wood, coal, gasoline, or natural gas. As incomplete combustion in a fire accident or in an engine, massively produced CO leads to a serious life threat because CO competes with oxygen (O(2)) binding to hemoglobin and makes people suffer from hypoxia. Although there is hyperbaric O(2) therapy for patients with CO poisoning, the nanoscale mechanism of CO dissociation in the O(2)-rich environment is not completely understood. In this study, we construct the classical force field parameters compatible with the CHARMM for simulating the coordination interactions between hemoglobin, CO, and O(2), and use the force field to reveal the impact of O(2) on the binding strength between hemoglobin and CO. Density functional theory and Car-Parrinello molecular dynamics simulations are used to obtain the bond energy and equilibrium geometry, and we used machine learning enabled via a feedforward neural network model to obtain the classical force field parameters. We used steered molecular dynamics simulations with a force field to characterize the mechanical strength of the hemoglobin-CO bond before rupture under different simulated O(2)-rich environments. The results show that as O(2) approaches the Fe(2+) of heme at a distance smaller than ∼2.8 Å, the coordination bond between CO and Fe(2+) is reduced to 50% bond strength in terms of the peak force observed in the rupture process. This weakening effect is also shown by the free energy landscape measured by our metadynamics simulation. Our work suggests that the O(2)-rich environment around the hemoglobin-CO bond effectively weakens the bonding, so that designing of O(2) delivery vector to the site is helpful for alleviating CO binding, which may shed light on de novo drug design for CO poisoning.