Structure and electrochemical properties for complexes of nitrocompounds with inorganic ions: A theoretical approach.

Reduction and oxidation (redox) reactions are widely used for removal of nitrocompounds from contaminated soil and water. Structures and redox properties for complexes of nitrocompounds, such as 2,4,6-trinitrotoluene (TNT), 2,4-dinitrotoluene (DNT), 2,4-dinitroanisole (DNAN), and 5-nitro-2,4-dihydro-3H-1,2,4-triazol-3-one (NTO), with common inorganic ions (Na(+) , Cl(-) , NO3-) were investigated at the SMD(Pauling)/PCM(Pauling)/MPWB1K/TZVP level of theory. Atoms in molecules (AIM) theory was applied to analyze the topological properties of the bond critical points involved in the interactions between the nitrocompounds and the ions. Topological analyses show that intermolecular interactions of the types O(N)…Na(+) , C-H…Cl(-) ( ONO2-), and C…Cl(-) ( ONO2-) may be discussed as noncovalent closed-shell interactions, while N-H···Cl(-) ( ONO2-) hydrogen bonds are partially covalent in nature. Complexation causes significant decrease of redox activity of the nitrocompounds. Analysis of the reduction potentials of the complexes obtained through application of the Pourbaix diagram of an iron/water system revealed that sodium complexes of NTO might be reduced by metallic iron. © 2016 Wiley Periodicals, Inc.

In Silico Alkaline Hydrolysis of Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine: Density Functional Theory Investigation.

HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine), an energetic material used in military applications, may be released to the environment during manufacturing, transportation, storage, training, and disposal. A detailed investigation of a possible mechanism of alkaline hydrolysis, as one of the most promising methods for HMX remediation, was performed by computational study at PCM(Pauling)/M06-2X/6-311++G(d,p) level. Obtained results suggest that HMX hydrolysis at pH 10 represents a highly exothermic multistep process involving initial deprotonation and nitrite elimination, hydroxide attachment accompanied by cycle cleavage, and further decomposition of cycle-opened intermediate to the products caused by a series of C-N bond ruptures, hydroxide attachments, and proton transfers. Computationally predicted products of HMX hydrolysis such as nitrite, 4-nitro-2,4-diazabutanal, formaldehyde, nitrous oxide, formate, and ammonia correspond to experimentally observed species. Based on computed reaction pathways for HMX decomposition by alkaline hydrolysis, the kinetics of the entire process was modeled. Very low efficiency of this reaction at pH 10 was observed. Computations predict significant increases (orders of magnitude) of the hydrolysis rate for hydrolysis reactions undertaken at pH 11, 12, and 13.

Are the reduction and oxidation properties of nitrocompounds dissolved in water different from those produced when adsorbed on a silica surface? A DFT M05-2X computational study.

The reduction and oxidation properties of four nitrocompounds (trinitrotoluene [TNT], 2,4-dinitrotoluene, 2,4-dinitroanisole, and 5-nitro-2,4-dihydro-3H-1,2,4-triazol-3-one [NTO]) dissolved in water as compared with the same properties for compounds adsorbed on a silica surface were studied. To consider the influence of adsorption, cluster models were developed at the M05/tzvp level. A hydroxylated silica (001) surface was chosen to represent a key component of soil. The PCM(Pauling) and SMD solvation models were used to model water bulk influence. The following properties were analyzed: electron affinity, ionization potential, reduction Gibbs free energy, oxidation Gibbs free energy, and reduction and oxidation potentials. It was found that adsorption and solvation decrease gas phase electron affinity, ionization potential, and Gibbs free energy of reduction and oxidation, and thus, promote redox transformation of nitrocompounds. However, in case of solvation, the changes are more significant than for adsorption. This means that nitrocompounds dissolved in water are easier to transform by reduction or oxidation than adsorbed ones. Among the considered compounds, TNT was found to be the most reactive in an electron attachment process and the least reactive for an electron detachment transformation. During ionization, a deprotonation of adsorbed NTO was found to occur.

Structure and Redox Properties of 5-Amino-3-nitro-1H-1,2,4-triazole (ANTA) Adsorbed on a Silica Surface: A DFT M05 Computational Study.

A cluster approximation was applied at the M05/tzvp level to model an adsorption of 5-amino-3-nitro-1H-1,2,4-triazole (ANTA) on the (001) surface of α-quartz. Structures of the obtained ANTA-silica complexes confirm a nearly parallel orientation of the nitro compound toward the surface. The atoms in molecules (AIM) method was applied to analyze binding between ANTA and the silica surface. Attachment or loss of an electron was found to lead to a significant deviation from coplanarity in the complexes and to a strengthening of a hydrogen bonding. Redox properties of the adsorbed ANTA were compared with those of gas-phase and hydrated species by calculation of the ionization potential, electron affinity, oxidation and reduction Gibbs free energies, and oxidation and reduction potentials. It was shown that the adsorbed ANTA has a lower ability to undergo redox transformations as compared to that of the hydrated one.

Comprehensive investigations of kinetics of alkaline hydrolysis of TNT (2,4,6-trinitrotoluene), DNT (2,4-dinitrotoluene), and DNAN (2,4-dinitroanisole).

Combined experimental and computational techniques were used to analyze multistep chemical reactions in the alkaline hydrolysis of three nitroaromatic compounds: 2,4,6-trinitrotoluene (TNT), 2,4-dinitrotoluene (DNT), and 2,4-dinitroanisole (DNAN). The study reveals common features and differences in the kinetic behavior of these compounds. The analysis of the predicted pathways includes modeling of the reactions, along with simulation of UV-vis spectra, experimental monitoring of reactions using LC/MS techniques, development of the kinetic model by designing and solving the system of differential equations, and obtaining computationally predicted kinetics for decay and accumulation of reactants and products. Obtained results suggest that DNT and DNAN are more resistant to alkaline hydrolysis than TNT. The direct substitution of a nitro group by a hydroxide represents the most favorable pathway for all considered compounds. The formation of Meisenheimer complexes leads to the kinetic first-step intermediates in the hydrolysis of TNT. Janovsky complexes can also be formed during hydrolysis of TNT and DNT but in small quantities. Methyl group abstraction is one of the suggested pathways of DNAN transformation during alkaline hydrolysis.

Theoretical study of ionization and one-electron oxidation potentials of N-heterocyclic compounds.

A number of density functionals was utilized to predict gas-phase adiabatic ionization potentials (IPs) for nitrogen-rich heterocyclic compounds. Various solvation models were applied to the calculation of difference in free energies of solvation of oxidized and reduced forms of heterocyclic compounds in acetonitrile (AN) for correct reproduction of their standard oxidation potentials. We developed generally applicable protocols that could successfully predict the gas-phase adiabatic ionization potentials of nitrogen-rich heterocyclic compounds and their standard oxidation potentials in AN. This approach is supported by a MPW1K/6-31+G(d) level of theory which uses SMD(UA0) approximation for estimation of solvation energy of neutral molecules and PCM(UA0) model for ionized ones. The mean absolute derivation (MAD) and root mean square error (RMSE) of the current theoretical models for IP are equal to 0.22 V and 0.26, respectively, and for oxidation potentials MAD = 0.13 V and RMSE = 0.17.

Mineral-organic interfacial processes: potential roles in the origins of life.

Life is believed to have originated on Earth ∼4.4-3.5 Ga ago, via processes in which organic compounds supplied by the environment self-organized, in some geochemical environmental niches, into systems capable of replication with hereditary mutation. This process is generally supposed to have occurred in an aqueous environment and, likely, in the presence of minerals. Mineral surfaces present rich opportunities for heterogeneous catalysis and concentration which may have significantly altered and directed the process of prebiotic organic complexification leading to life. We review here general concepts in prebiotic mineral-organic interfacial processes, as well as recent advances in the study of mineral surface-organic interactions of potential relevance to understanding the origin of life.

In silico structure-function analysis of E. cloacae nitroreductase.

Reduction, catalyzed by the bacterial nitroreductases, is the quintessential first step in the biodegradation of a variety of nitroaromatic compounds from contaminated waters and soil. The Enterobacter cloacae nitroreductase (EcNR) enzyme is considered as a prospective biotechnological tool for bioremediation of hazardous nitroaromatic compounds. Using diverse computational methods, we obtain insights into the structural basis of activity and mechanism of its function. We have performed molecular dynamics simulation of EcNR in three different states (free EcNR in oxidized form, fully reduced EcNR with benzoate inhibitor and fully reduced EcNR with nitrobenzene) in explicit solvent and with full electrostatics. Principal Component Analysis (PCA) of the variance-covariance matrix showed that the complexed nitroreductase becomes more flexible overall upon complexation, particularly helix H6, in the vicinity of the binding site. A multiple sequence alignment was also constructed in order to examine positional constraints on substitution in EcNR. Five regions which are highly conserved within the flavin mononucleotide (FMN) binding site were identified. Obtained results and their implications for EcNR functioning are discussed, and new plausible mechanism has been proposed.

Adsorption of nitrogen-containing compounds on hydroxylated α-quartz surfaces.

Adsorption energies of various nitrogen-containing compounds (specifically, 2,4,6-trinitrotoluene (TNT), 2,4-dinitrotoluene (DNT), 2,4-dinitroanisole (DNAn), and 3-nitro-1,2,4-triazole-5-one (NTO)) on the hydroxylated (001) and (100) α-quartz surfaces are computed. Different density functionals are utilized and both periodic as well as cluster approaches are applied. From the adsorption energies, partition coefficients on the considered α-quartz surfaces are derived. While TNT and DNT are preferably adsorbed on the (001) surface of α-quartz, NTO is rather located on both α-quartz surfaces.

Experimental and computational study of membrane affinity for selected energetic compounds.

The affinity of various energetic compounds for a biological membrane was investigated using experimental and computational techniques. We measured octanol-water (log(Kow)) and liposome-water (log(Klipw)) partition coefficients for the following chemicals: trinitrotoluene (TNT), 2,4-dinitrotoluene (2,4-DNT), 2,6-dinitrotoluene (2,6-DNT), 2,4-dinitroanisole (DNAN), 2methoxy-5-nitrophenol (2M5NP), 2,4,6-trinitrobenzene (TNB), and 2,4-dinitrophenol (2,4-DNP). In order to determine log(Klipw) experimentally, we used artificial solid supported lipid liposomes produced under trade mark TRANSIL. Log(Kow) value were predicted with several program packages including the COSMOthermX software. Log(Klipw) were estimated with COSMOmic as implemented in the COSMOthermX program package. In order to verify accuracy of our experimentally obtained results, we performed basic statistical analysis of data taken from the literature. We concluded that compounds considered in this study possess a moderate ability to penetrate into membranes. Comparison of both coefficients has shown that in general, the log(Kow) values are slightly smaller than log(Klipw).

Alkaline hydrolysis of hexahydro-1,3,5-trinitro-1,3,5-triazine: M06-2X investigation.

Alkaline hydrolysis mechanism of possible environmental contaminant RDX (hexahydro-1,3,5-trinitro-1,3,5-triazine) was investigated computationally at the PCM(Pauling)/M06-2X/6-311++G(d,p) level of theory. Results obtained show that the initial deprotonation of RDX by hydroxide leads to nitrite elimination and formation of a denitrated cyclohexene intermediate. Further nucleophilic attack by hydroxide onto cyclic CN double bond results in ring opening. It was shown that the presence of hydroxide is crucial for this stage of the reaction. The dominant decomposition pathway leading to a ring-opened intermediate was found to be formation of 4-nitro-2,4-diazabutanal. Hydrolytic transformation of its byproduct (methylene nitramine) leads to end products such as formaldehyde and nitrous oxide. Computational results are in a good agreement with experimental data on hydrolysis of RDX, suggesting that 4-nitro-2,4-diazabutanal, nitrite, formaldehyde, and nitrous oxide are main products for early stages of RDX decomposition under alkaline conditions.

Theoretical predictions of thermodynamic parameters of adsorption of nitrogen containing environmental contaminants on kaolinite.

In this study thermodynamic parameters of adsorption of nitrogen containing environmental contaminants (NCCs, 2,4,6, trinitrotoluene (TNT), 2,4-dinitrotoluene (DNT), 2,4-dinitroanisole (DNAN), and 3-one-1,2,4-triazol-5-one (NTO)) interacting with the tetrahedral and octahedral surfaces of kaolinite were predicted. Adsorption complexes were investigated using a density functional theory and both periodic and cluster approach. The complexes, modeled using the periodic boundary conditions approach, were fully optimized at the BLYP-D2 level to obtain the structures and adsorption energies. The relaxed kaolinite-NCCs structures were used to prepare cluster models to calculate thermodynamic parameters and partition coefficients at the M06-2X-D3 and BLYP-D2 levels from the gas phase. The entropy effect on the Gibbs free energies of adsorption of NCCS on kaolinite was also studied and compared with available experimental data. The results showed that in all calculated models, the NCCs molecules are physisorbed and they favor a parallel orientation toward both kaolinite surfaces. It was found that all calculated NCCs compounds are more stable on the octahedral than on the tetrahedral surface of kaolinite. The Gibbs free energies and partition coefficients were also predicted for interactions of NCCs with Na-kaolinite from aqueous solution. Calculations revealed adsorption of NCCs is effective from the gas phase on both cation free kaolinite surfaces and on Na-kaolinite from aqueous solution at room temperature. Theoretical data were validated against experimental results, and the reasons for small differences between calculated and measured partition coefficients are discussed.

Theoretical study of adsorption of nitrogen-containing environmental contaminants on kaolinite surfaces.

The adsorption of nitrogen-containing compounds (NCCs) including 2,4,6-trinitrotoluene (TNT), 2,4-dinitrotoluene (DNT), 2,4-dinitroanisole (DNAN), and 3-nitro-1,2,4-triazol-5-one (NTO) on kaolinite surfaces was investigated. The M06-2X and M06-2X-D3 density functionals were applied with the cluster approximation. Several different positions of NCCs relative to the adsorption sites of kaolinite were examined, including NCCs in perpendicular and parallel orientation toward both surface models of kaolinite. The binding between the target molecules and kaolinite surfaces was analyzed and bond energies were calculated applying the atoms in molecules (AIM) method. All NCCs were found to prefer a parallel orientation toward both kaolinite surfaces, and were bound more strongly to the octahedral than to the tetrahedral site. TNT exhibited the strongest interaction with the octahedral surface and DNAN with the tetrahedral surface of kaolinite. Hydrogen bonding was shown to be the dominant non-covalent interaction for NCCs interacting with the octahedral surface of kaolinite with a small stabilizing effect of dispersion interactions. In the case of adsorption on the tetrahedral surface, kaolonite-NCC binding was shown to be governed by the balance between hydrogen bonds and dispersion forces. The presence of water as a solvent leads to a significant decrease in the adsorption strength for all studied NCCs interacting with both kaolinite surfaces.

QSPR Prediction of Lipophilicity for Organic Compounds Using Random Forest Technique on the Basis of Simplex Representation of Molecular Structure.

The relationship between the octanol-water partition coefficient for more than twelve thousand organic compounds and their structures was investigated using a QSPR approach based on Simplex Representation of Molecular Structure (SiRMS). The dataset used in our study included 10973 compounds with experimental values of lipophilicity (LogKow ) for different chemical compounds. Random Forest (RF) method was used for statistical modeling at the 2D level of representation of molecular structure. Developed models are adequate and successfully validated with external test sets. Proposed models have clear interpretation due to the use of simplex representation of molecular structure and predict the LogKow values with the accuracy of the best modern models. Thus QSPR models proposed in this study represent powerful and easy-to use virtual screening tool that can be recommended for prediction of octanol-water partition coefficient.

DFT M06-2X investigation of alkaline hydrolysis of nitroaromatic compounds.

The nitroaromatic compounds 2,4,6-trinitrotoluene (TNT), 2,4-dinitrotoluene (DNT) and 2,4-dinitroanisole (DNAN) are potential environmental contaminants and their transformations under a variety of environmental conditions are consequently of great interest. One possible method to safely degrade these nitrocompounds is alkaline hydrolysis. A mechanism of the initial stages of this reaction was investigated computationally. Simulations of UV-VIS and NMR spectra for this mechanism were also produced. The results obtained were compared to available experimental data on the alkaline hydrolysis of TNT and suggest that the formation of Meisenheimer complexes and an anion of TNT are potential first-step intermediates in the reaction path. As the reaction proceeds, computational results indicate that polynegative complexes dominate the degradation pathway, followed by cycles of carbon chain opening and breaking. A second possible pathway was identified that leads to polymeric products through Janovsky complex formation. Results from this study indicate that the order of increasing resistance to alkaline hydrolysis is TNT, DNT and DNAN.

Molecular simulations of adsorption of RDX and TATP on IRMOF-1(Be).

The influence of different sorption sites of isoreticular metal-organic frameworks (IRMOFs) on interactions with explosive molecules is investigated. Different connector effects are taken into account by choosing IRMOF-1(Be) (IRMOF-1 with Zn replaced by Be), and two high explosive molecules: 1,3,5-trinitro-s-triazine (RDX) and triacetone triperoxide (TATP). The key interaction features (structural, electronic and energetic) of selected contaminants were analyzed by means of density functional calculations. The interaction of RDX and TATP with different IRMOF-1(Be) fragments is studied. The results show that physisorption is favored and occurs due to hydrogen bonding, which involves the C-H groups of both molecules and the carbonyl oxygen atoms of IRMOF-1(Be). Additional stabilization of RDX and TATP arises from weak electrostatic interactions. Interaction with IRMOF-1(Be) fragments leads to polarization of the target molecules. Of the molecular configurations we have studied, the Be-O-C cluster connected with six benzene linkers (1,4-benzenedicarboxylate, BDC), possesses the highest binding energy for the studied explosives (-16.4 kcal mol(-1) for RDX and -12.9 kcal mol(-1) for TATP). The main difference was discovered to be in the preferable adsorption site for adsorbates (RDX above the small and TATP placed above the big cage). Based on these results, IRMOF-1 can be suggested as an effective material for storage and also for separation of similar explosives. Hydration destabilizes most of the studied adsorption systems by 1-3 kcal mol(-1) but it leads to the same trend in the binding strength as found for the non-hydrated complexes.

Evaluation of the dependence of aqueous solubility of nitro compounds on temperature and salinity: a COSMO-RS simulation.

The solubility in pure and saline water at various temperatures was calculated for selected nitro compounds (nitrobenzene, 1,3,5-trinitrobenzene, 2-nitrotoluene, 3-nitrotoluene, 4-nitrotoluene, 2,4-dinitrotoluene, 2,6-dinitrotoluene, 2,3-dinitrotoluene, 3,4-dinitrotoluene, 2,4,6-trinitrotoluene) using the Conductor-like Screening model for Real Solvents (COSMO-RS). The results obtained were compared with experimental values. The COSMO-RS predictions have shown high accuracy in reproducing the trends of aqueous solubilities for both temperature and salinity. The proposed methodology was then applied to predict the aqueous solubilities of 19 nitro compounds in the temperature range of 5-50°C in saline solutions. The salting-out parameters of the Setschenow equation were also calculated. The predicted salting-out parameters were overestimated when compared to the measured values, but these parameters can still be used for qualitative estimation of the trends.

Application of quantum chemical approximations to environmental problems: prediction of water solubility for nitro compounds.

Water solubility values for 27 nitro compounds with experimentally measured values were computed using the conductor-like screening model for real solvent (COSMO-RS) based on the density functional theory and COSMO technique. We have found that the accuracy of the COSMO-RS approach for prediction of water solubility of liquid nitro compounds is impressively high (the errors are lower than 0.1 LU). However, for some solid nitro compounds, especially nitramines, there is sufficient disagreement between calculated and experimental values. In order to increase the accuracy of predictions the quantitative structure-property relationship (QSPR) part of the COSMO-RS approach has been modified. The solubility values calculated by the modified COSMO-RS method have shown much better agreement with the experimental values (the mean absolute errors are lower than 0.5 LU). Furthermore, this technique has been used for prediction of water solubility for an expanded set of 23 nitro compounds including nitroaromatic, nitramines, nitroanisoles, nitrogen rich compounds, and some their nitroso and amino derivatives with unknown experimental values. The solubility values predicted using the proposed computational technique could be useful for the determination of the environmental fate of military and industrial wastes and the development of remediation strategies for contaminated soils and waters. This predictive capability is especially important for unstable compounds and for compounds that have yet to be synthesized.