Yana Kholod

Ab Initio Kinetic Simulation of Gas-Phase Experiments: Tautomerization of Cytosine and Guanine

A novel kinetic approach based on ab initio calculated rate constants has been developed and implemented in the kTSim program. The proposed approach allows prediction of the distribution of reactant and product concentrations over time, based exclusively on computationally obtained rate constants. The newly developed methodology was used to simulate the process of evaporation and tautomerization of guanine and cytosine under thermal (T = 490 K, cytosine; T = 620 K, guanine) and laser (T = 1000 K, 24 ns laser pulse) desorption conditions. Both monomolecular and bimolecular mechanisms of the tautomerization were considered simultaneously. The rates of the reactions were estimated using the values of Gibbs free energies calculated at the MPWB1K/aug-cc-pVDZ level and specified in a kTSim input. We expect that the proposed approach can also be used for accurate kinetic simulation of a wide range of processes.

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.

Prediction of CL-20 chemical degradation pathways, theoretical and experimental evidence for dependence on competing modes of reaction

Highest occupied and lowest unoccupied molecular orbital energies, formation energies, bond lengths and FTIR spectra all suggest competing CL-20 degradation mechanisms. This second of two studies investigates recalcitrant, toxic, aromatic CL-20 intermediates that absorb from 370 to 430 nm. Our earlier study (Struct. Chem., 15, 2004) revealed that these intermediates were formed at high OH− concentrations via the chemically preferred pathway of breaking the C–C bond between the two cyclopentanes, thereby eliminating nitro groups, forming conjugated π bonds, and resulting in a pyrazine three-ring aromatic intermediate. In attempting to find and make dominant a more benign CL-20 transformation pathway, this current research validates hydroxylation results from both studies and examines CL-20 transformations via photo-induced free radical reactions. This article discusses CL-20 competing modes of degradation revealed through: computational calculation; UV/VIS and SF spectroscopy following alkaline hydrolysis; and photochemical irradiation to degrade CL-20 and its byproducts at their respective wavelengths of maximum absorption.

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.

Ab Initio Study of Molecular Interactions in Cellulose Iα

Biomass recalcitrance, the resistance of cellulosic biomass to degradation, is due in part to the stability of the hydrogen bond network and stacking forces between the polysaccharide chains in cellulose microfibers. The fragment molecular orbital (FMO) method at the correlated Møller-Plesset second order perturbation level of theory was used on a model of the crystalline cellulose Iα core with a total of 144 glucose units. These computations show that the intersheet chain interactions are stronger than the intrasheet chain interactions for the crystalline structure, while they are more similar to each other for a relaxed structure. An FMO chain pair interaction energy decomposition analysis for both the crystal and relaxed structures reveals an intricate interplay between electrostatic, dispersion, charge transfer, and exchange repulsion effects. The role of the primary alcohol groups in stabilizing the interchain hydrogen bond network in the inner sheet of the crystal and relaxed structures of cellulose Iα, where edge effects are absent, was analyzed. The maximum attractive intrasheet interaction is observed for the GT-TG residue pair with one intrasheet hydrogen bond, suggesting that the relative orientation of the residues is as important as the hydrogen bond network in strengthening the interaction between the residues.

A coarse-grained model for β-d-glucose based on force matching

Cellulosic ethanol production is a two-stage process that involves the hydrolysis of cellulose to form simple sugars and the fermentation of these sugars to ethanol. Hydrolysis of cellulose is the rate-limiting step, and there is a great need to characterize the process with numerical simulations to better understand the complex mechanisms involved. The ultimate goal is to generate accurate coarse-grained molecular models that are capable of predicting the structure of lignocellulose before and after pretreatment so that subsequent ab initio calculations can be performed to probe the degradation pathways. As a first step toward that goal, the force-matching method is used to derive coarse-grained models for β-d-glucose molecules in aqueous solution. Using the same reference, an all-atom molecular dynamics simulation trajectory, two sets of three- and six-site coarse-grained models of β-d-glucose are developed using two definitions of the coarse-grained center site location: center of mass (CG-CM) and geometric center (CG-GC). The performance of these coarse-grained models is evaluated by comparing the coarse-grained predictions for bond-length distributions and radial distribution functions to those obtained from the all-atom reference simulation. The six-site coarse-grained models retain more structural details than the three-site coarse-grained models. Comparison between center site definitions shows that CG-CM models generally predict local ordering better, while CG-GC models predict long-range structure better.

Progress in Predictions of Environmentally Important Physicochemical Properties of Energetic Materials: Applications of Quantum-Chemical Calculations

The review describes the advances of quantum-chemically based approximations (namely, COSMO-RS) in the prediction of several environmentally important physicochemical properties of energetic materials: vapor pressure, Henry’s law constants, water solubility, and octanol–water partition coefficients. It includes introduction, the section that briefly discusses COSMO-RS – the most popular quantum-chemistry-based statistical thermodynamics approximation, and the references on similar quantum-chemical approaches. Since the solubility, probably, plays the most important role in many environmental characteristics of energetic materials, the major section of the review describes the current status of the quantum-chemically based predictions of this property. Also, the description of a modeling of salinity effects is discussed. Then subsequent few sections review the current advancements of the calculations of other environmentally important physical properties of energetic compounds.

Excitation energy transfer pathways in light-harvesting proteins: Modeling with PyFREC

Excitation energy transfer (EET) determines the fate of sunlight energy absorbed by light‐harvesting proteins in natural photosynthetic systems and photovoltaic cells. As previously reported (D. Kosenkov, J. Comput. Chem. 2016, 37(19), 1847), PyFREC software enables computation of electronic couplings between organic molecules with a molecular fragmentation approach. The present work reports implementation of direct fragmentation‐based computation of the electronic couplings and EET rates in pigment–protein complexes within the Förster theory in PyFREC. The new feature enables assessment of EET pathways in a wide range of photosynthetic complexes, as well as artificial molecular architectures that include light‐harvesting proteins or tagged fluorescent biomolecules. The developed methodology has been tested analyzing EET in the Fenna–Matthews–Olson (FMO) pigment–protein complex. The pathways of excitation energy transfer in FMO have been identified based on the kinetics studies.

Evaluation of Proton Transfer in DNA Constituents: Development and Application of Ab Initio Based Reaction Kinetics

The kinetics of chemical reactions characterizes the rates of chemical processes, i.e. distribution of all reactants, intermediates and products over time. This information is of vital importance for all areas of chemistry: chemical technology to control organic or inorganic syntheses, chemical construction of nanomaterials, as well as for the investigation of biochemical processes. The chemical kinetics data provide a possibility to investigate the effect of different chemical, physical and environmental factors on the rate of a reaction, final products and by-products distribution, and even the direction of a chemical process. In the first part of the chapter the general introduction to the kinetics of chemical reactions is given. The classical kinetics of chemical reactions uses the outcome from experimental measurement of reaction rates. However, currently available reliable computational ab initio methods provide an alternative efficient way for estimation of the rate constants even for stepwise and multidirectional reactions. Another benefit of the computational investigations is the possibility to simulate a wide range of processes with duration from picoseconds to hours, days, or even for much longer time scales. Contemporary ab initio methods have been used for estimation and prediction of reaction rates for a number of different chemical reactions. Until recently most of the theoretical studies on kinetic parameters have not been extended beyond the calculations of the rate constants of chemical reactions. In the present review we describe the simulation of the chemical kinetics of proton transfer (tautomerization) in nucleic acid bases and their complexes with metal ions, also in the presence of water molecules. The considered models are based on the ab initio calculated rate constants of chemical reactions. Then, such predicted rate constants are used for further kinetic simulations. Biological consequences of investigated processes are also discussed.