On the nonpolar hydration free energy of proteins: surface area and continuum solvent models for the solute-solvent interaction energy

Implicit solvent hydration free energy models are an important component of most modern computational methods aimed at protein structure prediction, binding affinity prediction, and modeling of conformational equilibria. The nonpolar component of the hydration free energy, consisting of a repulsive cavity term and an attractive van der Waals solute-solvent interaction term, is often modeled using estimators based on the solvent exposed solute surface area. In this paper, we analyze the accuracy of linear surface area models for predicting the van der Waals solute-solvent interaction energies of native and non-native protein conformations, peptides and small molecules, and the desolvation penalty of protein-protein and protein-ligand binding complexes. The target values are obtained from explicit solvent simulations and from a continuum solvent van der Waals interaction energy model. The results indicate that the standard surface area model, while useful on a coarse-grained scale, may not be accurate or transferable enough for high resolution modeling studies of protein folding and binding. The continuum model constructed in the course of this study provides one path for the development of a computationally efficient implicit solvent nonpolar hydration free energy estimator suitable for high-resolution structural and thermodynamic modeling of biological macromolecules.     

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.     

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.     

On the molecular mechanism of stabilization of proteins by cosolvents: role of Lifshitz electrodynamic forces

Several ionic and nonionic additives are known to affect structural stability of proteins in aqueous solutions. At a fundamental level, the mechanism of stabilization or destabilization of proteins by cosolvents must be related to three-body interactions between the protein, additive, and the water medium. In this study, the role of the Lifshitz-van der Waals electrodynamic interaction between various additives (sucrose, glycerol, urea, poly(ethylene glycol)-200, betaine, taurine, proline, and valine) and bovine serum albumin (BSA) in water medium was examined. The electrodynamic interaction energy was attractive for all of the additives studied here when both far ultraviolet and infrared relaxations of the additives were included in their dielectric susceptibility representations. However, when only the infrared contribution was included for structure stabilizers and both far ultraviolet and infrared contributions for the structure destabilizers, the resulting electrodynamic interaction energy (E/kT) followed the structure stabilizing and/or destabilizing behavior of the additives; that is, the interaction was attractive for urea and PEG200 (structure destabilizers), whereas it was repulsive for sucrose, glycerol, betaine, taurine, alanine, valine, and proline (structure stabilizers). The electrodynamic interaction energy E/kT at any given surface-to-surface separation distance between the additives and BSA was positively correlated (r(2) = 0.92) with the experimental thermal denaturation temperature (T(d)) of BSA in 1 M solutions of the additives. These analyses provided a mechanistic basis for the experimental observations of exclusion of the structure-stabilizing additives from the protein-water interface and binding of the structure-destabilizing additives to the protein surface. The role of water structure in the three-body electrodynamic interaction is discussed. It is hypothesized that in the case of additives that enhance water structure the hydration shells formed around the additives effectively dampen the contribution of ultraviolet frequencies to the dielectric susceptibility of the additives and thus impart repulsive electrodyanamic interaction between the additive and the protein, whereas the opposite occurs in the case of additives that breakdown the hydrogen-bonded structure of water.     

温度对二甲基亚砜水溶液的结构和热力学性质的影响

应用参考作用格位模型理论计算了二甲基亚砜(DMSO)摩尔分数为0.002时不同温度下溶液的微观结构和热力学性质. 计算结果表明, DMSO加入到水中能够增强溶液的分子网络结构. 温度升高, 配位数减小, 溶液中分子排布趋向无序. 平均力势的波动增大表明分子间的诱导力表现为斥力. 计算得到的各种热力学性质显示: 温度升高, 溶液的熵和溶剂化自由能增加, 相互作用能和过剩化学位也增加, 即高温下溶液越来越偏离理想溶液; 空位形成能降低表明溶液分子结构在高温下更容易重组.     

Theoretical analysis on changes in thermodynamic quantities upon protein folding: essential role of hydration

The free energy change associated with the coil-to-native structural transition of protein G in aqueous solution is calculated by using the molecular theory of solvation, also known as the three-dimensional reference interaction site model theory, to uncover the molecular mechanism of protein folding. The free energy is decomposed into the protein intramolecular energy, the hydration energy, and the hydration entropy. The folding is accompanied with a large gain in the protein intramolecular energy. However, it is almost canceled by the correspondingly large loss in the hydration energy due to the dehydration, resulting in the total energy gain about an order of magnitude smaller than might occur in vacuum. The hydration entropy gain is found to be a substantial driving force in protein folding. It is comparable with or even larger than the total energy gain. The total energy gain coupled with the hydration entropy gain is capable of suppressing the conformational entropy loss in the folding. Based on careful analysis of the theoretical results, the authors present a challenging physical picture of protein folding where the overall folding process is driven by the water entropy effect.     

Protein-protein interactions in complex cosolvent solutions.

The effects of various kosmotropic and chaotropic cosolvents and salts on the intermolecular interaction potential of positively charged lysozyme is evaluated at varying protein concentrations by using synchrotron small-angle X-ray scattering in combination with liquid-state theoretical approaches. The experimentally derived static structure factors S(Q) obtained without and with added cosolvents and salts are analysed with a statistical mechanical model based on the Derjaguin-Landau-Verwey-Overbeek (DLVO) potential, which accounts for repulsive and attractive interactions between the protein molecules. Different cosolvents and salts influence the interactions between protein molecules differently as a result of changes in the hydration level or solvation, in charge screening, specific adsorption of the additives at the protein surface, or increased hydrophobic interactions. Intermolecular interaction effects are significant above protein concentrations of 1 wt %, and with increasing protein concentration, the repulsive nature of the intermolecular pair potential V(r) increases markedly. Kosmotropic cosolvents like glycerol and sucrose exhibit strong concentration-dependent effects on the interaction potential, leading to an increase of repulsive forces between the protein molecules at low to medium high osmolyte concentrations. Addition of trifluoroethanol exhibits a multiphasic effect on V(r) when changing its concentration. Salts like sodium chloride and potassium sulfate exhibit strong concentration-dependent changes of the interaction potential due to charge screening of the positively charged protein molecules. Guanidinium chloride (GdmCl) at low concentrations exhibits a similar charge-screening effect, resulting in increased attractive interactions between the protein molecules. At higher GdmCl concentrations, V(r) becomes more repulsive in nature due to the presence of high concentrations of Gdm(+) ions binding to the protein molecules. Our findings also imply that in calculations of thermodynamic properties of proteins in solution and cosolvent mixtures, activity coefficients may not generally be neglected in the concentration range above 1 wt % protein.     

Protein-water electrostatics and principles of bioenergetics

Despite its diversity, life universally relies on a simple basic mechanism of energy transfer in its energy chains-hopping electron transport between centers of electron localization on hydrated proteins and redox cofactors. Since many such hops connect the point of energy input with a catalytic site where energy is stored in chemical bonds, the question of energy losses in (nearly activationless) electron hops, i.e., energetic efficiency, becomes central for the understanding of the energetics of life. We show here that standard considerations based on rules of Gibbs thermodynamics are not sufficient, and the dynamics of the protein and the protein-water interface need to be involved. The rate of electronic transitions is primarily sensitive to the electrostatic potential at the center of electron localization. Numerical simulations show that the statistics of the electrostatic potential produced by hydration water are strongly non-Gaussian, with the breadth of the electrostatic noise far exceeding the expectations of the linear response. This phenomenon, which dramatically alters the energetic balance of a charge-transfer chain, is attributed to the formation of ferroelectric domains in the protein's hydration shell. These dynamically emerging and dissipating domains make the shell enveloping the protein highly polar, as gauged by the variance of the shell dipole which correlates with the variance of the protein dipole. The Stokes-shift dynamics of redox-active proteins are dominated by a slow component with the relaxation time of 100-500 ps. This slow relaxation mode is frozen on the time-scale of fast reactions, such as bacterial charge separation, resulting in a dramatically reduced reorganization free energy of fast electronic transitions. The electron transfer activation barrier becomes a function of the corresponding rate, self-consistently calculated from a non-ergodic version of the transition-state theory. The peculiar structure of the protein-water interface thus provides natural systems with two "non's"-non-Gaussian statistics and non-ergodic kinetics-to tune the efficiency of the redox energy transfer. Both act to reduce the amount of free energy released as heat in electronic transitions. These mechanisms are shown to increase the energetic efficiency of protein electron transfer by up to an order of magnitude compared to the "standard picture" based on canonical free energies and the linear response approximation. In other words, the protein-water tandem allows both the formation of a ferroelectric mesophase in the hydration shell and an efficient control of the energetics by manipulating the relaxation times.     

Secondary structure sensitivity of hydrogen bond lifetime dynamics in the protein hydration layer

The heterogeneous nature of a protein surface plays an essential role in its biological activity and molecular recognition, and this role is mediated at least partly through the surrounding water molecules. We have performed atomistic molecular dynamics simulations of an aqueous solution of HP-36 to investigate the correlation between the dynamics of the hydration layer water molecules and the lifetimes of protein-water hydrogen bonds. The nonexponential hydrogen bond lifetime correlation functions have been analyzed by using the formalism of Luzar and Chandler, which allowed identification of the quasi-bound states in the surface and quantification of the dynamic equilibrium between quasi-bound and free water molecules in terms of time-dependent rate of interconversion. It is noticed that, irrespective of the structural heterogeneity of different segments of the protein, namely the three alpha-helices, the positively charged amino acid residues form longer-lived hydrogen bonds with water. The overall relaxation behavior of protein-water hydrogen bonds is found to differ significantly among the three helices of the protein. Study of water number density fluctuation reveals that the hydration layer of helix-3 is much less rigid, which can be correlated with faster structural relaxation of the hydrogen bonds between its residues and water. This also agrees excellently with faster translational and rotational motions of water near helix-3, and hence the lower rigidity of its hydration layer. The lower rigidity of the helix-3 hydration layer also correlates well with the biological activity of the protein, as several of the active-site residues of HP-36 are located in helix-3.     

Free-energy analysis of the electron-density fluctuation in the quantum-mechanical/molecular-mechanical simulation combined with the theory of energy representation

The relationship is investigated for QM/MM (quantum-mechanical/molecular-mechanical) systems between the fluctuations of the electronic state of the QM subsystem and of the solvation effect due to the QM-MM interaction. The free-energy change due to the electron-density fluctuation around its average is highlighted, and is evaluated through an approximate functional formulated in terms of distribution functions of the many-body coupling (pairwise non-additive) part of the QM-MM interaction energy. A set of QM/MM simulations are conducted in MM water solvent for QM water solute in ambient and supercritical conditions and for QM glycine solute in the neutral and zwitterionic forms. The variation of the electronic distortion energy of the QM solute in the course of QM/MM simulation is then shown to be compensated by the corresponding variation of the free energy of solvation. The solvation free energy conditioned by the electronic distortion energy is further analyzed with its components. It is found that the many-body contribution is essentially equal between the free energy and the average sum of solute-solvent interaction energy.     

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.     

Free-energy analysis of the molecular binding into lipid membrane with the method of energy representation

A statistical-mechanical treatment of the molecular binding into lipid membrane is presented in combination with molecular simulation. The membrane solution is viewed as an inhomogeneous, mixed solvent system, and the free energy of solvation of a solute in membrane is computed with a realistic set of potential functions by the method of energy representation. Carbon monoxide, carbon dioxide, benzene, and ethylbenzene are adopted as model solutes to analyze the binding into 1,2-dimyristoyl-sn-glycero-3-phosphatidylcholine (DMPC) membrane. It is shown that the membrane inside is more favorable than bulk water and that the solute distribution is diffuse throughout the membrane inside. The membrane-water partition coefficient is then constructed with the help of the Kirkwood-Buff theory from the solvation free energy obtained separately in the hydrophobic, glycerol, headgroup, and aqueous regions. To discuss the role of repulsive and attractive interactions, the solvation free energy is partitioned into the DMPC and water contributions and the effect of water to stabilize the benzene and ethylbenzene solutes within the membrane is pointed out.     

Protein interactions studied by SAXS: effect of ionic strength and protein concentration for BSA in aqueous solutions

We have studied a series of samples of bovine serum albumin (BSA) solutions with protein concentration, c, ranging from 2 to 500 mg/mL and ionic strength, I, from 0 to 2 M by small-angle X-ray scattering (SAXS). The scattering intensity distribution was compared to simulations using an oblate ellipsoid form factor with radii of 17 x 42 x 42 A, combined with either a screened Coulomb, repulsive structure factor, SSC(q), or an attractive square-well structure factor, SSW(q). At pH = 7, BSA is negatively charged. At low ionic strength, I < 0.3 M, the total interaction exhibits a decrease of the repulsive interaction when compared to the salt-free solution, as the net surface charge is screened, and the data can be fitted by assuming an ellipsoid form factor and screened Coulomb interaction. At moderate ionic strength (0.3-0.5 M), the interaction is rather weak, and a hard-sphere structure factor has been used to simulate the data with a higher volume fraction. Upon further increase of the ionic strength (I >or= 1.0 M), the overall interaction potential was dominated by an additional attractive potential, and the data could be successfully fitted by an ellipsoid form factor and a square-well potential model. The fit parameters, well depth and well width, indicate that the attractive potential caused by a high salt concentration is weak and long-ranged. Although the long-range, attractive potential dominated the protein interaction, no gelation or precipitation was observed in any of the samples. This is explained by the increase of a short-range, repulsive interaction between protein molecules by forming a hydration layer with increasing salt concentration. The competition between long-range, attractive and short-range, repulsive interactions accounted for the stability of concentrated BSA solution at high ionic strength.     

Mapping Hydration Water around Alcohol Chains by THz Calorimetry

THz spectroscopy was used to probe changes that occur in the dynamics of the hydrogen bond network upon solvation of alcohol chains. The THz spectra can be decomposed into the spectrum of bulk water, tetrahedral hydration water, and more disordered (or interstitial) hydration water. The tetrahedrally ordered hydration water exhibits a band at 195 cm⁻¹ and is localized around the hydrophobic moiety of the alcohol. The interstitial component yields a band at 164 cm⁻¹ which is associated with hydration water in the first hydration shell. These temperature‐dependent changes in the low‐frequency spectrum of solvated alcohol chains can be correlated with changes of heat capacity, entropy, and free energy upon solvation. Surprisingly, not the tetrahedrally ordered component but the interstitial hydration water is found to be mainly responsible for the temperature‐dependent change in ΔC_p and ΔG. The solute‐specific offset in free energy is attributed to void formation and scales linearly with the chain length.     

Femtosecond studies of tryptophan fluorescence dynamics in proteins: local solvation and electronic quenching

We report our systematic examination of tryptophan fluorescence dynamics in proteins with femtosecond resolution. Distinct patterns of femtosecond-resolved fluorescence transients from the blue to the red side of emission have been characterized to distinguish local ultrafast solvation and electronic quenching. It is shown that tryptophan is an ideal local optical probe for hydration dynamics and protein-water interactions as well as an excellent local molecular reporter for ultrafast electron transfer in proteins, as demonstrated by a series of biological systems, here in melittin, human serum albumin, and human thioredoxin, and at lipid interfaces. These studies clarify the assignments in the literature about the ultrafast solvation or quenching dynamics of tryptophan in proteins. We also report a new observation of solvation dynamics at far red-side emission when the relaxation of the local environment is slower than 1 ps. These results provide a molecular basis for using tryptophan as a local molecular probe for ultrafast protein dynamics in general.     

Simple and exact approach to the electronic polarization effect on the solvation free energy: formulation for quantum-mechanical/molecular-mechanical system and its applications to aqueous solutions

We present a simple and exact numerical approach to compute the free energy contribution δμ in solvation due to the electron density polarization and fluctuation of a quantum-mechanical solute in the quantum-mechanical/molecular-mechanical (QM/MM) simulation combined with the theory of the energy representation (QM/MM-ER). Since the electron density fluctuation is responsible for the many-body QM-MM interactions, the standard version of the energy representation method cannot be applied directly. Instead of decomposing the QM-MM polarization energy into the pairwise additive and non-additive contributions, we take sum of the polarization energies in the QM-MM interaction and adopt it as a new energy coordinate for the method of energy representation. Then, it is demonstrated that the free energy δμ can be exactly formulated in terms of the energy distribution functions for the solution and reference systems with respect to this energy coordinate. The benchmark tests were performed to examine the numerical efficiency of the method with respect to the changes in the individual properties of the solvent and the solute. Explicitly, we computed the solvation free energy of a QM water molecule in ambient and supercritical water, and also the free-energy change associated with the isomerization reaction of glycine from neutral to zwitterionic structure in aqueous solution. In all the systems examined, it was demonstrated that the computed free energy δμ agrees with the experimental value, irrespective of the choice of the reference electron density of the QM solute. The present method was also applied to a prototype reaction of adenosine 5'-triphosphate hydrolysis where the effect of the electron density fluctuation is substantial due to the excess charge. It was demonstrated that the experimental free energy of the reaction has been accurately reproduced with the present approach.     

Comparison of methods to estimate the free energy of solvation:Importance in the modulation of the affinity of 3-benzazepines for the D1receptor

We computed the free energy of solvation for a series of ions and neutral molecules using two different continuum approaches. First, we used the AM1–SM1 technique, where the AM1 Fock matrix is modified to include a generalized Born contribution. Second, we applied the DelPhi approach, where the electrostatic component of the free energy of solvation is evaluated by resolving the Poisson–Boltzman equation by a finite difference method. Both methods appear equally reliable for ionic systems. For neutral compounds, AM1–SM1 performs better than DelPhi; however, the differences become less pronounced for compounds with larger free energies of solvation. In parallel, both methods were applied to study the influence of the solvation process in the overall drug receptor interaction for a series of closely related ligands for the D₁ dopamine receptor. An inverse linear relationship was found between the free energy of solvation and the logarithm of the affinity of the ligands; nevertheless, electrostatic properties are likely to modulate affinity as well. © 1993 John Wiley & Sons, Inc.     

Filled and empty states of carbon nanotubes in water: Dependence on nanotube diameter, wall thickness and dispersion interactions

We have carried out a series of molecular dynamics simulations of water containing a narrow carbon nanotube as a solute to investigate the filling and emptying of the nanotube and also the modifications of the density and hydrogen bond distributions of water inside and also in the vicinity of the outer surfaces of the nanotube. Our primary goal is to look at the effects of varying nanotube diameter, wall thickness and also solute-solvent interactions on the solvent structure in the confined region also near the outer surfaces of the solute. The thickness of the walls is varied by considering single and multi-walled nanotubes and the interaction potential is varied by tuning the attractive strength of the 12–6 pair interaction potential between a carbon atom of the nanotubes and a water molecule. The calculations are done for many different values of the tuning parameter ranging from fully Lennard-Jones to pure repulsive pair interactions. It is found that both the solvation characteristics and hydrogen bond distributions can depend rather strongly on the strength of the attractive part of the solute-water interaction potential. The thickness of the nanotube wall, however, is found to have only minor effects on the density profiles, hydrogen bond network and the wetting characteristics. This indicates that the long range electrostatic interactions between water molecules inside and on the outer side of the nanotube do not make any significant contribution to the overall solvation structure of these hydrophobic solutes. The solvation characteristics are primarily determined by the balance between the loss of energy due to hydrogen bond network disruption, cavity repulsion potential and offset of the same by attractive component of the solute-water interactions. Our studies with different system sizes show that the essential features of wetting and dewetting characteristics of narrow nanotubes for different diameter and interaction potentials are also present in relatively smaller systems consisting of about five hundred molecules. We dedicate this work to Professor Debashis Mukherjee on his 60th Birthday.     

Estimating the Gibbs energy of hydration from molecular dynamics trajectories obtained by integral equations of the theory of liquids in the RISM approximation

A method of integral equations of the theory of liquids in the reference interaction site model (RISM) approximation is used to estimate the Gibbs energy averaged over equilibrium trajectories computed by molecular mechanics. Peptide oxytocin is selected as the object of interest. The Gibbs energy is calculated using all chemical potential formulas introduced in the RISM approach for the excess chemical potential of solvation and is compared with estimates by the generalized Born model. Some formulas are shown to give the wrong sign of Gibbs energy changes when peptide passes from the gas phase into water environment; the other formulas give overestimated Gibbs energy changes with the right sign. Note that allowance for the repulsive correction in the approximate analytical expressions for the Gibbs energy derived by thermodynamic perturbation theory is not a remedy.