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. © 2017 Wiley Periodicals, Inc.

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.

Effect of a pH change on the conformational stability of the modified nucleotide queuosine monophosphate.

The naturally occurring modified nucleotide queuosine 5'-monophosphate (QMP) related to biochemical regulatory pathways in the cell was investigated using quantum chemical approaches. The relative stability of biologically relevant conformations of QMP in solvent under a pH change was predicted at the BVP86/TZVP and MP2/TZVP levels of theory. Hydrogen bonding in QMP was studied using Bader's approach. The acidity constants of QMP were estimated using the COSMO-RS theory. It has been found that the neutral and anionic forms of QMP are the most stable in the physiological pH range. These forms correspond to the anti/north conformation and exist as zwitterionic tautomers having a negatively charged phosphate group (-1 for neutral and -2 for anionic) and a positively charged secondary amine group in the side chain. It was also found that QMP possesses the syn conformation in the cationic state at pH < 5.0 and undergoes syn to anti conformation transition when the pH increases, remaining in the anti conformation at the higher pH values. The marker IR bands specific for the anionic and neutral QMP forms in the 2300-2700 cm(-1) region were assigned to H-bonded NH groups of the QMP side chain. The bands between 800 and 1300 cm(-1) of the "fingerprint" (400-1500 cm(-1)) region were assigned to the vibrations of the ribose ring, the phosphate group and the side chain of QMP. The predicted IR spectra can be useful for the assignment of vibration bands in the experiential spectra of QMP or identification of the QMP forms. The revealed peculiarities of the QMP conformation sensitivity to a pH change as well as additional formed H-bonds could be responsible for specific nucleotide interactions with enzymes.

CL-20 photodecomposition: ab initio foundations for identification of products.

1,5-Dihydrodiimidazo[4,5-b:4'5'e]pyrazine, 1H-imidazo[4,5-b]pyrazine, and 1H-imidazole were considered as possible products of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) photodecomposition. Since we took as a reference the product obtained after CL-20 irradiation in methanol solution, the nature of intermolecular bonds between heterocycles under study and methanol molecules was analyzed in detail. Existing hydrogen bonds were found to be quite strong, so dependence of calculations results on an influence of solvent was taken into account using both the polarizable continuum model (PCM) and the supermolecular approach. Electronic spectra of 1,5-dihydrodiimidazo[4,5-b:4'5'e]pyrazine, 1H-imidazo[4,5-b]pyrazine and 1H-imidazole were simulated using time dependent density functional theory (TD-DFT) and single-excitation configuration interaction (CIS) method. We observed that TD-DFT excitation energies are lower if compared to corresponding values obtained by the CIS method. Results of calculations with PCM and the supermolecular approach are very close. It was found that differences between calculated gas phase excitation energies and those values obtained by applying solvent models increases when the number of conjugated bonds in a molecule increases. Oscillator strengths of UV bands of the considered molecules are higher in the gas phase than in modeled methanol solutions. We found that the predicted spectrum of 1H-imidazole is in close agreement with the experimental UV spectrum of the CL-20 photolysis product.

Application of quantum-chemical approximations to environmental problems: Prediction of physical and chemical properties of TNT and related species.

This paper presents accurate predictions of ecologically important properties of nitroaromatic compounds and their derivatives, including vapor pressure, Henry's law constants, water solubility, octanol/water partition coefficients, heats of formation and ionization potentials. The proposed technique of calculations was based on quantum-chemical methods. The relationship between the chemical structure and mentioned physico-chemical parameters of such widespread military produced contaminants as trinitrotoluene and its derivatives was considered. We revealed that the DFT level of theory combined with the COSMO-RS technique is able to predict the studied parameters with an accuracy that results in error bars of less then one logarithmic unit.

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.

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.

Predicting water solubility of congeners: chloronaphthalenes--a case study.

Since the important physicochemical data for chloronaphtalenes (PCNs) are still scarce, we have predicted water solubility (logS) of all 75 congeners with the Quantitative Structure-Property Relationship (QSPR) scheme. The values of logS, predicted by the most efficient model, varied from 0.01 to 1660 microg dm(-3) (2.85 x 10(-11)-1.02 x 10(-5) mol dm(-3)), depending on the number of chlorine atoms present in the molecule and the substitution pattern. We found that the main factor determining relative differences in solubility between the congeners is the solvent accessible volume related to the cavitation process occurring in the solvent. The results are presented as a case study of QSPR modeling for those Persistent Organic Pollutants (POPs) that exist as families of congeners. By investigating the impact of (i) the way of the molecular descriptors' calculation, (ii) the size of applied database and (iii) chemometric method of modeling (Multiple Linear Regression, MLR, and/or Partial Least Squares regression, PLS) on the quality of the models we proposed general recommendations for dealing with congeners. We found that the combination of the B3LYP functional with 6-311++G(d,p) basis set was the most optimal technique of the molecular descriptors' calculation for congeners when comparing with semi-empirical PM3, ab initio Hartee-Fock (HF), and Møller-Pleset 2 (MP2) method carried out with different-size basis sets. Moreover, the model developed with a larger and more general database that includes chloronaphthalenes, polychlorinated dibezno-p-dioxins, furans and biphenyls predicted the values of logS for PCNs noticeable worse than the model calibrated only on PCNs. In the later case it was possible to obtain satisfactory results by employing even the simplest MLR method and only one molecular descriptor. The values of logS were also calculated with the WSKOWIN and COSMO-RS models as the reference techniques and then compared to our results.

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.

The mechanism of unimolecular decomposition of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane. A computational DFT study.

By using the B3LYP level of density functional theory, possible decomposition reaction pathways of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) in the gas phase have been investigated. We have found several types of reactions for this process: homolytic cleavage of an N-N bond to form the NO2* group; HONO elimination; C-C and C-N bonds breaking leading to ring opening; and H-migration. On the basis of the results of computation scanning of the potential energy surface, the most favorite pathway of CL-20 unimolecular decomposition that results in the formation of the stable aromatic compound 1,5-dihydrodiimidazo[4,5-b:4',5'-e]pyrazine has been proposed.