Local effects in hydrophobicity. Electrostatic interactions in the restructuring of highly polar solvent around less polar solute

Two simplifying assumptions are frequently used in the biophysical chemistry of aqueous solutions: (i) a dielectric mediates the interactions of polar and ionic molecules in aqueous phases and (ii) the dielectric constant of this medium is high and uniform up to molecular surfaces. Because of their great utility in rationalizing simple electrostatic and dielectric effects in such polar systems, it is important to examine whether these assumptions also lead to deductions that are locally consistent with the solvent restructuring observed in hydrophobic phenomena. In this paper, using a model polar fluid system, these macroscopic assumptions are applied to the rigorous, microscopic nonlinear integral equation for Wki, the potential of mean force between two adjacent polar molecules. In systems of high dielectric constant, linearization of Boltzmann exponentials and approximation of three-molecule potentials of mean force by superposition of two-molecule potentials permit reduction to a linear integral equation for Wki. It is shown that the strictly local electrostatic contributions to Wki exert an effect that is qualitatively similar to the global screening effect of a dielectric medium. Through the relation between Wki and configurational probabilities, it is further found that reducing the polarity of a molecule in a polar fluid shifts local pair probability density from energetically unfavorable to energetically favorable two-molecule configurations. This general effect, which clearly promotes local structure, would augment more specific hydrophobic mechanisms in aqueous systems. Thus, the assumptions upon which the highly successful Debye-Hückel and Onsager models are supported lead also to deductions about local structure that are consistent with hydrophobic structure enhancement.     

Fifty years of solvent denaturation

The author has worked in the area of solvent denaturation and stabilization, off and on, for approximately 50 years. This paper is a personal view of the progress which has been made since 1950. The topic is limited to the development of thermodynamic molecular models for the interpretation of the unfolding and stabilization of protein structures. The story starts with the work in Kauzmann's laboratory shortly after World War II and proceeds through models for multisite binding, the linear denaturation curve, the special considerations required for understanding weak solvent exchange and a new model for the balance of solvent contact interactions and excluded volume.     

Energy patterns in twist-opening models of DNA with solvent interactions

Energy localization, via modulation instability, is addressed in a modified twist-opening model of DNA with solvent interactions. The Fourier expansion method is used to reduce the complex roto-torsional equations of the system to a set of discrete coupled nonlinear Schrödinger equations, which are used to perform the analytical investigation of modulation instability. We find that the instability criterion is highly influenced by the solvent parameters. Direct numerical simulations, performed on the generic model, further confirm our analytical predictions, as solvent interactions bring about highly localized energy patterns. These patterns are also shown to be robust under thermal fluctuations.     

Trifluoroethanol may form a solvent matrix for assisted hydrophobic interactions between peptide side chains

Several models for interactions between trifluoroethanol (TFE) and peptides and proteins have recently been proposed, but none have been able to rationalize the puzzling observations that on the one hand TFE can stabilize some hydrophobic interactions in secondary structures, but on the other can also melt the hydrophobic cores of globular proteins. The former is illustrated in this paper by the effect of TFE on a short elastin peptide, GVG(VPGVG)(3), which forms type II beta-turns stabilized by hydrophobic interactions between two intra-turn valine side chains. This folding, driven by increasing the entropy of bulk water, is stimulated in TFE-water mixtures and/or by raising the temperature. To explain these apparently contradictory observations, we propose a model in which TFE clusters locally assist the folding of secondary structures by first breaking down interfacial water molecules on the peptide and then providing a solvent matrix for further side chain--side chain interactions. This model also provides an explanation for TFE-induced transitions between secondary structures, in which the TFE clusters may redirect non-local to local interactions.     

Role of solvent properties of aqueous media in macromolecular crowding effects

Analysis of the macromolecular crowding effects in polymer solutions show that the excluded volume effect is not the only factor affecting the behavior of biomolecules in a crowded environment. The observed inconsistencies are commonly explained by the so-called soft interactions, such as electrostatic, hydrophobic, and van der Waals interactions, between the crowding agent and the protein, in addition to the hard nonspecific steric interactions. We suggest that the changes in the solvent properties of aqueous media induced by the crowding agents may be the root of these “soft” interactions. To check this hypothesis, the solvatochromic comparison method was used to determine the solvent dipolarity/polarizability, hydrogen-bond donor acidity, and hydrogen-bond acceptor basicity of aqueous solutions of different polymers (dextran, poly(ethylene glycol), Ficoll, Ucon, and polyvinylpyrrolidone) with the polymer concentration up to 40% typically used as crowding agents. Polymer-induced changes in these features were found to be polymer type and concentration specific, and, in case of polyethylene glycol (PEG), molecular mass specific. Similarly sized polymers PEG and Ucon producing different changes in the solvent properties of water in their solutions induced morphologically different α-synuclein aggregates. It is shown that the crowding effects of some polymers on protein refolding and stability reported in the literature can be quantitatively described in terms of the established solvent features of the media in these polymers solutions. These results indicate that the crowding agents do induce changes in solvent properties of aqueous media in crowded environment. Therefore, these changes should be taken into account for crowding effect analysis.     

Analysis of the impact of solvent on contacts prediction in proteins

Background  

The correlated mutations concept is based on the assumption that interacting protein residues coevolve, so that a mutation in one of the interacting counterparts is compensated by a mutation in the other. Approaches based on this concept have been widely used for protein contacts prediction since the 90s. Previously, we have shown that water-mediated interactions play an important role in protein interfaces. We have observed that current "dry" correlated mutations approaches might not properly predict certain interactions in protein interfaces due to the fact that they are water-mediated.     

Halophilic proteins and the influence of solvent on protein stabilization

Competition between protein-solvent and protein-protein interactions is arguably the most important contributing factor to polypeptide folding in general. A study of halophilic proteins, correlating their stability and solution structures in different conditions, focuses on the effects of a high salt solvent. A mechanism is proposed to explain how these proteins have adapted to such an extreme environment.     

Polyethylene glycol behaves like weak organic solvent

Effect of polyethylene glycol (PEG) on protein solubility has been primarily ascribed to its large hydrodynamic size and thereby molecular crowding effect. However, PEG also shows characteristics of organic solvents. Here, we have examined the solubility of glycine and aliphatic and aromatic amino acids in PEG solutions. PEG400, PEG4000, and PEG20000 decreased the solubility of glycine, though to a much smaller magnitude than the level achieved by typical organic solvents, including ethanol and dimethyl sulfoxide. PEG4000 showed varying degree of interactions with amino acid side chains. The free energy of aliphatic side chains marginally increased by the addition of PEG4000, indicating their weak unfavorable interactions. However, it significantly decreased the free energy of the aromatic side chains and hence stabilized them. Thus, it was concluded that PEG behaves like weak organic solvents; namely PEG destabilized (interacted unfavorably with) polar and charged groups and stabilized (interacted favorably with) aromatic groups. Interestingly, the interaction of PEG20000, but neither PEG400 nor PEG4000, with glycine resulted in phase separation under the saturated concentration of glycine.     

Modeling loop reorganization free energies of acetylcholinesterase: a comparison of explicit and implicit solvent models

The treatment of hydration effects in protein dynamics simulations varies in model complexity and spans the range from the computationally intensive microscopic evaluation to simple dielectric screening of charge-charge interactions. This paper compares different solvent models applied to the problem of estimating the free-energy difference between two loop conformations in acetylcholinesterase. Molecular dynamics (MD) simulations were used to sample potential energy surfaces of the two basins with solvent treated by means of explicit and implicit methods. Implicit solvent methods studied include the generalized Born (GB) model, atomic solvation potential (ASP), and the distance-dependent dieletric constant. By using the linear response approximation (LRA), the explicit solvent calculations determined a free-energy difference that is in excellent agreement with the experimental estimate, while rescoring the protein conformations with GB or the Poisson equation showed inconsistent and inferior results. While the approach of rescoring conformations from explicit water simulations with implicit solvent models is popular among many applications, it perturbs the energy landscape by changing the solvent contribution to microstates without conformational relaxation, thus leading to non-optimal solvation free energies. Calculations applying MD with a GB solvent model produced results of comparable accuracy as observed with LRA, yet the electrostatic free-energy terms were significantly different due to optimization on a potential energy surface favored by an implicit solvent reaction field. The simpler methods of ASP and the distance-dependent scaling of the dielectric constant both produced considerable distortions in the protein internal free-energy terms and are consequently unreliable.     

Temperature and solvent effects on reaction centers from Chloroflexus aurantiacus and Chromatium tepidum.

Temperature and solvent effects on reaction center structures were examined in two thermophilic photosynthetic bacteria, Chloroflexus aurantiacus and Chromatium tepidum, in order to gain insight into the interactions among the reaction center proteins and pigment systems. Thermal stability of the reaction centers was found to be proportional to the optimum growth temperature. Circular dichroism (CD) spectra in the 250-300 nm region indicated that thermal denaturation destroyed tertiary structures (helix-to-helix interactions or amino acid residue conformation) in the native reaction center, keeping helical structures intact. Absorption and circular dichroism spectral changes showed that alcohol denatured the so-called special pair and the accessory BChl a independently. The alcohol denaturation further indicates that the coordination between BChl a and amino acid residue in the protein is one of the important interactions maintaining the pigment organization of the reaction centers.     

Energetics of solvent and ligand-induced conformational changes in alpha-lactalbumin

The energetics of structural changes in the holo and apo forms of a-lactalbumin and the transition between their native and denatured states induced by binding Ca2+ and Na+ have been studied by differential scanning and isothermal titration microcalorimetry and circular dichroism spectroscopy under various solvent conditions. Removal of Ca2+ from the protein enhances its sensitivity to pH and ionic conditions due to noncompensated negative charge-charge interactions at the cation binding site, which significantly reduces its overall stability. At neutral pH and low ionic strength, the native structure of apo-alpha-lactalbumin is stable below 14 C and undergoes a conformational change to a native-like molten globule intermediate at temperatures above 25 degrees C. The denaturation of either holo- or apo-alpha-lactalbumin is a highly cooperative process that is characterized by an enthalpy of similar magnitude when calculated at the same temperature. Measured by direct calorimetric titration, the enthalpy of Ca2+-binding to apo-LA at pH 7.5 is -7.1 kJ mol(-1) at 5.0 degrees C. which is essentially invariant to protonation effects. This small enthalpy effect infers that stabilization of alpha-lactalbumin by Ca2+ is primarily an entropy driven process, presumably arising from electrostatic interactions and the hydration effect. In contrast to the binding of calcium, the interaction of sodium with apo-LA does not produce a noticeable heat effect and is characterized by its ionic nature rather than specific binding to the metal-binding site. Characterization of the conformational stability and ligand binding energetics of alpha-lactalbumin as a function of solvent conditions furnishes significant insight regarding the molecular flexibility and regulatory mechanism mediated by this protein.     

A study of the solvent effect on the morphology of RDX crystal by molecular modeling method

Molecular dynamics simulations have been performed to investigate the effect of acetone solvent on the crystal morphology of RDX. The results show that the growth morphology of RDX crystal in vacuum is dominated by the (111), (020), (200), (002), and (210) faces using the BFDH laws, and (111) face is morphologically the most important. The analysis of surface structures of RDX crystal indicates that (020) face is non-polar, while (210), (111), (002), and (200) faces are polar among which (210) face has the strongest polarity. The interaction between acetone solvent and each RDX crystal face is different, and the order of binding energy on these surfaces is (210)?>?(111)?>?(002)?>?(200)?>?(020). The analysis of interactions among RDX and acetone molecules reveal that the system nonbond interactions are primary strong van der Waals and electrostatic interactions containing π-hole interactions, the weak hydrogen bond interactions are also existent. The effect of acetone on the growth of RDX crystal can be evaluated by comparing the binding energies of RDX crystalline faces. It can be predicted that compared to that in vacuum, in the process of RDX crystallization from acetone, the morphological importance of (210) face is increased more and (111) face is not the most important among RDX polar surfaces, while the non-polar (020) face probably disappears. The experimentally obtained RDX morphology grown from acetone is in agreement with the theoretical prediction.     

Effect of solvent diffusion on the apomyoglobin-water interface

Few techniques can identify interactions between proteins and individual water molecules when the protein is in solution. The present work has sought to bridge the gap between the molecular level studies and the search for a physical property of the solution (bathing the proteins) that would regulate the protein hydration level. The properties of the solution were varied by adding nondenaturing solutes and solvents to the protein solutions and then studying their effect on the intrinsic fluorescence of apomyoglobin. The resolution of the tryptophan emission into the two component spectra corresponding to tryptophans W7 (accessible to the solvent) and W14 (buried in the protein matrix) has allowed us to probe two specific parts of the protein. Whereas W14 is not affected when the medium is altered, the analysis of W7 fluorescence has shown that cosolvent diffusion plays a dominant role in the mobility of water molecules near the protein surface.     

Protein stability and electrostatic interactions between solvent exposed charged side chains

To investigate the contribution to protein stability of electrostatic interactions between charged surface residues, we have studied the effect of substituting three negatively charged solvent exposed residues with their side-chain amide analogs in bovine calbindin D9k--a small (Mr 8,500) globular protein of the calmodulin superfamily. The free energy of urea-induced unfolding for the wild-type and seven mutant proteins has been measured. The mutant proteins have increased stability towards unfolding relative to the wild-type. The experimental results correlate reasonably well with theoretically calculated relative free energies of unfolding and show that electrostatic interactions between charges on the surface of a protein can have significant effects on protein stability.     

Computational protein design with a generalized Born solvent model: application to Asparaginyl-tRNA synthetase

Computational Protein Design (CPD) is a promising method for high throughput protein and ligand mutagenesis. Recently, we developed a CPD method that used a polar-hydrogen energy function for protein interactions and a Coulomb/Accessible Surface Area (CASA) model for solvent effects. We applied this method to engineer aspartyl-adenylate (AspAMP) specificity into Asparaginyl-tRNA synthetase (AsnRS), whose substrate is asparaginyl-adenylate (AsnAMP). Here, we implement a more accurate function, with an all-atom energy for protein interactions and a residue-pairwise generalized Born model for solvent effects. As a first test, we compute aminoacid affinities for several point mutants of Aspartyl-tRNA synthetase (AspRS) and Tyrosyl-tRNA synthetase and stability changes for three helical peptides and compare with experiment. As a second test, we readdress the problem of AsnRS aminoacid engineering. We compare three design criteria, which optimize the folding free-energy, the absolute AspAMP affinity, and the relative (AspAMP-AsnAMP) affinity. The sequences and conformations are improved with respect to our previous, polar-hydrogen/CASA study: For several designed complexes, the AspAMP carboxylate forms three interactions with a conserved arginine and a designed lysine, as in the active site of the AspRS:AspAMP complex. The conformations and interactions are well maintained in molecular dynamics simulations and the sequences have an inverted specificity, favoring AspAMP over AsnAMP. The method is not fully successful, since experimental measurements with the seven most promising sequences show that they do not catalyze at a detectable level the adenylation of Asp (or Asn) with ATP. This may be due to weak AspAMP binding and/or disruption of transition-state stabilization.     

Excluded volume approximation to protein-solvent interaction. The solvent contact model.

Important properties of globular proteins, such as the stability of its folded state, depend sensitively on interactions with solvent molecules. Existing methods for estimating these interactions, such as the geometrical surface model, are either physically misleading or too time consuming to be applied routinely in energy calculations. As an alternative, we derive here a simple model for the interactions between protein atoms and solvent atoms in the first hydration layer, the solvent contact model, based on the conservation of the total number of atomic contacts, a consequence of the excluded-volume effect. The model has the conceptual advantage that protein-protein contacts and protein-solvent contacts are treated in the same language and the technical advantage that the solvent term becomes a particularly simple function of interatomic distances. The model allows rapid calculation of any physical property that depends only on the number and type of protein-solvent nearest-neighbor contacts. We propose use of the method in the calculation of protein solvation energies, conformational energy calculations, and molecular dynamics simulations.     

Identification of substrate binding sites in enzymes by computational solvent mapping

Enzyme structures determined in organic solvents show that most organic molecules cluster in the active site, delineating the binding pocket. We have developed algorithms to perform solvent mapping computationally, rather than experimentally, by placing molecular probes (small molecules or functional groups) on a protein surface, and finding the regions with the most favorable binding free energy. The method then finds the consensus site that binds the highest number of different probes. The probe-protein interactions at this site are compared to the intermolecular interactions seen in the known complexes of the enzyme with various ligands (substrate analogs, products, and inhibitors). We have mapped thermolysin, for which experimental mapping results are also available, and six further enzymes that have no experimental mapping data, but whose binding sites are well characterized. With the exception of haloalkane dehalogenase, which binds very small substrates in a narrow channel, the consensus site found by the mapping is always a major subsite of the substrate-binding site. Furthermore, the probes at this location form hydrogen bonds and non-bonded interactions with the same residues that interact with the specific ligands of the enzyme. Thus, once the structure of an enzyme is known, computational solvent mapping can provide detailed and reliable information on its substrate-binding site. Calculations on ligand-bound and apo structures of enzymes show that the mapping results are not very sensitive to moderate variations in the protein coordinates.