Predicting Cytotoxicity of Metal Oxide Nanoparticles using Isalos Analytics Platform.

A literature curated dataset containing 24 distinct metal oxide (MexOy) nanoparticles (NPs), including 15 physicochemical, structural and assay-related descriptors, was enriched with 62 atomistic computational descriptors and exploited to produce a robust and validated in silico model for prediction of NP cytotoxicity. The model can be used to predict the cytotoxicity (cell viability) of MexOy NPs based on the colorimetric lactate dehydrogenase (LDH) assay and the luminometric adenosine triphosphate (ATP) assay, both of which quantify irreversible cell membrane damage. Out of the 77 total descriptors used, 7 were identified as being significant for induction of cytotoxicity by MexOy NPs. These were NP core size, hydrodynamic size, assay type, exposure dose, the energy of the MexOy conduction band (EC), the coordination number of the metal atoms on the NP surface (Avg. C.N. Me atoms surface) and the average force vector surface normal component of all metal atoms (v⟂ Me atoms surface). The significance and effect of these descriptors is discussed to demonstrate their direct correlation with cytotoxicity. The produced model has been made publicly available by the Horizon 2020 (H2020) NanoSolveIT project and will be added to the project's Integrated Approach to Testing and Assessment (IATA).

Fe-Doped ZnO nanoparticle toxicity: assessment by a new generation of nanodescriptors.

In the search for novel tools to combat cancer, nanoparticles (NPs) have attracted a lot of attention. Recently, the controlled release of cancer-cell-killing metal ions from doped NPs has shown promise, but fine tuning of dissolution kinetics is required to ensure specificity and minimize undesirable toxic side-effects. Theoretical tools to help in reaching a proper understanding and finally be able to control the dissolution kinetics by NP design have not been available until now. Here, we present a novel set of true nanodescriptors to analyze the charge distribution, the effect of doping and surface coating of whole metal oxide NP structures. The polarizable model of oxygen atoms enables light to be shed on the charge distribution on the NP surface, allowing the in detail study of the factors influencing the release of metal ions from NPs. The descriptors and their capabilities are demonstrated on a Fe-doped ZnO nanoparticle system, a system with practical outlook and available experimental data.

Theoretical modeling of HPV: QSAR and novodesign with fragment approach.

Structure-activity relationships in a data set of HPV6-E1 helicase ATPase inhibitors were investigated based on two different sets of descriptors. Statistically significant four parameter Quantitative Structure-Activity Relationships (QSAR) models were constructed and validated in both cases (R(2)=0.849; R(2) cv=0.811; F=52.20; s(2)=0.25; N=42). A Fragment based QSAR (FQSAR) approach was applied for developing a fragment-QSAR equation, which enabled the construction of virtual structures for novel ATPase inhibitors with desired or pre-defined activity.

Gas-phase lithium cation basicity: revisiting the high basicity range by experiment and theory.

According to high level calculations, the upper part of the previously published FT-ICR lithium cation basicity (LiCB at 373 K) scale appeared to be biased by a systematic downward shift. The purpose of this work was to determine the source of this systematic difference. New experimental LiCB values at 373 K have been measured for 31 ligands by proton-transfer equilibrium techniques, ranging from tetrahydrofuran (137.2 kJ mol(-1)) to 1,2-dimethoxyethane (202.7 kJ mol(-1)). The relative basicities (ΔLiCB) were included in a single self-consistent ladder anchored to the absolute LiCB value of pyridine (146.7 kJ mol(-1)). This new LiCB scale exhibits a good agreement with theoretical values obtained at G2(MP2) level. By means of kinetic modeling, it was also shown that equilibrium measurements can be performed in spite of the formation of Li(+) bound dimers. The key feature for achieving accurate equilibrium measurements is the ion trapping time. The potential causes of discrepancies between the new data and previous experimental measurements were analyzed. It was concluded that the disagreement essentially finds its origin in the estimation of temperature and the calibration of Cook's kinetic method.

Catalytic effect of cesium cation adduct formation on the decarboxylation of carboxylate ions in the gas phase.

The effect of Cs(+) ligation on the decarboxylation of malonic acids (unsubstituted and methyl-, dimethyl-, ethyl-, and phenyl-substituted) in their carboxylate form was studied in the gas phase using tandem mass spectrometry. The study is based on the comparison of the decarboxylation of the bare monoanion (hydrogen malonates) and of the cesium adduct of the cesium salt (Cs(+) [cesium hydrogen malonates]) under collisional activation. Energy-resolved dissociation curves of the negative and positive ions exhibit major differences. Decarboxylation of the cationic adducts of substituted malonic acid salts occurs at significantly lower collisional activation than for the corresponding bare hydrogen malonate anions. The conclusions from these experiments are supported by DFT calculations. The calculated activation parameters (enthalpy and Gibbs energy) confirm that the cesium cation coordination assists the decarboxylation of the carboxylate form.

Gas phase acidities of N-substituted amine-boranes.

Complexation energies and acidities of 19 primary, secondary and tertiary amine-boranes were investigated using MP2/6-311+G(d,p) and B3LYP/6-311+G(d,p) methods. Gas phase acidities for free amines were also calculated. Acidity values for studied complexes range from 327.3 to 349.1 kcal mol(-1) and the most acidic are the ones with direct connection between deprotonation center and a π-system. Results obtained by both computational methods are in good agreement with each other and with known experimental data. Addition of BH3 increases the acidity of amines by 30 to 50 kcal mol(-1). This enhancement effect was compared to the respective effect witnessed in phosphine-boranes and traced back to changes of charge delocalization on nitrogen. A question about the structural stability of several deprotonated amine-borane anions in the gas phase was also raised.

Critical Test of Some Computational Chemistry Methods for Prediction of Gas-Phase Acidities and Basicities.

Gas-phase acidities and basicities were calculated for 64 neutral bases (covering the scale from 139.9 kcal/mol to 251.9 kcal/mol) and 53 neutral acids (covering the scale from 299.5 kcal/mol to 411.7 kcal/mol). The following methods were used: AM1, PM3, PM6, PDDG, G2, G2MP2, G3, G3MP2, G4, G4MP2, CBS-QB3, B1B95, B2PLYP, B2PLYPD, B3LYP, B3PW91, B97D, B98, BLYP, BMK, BP86, CAM-B3LYP, HSEh1PBE, M06, M062X, M06HF, M06L, mPW2PLYP, mPW2PLYPD, O3LYP, OLYP, PBE1PBE, PBEPBE, tHCTHhyb, TPSSh, VSXC, X3LYP. The addition of the Grimmes empirical dispersion correction (D) to B2PLYP and mPW2PLYP was evaluated, and it was found that adding this correction gave more-accurate results when considering acidities. Calculations with B3LYP, B97D, BLYP, B2PLYPD, and PBE1PBE methods were carried out with five basis sets (6-311G**, 6-311+G**, TZVP, cc-pVTZ, and aug-cc-pVTZ) to evaluate the effect of basis sets on the accuracy of calculations. It was found that the best basis sets when considering accuracy of results and needed time were 6-311+G** and TZVP. Among semiempirical methods AM1 had the best ability to reproduce experimental acidities and basicities (the mean absolute error (mae) was 7.3 kcal/mol). Among DFT methods the best method considering accuracy, robustness, and computation time was PBE1PBE/6-311+G** (mae = 2.7 kcal/mol). Four Gaussian-type methods (G2, G2MP2, G4, and G4MP2) gave similar results to each other (mae = 2.3 kcal/mol). Gaussian-type methods are quite accurate, but their downside is the relatively long computational time.

Computational study of the Sonogashira cross-coupling reaction in the gas phase and in dichloromethane solution.

The Sonogashira cross-couplig reaction, consisting of oxidative addition, cis-trans isomerization, transmetalation, and reductive elimination, was computationally modeled using the DFT B3LYP/cc-pVDZ method for reaction between bromobenzene and phenylacetylene. Palladium diphosphane was used as a catalyst, copper(I) bromide as a co-catalyst and trimethylamine as a base. The reaction mechanism was studied both in the gas phase and in dichloromethane solution using PCM method. The complete catalytic cycle is thermodynamically strongly shifted toward products (diphenylacetylene and regenerated palladium catalyst) and is exothermic being in accordance with experimental data. The rate-determining step is the oxidative addition, since the highest point on the Gibbs energy graph of the complete reaction is the transition state of this step. This conclusion is also supported by recent experimental data. The computed energy profile suggests that the transmetalation step is initiated by the dissociation of neutral ligand, while the activation Gibbs energy of this step is 0.1 kcal mol(-1) in the gas phase.

Gas-phase basicities around and below water revisited.

This work employs Fourier transform ion cyclotron resonance (FT-ICR) and the Gaussian quantum chemistry composite methods W1 and G2 to experimentally and computationally analyze gas-phase basicities (GB) for a series of weak bases in the basicity region around and below water. The study aims to clarify the long-standing discrepancy between reported GB values for weak bases obtained via high-pressure mass spectrometry (HPMS) and ICR; the ICR scale is observed to be more than 2 times contracted compared to the HPMS scale. The computational results of this work support published HPMS data. This agreement improves with increasing sophistication of the computational method and is excellent at the W1 level. Several equilibria were also re-examined experimentally using FT-ICR. In the experiments with some polyfluorinated weak bases (hexafluoro-2-propanol and nonafluoro-2-methyl-2-propanol), it was found that two protonation processes compete in the gas phase: protonation on oxygen and protonation on fluorine. In these species, protonation on fluorine proceeds faster and is statistically favored over protonation on oxygen but leads to cations that are thermodynamically less stable than oxygen-protonated cations. The process may also lead to the irreversible loss of HF. The rearrangement of fluorine-protonated cations to oxygen-protonated cations is very slow and is further suppressed by the process of HF abstraction. These results at least partially explain the discrepancy between published HPMS data and earlier FT-ICR findings and call for the utmost care in using FT-ICR for gas-phase basicity measurements of heavily fluorinated compounds. The narrower dynamic range of ICR necessitates the measurement of several problematic bases and produces some differences between the ICR results in the present work and the published HPMS data; the wider dynamic range allows HPMS to overcome these difficulties in connecting the ladder.

Interaction of the cesium cation with mono-, di-, and tricarboxylic acids in the gas phase. A Cs+ affinity scale for cesium carboxylates ion pairs.

Humic substances (HS), including humic and fulvic acids, play a significant role in the fate of metals in soils. The interaction of metal cations with HS occurs predominantly through the ionized (anionic) acidic functions. In the context of the effect of HS on transport of radioactive cesium isotopes in soils, a study of the interaction between the cesium cation and model carboxylic acids was undertaken. Structure and energetics of the adducts formed between Cs+ and cesium carboxylate salts [Cs+RCOO-] were studied by the kinetic method and density functional theory (DFT). Clusters generated by electrospray ionization mass spectrometry from mixtures of a cesium salt (nitrate, iodide, trifluoroacetate) and carboxylic acids were quantitatively studied by CID. By combining the results of the kinetic method and the energetic data from DFT calculations, a scale of cesium cation affinity, CsCA, was built for 33 cesium carboxylates representing the first scale of cation affinity of molecular salts. The structural effects on the CsCA values are discussed.