Preferential hydration of lysozyme in water/glycerol mixtures: a small-angle neutron scattering study

In solution small-angle neutron scattering has been used to study the solvation properties of lysozyme dissolved in water/glycerol mixtures. To detect the characteristics of the protein-solvent interface, 35 different experimental conditions (i.e., protein concentration, water/glycerol fraction in the solvent, content of deuterated compounds) have been considered and a suitable software has been developed to fit simultaneously the whole set of scattering data. The average composition of the solvent in the close vicinity of the protein surface at each experimental condition has been derived. In all the investigated conditions, glycerol resulted especially excluded from the protein surface, confirming that lysozyme is preferentially hydrated. By considering a thermodynamic hydration model based on an equilibrium exchange between water and glycerol from the solvation layer to the bulk, the preferential binding coefficient and the excess solvation number have been estimated. Results were compared with data previously derived for ribonuclease A in the same mixed solvent: even if the investigated solvent compositions were very different, the agreement between data is noticeable, suggesting that a unique mechanism presides over the preferential hydration process. Moreover, the curve describing the excess solvation number as a function of the solvent composition shows the occurrence of a region of maximal hydration, which probably accounts for the changes in protein stability detected in the presence of cosolvents.     

Solvation of the Whelk-O1 chiral stationary phase: a molecular dynamics study

Density functional theory calculations and molecular dynamics simulations are employed to explore the solvation of the Whelk-O1 chiral stationary phase. First, a semi-flexible representation of the Whelk-O1 selective molecule is extracted from an extensive series of B3LYP/6-311+ G(2d,p) calculations. The resulting model is used to build a chiral surface, including end-caps, for molecular dynamics study of the interface between solvent and Whelk-O1. Three solvent environments in common use for Whelk-O1 HPLC have been examined: a normal-phase solvent of n-hexane/2-propanol; a reversed-phase solvent of water/methanol; and a supercritical solvent of CO(2) and methanol. In each case, we analyze the interface with an emphasis on solvent composition and solvent hydrogen bonding to the Whelk-O1 selector.     

Microcalorimetric study of thermal unfolding of lysozyme in water/glycerol mixtures: an analysis by solvent exchange model

Folded protein stabilization or destabilization induced by cosolvent in mixed aqueous solutions has been studied by differential scanning microcalorimetry and related to difference in preferential solvation of native and denatured states. In particular, the thermal denaturation of a model system formed by lysozyme dissolved in water in the presence of the stabilizing cosolvent glycerol has been considered. Transition temperatures and enthalpies, heat capacity, and standard free energy changes have been determined when applying a two-state denaturation model to microcalorimetric data. Thermodynamic parameters show an unexpected, not linear, trend as a function of solvent composition; in particular, the lysozyme thermodynamic stability shows a maximum centered at water molar fraction of about 0.6. Using a thermodynamic hydration model based on the exchange equilibrium between glycerol and water molecules from the protein solvation layer to the bulk, the contribution of protein-solvent interactions to the unfolding free energy and the changes of this contribution with solvent composition have been derived. The preferential solvation data indicate that lysozyme unfolding involves an increase in the solvation surface, with a small reduction of the protein-preferential hydration. Moreover, the derived changes in the excess solvation numbers at denaturation show that only few solvent molecules are responsible for the variation of lysozyme stability in relation to the solvent composition.     

Shells theory of solvation and the long-range Born correction

A new solvation model, named shells theory of solvation, is proposed. In this approach, the solvent is divided in two regions, the S₁ shell, close to the solute and describing specific solute–solvent interactions, and the S₂ shell, representing the remain solvent and accounting for the long-range interaction contribution. A simple theoretical equation can be derived which allows the computation of the solvation free energy using two-point thermodynamic integration and configurations generated from molecular dynamics simulation. The discrete/continuum version of this theory provides rigorous theoretical foundations for the popular long-range Born correction and presents a new reliable expression for including this contribution. Further, it converges to the full discrete representation of the solvent when the number of solvent molecules goes to infinity. The method can be easily applied when the solute–solvent interaction (S₁ shell) is treated by full quantum mechanics, while the S₂ shell is described by a dielectric continuum solvation method. A simple test of the theory was done for solvation of fluoride ion in benzene solution. The S₁ shell was composed of the fluoride ion plus 32 benzene molecules, and the interaction with the S₂ shell was calculated at Hartree–Fock level with the MINI basis set and using the polarizable continuum model.     

Design of surface active soluble peptide molecules at the air/water interface

Surface active molecules collect at interfaces and have the potential to be used for water evaporation reduction. The objective of this work is to design surface active soluble peptides that collect at the air/water interface using molecular simulations. Rotational isomeric state Monte Carlo (RISMC) sampling together with a solvation model that we recently invented, the AAD solvation model [Gu, C.; Lustig, S.; Trout, B. J. Phys. Chem. B 2006, 110 (3), 1476-1484] was applied to calculate the adsorption free energy of the peptide molecule at the air/water interface. The results were validated by both molecular dynamics simulations with an explicit solvent model and surface tension measurements on synthesized peptides. It was demonstrated that this approach is able to give a reasonable prediction of surface activity with an approximately 50% hit rate in terms of designed surface active molecules actually being surface active. The relationship between the chemical composition and the surface morphology is also discussed.     

Molecular probes of solvation phenomena

The properties of the molecules present in any chemical or biological system are dependent on interactions with the environment, and a quantitative understanding of solvation phenomena remains a major challenge. Molecular recognition probes provide a new approach to quantitatively measure the properties of solvents. Traditionally, solvent polarity scales have been based on spectroscopic probes that provide insight into the nature of solvent-solute interactions. This review compares the solvent polarity parameters obtained from the wavelengths of UV/Visible absorption maxima with solute H-bond parameters obtained from the free energies of solution equilibria. The similarity of the solvent and solute H-bond scales leads to a general H-bond scale that uses the same parameters to describe both solvent and solute. The general H-bond scale provides a framework for understanding the relationship between local intermolecular interactions and the properties of the bulk medium. Intermolecular interactions are sensitive to solvation equilibria, so molecular recognition probes provide fundamentally different information from spectroscopic probes that are sensitive to the populations of different solvation states of the solute. Studies of mixed solvents demonstrate the potential of molecular recognition probes for providing new insights into solvation phenomena.     

Impact of the solvent on the conformational isomerism of calix[4]arenes: a study based on continuum solvation models

The influence of solvation on the conformational isomerism of calix[4]arene and p-tert-butylcalix[4]arene has been investigated by using the continuum model reported by Miertus, Scrocco, and Tomasi (MST). The quantum mechanical (QM) and semiclassical (SC) formalisms of the MST model have been considered for two different solvents (chloroform and water). The suitability of the QM-MST and SC-MST methods has been examined by comparison with previous results derived from classical molecular dynamics (MD) simulations with explicit solvent molecules. The application of the continuum model to the solute configurations generated by using in vacuo classical MD simulations provides a fast strategy to evaluate the effects of the solvent on the conformational preferences of calixarenes. These encouraging results allow us to propose the use of continuum models to solutes with complex molecular structures, which are traditionally studied by MD simulations.     

Partition coefficients of methylated DNA bases obtained from free energy calculations with molecular electron density derived atomic charges

Partition coefficients serve in various areas as pharmacology and environmental sciences to predict the hydrophobicity of different substances. Recently, they have also been used to address the accuracy of force fields for various organic compounds and specifically the methylated DNA bases. In this study, atomic charges were derived by different partitioning methods (Hirshfeld and Minimal Basis Iterative Stockholder) directly from the electron density obtained by electronic structure calculations in a vacuum, with an implicit solvation model or with explicit solvation taking the dynamics of the solute and the solvent into account. To test the ability of these charges to describe electrostatic interactions in force fields for condensed phases, the original atomic charges of the AMBER99 force field were replaced with the new atomic charges and combined with different solvent models to obtain the hydration and chloroform solvation free energies by molecular dynamics simulations. Chloroform–water partition coefficients derived from the obtained free energies were compared to experimental and previously reported values obtained with the GAFF or the AMBER‐99 force field. The results show that good agreement with experimental data is obtained when the polarization of the electron density by the solvent has been taken into account, and when the energy needed to polarize the electron density of the solute has been considered in the transfer free energy. These results were further confirmed by hydration free energies of polar and aromatic amino acid side chain analogs. Comparison of the two partitioning methods, Hirshfeld‐I and Minimal Basis Iterative Stockholder (MBIS), revealed some deficiencies in the Hirshfeld‐I method related to the unstable isolated anionic nitrogen pro‐atom used in the method. Hydration free energies and partitioning coefficients obtained with atomic charges from the MBIS partitioning method accounting for polarization by the implicit solvation model are in good agreement with the experimental values. © 2018 Wiley Periodicals, Inc.     

Spectral shift of the n → π* transition for acetone and formic acid with an explicit solvent model

In our recent work, a new form of the electrostatic solvation energy for the nonequilibrium polarization has been derived by introducing the method of constrained equilibrium state in the framework of continuous medium theory. Up until now, the idea of the constrained equilibrium state method has not been introduced into the explicit solvent model by others; therefore this nonequilibrium energy form was further equivalently extended to the explicit solvent model in this work based on the discrete representation of the solvent permanent charges and induced dipoles. Making use of this expression in explicit solvent model, we modified the nonequilibrium module in the averaged solvent electrostatic potential/molecular dynamics program to implement numerical calculations. Subsequently, the new codes were applied to study the solvatochromic shifts of the n → π* absorption spectra for acetone and trans-formic acid in aqueous solution. The calculation results show a good agreement with the experimental observations. When our results of spectral shift are compared with those achieved directly from the continuum model, it can be seen that both the explicit solvent model and continuum model derived based on the constrained equilibrium approach can give reasonable predictions. The hydrogen bond effect was also discussed and deemed to be a dominant contribution to the spectral shift by calculating the n → π* absorption spectra of acetone-water complexes.     

Solvation thermodynamics of neutral and charged solutes using the solvation-layer interface condition continuum dielectric model

We demonstrate that the solvation-layer interface condition (SLIC) continuum dielectric model for molecular electrostatics, combined with a simple solvent-accessible-surface-area (SASA)-proportional model for nonpolar solvent effects, accurately predicts solvation entropies of neutral and charged small molecules. The SLIC/SASA model has only seven fitting parameters in total and achieves this accuracy using a training set with only 20 compounds. Despite this simplicity, solvation free energies and entropies are nearly as accurate as those predicted by the more sophisticated Langevin dipoles solvation model. Surprisingly, the model automatically reproduces the negligible contribution of electrostatics to the solvation of hydrophobic compounds. Opportunities for improvement include nonpolar solvation, anion solvation entropies, and heat capacities. More molecular realism may be needed for these quantities. To enable a future, explicit-solvent-based assessment of the SLIC/SASA implicit-solvent model, we predict solvation entropies for the Mobley test set, which are available as Supporting Information.     

Implicit nonpolar solvent models

We have systematically analyzed a new nonpolar solvent model that separates nonpolar solvation free energy into repulsive and attractive components. Our analysis shows that either molecular surfaces or volumes can be used to correlate with repulsive free energies of tested molecules in explicit solvent with correlation coefficients higher than 0.99. In addition, the attractive free energies in explicit solvent can also be reproduced with the new model with a correlation coefficient higher than 0.999. Given each component optimized, the new nonpolar solvent model is found to reproduce monomer nonpolar solvation free energies in explicit solvent very well. However, the overall accuracy of the nonpolar solvation free energies is lower than that of each component. In the more challenging dimer test cases, the agreement of the new model with explicit solvent is less impressive. Nevertheless, it is found that the new model works reasonably well for reproducing the relative nonpolar free energy landscapes near the global minimum of the dimer complexes.     

Hybrid supermolecule-polarizable continuum approach to solvation: Application to the mechanism of the Stevens rearrangement

Semiempirical molecular orbital theory has been used to study the effects of solvation by acetonitrile on the Stevens rearrangement of methylammonium formylmethylide to 2-aminopropanal. Three methods of solvation have been used to investigate both the electrostatic and specific solvent–solute effects of solvation: a supermolecule calculation involving the complete geometry optimization of up to six solvent molecules about the solute, the conductor-like screening model (COSMO) polarizable continuum method which allows for geometry optimization of the solute in a solvent defined by its dielectric constant, and a hybrid method in which up to five solvent molecules are incorporated inside the solute cavity and complete geometry optimization of the complex is carried out within the polarizable continuum. A comparison of the calculated geometries, rearrangement activation energies, and enthalpies of solvation from these approaches is presented, and the explicit versus bulk solvation effects are discussed. The overall effect of all methods for incorporating solvation effects is that the radical pair pathway is perferred over the concerted mechanism. © 1996 by John Wiley & Sons, Inc.     

Redox entropy of plastocyanin: developing a microscopic view of mesoscopic polar solvation

We report applications of analytical formalisms and molecular dynamics (MD) simulations to the calculation of redox entropy of plastocyanin metalloprotein in aqueous solution. The goal of our analysis is to establish critical components of the theory required to describe polar solvation at the mesoscopic scale. The analytical techniques include a microscopic formalism based on structure factors of the solvent dipolar orientations and density and continuum dielectric theories. The microscopic theory employs the atomistic structure of the protein with force-field atomic charges and solvent structure factors obtained from separate MD simulations of the homogeneous solvent. The MD simulations provide linear response solvation free energies and reorganization energies of electron transfer in the temperature range of 280-310 K. We found that continuum models universally underestimate solvation entropies, and a more favorable agreement is reported between the microscopic calculations and MD simulations. The analysis of simulations also suggests that difficulties of extending standard formalisms to protein solvation are related to the inhomogeneous structure of the solvation shell at the protein-water interface combining islands of highly structured water around ionized residues along with partial dewetting of hydrophobic patches. Quantitative theories of electrostatic protein hydration need to incorporate realistic density profile of water at the protein-water interface.     

Solvation dynamics in acetonitrile: a study incorporating solute electronic response and nuclear relaxation

The solvent reorganization process after electronic excitation of a polar solute in a polar solvent such as acetonitrile is related mainly to the time evolution of the solute-solvent electrostatic interaction. Modern laser-based techniques have sufficient time resolution to follow this decay in real time, providing information to be confirmed and interpreted by theories and models. We present here a study aimed at the investigation of the different steps involved in the process taking place after a vertical S(0) --> S(1) excitation of a large size chromophore, coumarin 153 (C153), in acetonitrile, from both the solute and the solvent points of view. To do this, we use accurate quantum mechanical calculations for the solute properties within the polarizable continuum model (PCM) and classical molecular dynamics (MD) simulations, both equilibrium and nonequilibrium, for C153 in the presence of the solvent. The geometry of the solute is allowed to change in order to study the role of internal motions in the time-dependent solvation process. The solvent response function has been obtained from the simulation data and compared to experiment, while the comparison between equilibrium and nonequilibrium MD results for the solvation response confirms the validity of the linear response approximation in the C153-acetonitrile system. The MD trajectories have also been used to monitor the structure of the solvation shell and to determine its change in response to the change in the solute partial charges.