Electron solvation by highly polar molecules: density functional theory study of atomic sodium interaction with water, ammonia, and methanol

This study further extends the scope of a previous paper [Y. Ferro and A. Allouche, J. Chem. Phys. 118, 10461 (2003)] on the reactivity of atomic Na with water to some other highly polar molecules known for their solvation properties connected to efficient hydrogen bonding. The solvation mechanisms of ammonia and methanol are compared to the hydration mechanism. It is shown that in the case of ammonia, the stability of the solvated system is only ensured by electrostatic interactions, whereas the methanol action is more similar to that of water. More specific attention is given to the solvation process of the valence 3s Na electron. The consequences on the chemical reactivity are analyzed: Whereas ammonia is nonreactive when interacting with atomic sodium, two chemical reactions are proposed for methanol. The first process is dehydrogenation and yields methoxy species and hydrogen. The other one is dehydration and the final products are methoxy species, but also methyl radical and water. The respective roles of electron solvation and hydrogen bonds network are analyzed in detail in view of the density of states of the reactive systems.     

Preferential Solvation of Some Sulfonamides in Propylene Glycol + Water Solvent Mixtures According to the IKBI and QLQC Methods

The preferential solvation parameters, which represent differences between the local and bulk mole fractions of the solvents near to the solute, in solutions of some sulfonamides in propylene glycol + water binary mixtures are derived from their thermodynamic properties by means of the inverse Kirkwood?Buff integrals (IKBI) and the Quasi-Lattice Quasi-Chemical (QLQC) method. From solvent effect studies, it is found that sulfonamides are sensitive to solvation effects; the preferential solvation parameter, δx _PG,S, is negative in water-rich mixtures but positive in compositions from 0.20 to 1.00 in mole fraction of propylene glycol according to IKBI method and positive in all co-solvent compositions if the QLQC method is considered. It is conjecturable that in water-rich mixtures, hydrophobic hydration around the aromatic ring and/or other non-polar groups plays a relevant role in the solvation. The greater solvation by propylene glycol mixtures of similar solvent compositions and in co-solvent-rich mixtures could be due mainly to polarity effects and acidic behavior of the sulfonamides, in contrast to the more basic solvent propylene glycol.     

Communication: Free-energy analysis of hydration effect on protein with explicit solvent: equilibrium fluctuation of cytochrome c

The relationship between the protein conformation and the hydration effect is investigated for the equilibrium fluctuation of cytochrome c. To elucidate the hydration effect with explicit solvent, the solvation free energy of the protein immersed in water was calculated using the molecular dynamics simulation coupled with the method of energy representation. The variations of the protein intramolecular energy and the solvation free energy are found to compensate each other in the course of equilibrium structural fluctuation. The roles of the attractive and repulsive components in the protein-water interaction are further examined for the solvation free energy. The attractive component represented as the average sum of protein-water interaction energy is dominated by the electrostatic effect and is correlated to the solvation free energy through the linear-response-type relationship. No correlation with the (total) solvation free energy is seen, on the other hand, for the repulsive component expressed as the excluded-volume effect.     

PM3-SM3: A general parameterization for including aqueous solvation effects in the PM3 molecular orbital model

Our recently proposed scheme for including aqueous solvation free energies in parameterized NDDO SCF models is extended to the Parameterized Model 3 semiempirical Hamiltonian. The solvation model takes accurate account of the hydrophobic effect for hydrocarbons, as well as electric polarization of the solvent, the free energy of cavitation, and dispersion interactions. Eight heteroatoms are included (along with H and C), and the new model is parameterized accurately for the water molecule itself, which allows meaningful treatments of specifically hydrogen bonded water molecules. The unphysical partial charges on nitrogen atoms predicted by the Parameterized Model 3 Hamiltonian limit the accuracy of the predicted solvation energies for some compounds containing nitrogen, but the model may be very useful for other systems, especially those for which PM3 is preferred over AM1 for the solute properties of the particular system under study. © 1992 by John Wiley & Sons, Inc.     

Sensitivity of polar solvation dynamics to the secondary structures of aqueous proteins and the role of surface exposure of the probe

The structure and dynamics of water around a protein is expected to be sensitive to the details of the adjacent secondary structure of the protein. In this article, we explore this sensitivity by calculating both the orientational dynamics of the surface water molecules and the equilibrium solvation time correlation function of the polar amino acid residues in each of the three helical segments of the protein HP-36, using atomistic molecular dynamics simulations. The solvation dynamics of polar amino acid residues in helix-2 is found to be faster than that of the other two helices (the average time constant is smaller by a factor of 2), although the interfacial water molecules around helix-2 exhibit much slower orientational dynamics than that around the other two helices. A careful analysis shows that the origin of such a counterintuitive behavior lies in the dependence of the solvation time correlation function on the surface exposure of the probe-the more exposed is the probe, the faster the solvation dynamics. We discuss that these results are useful in explaining recent solvation dynamics experiments.     

A method for calculating the Gibbs energy of nonspecific solvation

A method for calculating the Gibbs energy of nonspecific solvation of nonelectrolytes was suggested. The new equation for the Gibbs energy of nonspecific solvation contains one solvent parameter that characterize nonspecific solvent-solute interactions and two experimental Gibbs energies of solvation in two standard solvents. The method is applicable to a wide range of solutes and solvents. It was successfully used to describe some 800 Gibbs energies of solvation for systems without specific solvent-solute interactions.     

The impact of surface area,volume, curvature,and Lennard–Jones potential to solvation modeling

This article explores the impact of surface area, volume, curvature, and Lennard–Jones (LJ) potential on solvation free energy predictions. Rigidity surfaces are utilized to generate robust analytical expressions for maximum, minimum, mean, and Gaussian curvatures of solvent–solute interfaces, and define a generalized Poisson–Boltzmann (GPB) equation with a smooth dielectric profile. Extensive correlation analysis is performed to examine the linear dependence of surface area, surface enclosed volume, maximum curvature, minimum curvature, mean curvature, and Gaussian curvature for solvation modeling. It is found that surface area and surfaces enclosed volumes are highly correlated to each other's, and poorly correlated to various curvatures for six test sets of molecules. Different curvatures are weakly correlated to each other for six test sets of molecules, but are strongly correlated to each other within each test set of molecules. Based on correlation analysis, we construct twenty six nontrivial nonpolar solvation models. Our numerical results reveal that the LJ potential plays a vital role in nonpolar solvation modeling, especially for molecules involving strong van der Waals interactions. It is found that curvatures are at least as important as surface area or surface enclosed volume in nonpolar solvation modeling. In conjugation with the GPB model, various curvature‐based nonpolar solvation models are shown to offer some of the best solvation free energy predictions for a wide range of test sets. For example, root mean square errors from a model constituting surface area, volume, mean curvature, and LJ potential are less than 0.42 kcal/mol for all test sets. © 2016 Wiley Periodicals, Inc.     

Solvation model based on order parameters and a fast sampling method for the calculation of the solvation free energies of peptides

An analytical solvation model is proposed as a function of an order parameter, which represents the local arrangement of water molecules in the first solvation shell of peptide atoms. The model is combined with a fast sampling method, rotational isomeric state Monte Carlo, to sample efficiently the torsional degrees of freedom on a peptide backbone. This order parameter solvation model is shown to reproduce without ad hoc fitting parameters the solvation free energies of single amino acids and tripeptides with slightly better accuracy than the generalized Born model but with several orders of magnitude improvement in efficiency. This method is a potential candidate for efficiently and accurately tackling some important issues in biophysical chemistry that are related to solvation, for example, protein folding, ligand binding, etc. Our results also present fundamental new insights into solvation. Specifically, the local water geometry, represented in this work by a properly defined order parameter, carries the majority, if not all, of the energetic information of solvation, including solute-solvent interactions and solvent reorganization in the presence of the solute.     

Characterization of metal ions in coordinating solvent mixtures by means of Raman spectroscopy.

Titration Raman spectroscopy has been developed for studying the solvation structure of metal ions in solution. The method affords us the solvation number, and the value thus obtained in neat solvents is in good agreement with that determined by EXAFS. The method is then applied to solvent mixtures, and the individual solvation number for each solvent is extracted. In a solvent mixture of N,N-dimethylformamide (DMF) and N,N,N',N'-tetramethylurea (TMU), the metal ion prefers DMF to TMU, which is ascribed to the solvation steric effect. The same applies also for the solvent mixture of N,N-dimethylpropionamide (DMPA) and DMF. However, unlike TMU, DMPA changes its conformation from the planar cis to non-planar staggered upon solvation to the metal ion. The enthalpy for the conformational change of DMPA is positive in the bulk, while it is significantly negative in the coordination sphere of the manganese(II) ion. Here, we briefly describe the procedure of measurements and analyses for the titration Raman spectroscopy, and review the solvation structure of the alkaline earth, first transition metal(II) and lanthanide(III) ions in some solvent mixtures in view of solvation steric effect.     

Adiabatic compressibility and solvation of drag reducing polymers in aqueous solutions

Adiabatic compressibilities of aqueous solutions of some drag reducing polymers have been evaluated from ultrasonic velocity and density measurements. The solvation numbers of the repeat units of the polymers have been further evaluated by using Passynsky equation. It is observed that the solvation number increases with the shear stability of drag reducing polymers in turbulent flow.     

One-dimensional model for water and aqueous solutions. II. Solvation of inert solutes in water

The two one-dimensional models introduced in Part I are used to study the thermodynamics of solvation of inert solutes in water. It is shown that the anomalously large Gibbs energy of solvation of inert solutes in water, on one hand, and the large negative entropy of solvation, on the other hand, arise from different molecular sources. While the primitive model can give rise to a large positive solvation Gibbs energy, it fails to show large negative entropy and enthalpy of solvation. It is remarkable that the primitive cluster model can show both the large positive Gibbs energy and enthalpy of solvation.     

Differential geometry based solvation model. III. Quantum formulation

Solvation is of fundamental importance to biomolecular systems. Implicit solvent models, particularly those based on the Poisson-Boltzmann equation for electrostatic analysis, are established approaches for solvation analysis. However, ad hoc solvent-solute interfaces are commonly used in the implicit solvent theory. Recently, we have introduced differential geometry based solvation models which allow the solvent-solute interface to be determined by the variation of a total free energy functional. Atomic fixed partial charges (point charges) are used in our earlier models, which depends on existing molecular mechanical force field software packages for partial charge assignments. As most force field models are parameterized for a certain class of molecules or materials, the use of partial charges limits the accuracy and applicability of our earlier models. Moreover, fixed partial charges do not account for the charge rearrangement during the solvation process. The present work proposes a differential geometry based multiscale solvation model which makes use of the electron density computed directly from the quantum mechanical principle. To this end, we construct a new multiscale total energy functional which consists of not only polar and nonpolar solvation contributions, but also the electronic kinetic and potential energies. By using the Euler-Lagrange variation, we derive a system of three coupled governing equations, i.e., the generalized Poisson-Boltzmann equation for the electrostatic potential, the generalized Laplace-Beltrami equation for the solvent-solute boundary, and the Kohn-Sham equations for the electronic structure. We develop an iterative procedure to solve three coupled equations and to minimize the solvation free energy. The present multiscale model is numerically validated for its stability, consistency and accuracy, and is applied to a few sets of molecules, including a case which is difficult for existing solvation models. Comparison is made to many other classic and quantum models. By using experimental data, we show that the present quantum formulation of our differential geometry based multiscale solvation model improves the prediction of our earlier models, and outperforms some explicit solvation model.     

Ion solvation thermodynamics from simulation with a polarizable force field

Thermodynamic measurements of the solvation of salts and electrolytes are relatively straightforward, but it is not possible to separate total solvation free energies into distinct cation and anion contributions without reference to an additional extrathermodynamic assumption. The present work attempts to resolve this difficulty using molecular dynamics simulations with the AMOEBA polarizable force field and perturbation techniques to directly compute absolute solvation free energies for potassium, sodium, and chloride ions in liquid water and formamide. Corresponding calculations are also performed with two widely used nonpolarizable force fields. The simulations with the polarizable force field accurately reproduce in vacuo quantum mechanical results, experimental ion-cluster solvation enthalpies, and experimental solvation free energies for whole salts, while the other force fields do not. The results indicate that calculations with a polarizable force field can capture the thermodynamics of ion solvation and that the solvation free energies of the individual ions differ by several kilocalories from commonly cited values.     

Site-specific solvation determined by intermolecular nuclear Overhauser effect--measurements and molecular dynamics

Site-specific solvation has been determined by intermolecular NOE measurements between solvent and solute. The experimental effect is shown on the four compounds 2-butanol, L-alanyl-L-tryptophan (Ala-Trp), adenosine and the disodium salt of adenosine 5'-monophosphate (5'-AMP) in the two solvents water and dimethyl sulfoxide (DMSO). The strength of NOE transfer correlates with the average distribution of solvent molecules around the corresponding solvation sites represented by the number of solvent molecules in a first solvation sphere, which can be obtained from molecular dynamics simulations in water. Saturation transfer between exchanging protons explains some deviations from this correlation. The NOE transfer measurements provide information on specific solute-solvent interactions and contribute to a better understanding of solvation phenomena. On the basis of a distinct relationship between steric solvation hindrance and the strength of NOE transfer, the application of such measurements for conformational analysis has been demonstrated for the first time.     

Contribution of Solvation Energy in Protein‐Peptide Recognition Systems

The contribution of solvation energy to guiding molecular recognition for six rigid protein‐peptide systems bad been evaluated by the variation in the number of the identified native‐like configurations and in the driving force of specific interaction resulting from the addition of the explicit solvation term in the force field function. The AMBER force field energy and the total energy including the force field energy and the WZS solvation energy were calculated for sampled configurations. The results obtained by the calculations of both force field and total energies were compared with each other. It suggests that the contribution of solvation energy is important to guiding the specific recognition of the systems in which the ligands possess larger hydrophobic or aromatic residues while the protein receptors provide the active surfaces with hydrophobic property.     

Ionic and dipolar solvation dynamics in liquid water

The solvation time correlation function for solvation in liquid water was measured recently. The solvation was found to be very fast, with a time constant equal to 55 fs. In this article we present theoretical studies on solvation dynamics of ionic and dipolar solutes in liquid water, based on the molecular hydrodynamic approach developed earlier. The molecular hydrodynamic theory can successfully predict the ultrafast dynamics of solvation in liquid water as observed from recent experiments. The present study also reveals some interesting aspects of dipolar solvation dynamics, which differs significantly from that of ionic solvation. Dedicated to Prof. C N R Rao on his 60th birthday     

Cavity shape effect on pKa prediction of small amines

Cavity shape affects the free energy of solvation by which pK_a changes. In this work, we have used different level of theories and basis sets combine with the PCM model of solvation to predict acidity constants of small ammonium ions during their protonation in aqueous solution. Three different kinds of models including UAHF, Bondi, and Pauling have been chosen to study that how a cavity shape changes acidity constant. Contrary to the previous results of literature for bulky ammonium ions, the difference between our results is not very much. Predictions show that in some levels and basis sets UAHF model is the best, while in some other cases it is the worst.     

Coupled reference interaction site model/simulation approach for thermochemistry of solvation: theory and prospects

We present a new methodology for computing solvation free energy, which is based upon the reference interaction site model (RISM)/hypernetted chain (HNC) solvation free energy expression, but which substitutes radial distribution functions taken from simulations for those calculated by simultaneous solution of the RISM and HNC equations. Consequently, solvation free energy can be obtained from a single molecular dynamics or Monte Carlo simulation. Here we describe in detail the coupled RISM/simulation approach, and offer some error analysis. Finally we give the results of its application to a set of small test molecules in aqueous solution. The success shown in some of our results demonstrates that the coupled RISM/simulation approach is worth considering further as a potentially useful tool in studies of solvated systems, such as aqueous molecular biosystems.     

A comparative theoretical study of dipeptide solvation in water

Molecular dynamics studies have been performed on the zwitterionic form of the dipeptide glycine-alanine in water, with focus on the solvation and electrostatic properties using a range of theoretical methods, from purely classical force fields, through mixed quantum mechanical/molecular mechanical simulations, to fully quantum mechanical Car-Parrinello calculations. The results of these studies show that the solvation pattern is similar for all methods used for most atoms in the dipeptide, but can differ substantially for some groups; namely the carboxy and aminoterminii, and the backbone amid NH group. This might have implications in other theoretical studies of peptides and proteins with charged -NH(3) (+) and -CO(2) (-) side chains solvated in water. Hybrid quantum mechanical/molecular mechanical simulations successfully reproduce the solvation patterns from the fully quantum mechanical simulations (PACS numbers: 87.14.Ee, 87.15.Aa, 87.15.He, 71.15.Pd).