Physics and biophysics of solvent induced forces: hydrophobic interactions and context-dependent hydration

Solvent induced forces (SIFs) among solutes derive from solvent structural modification due to solutes, and consequent thermodynamic drive towards minimization of related free energy costs. The role of SIFs in biomolecular conformation and function is appreciated by observing that typical SIF values fall within the 20–200 pN interval, and that proteins are stable by only a few kcal mol^–1 (1 kcal mol^–1 corresponds to 70 pN Å). Here we study SIFs, in systems of increasing complexity, using Molecular Dynamics (MD) simulations which give time- and space-resolved details on the biologically significant scale of single protein residues and sidechains. Of particular biological relevance among our results are a strong modulability of hydrophobic SIFs by electric charges and the dependence of this modulability upon charge sign. More generally, the present results extend our understanding of the recently reported strong context-dependence of SIFs and the related potential of mean force (PMF). This context-dependence can be strong enough to propagate (by relay action) along a composite solute, and to reverse SIFs acting on a given element, relative to expectations based on its specific character (hydrophobic/ philic, charged). High specificity such as that of SIFs highlighted by the present results is of course central to biological function. Biological implications of the present results cover issues such as biomolecular functional interactions and folding (including chaperoning and pathological conformational changes), coagulation, molecular recognition, effects of phosphorylation and more.     

Influence of medium and long range interactions in protein folding.

Protein structures are stabilized by both local and long range interactions. In this work, we analyze the residue-residue contacts and the role of medium- and long-range interactions in globular proteins belonging to different structural classes. The results show that while medium range interactions predominate in all-alpha class proteins, long-range interactions predominate in all-beta class. Based on this, we analyze the performance of several structure prediction methods in different structural classes of globular proteins and found that all the methods predict the secondary structures of all-alpha proteins more accurately than other classes. Also, we observed that the residues occurring in the range of 21-30 residues apart contributes more towards long-range contacts and about 85% of residues are involved in long-range contacts. Further, the preference of residue pairs to the folding and stability of globular proteins is discussed.     

Species-specific long range interactions between receptor/ligand pairs.

Total internal reflection microscopy (TIRM) monitors Brownian fluctuations in elevation as small as 1 nm by measuring the scattering of a single sphere illuminated by an evanescent wave when the sphere is levitated by colloidal forces such as electrostatic double-layer repulsion. From the Boltzmann distribution of elevations sampled by the sphere over time, the potential energy profile can be determined with a resolution of approximately 0.1 of the thermal energy kT. Thus, the interaction between a receptor-coated (goat, horse, or rabbit immunoglobulin G (IgG)) latex sphere and a protein A (SpA)-coated glass microscope slide was studied. A typical TIRM potential energy profile measured between a bare sphere and a bare glass plate, where the sphere fluctuates around the secondary potential energy minimum formed between double-layer repulsion and gravitational attraction, agrees well with DLVO theory. The interactions measured between IgG-coated spheres and SpA-coated slides, on the other hand, displayed a weaker repulsion compared with that observed between bare surfaces under the same conditions. Analysis of the results obtained between the coated surfaces suggests an additional attractive force. The decay length of this attraction correlates with the known dissociation constants for the binding of IgG with SpA in free solution.     

Model of membrane fusion: Continuous transition to fusion pore with regard of hydrophobic and hydration interactions

We consider the process of fusion of lipid membranes from the stage of stalk with minimal radius to the stage of fusion pore. We assume that stalk directly developed into the fusion pore, omitting the stage of hemifusion diaphragm. Energy of intermediate stages is calculated on the basis of the classical elasticity theory of liquid crystals adapted for lipid membranes. The trajectory of transition from stalk to pore is obtained with regard to hydrophobic and hydration interactions. Continuous change of orientation of lipids in distal monolayers occurs along the trajectory. The orientation changes from the direction along rotational axis of the system specific to stalk to the direction corresponding to the fusion pore. Dependence of energy of intermediate stages on the value of spontaneous curvature of distal monolayers of the fusing membranes is obtained. We demonstrate that the energy barrier of the stalk-to-pore transition decreases when distal monolayers have positive spontaneous curvature, which is in accordance with available experimental data.     

An alanine-zipper structure determined by long range intermolecular interactions

A major challenge in protein folding is to identify and quantify specific structural determinants that allow native proteins to acquire their unique folded structures. Here we report the engineering of a 52-residue protein (Ala-14) that contains exclusively alanine residues at the hydrophobic a and d positions of a natural heptad-repeat sequence. Ala-14 is unfolded under normal solution conditions yet forms a parallel three-stranded alpha-helical coiled coil in crystals. Ala-14 trimers in the solid state associate with each other through the pairing of polar side chains and formation of an extended network of water-mediated hydrogen bonds. In contrast to the classical view that local intramolecular tertiary interactions dictate the three-dimensional structure of small single-domain proteins, Ala-14 shows that long range intermolecular interactions can be essential in determining the metastable alanine-zipper structure. A similar interplay between short range local and longer range global forces may underlie the conformational properties of the growing class of natively unstructured proteins in biological processes.     

Glycan specificity of a testis-specific lectin chaperone calmegin and effects of hydrophobic interactions

Background

Testis-specific chaperone calmegin is required for the generation of normal spermatozoa. Calmegin is known to be a homologue of endoplasmic reticulum (ER) residing lectin chaperone calnexin. Although functional similarity between calnexin and calmegin has been predicted, detailed information concerned with substrate recognition by calmegin, such as glycan specificity, chaperone function and binding affinity, are obscure.

Methods

In this study, biochemical properties of calmegin and calnexin were compared using synthetic glycans and glycosylated or non-glycosylated proteins as substrates.

Results

Whereas their amino acid sequences are quite similar to each other, a certain difference in secondary structures was indicated by circular dichroism (CD) spectrum. While both of them inhibited protein heat-aggregation to a similar extent, calnexin exhibited a higher ability to facilitate protein folding. Similarly to calnexin, calmegin preferentially recognizes monoglucosylated glycans such as Glc₁Man₉GlcNAc₂ (G1M9). While the surface hydrophobicity of calmegin was higher than that of calnexin, calnexin showed stronger binding to substrate. We reasoned that lectin activity, in addition to hydrophobic interaction, contributes to this strong affinity between calnexin and substrate.

Conclusions

Although their similarity in carbohydrate binding specificities is high, there seems to be some differences in the mode of substrate recognition between calmegin and calnexin.

General significance

Properties of calmegin as a lectin-chaperone were revealed in comparison with calnexin.     

CAP and RNA polymerase interactions with the lac promoter: binding stoichiometry and long range effects.

The binding stoichiometries of the complexes formed when the E. coli cyclic AMP receptor protein (CAP) binds to 203 bp lac promoter-operator restriction fragments have been determined. Under quantitative binding conditions, a single dimer of CAP occupies each of two sites in the promoter. Different electrophoretic mobilities are observed for 1:1 complexes formed with L8-UV5 mutant, L305 mutant, and wild type promoter fragments, indicating sequence-specific structural differences between the complexes. The differences in gel mobility between L8-UV5 and wild type complexes disappear when the promoter fragments are cleaved with Hpa II restriction endonuclease. Models in which CAP alters DNA conformation or in which CAP forms a transient intramolecular bridge between two domains of a DNA molecule could account for these observations. The selective binding of RNA polymerase to CAP-promoter complexes is demonstrated: the binding of a single CAP dimer to the promoter is sufficient to stimulate subsequent polymerase binding. Functional CAP molecules are not released from the promoter on polymerase binding.     

Formins regulate actin filament flexibility through long range allosteric interactions

The members of the formin family nucleate actin polymerization and play essential roles in the regulation of the actin cytoskeleton during a wide range of cellular and developmental processes. In the present work, we describe the effects of mDia1-FH2 on the conformation of actin filaments by using a temperature-dependent fluorescence resonance energy transfer method. Our results revealed that actin filaments were more flexible in the presence than in the absence of formin. The effect strongly depends on the mDia1-FH2 concentration in a way that indicates that more than one mechanism is responsible for the formin effect. In accordance with the more flexible filament structure, the thermal stability of actin decreased and the rate of phosphate dissociation from actin filaments increased in the presence of formin. The interpretation of the results supports a model in which formin binding to barbed ends makes filaments more flexible through long range allosteric interactions, whereas binding of formin to the sides of the filaments stabilizes the protomer-protomer interactions. These results suggest that formins can regulate the conformation of actin filaments and may thus also modulate the affinity of actin-binding proteins to filaments nucleated/capped by formins.     

Atomic-scale analysis of the solvation thermodynamics of hydrophobic hydration.

Molecular dynamics simulations are used to model the transfer thermodynamics of krypton from the gas phase into water. Extra long, nanosecond simulations are required to reduce the statistical uncertainty of the calculated "solvation" enthalpy to an acceptable level. Thermodynamic integration is used to calculate the "solvation" free energy, which together with the enthalpy is used to calculate the "solvation" entropy. A comparison series of simulations are conducted using a single Lennard-Jones sphere model of water to identify the contribution of hydrogen bonding to the thermodynamic quantities. In contrast to the classical "iceberg" model of hydrophobic hydration, the favorable enthalpy change for the transfer process at room temperature is found to be due primarily to the strong van der Waals interaction between the solute and solvent. Although some stabilization of hydrogen bonding does occur in the solvation shell, this is overshadowed by a destabilization due to packing constraints. Similarly, whereas some of the unfavorable change in entropy is attributed to the reduced rotational motion of the solvation shell waters, the major component is due to a decrease in the number of positional arrangements associated with the translational motions.     

Differences in hydration structure near hydrophobic and hydrophilic amino acids.

We use molecular dynamics to simulate recent neutron scattering experiments on aqueous solutions of N-acetyl-leucine-amide and N-acetyl-glutamine-amide, and break down the total scattering function into contributions from solute-solute, solute-water, water-water, and intramolecular correlations. We show that the shift of the main diffraction peak to smaller angle that is observed for leucine, but not for glutamine, is attributable primarily to alterations in water-water correlations relative to bulk. The perturbation of the water hydrogen-bonded network extends roughly two solvation layers from the hydrophobic side chain surface, and is characterized by a distribution of hydrogen bonded ring sizes that are more planar and are dominated by pentagons in particular than those near the hydrophilic side chain. The different structural organization of water near the hydrophobic solute that gives rise to the inward shift in the main neutron diffraction peak under ambient conditions may also provide insight into the same directional shift for pure liquid water as it is cooled and supercooled.     

Direct measurement of the intermolecular forces between counterion-condensed DNA double helices. Evidence for long range attractive hydration forces.

Rather than acting by modifying van der Waals or electrostatic double layer interactions or by directly bridging neighboring molecules, polyvalent ligands bound to DNA double helices appear to act by reconfiguring the water between macromolecular surfaces to create attractive long range hydration forces. We have reached this conclusion by directly measuring the repulsive forces between parallel B-form DNA double helices pushed together from the separations at which they have self organized into hexagonal arrays of parallel rods. For all of the wide variety of condensing agents from divalent Mn to polymeric protamines, the resulting intermolecular force varies exponentially with a decay rate of 1.4-1.5 A, exactly one-half that seen previously for hydration repulsion. Such behavior qualitatively contradicts the predictions of all electrostatic double layer and van der Waals force potentials previously suggested. It fits remarkably well with the idea, developed and tested here, that multivalent counterion adsorption reorganizes the water at discrete sites complementary to unadsorbed sites on the apposing surface. The measured strength and range of these attractive forces together with their apparent specificity suggest the presence of a previously unexpected force in molecular organization.     

A ribonucleoprotein fragment of the 30 S ribosome of E. coli: evidence for long range RNA-RNA interactions.

A stable homogeneous ribonucleoprotein fragment of the 30 S ribosomal subunit of E. coli has been prepared by mild nuclease digestion and heating in a constant ionic environment. The fragment contains about half of the 16 S ribosomal RNa and six proteins: S4, S7, S9, S13, S16 and S19. The RNA moiety contains the reported binding sites of all six proteins. After deproteinization, 80% of the RNA migrated as two major electrophoretic bands, which were isolated and sequenced. Each band contained sequences from the 5' and 3' thirds of the 16 S RNA but none from the central third. That these two noncontiguous RNA domains migrated together electrophoretically in Mg++-containing gels after deproteinization constitutes direct evidence that the 16 S RNA is folded in the intact ribosome so as to bring the two domains close together and that there are RNA-RNA interactions between them in the presence of Mg++.     

The effect of vicinal polar and charged groups on hydrophobic hydration.

The use of a linear relationship between free energy of hydrophobic hydration and solvent-accessible apolar surface area has been helpful in interpreting the thermodynamics of biological macromolecules. However, a recent study (Y.-K. Cheng, P. J. Rossky, Nature 1998, Vol. 392, pp. 696-699) has established a substantial enthalpic dependence on biomolecular surface topography, originating from solvent hydrogen-bonding loss in a restrictive geometry. In this study, we use molecular dynamics simulations of 2-Zn insulin in water solvent to explore the further effect of vicinal polar or charged groups on hydrophobic hydration at a biomolecular surface. In contrast to the case for solvent more isolated from such polar solute influences, the binding energies of the water that is proximal to the hydrophobic dimeric interface of insulin and vicinal to polar and charged groups are comparable to the bulk solvent value, a result of compensating interaction primarily with the solute counterions. The results suggest a special importance for such polar/charged groups in biological processes involving hydrophobic surface regions of restricted geometry and also suggest a general route for tuning the hydrophobicity of interfaces.     

Lipids and topological rules of membrane protein assembly: balance between long and short range lipid-protein interactions

The N-terminal six-transmembrane domain (TM) bundle of lactose permease of Escherichia coli is uniformly inverted when assembled in membranes lacking phosphatidylethanolamine (PE). Inversion is dependent on the net charge of cytoplasmically exposed protein domains containing positive and negative residues, net charge of the membrane surface, and low hydrophobicity of TM VII acting as a molecular hinge between the two halves of lactose permease (Bogdanov, M., Xie, J., Heacock, P., and Dowhan, W. (2008) J. Cell Biol. 182, 925-935). Net neutral lipids suppress the membrane translocation potential of negatively charged amino acids, thus increasing the cytoplasmic retention potential of positively charged amino acids. Herein, TM organization of sucrose permease (CscB) and phenylalanine permease (PheP) as a function of membrane lipid composition was investigated to extend these principles to other proteins. For CscB, topological dependence on PE only becomes evident after a significant increase in the net negative charge of the cytoplasmic surface of the N-terminal TM bundle. High negative charge is required to overcome the thermodynamic block to inversion due to the high hydrophobicity of TM VII. Increasing the positive charge of the cytoplasmic surface of the N-terminal TM hairpin of PheP, which is misoriented in PE-lacking cells, favors native orientation in the absence of PE. PheP and CscB also display co-existing dual topologies dependent on changes in the charge balance between protein domains and the membrane lipids. Therefore, the topology of both permeases is dependent on PE. However, CscB topology is governed by thermodynamic balance between opposing lipid-dependent electrostatic and hydrophobic interactions.     

Voltage-dependent hydration and conduction properties of the hydrophobic pore of the mechanosensitive channel of small conductance

A detailed picture of water and ion properties in small pores is important for understanding the behavior of biological ion channels. Several recent modeling studies have shown that small, hydrophobic pores exclude water and ions even if they are physically large enough to accommodate them, a mechanism called hydrophobic gating. This mechanism has been implicated in the gating of several channels, including the mechanosensitive channel of small conductance (MscS). Although the pore in the crystal structure of MscS is wide and was initially hypothesized to be open, it is lined by hydrophobic residues and may represent a nonconducting state. Molecular dynamics simulations were performed on MscS to determine whether or not the structure can conduct ions. Unlike previous simulations of hydrophobic nanopores, electric fields were applied to this system to model the transmembrane potential, which proved to be important. Although simulations without a potential resulted in a dehydrated, occluded pore, the application of a potential increased the hydration of the pore and resulted in current flow through the channel. The calculated channel conductance was in good agreement with experiment. Therefore, it is likely that the MscS crystal structure is closer to a conducting than a nonconducting state.     

Protein stability and hydrophobic interactions

The paper summarizes results of calorimetric studies of protein denaturation and of dissolution of non-polar substances in water. The analysis of the available experimental data shows that the positive contribution of the hydrophobic interactions in stabilization of the protein compact state is due to van der Waals interactions between the protein non-polar groups, while the contribution of water solvation by these groups, in spite of the widely spread opinion, appears to be always negative. This destabilizing action of water solvation on the protein increases as the temperature decreases, and at a significantly low temperature causes unfolding of the compact structure of protein, i. e. cold denaturation.     

Influence of hydrophobic interactions on the structure and steroid-binding properties of glucocorticoid receptors

Glucocorticoid-receptor complexes in rat thymus cytosol were characterized by gel-filtration and ion-exchange chromatography and by other procedures. Two forms of non-transformed complex were identified at low ionic strength in the presence of molybdate, with Stokes radii of approx. 8 and 6 nm. The 8 nm molybdate-stabilized form could be converted to the 6 nm form by chromatography on Sephacryl S-300 or Lipidex 1000 or by incubation with charcoal or phospholipase C, but not by chromatography on Sephadex G-25. The dissociation rate of the complex was reduced by treatment with charcoal or Lipidex 1000, but none of the treatments caused transformation to a DNA-binding form. Transformation of the complex, by exposure to elevated temperature or ionic strength in the absence of molybdate, resulted in the appearance of a different 6 nm form, distinguished by an increased affinity for DNA-cellulose and a reduced affinity for DEAE-cellulose. These results suggest that receptor transformation is preceded by structural changes associated with the loss of a lipid factor from the complex. Non-polar steroid antagonists, and lipophilic compounds such as phenothiazines, were found to bind to secondary, hydrophobic sites on the receptor and to exert allosteric effects on the primary steroid-binding site; these and other observations emphasize the importance of hydrophobic interactions as determinants of the structure and properties of glucocorticoid receptors.     

Entropy convergence in hydrophobic hydration: a scaled particle theory analysis

The occurrence of entropy convergence in hydrophobic hydration is verified from available experimental thermodynamic data for both noble gases and hydrocarbons. The entropy convergence phenomenon can be reproduced by means of the scaled particle theory, provided that a temperature-dependent hard sphere diameter is used for water molecules. The calculated work of cavity creation shows a non-monotonic temperature dependence with a flat maximum slightly above 100 degrees C, irrespective of the cavity size. The corresponding cavity entropy changes converge approximately 100 degrees C, in qualitative agreement with the experimental finding.     

Volume and hydration changes of DNA-ligand interactions

We report the volumetric and other thermodynamic properties of ethidium bromide (EB), propidium iodide (PI) and daunomycin (DAU) intercalating with poly(dA).poly(dT), poly[d(A-T)].poly[d(A-T)], and poly[d(G-C)].poly[d(G-C)], respectively, as well as minor groove binder Hoechst 33258 binding with poly[d(A-T)].poly[d(A-T)]. The data were obtained using fluorescence titration and hydrostatic pressure measurements. Our thermodynamic data are combined with enthalpies from literature reports to analyze the thermodynamic characteristics of the different interactions. The differences are interpreted based on three processes related to hydration: I. burial of non-polar hydrophobic solvent accessible surface, II. burial of polar surface and formation of solute-solute H-bonds, and III. disruption of "structural" hydration. Sequence dependent conformational changes may also be important when comparing ligand binding to different DNA sequences. We conclude that a combination of different thermodynamic parameters, especially volume change, is essential in order to understand the role of hydration in the energetics of DNA-ligand interactions.