COSMO-RS based predictions for the SAMPL6 logP challenge.

Within the framework of the 6th physical property blind challenge (SAMPL6) the authors have participated in predicting the octanol-water partition coefficients (logP) for several small drug like molecules. Those logP values where experimentally known by the organizers but only revealed after the submissions of the predictions. Two different sets of predictions were submitted by the authors, both based on the COSMOtherm implementation of COSMO-RS theory. COSMOtherm predictions using the FINE parametrization level (hmz0n) obtained the highest accuracy among all submissions as measured by the root mean squared error. COSMOquick predictions using a fast algorithm to estimate σ-profiles and an a posterio machine learning correction on top of the COSMOtherm results (3vqbi) scored 3rd out of 91 submissions. Both results underline the high quality of COSMO-RS derived molecular free energies in solution.

Computational Screening of Drug Solvates.

PURPOSE: Solvates are mainly undesired by-products during the pharmaceutical development of new drugs. In addition, solvate formation may also distort solubility measurements. The presented study introduces a simple computational approach that allows for the identification of drug solvent pairs which most likely form crystalline solid phases. METHODS: The mixing enthalpy as a measure for drug-solvent complementarity is obtained by computational liquid phase thermodynamics (COSMO-RS theory). In addition a few other simple descriptors were taking into account describing the shape and topology of the drug and the solvent. Using an extensive dataset of drug solvent pairs a simple and statistically robust model is developed which allows for a rough assessment of a solvent's ability to form a solvate. RESULTS: Similar to the related issue of cocrystal screening, the mixing (or excess) enthalpy of the subcooled liquid mixture of the drug-solvent pair proves to be an important quantity controlling solvate formation. Due to the fact that many solvates form inclusion compounds, the solvent shape is another important factor influencing solvate formation. Solvates forming channel-like voids in the solid state are predicted less well. CONCLUSION: The approach ranks any drug-solvent pair that forms a solvate before any non-solvate by a probability of about 81% (AUC = 0.81), giving a significant advantage over any trial and error approach. Hence it can help to identify suitable solvent candidates early in the drug development process.

Solubility prediction, solvate and cocrystal screening as tools for rational crystal engineering.

OBJECTIVES: The fact that novel drug candidates are becoming increasingly insoluble is a major problem of current drug development. Computational tools may address this issue by screening for suitable solvents or by identifying potential novel cocrystal formers that increase bioavailability. In contrast to other more specialized methods, the fluid phase thermodynamics approach COSMO-RS (conductor-like screening model for real solvents) allows for a comprehensive treatment of drug solubility, solvate and cocrystal formation and many other thermodynamics properties in liquids. This article gives an overview of recent COSMO-RS developments that are of interest for drug development and contains several new application examples for solubility prediction and solvate/cocrystal screening. METHODS: For all property predictions COSMO-RS has been used. The basic concept of COSMO-RS consists of using the screening charge density as computed from first principles calculations in combination with fast statistical thermodynamics to compute the chemical potential of a compound in solution. KEY FINDING: The fast and accurate assessment of drug solubility and the identification of suitable solvents, solvate or cocrystal formers is nowadays possible and may be used to complement modern drug development. Efficiency is increased by avoiding costly quantum-chemical computations using a database of previously computed molecular fragments. SUMMARY: COSMO-RS theory can be applied to a range of physico-chemical properties, which are of interest in rational crystal engineering. Most notably, in combination with experimental reference data, accurate quantitative solubility predictions in any solvent or solvent mixture are possible. Additionally, COSMO-RS can be extended to the prediction of cocrystal formation, which results in considerable predictive accuracy concerning coformer screening. In a recent variant costly quantum chemical calculations are avoided resulting in a significant speed-up and ease-of-use.

Rational coformer or solvent selection for pharmaceutical cocrystallization or desolvation.

It is demonstrated that the fluid-phase thermodynamics theory conductor-like screening model for real solvents (COSMO-RS) as implemented in the COSMOtherm software can be used for accurate and efficient screening of coformers for active pharmaceutical ingredient (API) cocrystallization. The excess enthalpy, H(ex) , between an API-coformer mixture relative to the pure components reflects the tendency of those two compounds to cocrystallize. Thus, predictive calculations may be performed with decent effort on a large set of molecular data in order to identify potentially new cocrystal systems. In addition, it is demonstrated that COSMO-RS theory allows reasonable ranking of coformers for API solubility improvement. As a result, experiments may be focused on those coformers, which have an increased probability of cocrystallization, leading to the largest improvement of the API solubility. In a similar way as potential coformers are identified for cocrystallization, solvents that do not tend to form solvates may be determined based on the highest H(ex) s with the API. The approach was successfully tested on tyrosine kinase inhibitor axitinib, which has a propensity to form relatively stable solvated structures with the majority of common solvents, as well as on thiophanate-methyl and thiophanate-ethyl benzimidazole fungicides, which form channel solvates.

Performance of plane-wave-based LDA+U and GGA+U approaches to describe magnetic coupling in molecular systems.

This work explores the performance of periodic plane wave density functional theory calculations with an on-site Coulomb correction to the standard LDA and GGA exchange-correlation potential--commonly used to describe strongly correlated solids--in describing the magnetic coupling constant of a series of molecular compounds representative of dinuclear Cu complexes and of organic diradicals. The resulting LDA+U or GGA+U formalisms, lead to results comparable to experiment and to those obtained by means of standard hybrid functionals provided that the value of the U parameter is adequately chosen. Hence, these methods offer an alternative efficient computational scheme to correct LDA and GGA approaches to adequately describe the electronic structure and magnetic coupling in large molecular magnetic systems, although at the expenses of introducing an empirical (U) parameter. For all investigated copper dinuclear systems, the LDA+U and GGA+U approaches lead to an improvement in the description of magnetic properties over the original LDA and GGA schemes with an accuracy similar to that arising from the hybrid B3LYP functional, by increasing the on-site Coulomb repulsion with a moderate U value. Nevertheless, the introduction of an arbitrary U value in the 0-10 eV range most often provides the correct ground-state spin distribution and the correct sign of the magnetic coupling constant.

Density functional studies of model cerium oxide nanoparticles.

Density functional plane-wave calculations have been performed to investigate a series of ceria nanoparticles (CeO2-x)(n), n Ce3+ reduction have been accounted for through the use of an effective on-site Coulomb repulsive interaction within the so-called DFT+U approach. Twelve nanoparticles of up to 2 nm in diameter and of both cuboctahedral and octahedral forms are chosen as representative model systems. Energetic and structural effects of oxygen vacancy formation in these nanoparticles are discussed with respect to those in the bulk and on extended surfaces. We show that the average interatomic distances of the nanoparticles are most significantly affected by the creation of oxygen vacancies. The formation energies of non-stoichiometric nanoparticles (CeO2-x)(n) are found to scale linearly with the average coordination number of Ce atoms; where x < 0 species, containing partially reduced O atoms, are less stable. The stability of octahedral ceria particles at small sizes, and the predicted strong propensity of Ce cations to acquire a reduced state at lower coordinated sites, is supported by interatomic potential-based global optimisations probing the low energy isomers of the Ce19O32 nanoparticle.

Electronic structure of CO--an exercise in modern chemical bonding theory.

This paper discusses recent progress that has been made in the understanding of the electronic structure and bonding situation of carbon monoxide which was analyzed using modern quantum chemical methods. The new results are compared with standard models of chemical bonding. The electronic charge distribution and the dipole moment, the nature of the HOMO and the bond dissociation energy are discussed in detail.

A first structural and theoretical comparison of pyridinylidene-type rNHC (remote N-heterocyclic carbene) and NHC complexes of Ni(II) obtained by oxidative substitution.

A series of cationic pyridinylidene and quinolinylidene complexes of chlorobis(triphenylphosphine)-nickel(II) were prepared by oxidative substitution of Ni(PPh3)4 with methylated chloropyridines or chloroquinolines. NMR as well as X-ray crystallographic studies confirmed the trans arrangement of the two phosphines in the products. Calculations, using suitable model compounds at the BP86/TZVP level, clearly differentiate between a standard imidazolylidene complex and new complexes of the NHC-type on the one hand, and new complexes classified as rNHC-types-with the heteroatom distant from the carbene carbon-on the other. The latter form significantly stronger bonds-mainly of an electrostatic nature-with the metal.

Why are olefins oxidized by RuO4 under cleavage of the carbon-carbon bond whereas oxidation by OsO4 yields cis-diols?

Quantum chemical calculations using gradient-corrected (B3LYP) density functional theory have been carried out to investigate the mechanism of the oxidative cleavage of alkenes by ruthenium tetraoxide. The initial reaction of the tetraoxide with the olefin occurs via a [3+2] cycloaddition as in the case of osmium tetraoxide. The results clearly show that the bond cleavage does not take place at the primary adduct, but much later in the reaction path. After the formation of the ruthenium(VI)dioxo-2,5-dioxolane, the reaction proceeds with the addition of a second olefin to yield ruthenium(IV)-bis(2,5-dioxolane), which in turn becomes oxidized first to rutheniumoxo(VI)-bis(2,5-dioxolane) 6(Ru) and then to ruthenium(VIII)-dioxo-bis(2,5-dioxolane) 7(Ru). Only in complexes containing the metal center in the formal oxidation state +VIII are low activation barriers for C-C bond cleavage and exothermic formation of carbonyl compounds as products calculated. The lowest activation barrier, DeltaH(++) = 2.5 kcal/mol, is calculated for the C-C bond breaking reaction of 7(Ru) which is predicted as the pivotal intermediate of the oxidation reaction. The calculations of the oxidation reaction with OsO(4) show that those reactions where the oxidation state of the metal increases have larger activation barriers for M = Ru than for M = Os, while reactions which reduce the oxidation state have a lower activation barrier for ruthenium compounds. Also, reactions which increase the oxidation state of the metal are in the case of M = Os more exothermic than for M = Ru. In this work, all important points of the potential energy surface (PES) are reported, and the complete catalytic cycle for the oxidative cleavage of olefins by ruthenium tetraoxide is presented.

Quantum chemical investigations and bonding analysis of iron complexes with mixed cyano and carbonyl ligands.

The equilibrium structures and vibrational frequencies of the iron complexes [Fe(CN)(x)(CO)(y)](q) (x = 0-6 and y = 0-5) have been calculated at the BP86 level of theory. The nature of the Fe-CN and Fe-CO has been analyzed with an energy partitioning method. The calculated Fe-CO bond lengths are in good agreement with the results of X-ray structure analysis whereas the Fe-CN bonds are calculated somewhat longer than the experimental values. The theoretically predicted vibrational frequencies of the C-O stretching mode are always lower and the calculated CN(-) frequencies are higher than the observed fundamental modes. The results of the bonding analysis suggest that the Fe-CO binding interactions have approximately 55% electrostatic character and approximately 45% covalent character. There is a significant contribution of the pi orbital interaction to the Fe-CO covalent bonding which increases when the complexes become negatively charged. The strength of deltaE(pi) may even be larger than deltaE(sigma). The Fe-CN(-) bonds have much less pi character. The calculated binding energy of the Fe-CO pi-interactions correlates very well with the C-O stretching frequencies.