Environmental effects on the spectroscopic properties of gallic acid: a combined classical and quantum mechanical study

The solvation of gallic acid (in water and acetonitrile) is studied by means of its spectroscopic properties. IR, UV, and NMR spectra are predicted by using various solvation models obtained in terms of both purely classical and density functional approaches. Comparison with experiments is used to validate solvation models. Hydrogen-bond and long-range (or bulk) effects are evaluated by comparing different solvation models. A continuum-only approach, a purely discrete, and a mixed continuum/discrete approach based on quantum-mechanical and classical molecular-dynamics solute-solvent clusters are tested.     

High-throughput prediction of blood-brain partitioning: a thermodynamic approach

A high-throughput in silico screening tool for potentially CNS active compounds was developed on the basis of the correlation of solvation free energies and blood-brain partitioning (log(cbrain/cblood) = log BB) data available from experimental sources. Utilizing a thermodynamic approach, solvation free energies were calculated by the fast and efficient generalized Born/surface area continuum solvation model, which enabled us to evaluate more than 10 compounds/min. Our training set involved a structurally diverse set of 55 compounds and yielded a function of log BB = 0.035Gsolv + 0.2592 (r = 0.85, standard error 0.37). Calculation of solvation free energies for 8700 CNS active compounds (CIPSLINE database) revealed that Gsolv is higher than -50 kJ/mol for the 96% of these compounds which can be used as suitable criteria for the identification of compounds preferable for CNS penetration.     

A novel quantum mechanical/molecular mechanical approach to the free energy calculation for isomerization of glycine in aqueous solution

The free energy change associated with the isomerization reaction of glycine in water solution has been studied by a hybrid quantum mechanical/molecular mechanical (QM/MM) approach combined with the theory of energy representation (QM/MM-ER) recently developed. The solvation free energies for both neutral and zwitterionic form of glycine have been determined by means of the QM/MM-ER simulation. The contributions of the electronic polarization and the fluctuation of the QM solute to the solvation free energy have been investigated. It has been found that the contribution of the density fluctuation of the zwitterionic solute is estimated as -4.2 kcal/mol in the total solvation free energy of -46.1 kcal/mol, while that of the neutral form is computed as -3.0 kcal/mol in the solvation free energy of -15.6 kcal/mol. The resultant free energy change associated with the isomerization of glycine in water has been obtained as -7.8 kcal/mol, in excellent agreement with the experimental data of -7.3 or -7.7 kcal/mol, implying the accuracy of the QM/MM-ER approach. The results have also been compared with those computed by other methodologies such as the polarizable continuum model and the classical molecular simulation. The efficiency and advantage of the QM/MM-ER method has been discussed.     

Calculation of solvent shifts on electronic g-tensors with the conductor-like screening model (COSMO) and its self-consistent generalization to real solvents (direct COSMO-RS)

The conductor-like screening model (COSMO) was used to investigate the solvent influence on electronic g-values of organic radicals. The previously studied diphenyl nitric oxide and di-tert-butyl nitric oxide radicals were taken as test cases. The calculations employed spin-unrestricted density functional theory and the BP and B3LYP density functionals. The g-tensors were calculated as mixed second derivative properties with respect to the external magnetic field and the electron magnetic moment. The first-order response of the Kohn-Sham orbitals with respect to the external magnetic field was determined through the coupled-perturbed DFT approach. The spin-orbit coupling operator was treated using an accurate multicenter spin-orbit mean-field (SOMF) approach. Provided that important hydrogen bonds are explicitly modeled by a supermolecule approach and that the basis set is sufficiently saturated, the COSMO calculations lead to accurate predictions of isotropic g-shifts with deviations of not more than 100 ppm relative to experiment. Very accurate results were obtained by employing a recently developed self-consistent modification of the COSMO method to real solvents (COSMO-RS), which we briefly introduce in this paper as direct COSMO-RS (D-COSMO-RS). This model gives isotropic g-shifts of similar high accuracy for water without using the supermolecule approach. This is an important result because it solves many of the problems associated with the supermolecule approach such as local minima and the choice of a suitable model system. Thus, the self-consistent D-COSMO-RS incorporates some specific solvation effects into continuum models, in particular it appears to successfully model the effects of hydrogen bonding. Although not yet widely validated, this opens a novel approach for the calculation of properties which so far only could be calculated by the inclusion of explicit solvent molecules in continuum solvation methods.     

On the solvation of the Zn2+ ion in methanol: a combined quantum mechanics, molecular dynamics, and EXAFS approach

The solvation properties of the Zn(2+) ion in methanol solution have been investigated using a combined approach based on molecular dynamics (MD) simulations and extended X-ray absorption fine structure (EXAFS) experimental results. The quantum mechanical potential energy surface for the interaction of the Zn(2+) ion with a methanol molecule has been calculated taking into account the effect of bulk solvent by the polarizable continuum model (PCM). The effective Zn-methanol interactions have been fitted by suitable analytical potentials, and have been utilized in the MD simulation to obtain the structural properties of the solution. The reliability of the whole procedure has been assessed by comparing the theoretical structural results with the EXAFS experimental data. The structural parameters of the first solvation shells issuing from the MD simulations provide an effective complement to the EXAFS experiments.     

Peculiarities of the environmental influence on the optical properties of push-pull nonlinear optical molecules: a theoretical study

Gas-phase geometry optimization of NLO-active molecules is one of the standard approaches in the first principle computational methodology, whereas the important role of the environment is usually not considered during the evaluation of structural parameters. With a wide variety of environmentally influenced models in most cases only the high quality single point calculations are prepared. Among different approaches, the most used polarizable continuum model (PCM) seems to be promising. In this study, we have compared the electronic properties of gas-phase optimized geometries of imidazole-derived push-pull compounds with those optimized using PCM solvation approach including CH(2)Cl(2) and PMMA as media. We have focused particularly on the linear optical properties of investigated molecules, namely on the UV-vis absorption spectra. The analysis of presented results shows the applicability of the different quantum chemical (QC) methods for the UV-vis spectra calculations of linear NLO molecules. Herein we also present the need of molecule geometry optimization affected by the environment. Following the performed calculations, the electronic properties of gas-phase optimized molecules give conformable results with respect to those obtained by more time-consuming continuum optimizations. All computational data are supported by experimental investigations.     

Rapid prediction of solvation free energy. 3. Application to the SAMPL2 challenge

The SAMPL2 hydration free energy blind prediction challenge consisted of a data set of 41 molecules divided into three subsets: explanatory, obscure and investigatory, where experimental hydration free energies were given for the explanatory, withheld for the obscure, and not known for the investigatory molecules. We employed two solvation models for this challenge, a linear interaction energy (LIE) model based on explicit-water molecular dynamics simulations, and the first-shell hydration (FiSH) continuum model previously calibrated to mimic LIE data. On the 23 compounds from the obscure (blind) dataset, the prospectively submitted LIE and FiSH models provided predictions highly correlated with experimental hydration free energy data, with mean-unsigned-errors of 1.69 and 1.71 kcal/mol, respectively. We investigated several parameters that may affect the performance of these models, namely, the solute flexibility for the LIE explicit-solvent model, the solute partial charging method, and the incorporation of the difference in intramolecular energy between gas and solution phases for both models. We extended this analysis to the various chemical classes that can be formed within the SAMPL2 dataset. Our results strengthen previous findings on the excellent accuracy and transferability of the LIE explicit-solvent approach to predict transfer free energies across a wide spectrum of functional classes. Further, the current results on the SAMPL2 test dataset provide additional support for the FiSH continuum model as a fast yet accurate alternative to the LIE explicit-solvent model. Overall, both the LIE explicit-solvent model and the FiSH continuum solvation model show considerable improvement on the SAMPL2 data set over our previous continuum electrostatics-dispersion solvation model used in the SAMPL1 blind challenge.     

Solvation Dynamics of a Radical Ion Pair in Micro‐Heterogeneous Binary Solvents: A Semi‐Quantitative Study Utilizing MARY Line‐Broadening Experiments

This work aims at elucidating the mechanism of solvation of a radical ion pair (RIP) in a micro‐heterogeneous binary solvent mixture using magnetically affected reaction yield (MARY) spectroscopy. For the exciplex‐forming 9,10‐dimethylanthracene/N,N‐dimethylaniline system a comparative, composition‐dependent MARY line‐broadening study is undertaken in a heterogeneous (toluene/dimethylsulfoxide) and a quasi‐homogenous (propyl acetate/butyronitrile) solvent mixture. The half‐saturation field extrapolated to zero‐quencher concentration, B_1/2, and the self‐exchange rate constants are analyzed in the light of solvent dynamical properties of the mixtures and a dielectric continuum solvation model. The dependence of B_1/2 on the solvent composition is explained by cluster formation giving rise to shortened RIP lifetimes. The results are in qualitative agreement with the continuum solvation model suggesting that it could serve as a theoretical basis for quantitative modeling.     

Proton-donating power of 100% perchloric acid: An MP2 study

The ab initio quantum chemical method MP2 with a 6-311++G(d, p) basis set is used to calculate the energy of gas-phase solvation of the H₂Cl₄⁺ cation and Cl₄⁻ anion by one perchloric acid molecule. The energy of additional solvation of the resulting complexes by liquid perchloric acid is estimated within the continuum model of solvation by the PCM method, with the acid modeled as a continuum with a large dielectric constant of ɛ = 115. The calculated data have provided an almost quantitative estimate for the energy of selfionization of 100% liquid perchloric acid. The similarly calculated energy of solvation of protons by 100% perchloric acid is 30.7 kcal/mol lower than the heat of hydration of protons in aqueous solution. This result explains the fact that anhydrous perchloric acid exhibits the properties of a superacid.     

Smooth solvation method for d-orbital semiempirical calculations of biological reactions. 2. Application to transphosphorylation thio effects in solution

Density-functional and semiempirical quantum methods and continuum dielectric and explicit solvation models are applied to study the role of solvation on the stabilization of native and thio-substituted transphosphorylation reactions. Extensive comparison is made between results obtained from the different methods. For the semiempirical methods, explicit solvation was treated using a hybrid quantum mechanical/molecular mechanical (QM/MM) approach and the implicit solvation was treated using a recently developed smooth solvation model implemented into a d-orbital semiempirical framework (MNDO/d-SCOSMO) within CHARMM. The different quantum and solvation methods were applied to the transesterification of a 3'-ribose,5'-methyl phosphodiester that serves as a nonenzymatic model for the self-cleavage reaction catalyzed by the hammerhead and hairpin ribozymes. Thio effects were studied for a double sulfur substitution at the nonbridging phosphoryl oxygen positions. The reaction profiles of both the native and double sulfur-substituted reactions from the MNDO/d-SCOSMO calculations were similar to the QM/MM results and consistent with the experimentally observed trends. These results underscore the need for a d-orbital semiempirical representation for phosphorus and sulfur for the study of experimentally observed thio effects in enzymatic and nonenzymatic phosphoryl transfer reactions. One of the major advantages of the present approach is that it can be applied to model chemical reactions at a significantly lower computational cost than either the density-functional calculations with implicit solvation or the semiempirical QM/MM simulations with explicit solvent.     

A new quantum method for electrostatic solvation energy of protein

A new method that incorporates the conductorlike polarizable continuum model (CPCM) with the recently developed molecular fractionation with conjugate caps (MFCC) approach is developed for ab initio calculation of electrostatic solvation energy of protein. The application of the MFCC method makes it practical to apply CPCM to calculate electrostatic solvation energy of protein or other macromolecules in solution. In this MFCC-CPCM method, calculation of protein solvation is divided into calculations of individual solvation energies of fragments (residues) embedded in a common cavity defined with respect to the entire protein. Besides computational efficiency, the current approach also provides additional information about contribution to protein solvation from specific fragments. Numerical studies are carried out to calculate solvation energies for a variety of peptides including alpha helices and beta sheets. Excellent agreement between the MFCC-CPCM result and those from the standard full system CPCM calculation is obtained. Finally, the MFCC-CPCM calculation is applied to several real proteins and the results are compared to classical molecular mechanics Poisson-Boltzmann (MM/PB) and quantum Divid-and-Conque Poisson-Boltzmann (D&C-PB) calculations. Large wave function distortion energy (solute polarization energy) is obtained from the quantum calculation which is missing in the classical calculation. The present study demonstrates that the MFCC-CPCM method is readily applicable to studying solvation of proteins.     

New approach to free energy of solvation applying continuum models to molecular dynamics simulation

A new approach to the calculation of the free energy of solvation from trajectories obtained by molecular dynamics simulation is presented. The free energy of solvation is computed as the sum of three contributions originated at the cavitation of the solute by the solvent, the solute-solvent nonpolar (repulsion and dispersion) interactions, and the electrostatic solvation of the solute. The electrostatic term is calculated based on ideas developed for the broadly used continuum models, the cavitational contribution from the excluded volume by the Claverie-Pierotti model, and the Van der Waals term directly from the molecular dynamics simulation. The proposed model is tested for diluted aqueous solutions of simple molecules containing a variety of chemically important functions: methanol, methylamine, water, methanethiol, and dichloromethane. These solutions were treated by molecular dynamics simulations using SPC/E water and the OPLS force field for the organic molecules. Obtained free energies of solvation are in very good agreement with experimental data.     

Microsolvated transition state models for improved insight into chemical properties and reaction mechanisms

Over the years, several methods have been developed to effectively represent the chemical behavior of solutes in solvents. The environmental effects arising due to solvation can generally be achieved either through inclusion of discrete solvent molecules or by inscribing into a cavity in a homogeneous and continuum dielectric medium. In both these approaches of computational origin, the perturbations on the solute induced by the surrounding solvent are at the focus of the problem. While the rigor and method of inclusion of solvent effects vary, such solvation models have found widespread applications, as evident from modern chemical literature. A hybrid method, commonly referred to as cluster-continuum model (CCM), brings together the key advantages of discrete and continuum models. In this perspective, we intend to highlight the latent potential of CCM toward obtaining accurate estimates on a number of properties as well as reactions of contemporary significance. The objective has generally been achieved by choosing illustrative examples from the literature, besides expending efforts to bring out the complementary advantages of CCM as compared to continuum or discrete solvation models. The majority of examples emanate from the prevalent applications of CCM to organic reactions, although a handful of interesting organometallic reactions have also been discussed. In addition, increasingly accurate computations of properties like pK(a) and solvation of ions obtained using the CCM protocol are also presented.     

Dispersion terms and analysis of size- and charge dependence in an enhanced Poisson-Boltzmann approach

We implement a well-established concept to consider dispersion effects within a Poisson-Boltzmann approach of continuum solvation of proteins. The theoretical framework is particularly suited for boundary element methods. Free parameters are determined by comparison to experimental data as well as high-level quantum mechanical reference calculations. The method is general and can be easily extended in several directions. The model is tested on various chemical substances and found to yield good-quality estimates of the solvation free energy without obvious indication of any introduced bias. Once optimized, the model is applied to a series of proteins, and factors such as protein size or partial charge assignments are studied.     

Elucidation of the thermochemical properties of triphenyl- or tributyl-substituted Si-, Ge-, and Sn-centered radicals by means of electrochemical approaches and computations

Redox potentials of a number of triphenyl- or tributyl-substituted Si-, Ge-, or Sn-centered radicals, R(3)M(*), have been measured in acetonitrile, tetrahydrofuran, or dimethyl sulfoxide by photomodulated voltammetry or through a study of the oxidation process of the corresponding anions in linear sweep voltammetry. For the results pertaining to the Ph(3)M(*) series (including literature data for M = C), the order of reduction potentials follows Sn > Ge > C > Si, while for the two oxidation potentials, it is C > Si. The effect of the R group on the redox properties of R(3)Sn(*) is pronounced in that the reduction potential is more negative by 490 mV in tetrahydrofuran (390 mV in dimethyl sulfoxide) when R is a butyl rather than a phenyl group. The experimental trends have been substantiated through quantum chemical calculations, and they can be explained qualitatively by considering a combination of effects, such as charge capacity being most pronounced for the heavier elements, resonance stabilization present for the planar Ph(3)C(*) and all R(3)M(+)(), and finally a contribution from solvation. The solvation of R(3)M(-) is observed to be relatively strong because of a rather localized negative charge in the pyramidal geometry. However, there is no evidence in the calculations to support the existence of covalent interactions between solvent and anions. The solvation of R(3)M(+)() is relatively weak, which may be attributed to the planar geometry around the center atom, leading to more spread out charge than that for a pyramidal geometry. Although the calculated solvation energies based on the polarizable continuum model approach exhibit the expected trends, they are not able to reproduce the experimentally derived values on a detailed level for these types of ions. An evaluation of the general performance of the continuum model is provided on the basis of present and previous studies.     

基于分子动力学模拟和连续介质模型的自由能计算方法*

近些年,基于分子动力学模拟和连续介质模型的自由能计算方法受到了越来越多的关注,其中MM/PBSA就是最具代表性的方法.在MM/PBSA中,体系的焓变采用分子力学(MM)的方法计算得到;溶剂效应中极性部分对自由能的贡献通过解Poisson-Boltzmann(PB)方程的方法计算得到;溶液效应中非极性部分对自由能的贡献则通过分子表面积(SA)计算得到.本文结合我们科研组的工作,就近几年MM/PBSA方法的最新进展做了较为详细的阐述,同时对MM/PBSA的发展前景进行了展望.     

Thermochemistry of arylselanyl radicals and the pertinent ions in acetonitrile

Reduction and oxidation potentials of a series of parasubstituted phenylselanyl radicals, XC(6)H(4)Se(*), have been measured using photomodulated voltammetry in acetonitrile. The thermodynamic significance of these data was substantiated through a study of the oxidation process of the pertinent selenolates in linear sweep voltammetry. Both the reduction and the oxidation potentials correlate linearly with the Hammett substituent coefficients sigma and sigma(+) leading in the latter case to slopes, rho(+), of 2.5 and 3.8, respectively. Through comparison of these slopes with those published previously for the O- and S-centered analogues, it is revealed that the pi-interaction becomes progressively smaller as the size of the radical center increases in the order O, S, and Se. Solvation energies of the pertinent selenolates and selanylium ions have been extracted from thermochemical cycles incorporating the measured electrode potentials for XC(6)H(4)Se(*) as well as electron affinities and ionization potentials obtained from theoretical calculations at the B3LYP/6-31+G(d) level. The extracted data show the expected overall substituent dependency for both kinds of ions; that is, the absolute value of the solvation energy decreases as the charge becomes more delocalized. The data have also been compared with solvation energies computed using the polarizable continuum model (PCM). Interestingly, we find that, while the model seems to work well for selenolates, it underestimates the solvation of selanylium ions in acetonitrile by as much as 25 kcal mol(-)(1). These large deviations are ascribed to the fact that the PCM method does not take specific solvent effects into account as it treats the solvent as a continuum described solely by its dielectric constant. Gas-phase calculations show that the arylselanylium ions can coordinate covalently to one or two molecules of acetonitrile in strong Ritter-type adducts. When this strong interaction is included in the solvation energy calculations by means of a combined supermolecule and PCM approach, the experimental data are reproduced within a few kcal mol(-)(1). Although the energy difference of the singlet and triplet spin states of the arylselanylium ions is small for the gas-phase structures, the singlet cation is undoubtedly the dominating species in solution because the triplet cation lacks the ability to form covalent bonds.     

Solubility of naproxen in several organic solvents at different temperatures

By using the van’t Hoff and Gibbs equations the thermodynamic functions Gibbs free energy, enthalpy, and entropy of solution, were evaluated from solubility data of naproxen (NAP) determined at several temperatures in octanol, isopropyl myristate, chloroform, and cyclohexane, as pure solvents. The water-saturated organic solvents also were studied except cyclohexane. The excess free energy and the activity coefficients of the solutes, and the mixing and solvation thermodynamic quantities were also determined. The NAP solubilities were higher in chloroform and octanol with respect to those obtained in cyclohexane. In addition, by using literature values for NAP aqueous solubility, the thermodynamic functions relative to transfer of this drug from water to organic solvents were also estimated.