Amino acid side chain interactions in the presence of salts

The effects of salt on the intermolecular interactions between polar/charged amino acids are investigated through molecular dynamics simulations. The mean forces and associated potentials are calculated for NaCl salt in the 0-2 M concentration range at 298 K. It is found that the addition of salt may stabilize or destabilize the interactions, depending on the nature of the interacting molecules. The degree of (de)stabilization is quantified, and the origin of the salt-dependent modulation is discussed based upon an analysis of solvent density profiles. To gain insight into the molecular origin of the salt modulation, spatial distribution functions (sdf's) are calculated, revealing a high degree of solvent structuredness in all cases. The peaks in the sdf's are consistent with long-range hydrogen-bonding networks connecting the solute hydrophilic groups, and that contribute to their intermolecular solvent-induced forces. The restructuring of water around the solutes as they dissociate from close contact is analyzed. This analysis offers clues on how the solvent structure modulates the effective intermolecular interactions in complex solutes. This modulation results from a critical balance between bulk electrostatic forces and those exerted by (i) the water molecules in the structured region between the monomers, which is disrupted by ions that transiently enter the hydration shells, and (ii) the ions in the hydration shells in direct interactions with the solutes. The implications of these findings in protein/ligand (noncovalent) association/dissociation mechanisms are briefly discussed.     

On the influence of solute polarizability on the hydrophobic interaction

The authors have performed molecular dynamics simulations of polarizable solutes in water to investigate how solute polarizability affects solute-solute hydrophophic interactions. A degree of polarization similar to the one expected in biomolecules, corresponding to a dielectric response of epsilon=2-20, results in dramatic changes in the hydrophobic forces. They find that this degree of polarizability is enough to inhibit drying between hydrophobic solutes and to stabilize a reduced water density phase whose density is smaller than the bulk water density. The hydrophobic forces associated with such reduced density states are still very significant with values of the order of several tens of piconewtons. Their results suggest that polarizability plays an important role in determining the hydrophobic force acting between weakly polar surfaces.     

Combined depletion and electrostatic forces in polymer-induced membrane adhesion: a theoretical model

We develop a semi-quantitative analytical theory to describe adhesion between two identical planar charged surfaces embedded in a polymer-containing electrolyte solution. Polymer chains are uncharged and differ from the solvent by their lower dielectric permittivity. The solution mimics physiological fluids: It contains 0.1 M of monovalent ions and a small number of divalent cations that form tight bonds with the headgroups of charged lipids. The components have heterogeneous spatial distributions. The model was derived self-consistently by combining: (a) a Poisson-Boltzmann like equation for the charge densities, (b) a continuum mean-field theory for the polymer profile, (c) a solvation energy forcing the ions toward the polymer-poor regions, and (d) surface interactions of polymers and electrolytes. We validated the theory via extensive coarse-grained Molecular Dynamics (MD) simulations. The results confirm our analytical model and reveal interesting details not detected by the theory. At high surface charges, polymer chains are mainly excluded from the gap region, while the concentration of ions increases. The model shows a strong coupling between osmotic forces, surface potential and salting-out effects of the slightly polar polymer chains. It highlights some of the key differences in the behaviour of monomeric and polymeric mixed solvents and their responses to Coulomb interactions. Our main findings are: (a) the onset of long-ranged ion-induced polymer depletion force that increases with surface charge density and (b) a polymer-modified repulsive Coulomb force that increases with surface charge density. Overall, the system exhibits homeostatic behaviour, resulting in robustness against variations in the amount of charges. Applications and extensions of the model are briefly discussed.     

Accounting for electronic polarization in non-polarizable force fields

The issues of electronic polarizability in molecular dynamics simulations are discussed. We argue that the charges of ionized groups in proteins, and charges of ions in conventional non-polarizable force fields such as CHARMM, AMBER, GROMOS, etc should be scaled by a factor about 0.7. Our model explains why a neglect of electronic solvation energy, which typically amounts to about a half of total solvation energy, in non-polarizable simulations with un-scaled charges can produce a correct result; however, the correct solvation energy of ions does not guarantee the correctness of ion-ion pair interactions in many non-polarizable simulations. The inclusion of electronic screening for charged moieties is shown to result in significant changes in protein dynamics and can give rise to new qualitative results compared with the traditional non-polarizable force field simulations. The model also explains the striking difference between the value of water dipole μ～ 3D reported in recent ab initio and experimental studies with the value μ(eff)～ 2.3D typically used in the empirical potentials, such as TIP3P or SPC/E. It is shown that the effective dipole of water can be understood as a scaled value μ(eff) = μ/√ε(el), where ε(el) = 1.78 is the electronic (high-frequency) dielectric constant of water. This simple theoretical framework provides important insights into the nature of the effective parameters, which is crucial when the computational models of liquid water are used for simulations in different environments, such as proteins, or for interaction with solutes.     

Competition of hydrophobic and Coulombic interactions between nanosized solutes

The solvation of charged, nanometer-sized spherical solutes in water, and the effective, solvent-induced force between two such solutes are investigated by constant temperature and pressure molecular dynamics simulations of model solutes carrying various charge patterns. The results for neutral solutes agree well with earlier findings, and with predictions of simple macroscopic considerations: substantial hydrophobic attraction may be traced back to strong depletion ("drying") of the solvent between the solutes. This hydrophobic attraction is strongly reduced when the solutes are uniformly charged, and the total force becomes repulsive at sufficiently high charge; there is a significant asymmetry between anionic and cationic solute pairs, the latter experiencing a lesser hydrophobic attraction. The situation becomes more complex when the solutes carry discrete (rather than uniform) charge patterns. Due to antagonistic effects of the resulting hydrophilic and hydrophobic "patches" on the solvent molecules, water is once more significantly depleted around the solutes, and the effective interaction reverts to being mainly attractive, despite the direct electrostatic repulsion between solutes. Examination of a highly coarse-grained configurational probability density shows that the relative orientation of the two solutes is very different in explicit solvent, compared to the prediction of the crude implicit solvent representation. The present study strongly suggests that a realistic modeling of the charge distribution on the surface of globular proteins, as well as the molecular treatment of water, are essential prerequisites for any reliable study of protein aggregation.     

Brownian dynamics simulation of a model simple electrolyte in solvents of low dielectric constant

Brownian dynamics simulations are performed to investigate the ionic transport of model simple electrolytes, in which ions are interacting with each other through the repulsive core and Coulombic interactions. The equivalent conductivity and self-diffusion coefficient show minima as the function of the number density of ions when the dielectric constant of the solvent is low. Although the minimum of the former is in harmony with various experiments, no experiment has ever been reported on that of the latter. The analysis of time-dependent transport coefficients reveals that the presence of the minima is ascribed to the slow dynamics, rather than to static association models. The inclusion of a model function that resembles the short-range part of the potential of mean force induced by solvent affects the transport coefficients qualitatively, which suggests the importance of solvent-induced potential of mean force in the conduction mechanism of electrolytes in solvents of low dielectric constant.     

Density functional theory for efficient ab initio molecular dynamics simulations in solution

We present a density functional for first-principles molecular dynamics simulations that includes the electrostatic effects of a continuous dielectric medium. It allows for numerical simulations of molecules in solution in a model polar solvent. We propose a smooth dielectric model function to model solvation into water and demonstrate its good numerical properties for total energy calculations and constant energy molecular dynamics.     

Molecular density functional theory of solvation: from polar solvents to water

A classical density functional theory approach to solvation in molecular solvent is presented. The solvation properties of an arbitrary solute in a given solvent, both described by a molecular force field, can be obtained by minimization of a position and orientation-dependent free-energy density functional. In the homogeneous reference fluid approximation, limited to two-body correlations, the unknown excess term of the functional approximated by the angular-dependent direct correlation function of the pure solvent. We show that this function can be extracted from a preliminary MD simulation of the pure solvent by computing the angular-dependent pair distribution function and solving subsequently the molecular Ornstein-Zernike equation using a discrete angular representation. The corresponding functional can then be minimized in the presence of an arbitrary solute on a three-dimensional cubic grid for positions and Gauss-Legendre angular grid for orientations to provide the solvation structure and free-energy. This two-step procedure is proved to be much more efficient than direct molecular dynamics simulations combined to thermodynamic integration schemes. The approach is shown to be relevant and accurate for prototype polar solvents such as the Stockmayer solvent or acetonitrile. For water, although correct for neutral or moderately charged solute, it tends to underestimate the tetrahedral solvation structure around H-bonded solutes, such as spherical ions. This can be corrected by introducing suitable three-body correlation terms that restore both an accurate hydration structure and a satisfactory energetics.     

Incorporating solvent and ion screening into molecular dynamics using the finite-difference Poisson–Boltzmann method

The finite difference method for solving the Poisson–Boltzmann equation is used to calculate the reaction field acting on a macromolecular solute due to the surrounding water and ions. Comparisons with analytical test cases indicate that the solvation forces can be calculated rapidly and accurately with this method. These forces act to move charged solute atoms towards the solvent where they are better solvated, and to screen interactions between charges. A way of combining such calculations with conventional molecular dynamics force fields is proposed which requires little modification of existing molecular dynamics programs. Simulations on the alanine dipeptide show that solvent forces affect the conformational dynamics by reducing the preference for internal H-bonding forms, increasing the R-alpha helix preference and reducing transition barriers. These solvent effects are similar to previous explicit solvent simulations, but require little more computation than vacuum simulations. The method should scale up with little increase in computational cost to larger molecules such as proteins and nucleic acids.     

Charge asymmetries in hydration of polar solutes

We study the solvation of polar molecules in water. The center of water's dipole moment is offset from its steric center. In common water models, the Lennard-Jones center is closer to the negatively charged oxygen than to the positively charged hydrogens. This asymmetry of water's charge sites leads to different hydration free energies of positive versus negative ions of the same size. Here, we explore these hydration effects for some hypothetical neutral solutes, and two real solutes, with molecular dynamics simulations using several different water models. We find that, like ions, polar solutes are solvated differently in water depending on the sign of the partial charges. Solutes having a large negative charge balancing diffuse positive charges are preferentially solvated relative to those having a large positive charge balancing diffuse negative charges. Asymmetries in hydration free energies can be as large as 10 kcal/mol for neutral benzene-sized solutes. These asymmetries are mainly enthalpic, arising primarily from the first solvation shell water structure. Such effects are not readily captured by implicit solvent models, which respond symmetrically with respect to charge.     

Contact mechanics of nanometer-scale molecular contacts: correlation between adhesion, friction, and hydrogen bond thermodynamics

Using a scanning force microscope, adhesion forces have been measured between carboxylic acid terminated self-assembled monolayers in different nonpolar solvents or in two-component liquid mixtures consisting of a polar solvent (ethyl acetate or acetone) in heptane. The adhesion forces measured in pure acetone and ethyl acetate were small (0.24 nN) but increased logarithmically as the concentration of the polar solvent decreased to reach a maximum value (2.77 nN), equal to that measured in pure heptane, and for lower concentrations of polar solvent, the adhesion force remained constant. This behavior is identical to that observed for association constants measured for the formation of 1:1 H-bonded complexes between dilute solutes in solvent mixtures. The transition between the solvent-dependent and -independent regimes occurs at a polar solvent concentration corresponding to 1/K(S), where K(S) is the equilibrium constant for solvation of a carboxylic acid by the polar solvent in heptane. A simple model, in which the solvation of the carboxylic acid groups may be estimated by considering the concentration and polarity of functional groups in the liquid, accurately predicts values of K(S) that were found to correlate very well with the observed solvent-dependence of the adhesion force. Friction-load relationships were measured using friction-force microscopy. In pure acetone and ethyl acetate, a linear friction-load relationship was observed, in agreement with Amontons' law. However, as the concentration of polar solvent was reduced, a nonlinear relationship was observed and the friction-load relationship was found to fit the Derjaguin-Müller-Toporov (DMT) model for single asperity contacts. For pure heptane and a range of other nonpolar liquids with identical dielectric constants, the friction-load relationship was described by DMT mechanics. Exceptionally, for perfluorodecalin, Johnson-Kendall-Roberts mechanics was observed. These observations may be rationalized by treating the friction force as the sum of load-dependent and shear contributions. Under conditions of low adhesion, where the carboxylic acid surface is solvated by polar solvent molecules, the shear term is negligible and the sliding interaction is dominated by load-dependent friction. As the degree of solvation of the carboxylic acid groups decreases and the adhesion force increases, the shear friction contribution increases, dominating the interaction for media in which the adhesion force is greater than ca. 0.6 nN.     

Effects of solvent density on retention in gas-liquid chromatography. II. Polar solutes in poly(ethylene glycol) stationary phases

Effects of solvent density on the solubility of polar probes which undergo specific interactions with poly(oxyethylene) are studied. The analysis of retention data on capillary columns coated with oligomeric poly(oxyethylene) stationary phases shows that, within the experimental error, the enthalpic contribution to the solubility is practically independent of variations in the solvent density. Average values of enthalpies of solute transfer are reported for different probes and temperatures. The observed systematic decrease of solubility with the increasing density is due to a change of entropy. Some thermodynamic consequences inferred from these general results are discussed. One relevant observation is that the influence of solvent's final groups must be negligible. This is even the case for oligomers with number-average degrees of polymerization as low as 13, hosting solutes capable of strong interactions with the end hydroxyl groups of linear poly(ethylene glycols). Possible explanations for this behavior are explored through molecular dynamics simulations of the liquid solvent.     

Dielectric saturation of the ion hydration shell and interaction between two double helices of DNA in mono- and multivalent electrolyte solutions: foundations of the epsilon-modified Poisson-Boltzmann theory

Potentials of mean force between single Na+, Ca2+, and Mg2+ cations and a highly charged spherical macroion in SPC/E water have been determined using molecular dynamics simulations. Results are compared to the electrostatic energy calculations for the primitive polarization model (PPM) of hydrated cations describing the ion hydration shell as a dielectric sphere of low permittivity (Gavryushov, S.; Linse, P. J. Phys. Chem. B 2003, 107, 7135). Parameters of the ion dielectric sphere and radius of the macroion/water dielectric boundary were extracted by means of this comparison to approximate the short-range repulsion of ions near the interface. To explore the counterion distributions around a simplified model of DNA, the obtained PPM parameters for Na+ and Ca2+ have been substituted into the modified Poisson-Boltzmann (MPB) equations derived for the PPM and named the epsilon-MPB (epsilon-MPB) theory. epsilon-MPB results for DNA suggest that such polarization effects are important in the case of 2:1 electrolyte and highly charged macromolecules. The three-dimensional implementation of the epsilon-MPB theory was also applied to calculation of the energies of interaction between two parallel macromolecules of DNA in solutions of NaCl and CaCl2. Being compared to results of MPB calculations without the ion polarization effects, it suggests that the ion hydration shell polarization and inhomogeneous solvent permittivity might be essential factors in the experimentally known hydration forces acting between charged macromolecules and bilayers at separations of less than 20 A between their surfaces.     

Role of ions on structure and stability of a synthetic gramicidin ion channel in solution. A molecular dynamics study

We performed a molecular dynamics (MD) simulation to the investigate structure and stability of a synthetic gramicidin-like peptide in solution with and without ions. The starting structures of the MD simulations were taken from two recently solved NMR structures of this peptide in isotropic solution, which forms stable monomers or dimers in the presence or absence of ions, respectively. The monomeric structure is channel-like and is assumed to be stabilized by the presence of two Cs(+) ions bound in the channel, each one close to one channel entrance. In our MD simulations, we observed how the Cs(+) ions bind in the channel formed by the monomeric gramicidin-like peptide using implicit solvent and explicit ions with a concentration of 2 M. MD simulations were performed with and without explicit ions but with an implicit solvent model defined by the generalized Born approximation, which was used to mimic the dielectric properties of the solvent and to speed up the computations.     

Effective interaction between two ion-adsorbing plates: Hofmeister series and salting-in/salting-out phase diagrams from a global mean-field analysis

Interactions between ions and solutes are key to ion-specificity. A generic model in which ions interact via square well potentials of finite range with charged plates is solved analytically on the Poisson-Boltzmann level and analyzed globally for varying surface charge, salt concentration, and ion-surface affinity. Ion adsorption as well as depletion can lead to stably bound plates at finite separation, relevant for the equilibrium salting-out of small solutes such as proteins. The interplate pressure at large plate separation, relevant for aggregation kinetics of large solutes, exhibits direct as well as indirect Hofmeister ordering, depending on surface charge and salt concentration. A simple method for mapping explicit ion-surface potentials of mean force as obtained from solvent-explicit molecular dynamics simulations onto square-well potential parameters is demonstrated.     

Self-consistent treatment of the local dielectric permittivity and electrostatic potential in solution for polarizable macromolecular force fields

A self-consistent method is presented for the calculation of the local dielectric permittivity and electrostatic potential generated by a solute of arbitrary shape and charge distribution in a polar and polarizable liquid. The structure and dynamics behavior of the liquid at the solute∕liquid interface determine the spatial variations of the density and the dielectric response. Emphasis here is on the treatment of the interface. The method is an extension of conventional methods used in continuum protein electrostatics, and can be used to estimate changes in the static dielectric response of the liquid as it adapts to charge redistribution within the solute. This is most relevant in the context of polarizable force fields, during electron structure optimization in quantum chemical calculations, or upon charge transfer. The method is computationally efficient and well suited for code parallelization, and can be used for on-the-fly calculations of the local permittivity in dynamics simulations of systems with large and heterogeneous charge distributions, such as proteins, nucleic acids, and polyelectrolytes. Numerical calculation of the system free energy is discussed for the general case of a liquid with field-dependent dielectric response.     

Balancing solvation and intramolecular interactions: toward a consistent generalized Born force field

The efficient and accurate characterization of solvent effects is a key element in the theoretical and computational study of biological problems. Implicit solvent models, particularly generalized Born (GB) continuum electrostatics, have emerged as an attractive tool to study the structure and dynamics of biomolecules in various environments. Despite recent advances in this methodology, there remain limitations in the parametrization of many of these models. In the present work, we demonstrate that it is possible to achieve a balanced implicit solvent force field by further optimizing the input atomic radii in combination with adjusting the protein backbone torsional energetics. This parameter optimization is guided by the potentials of mean force (PMFs) between amino acid polar groups, calculated from explicit solvent free energy simulations, and by conformational equilibria of short peptides, obtained from extensive folding and unfolding replica exchange molecular dynamics (REX-MD) simulations. Through the application of this protocol, the delicate balance between the competing solvation forces and intramolecular forces appears to be better captured, and correct conformational equilibria for a range of both helical and beta-hairpin peptides are obtained. The same optimized force field also successfully folds both beta-hairpin trpzip2 and mini-protein Trp-Cage, indicating that it is quite robust. Such a balanced, physics-based force field will be highly applicable to a range of biological problems including protein folding and protein structural dynamics.     

Electrode screening by ionic liquids

In this work we are concerned with the short-range screening provided by the ionic liquid dimethylimidazolium chloride near a charged wall. We study the free energy profiles (or potentials of mean force) for charged and neutral solutes as a function of distance from a charged wall. Four different wall charge densities are used in addition to a wall with zero charge. The highest magnitude of the charge densities is ±1 e nm(-2) which is close to the maximum limit of charge densities accessible in experiments, while the intermediate charges ±0.5 e nm(-2) are in the range of densities typically used in most of the experimental studies. Positively and negatively charged solutes of approximately the size of a BF ion and a Cl(-) ion are used as probes. We find that the ionic liquid provides excellent electrostatic screening at a distance of 1-2 nm. The free energy profiles show minima which are due to layering in the ionic liquid near the electrodes. This indicates that the solute ions tend to displace ionic liquid ions in the layers when approaching the electrode. The important role of non-electrostatic forces is demonstrated by the oscillations in the free energy profiles of uncharged solutes as a function of distance from the wall.     

Interaction of urea with amino acids: implications for urea-induced protein denaturation

The molecular mechanism of urea-induced protein denaturation is not yet fully understood. Mainly two opposing mechanisms are controversially discussed, according to which either hydrophobic, or polar interactions are the dominant driving force. To resolve this question, we have investigated the interactions between urea and all 20 amino acids by comprehensive molecular dynamics simulations of 22 tripeptides. Calculation of atomic contact frequencies between the amino acids and solvent molecules revealed a clear profile of solvation preferences by either water or urea. Almost all amino acids showed preference for contacts with urea molecules, whereas charged and polar amino acids were found to have slight preferences for contact with water molecules. Particularly strong preference for contacts to urea were seen for aromatic and apolar side-chains, as well as for the protein backbone of all amino acids. Further, protein-urea hydrogen bonds were found to be significantly weaker than protein-water or water-water hydrogen bonds. Our results suggest that hydrophobic interactions are the dominant driving force, while hydrogen bonds between urea and the protein backbone contribute markedly to the overall energetics by avoiding unfavorable unsatisfied hydrogen bond sites on the backbone. In summary, we suggest a combined mechanism that unifies the two current and seemingly opposing views.