Molecular theory of solvation: Methodology summary and illustrations

Integral equation theory of molecular liquids based on statistical mechanics is quite promising as an essential part of multiscale methodology for chemical and biomolecular nanosystems in solution. Beginning with a molecular interaction potential force field, it uses diagrammatic analysis of the solvation free energy to derive integral equations for correlation functions between molecules in solution in the statistical-mechanical ensemble. The infinite chain of coupled integral equations for many-body correlation functions is reduced to a tractable form for 2- or 3-body correlations by applying the so-called closure relations. Solving these equations produces the solvation structure with accuracy comparable to molecular simulations that have converged but has a critical advantage of readily treating the effects and processes spanning over a large space and slow time scales, by far not feasible for explicit solvent molecular simulations. One of the versions of this formalism, the three-dimensional reference interaction site model (3D-RISM) integral equation complemented with the Kovalenko-Hirata (KH) closure approximation, yields the solvation structure in terms of 3D maps of correlation functions, including density distributions, of solvent interaction sites around a solute (supra)molecule with full consistent account for the effects of chemical functionalities of all species in the solution. The solvation free energy and the subsequent thermodynamics are then obtained at once as a simple integral of the 3D correlation functions by performing thermodynamic integration analytically.Comment: 24 pages, 10 figures, Revie

Computation of liquid-liquid equilibria and phase stabilities: implications for RH-dependent gas/particle partitioning of organic-inorganic aerosols

Semivolatile organic and inorganic aerosol species partition between the gas and aerosol particle phases to maintain thermodynamic equilibrium. Liquid-liquid phase separation into an organic-rich and an aqueous electrolyte phase can occur in the aerosol as a result of the salting-out effect. Such liquid-liquid equilibria (LLE) affect the gas/particle partitioning of the different semivolatile compounds and might significantly alter both particle mass and composition as compared to a one-phase particle. We present a new liquid-liquid equilibrium and gas/particle partitioning model, using as a basis the group-contribution model AIOMFAC (Zuend et al., 2008). This model allows the reliable computation of the liquid-liquid coexistence curve (binodal), corresponding tie-lines, the limit of stability/metastability (spinodal), and further thermodynamic properties of multicomponent systems. Calculations for ternary and multicomponent alcohol/polyol-water-salt mixtures suggest that LLE are a prevalent feature of organic-inorganic aerosol systems. A six-component polyol-water-ammonium sulphate system is used to simulate effects of relative humidity (RH) and the presence of liquid-liquid phase separation on the gas/particle partitioning. RH, salt concentration, and hydrophilicity (water-solubility) are identified as key features in defining the region of a miscibility gap and govern the extent to which compound partitioning is affected by changes in RH. The model predicts that liquid-liquid phase separation can lead to either an increase or decrease in total particulate mass, depending on the overall composition of a system and the particle water content, which is related to the hydrophilicity of the different organic and inorganic compounds. Neglecting non-ideality and liquid-liquid phase separations by assuming an ideal mixture leads to an overestimation of the total particulate mass by up to 30% for the composition and RH range considered in the six-component system simulation. For simplified partitioning parametrizations, we suggest a modified definition of the effective saturation concentration, C_j^*, by including water and other inorganics in the absorbing phase. Such a C_j^* definition reduces the RH-dependency of the gas/particle partitioning of semivolatile organics in organic-inorganic aerosols by an order of magnitude as compared to the currently accepted definition, which considers the organic species only

Prediction of Partition Coefficients in SDS Micelles by DFT Calculations

A computational methodology using Density-Functional Theory (DFT) calculations was developed to determine the partition coefficient of a compound in a solution of Sodium Dodecyl Sulfate (SDS) micelles. Different sets of DFT calculations were used to predict the free energy of a set of 63 molecules in 15 different solvents with the purpose of identifying the solvents with similar physicochemical characteristics to the studied micelles. Experimental partition coefficients were obtained from Micellar Electrokinetic Chromatography (MEKC) measurements. The experimental partition coefficient of these molecules was compared with the predicted partition coefficient in heptane/water, cyclohexane/water, N-dodecane/water, pyridine/water, acetic acid/water, decan-1- ol/water, octanol/water, propan-2-ol/water, acetone/water, propan-1-ol/water, methanol/water, 1,2-ethane diol/water, dimethyl sulfoxide/water, formic acid/water, and diethyl sulphide/water systems. It is observed that the combination of pronan-1-ol/water solvent was the most appropriated to estimate the partition coefficient for SDS micelles. This approach allowed us to estimate the partition coefficient orders of magnitude faster than the classical molecular dynamics simulations. The DFT calculations were carried out with the well-known exchange correlation functional B3LYP and with the global hybrid functional M06-2X from the Minnesota functionals with 6-31++G ** basis set. The effect of solvation was considered by the continuum model based on density (SMD). The proposed workflow for the prediction rate of the participation coefficient unveiled the symmetric balance between the experimental data and the computational methods

Protein partitioning in polymer/polymer aqueous two-phase systems

Tese de doutoramento. Engenharia Química e Biológica. Faculdade de Engenharia. Universidade do Porto. 200

Crossover Scales at the Critical Points of Fluids with Electrostatic Interactions

Criticality in a fluid of dielectric constant D that exhibits Ising-type behavior is studied as additional electrostatic (i.e., ionic) interactions are turned on. An exploratory perturbative calculation is performed for small ionicity as measured by the ratio of the electrostatic energy to the strength of the short-range nonionic (i.e., van der Waals) interactions in the uncharged fluid. With the aid of distinct transformations for the short-range and for the Coulombic interactions, an effective Hamiltonian with coefficients depending on the ionicity is derived at the Debye-Hueckel limiting-law level for a fully symmetric model. The crossover between classical (mean-field) and Ising behavior is then estimated using a Ginzburg criterion. This indicates that the reduced crossover temperature depends only weakly on the ionicity (and on the range of the nonionic potentials); however, the trends do correlate with the, much stronger, dependence observed experimentally.Comment: 25 pages, 4 figure; submitted to J. Chem. Phy

Counter-current chromatography for the separation of terpenoids: A comprehensive review with respect to the solvent systems employed

Copyright @ 2014 The Authors.This article is distributed under the terms of the Creative Commons Attribution License which permits any use, distribution, and reproduction in any medium, provided the original author(s) and the source are credited.Natural products extracts are commonly highly complex mixtures of active compounds and consequently their purification becomes a particularly challenging task. The development of a purification protocol to extract a single active component from the many hundreds that are often present in the mixture is something that can take months or even years to achieve, thus it is important for the natural product chemist to have, at their disposal, a broad range of diverse purification techniques. Counter-current chromatography (CCC) is one such separation technique utilising two immiscible phases, one as the stationary phase (retained in a spinning coil by centrifugal forces) and the second as the mobile phase. The method benefits from a number of advantages when compared with the more traditional liquid-solid separation methods, such as no irreversible adsorption, total recovery of the injected sample, minimal tailing of peaks, low risk of sample denaturation, the ability to accept particulates, and a low solvent consumption. The selection of an appropriate two-phase solvent system is critical to the running of CCC since this is both the mobile and the stationary phase of the system. However, this is also by far the most time consuming aspect of the technique and the one that most inhibits its general take-up. In recent years, numerous natural product purifications have been published using CCC from almost every country across the globe. Many of these papers are devoted to terpenoids-one of the most diverse groups. Naturally occurring terpenoids provide opportunities to discover new drugs but many of them are available at very low levels in nature and a huge number of them still remain unexplored. The collective knowledge on performing successful CCC separations of terpenoids has been gathered and reviewed by the authors, in order to create a comprehensive document that will be of great assistance in performing future purifications. © 2014 The Author(s)

Fast Computation of Solvation Free Energies with Molecular Density Functional Theory: Thermodynamic-Ensemble Partial Molar Volume Corrections

Molecular Density Functional Theory (MDFT) offers an efficient implicit- solvent method to estimate molecule solvation free-energies whereas conserving a fully molecular representation of the solvent. Even within a second order ap- proximation for the free-energy functional, the so-called homogeneous reference uid approximation, we show that the hydration free-energies computed for a dataset of 500 organic compounds are of similar quality as those obtained from molecular dynamics free-energy perturbation simulations, with a computer cost reduced by two to three orders of magnitude. This requires to introduce the proper partial volume correction to transform the results from the grand canoni- cal to the isobaric-isotherm ensemble that is pertinent to experiments. We show that this correction can be extended to 3D-RISM calculations, giving a sound theoretical justifcation to empirical partial molar volume corrections that have been proposed recently.Comment: Version with correct equation numbers is here: http://compchemmpi.wikispaces.com/file/view/sergiievskyi_et_al.pdf/513575848/sergiievskyi_et_al.pdf Supporting information available online at: http://compchemmpi.wikispaces.com/file/view/SuppInf_sergiievskyi_et_al_07-04-2014.pdf/513576008/SuppInf_sergiievskyi_et_al_07-04-2014.pd