Theoretical Modeling of Structure-Toxicity Relationship of Cyanides

Marcin Molski
DOI: https://doi.org/10.48550/arXiv.2104.10031
2021-04-15
Abstract:The global descriptors of chemical activity:ionization potential IP, electron affinity EA, chemical potential mu, absolute electronegativity chi, molecular hardness eta and softness S,electrophilicity index omega for cyanides X(CN)k with X=H, Na, K, Ag, Cu, Ca, Hg, Cd, Zn in the gas phase and water medium have been determined by taking advantage of the quantum-chemical <a class="link-external link-http" href="http://computations.To" rel="external noopener nofollow">this http URL</a> this aim, the HOMO and LUMO energy levels were calculated using DFT B3LYP method and QZVP (Valence Quadruple-Zeta Polarization) basis set, which enables precise calculations for hydrogen cyanide and its salts containing both light (H, Na, Ca) and heavy (K, Ag, Cu, Cd, Hg, Zn) <a class="link-external link-http" href="http://atoms.The" rel="external noopener nofollow">this http URL</a> results obtained indicate that while the EA-parameter roughly determines the LD50 values for the cyanides considered, the omega-descriptor is related rather to the product of cyanide LD50 and hydrolysis <a class="link-external link-http" href="http://n-degree.Hence" rel="external noopener nofollow">this http URL</a>,the theoretical LD50CN(omega) function proposed is interconnected with the nLD50 collective variable, whereas LD50(EA)directly approximates LD50 indicating that the toxicity of cyanides decreases with increasing EA, <a class="link-external link-http" href="http://omega-values.The" rel="external noopener nofollow">this http URL</a> calculations carried out suggest that some of the LD50(Exp)experimentally determined are incorrect (AgCN) or inaccurate (Cd(CN)2) and require <a class="link-external link-http" href="http://revision.Comparison" rel="external noopener nofollow">this http URL</a> of the theoretically estimated LD50(omega) with NOAEL toxicity parameters indicates that they are well correlated in contrast to LD50(EA) and LD50(Exp), exhibiting lower and marginal degrees of correlation, respectively.
Chemical Physics,Other Quantitative Biology
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to explore the relationship between the structure and toxicity of cyanide through theoretical calculation methods. The specific objectives include: 1. **Determination of global descriptors of cyanide**: Through quantum - chemical calculation methods (such as DFT B3LYP method and QZVP basis set), calculate the HOMO and LUMO energy levels of cyanide (such as HCN and its salts) in gas phase and aqueous medium, and then determine its global descriptors, including ionization potential (IP), electron affinity (EA), chemical potential (μ), absolute electronegativity (χ), molecular hardness (η), softness (S) and electrophilic index (ω). 2. **Establishment of a theoretical model for cyanide toxicity**: Study the relationship between these global descriptors and the median lethal dose (LD50) of cyanide, especially explore the influence of EA and ω descriptors on cyanide toxicity. The author proposes two theoretical functions \(T(EA)\) and \(T(ω)\) to describe the non - linear relationship between cyanide toxicity and EA and ω. 3. **Verification of the accuracy of experimental data**: By comparing the LD50 values obtained from theoretical calculations with experimental data, evaluate the accuracy of some experimentally determined LD50 values. For example, the paper points out that the experimental LD50 value of AgCN may be incorrect and needs to be re - evaluated. 4. **Exploration of more reliable toxicity parameters**: Compare the theoretically calculated LD50 values with more advanced toxicity parameters (such as NOAEL and AEGL) to evaluate the correlation between these parameters and explore their application value in cyanide toxicity prediction. ### Formula summary - **HOMO - LUMO energy gap**: \[ \Delta E=E_{\text{LUMO}}-E_{\text{HOMO}} \] - **Ionization potential**: \[ IP = - E_{\text{HOMO}} \] - **Electron affinity**: \[ EA=-E_{\text{LUMO}} \] - **Chemical hardness**: \[ \eta\approx\frac{1}{2}(E_{\text{LUMO}}-E_{\text{HOMO}}) \] - **Chemical softness**: \[ S\approx\frac{1}{E_{\text{LUMO}}-E_{\text{HOMO}}} \] - **Chemical potential**: \[ \mu\approx\frac{1}{2}(E_{\text{LUMO}}+E_{\text{HOMO}}) \] - **Electronegativity**: \[ \chi\approx\frac{1}{2}(E_{\text{LUMO}}+E_{\text{HOMO}})=\frac{1}{2}(IP - EA) \] - **Electrophilic index**: \[ \omega\approx\frac{\chi^{2}}{2\eta}=\frac{(E_{\text{LUMO}}+E_{\text{HOMO}})^{2}}{8(E_{\text{LUMO}}-E_{\text{HOMO}})} \] - **Toxicity function**: \[ T(EA)=T_{0}\exp(a(EA - EA_{0})) \] \[ T(\omega)=T_{0}\exp(a(\omega-\omega_{0})) \] Through these calculations and analyses, the paper aims to provide a reliable method for predicting the toxicity of cyanide and provide a theoretical basis for the accuracy of experimental data.