Collective properties of hydration: long range and specificity of hydrophobic interactions

Parallelization and improvements of the generalized born model with a simple sWitching function for modern graphics processors

Two fundamental challenges of simulating biologically relevant systems are the rapid calculation of the energy of solvation and the trajectory length of a given simulation. The Generalized Born model with a Simple sWitching function (GBSW) addresses these issues by using an efficient approximation of Poisson-Boltzmann (PB) theory to calculate each solute atom's free energy of solvation, the gradient of this potential, and the subsequent forces of solvation without the need for explicit solvent molecules. This study presents a parallel refactoring of the original GBSW algorithm and its implementation on newly available, low cost graphics chips with thousands of processing cores. Depending on the system size and nonbonded force cutoffs, the new GBSW algorithm offers speed increases of between one and two orders of magnitude over previous implementations while maintaining similar levels of accuracy. We find that much of the algorithm scales linearly with an increase of system size, which makes this water model cost effective for solvating large systems. Additionally, we utilize our GPU-accelerated GBSW model to fold the model system chignolin, and in doing so we demonstrate that these speed enhancements now make accessible folding studies of peptides and potentially small proteins.

The role of solvent in protein folding and in aggregation

We discuss features of the effect of solvent on protein folding andaggregation, highlighting the physics related to the particulate nature and the peculiar structure of the aqueous solvent, and the biological significance of interactions between solvent and proteins. To this purpose we use a generalized energy landscape of extended dimensionality. A closer look at the properties of solvent induced interactions and forces proves useful for understanding the physical grounds of `ad hoc' interactions and for devising realistic ways of accounting for solvent effects. The solvent has long been known to be a crucially important part of biological systems, and times appear mature for it to be adequately accounted for in the protein folding problem. Use of the extended dimensionality energy landscape helpseliciting the possibility of coupling among conformational changes and aggregation, such as proved by experimental data in the literature.

Biochemical and biophysical analysis of five disease-associated human adenylosuccinate lyase mutants

Adenylosuccinate lyase (ASL), a catalyst of key reactions in purine biosynthesis, is normally a homotetramer in which three subunits contribute to each of four active sites. Human ASL deficiency is an inherited metabolic disease associated with autism and mental retardation. We have characterized five disease-associated ASL mutants: R194C and K246E are located at subunit interfaces, L311V is in the central helical region away from the active site, and R396C and R396H are at the entrance to the active site. The V(max) (at 25 degrees C) for R194C is comparable to that of WT, while those of L311V, R396C, R396H, and K246E are considerably reduced and affinity for adenylosuccinate is retained. The mutant enzymes have decreased positive cooperativity as compared to WT. K246E exists mainly as dimer or monomer, accounting for its negligible activity, whereas the other mutant enzymes are similar to WT in the predominance of tetramer. At 37 degrees C, the specific activity of WT and these mutant enzymes slowly decreases 30-40% with time and reaches a limiting specific activity without changing significantly the amount of tetramer. Mutant R194C is unique in being rapidly inactivated at the harsher temperature of 60 degrees C, indicating that it is the least stable enzyme in vitro. Conformational changes in the mutant enzymes are evident from protein fluorescence intensity at 25 degrees C and after incubation at 37 degrees C, which correlates with the loss of enzymatic activity. Thus, these disease-associated single mutations can yield enzyme with reduced activity either by affecting the active site or by perturbing the enzyme's structure and/or native conformation which are required for catalytic function.

Solvent-induced free energy landscape and solute-solvent dynamic coupling in a multielement solute

Molecular dynamics simulations using a simple multielement model solute with internal degrees of freedom and accounting for solvent-induced interactions to all orders in explicit water are reported. The potential energy landscape of the solute is flat in vacuo. However, the sole untruncated solvent-induced interactions between apolar (hydrophobic) and charged elements generate a rich landscape of potential of mean force exhibiting typical features of protein landscapes. Despite the simplicity of our solute, the depth of minima in this landscape is not far in size from free energies that stabilize protein conformations. Dynamical coupling between configurational switching of the system and hydration reconfiguration is also elicited. Switching is seen to occur on a time scale two orders of magnitude longer than that of the reconfiguration time of the solute taken alone, or that of the unperturbed solvent. Qualitatively, these results are unaffected by a different choice of the water-water interaction potential. They show that already at an elementary level, solvent-induced interactions alone, when fully accounted for, can be responsible for configurational and dynamical features essential to protein folding and function.

Protein aggregation/crystallization and minor structural changes: universal versus specific aspects

Protein association covers wide interests in biophysics, protein science, and biotechnologies, and it is often viewed as governed by conformation details. More recently, the existence of a universal physical principle governing aggregation/crystallization processes has been suggested by a series of experiments and shown to be linked to the universal scaling properties of concentration fluctuations occurring in the proximity of a phase transition (spinodal demixing in the specific case). Such properties have provided a quantitative basis for capturing kinetic association data on a universal master curve, ruled by the normalized distance of the state of the system from its instability region. Here we report new data on lysozyme crystal nucleation. They strengthen the evidence in favor of universality and show that the system enters the region of universal behavior in a stepwise manner as a result of minor conformation changes. Results also show that the link between conformation details and universal behavior is actuated by interactions mediated by the solvent. Outside the region of universal behavior, nucleation rates become unpredictable and undetectably long.

Pressure and temperature dependence of hydrophobic hydration: Volumetric, compressibility, and thermodynamic signatures

The combined effect of pressure and temperature on hydrophobic hydration of a nonpolar methanelike solute is investigated by extensive simulations in the TIP4P model of water. Using test-particle insertion techniques, free energies of hydration under a range of pressures from 1 to 3000 atm are computed at eight temperatures ranging from 278.15 to 368.15 K. Corresponding enthalpy, entropy, and heat capacity accompanying the hydration process are estimated from the temperature dependence of the free energies. Partial molar and excess volumes calculated using pressure derivatives of the simulated free energies are consistent with those determined by direct volume simulations; but direct volume determination offers more reliable estimates for compressibility. At 298.15 K, partial molar and excess isothermal compressibilities of methane are negative at 1 atm. Partial molar and excess adiabatic (isentropic) compressibilities are estimated to be also negative under the same conditions. But partial molar and excess isothermal compressibilities are positive at high pressures, with a crossover from negative to positive compressibility at approximately 100-1000 atm. This trend is consistent with experiments on aliphatic amino acids and pressure-unfolded states of proteins. For the range of pressures simulated, hydration heat capacity exhibits little pressure dependence, also in apparent agreement with experiment. When pressure is raised at constant room temperature, hydration free energy increases while its entropic component remains essentially constant. Thus, the increasing unfavorability of hydration under raised pressure is seen as largely an enthalpic effect. Ramifications of the findings of the authors for biopolymer conformational transitions are discussed.

Effect of T-R conformational change on sickle-cell hemoglobin interactions and aggregation

We compare the role of a conformational switch and that of a point mutation in the thermodynamic stability of a protein solution and in the consequent propensity toward aggregation. We study sickle-cell hemoglobin (HbS), the beta6 Glu-Val point mutant of adult human hemoglobin (HbA), in its R (CO-liganded) conformation, and compare its aggregation properties to those of both HbS and HbA in their T (unliganded) conformation. Static and dynamic light scattering measurements performed for various hemoglobin concentrations showed critical divergences with mean field exponents as temperature was increased. This allowed determining spinodal data points T(S)(c) by extrapolation. These points were fitted to theoretical expressions of the T(S)(c) spinodal line, which delimits the region where the homogeneous solution becomes thermodynamically unstable against demixing in two sets of denser and dilute mesoscopic domains, while remaining still liquid. Fitting provided model-free numerical values of enthalpy and entropy parameters measuring the stability of solutions against demixing, namely, 93.2 kJ/mol and 314 J/ degrees K-mol, respectively. Aggregation was observed also for R-HbS, but in amorphous form and above physiological temperatures close to the spinodal, consistent with the role played in nucleation by anomalous fluctuations governed by the parameter epsilon = (T - T(S))/T(S). Fourier transform infrared (FTIR) and optical spectroscopy showed that aggregation is neither preceded nor followed by denaturation. Transient multiple interprotein contacts occur in the denser liquid domains for R-HbS, T-HbS, and T-HbA. The distinct effects of their specific nature and configurations, and those of desolvation on the demixing and aggregation thermodynamics, and on the aggregate structure are highlighted.

Internal dynamics and protein–matrix coupling in trehalose-coated proteins

We review recent studies on the role played by non-liquid, water-containing matrices on the dynamics and structure of embedded proteins. Two proteins were studied, in water-trehalose matrices: a water-soluble protein (carboxy derivative of horse heart myoglobin) and a membrane protein (reaction centre from Rhodobacter sphaeroides). Several experimental techniques were used: Mossbauer spectroscopy, elastic neutron scattering, FTIR spectroscopy, CO recombination after flash photolysis in carboxy-myoglobin, kinetic optical absorption spectroscopy following pulsed and continuous photoexcitation in Q(B) containing or Q(B) deprived reaction centre from R. sphaeroides. Experimental results, together with the outcome of molecular dynamics simulations, concurred to give a picture of how water-containing matrices control the internal dynamics of the embedded proteins. This occurs, in particular, via the formation of hydrogen bond networks that anchor the protein surface to the surrounding matrix, whose stiffness increases by lowering the sample water content. In the conclusion section, we also briefly speculate on how the protein-matrix interactions observed in our samples may shed light on the protein-solvent coupling also in liquid aqueous solutions.

Temperature Dependence of Three-Body Hydrophobic Interactions: Potential of Mean Force, Enthalpy, Entropy, Heat Capacity, and Nonadditivity

Temperature-dependent three-body hydrophobic interactions are investigated by extensive constant-pressure simulations of methane-like nonpolar solutes in TIP4P model water at six temperatures. A multiple-body hydrophobic interaction is considered to be (i) additive, (ii) cooperative, or (iii) anti-cooperative if its potential of mean force (PMF) is (i) equal to, (ii) smaller than, or (iii) larger than the corresponding pairwise sum of two-methane PMFs. We found that three-methane hydrophobic interactions at the desolvation barrier are anti-cooperative at low to intermediate T, and vary from essentially additive to slightly cooperative at high T. Interactions at the contact minimum are slightly anti-cooperative over a wider temperature range. Enthalpy, entropy, and heat capacity are estimated from the computed PMFs. Contrary to the common expectation that burial of solvent-accessible nonpolar surface area always leads to a decrease in heat capacity, the present results show that the change in heat capacity upon three-methane association is significantly positive at the desolvation barrier and slightly positive at the contact minimum. This suggests that the heat capacity signature of a hydrophobic polymer need not vary uniformly nor monotonically with conformational compactness. Ramifications for protein folding are discussed.

Dielectric modulation of biological water

We show that water constrained by vicinal hydrophobes undergoes a librational dynamics that lowers the dielectric susceptibility and induces a "redshift" of the relaxation frequency in the hydration shell. The results shed light on the way proteins enhance their intramolecular interactions as they fold or associate.

Time scale of protein aggregation dictated by liquid-liquid demixing

The growing impact of protein aggregation pathologies, together with the current high need for extensive information on protein structures are focusing much interest on the physics underlying the nucleation and growth of protein aggregates and crystals. Sickle Cell Hemoglobin (HbS), a point-mutant form of normal human Hemoglobin (HbA), is the first recognized and best-studied case of pathologically aggregating protein. Here we reanalyze kinetic data on nucleation of deoxy-HbS aggregates by referring them to the (concentration-dependent) temperature T(s) characterizing the occurrence of the phase transition of liquid-liquid demixing (LLD) of the solution. In this way, and by appropriate scaling of kinetic data at different concentrations, so as to normalize their spans, the apparently disparate sets of data are seen to fall on a master curve. Expressing the master curve vs. the parameter epsilon = (T - T(s)) / T(s), familiar from phase transition theory, allows eliciting the role of anomalously large concentration fluctuations associated with the LLD phase transition and also allows decoupling quantitatively the role of such fluctuations from that of microscopic, inter-protein interactions leading to nucleation. Referring to epsilon shows how in a narrow temperature span, that is at T - T(s), nucleation kinetics can undergo orders-of-magnitude changes, unexpected in terms of ordinary chemical kinetics. The same is true for similarly small changes of other parameters (pH, salts, precipitants), capable of altering T(s) and consequently epsilon. This offers the rationale for understanding how apparently minor changes of parameters can dramatically affect protein aggregation and related diseases.

Anti-cooperativity and cooperativity in hydrophobic interactions: Three-body free energy landscapes and comparison with implicit-solvent potential functions for proteins

Potentials of mean force (PMFs) of three-body hydrophobic association are investigated to gain insight into similar processes in protein folding. Free energy landscapes obtained from explicit simulations of three methanes in water are compared with that predicted by popular implicit-solvent effective potentials for the study of proteins. Explicit-water simulations show that for an extended range of three-methane configurations, hydrophobic association at 25 degrees C under atmospheric pressure is mostly anti-cooperative, that is, less favorable than if the interaction free energies were pairwise additive. Effects of free energy nonadditivity on the kinetic path of association and the temperature dependence of additivity are explored by using a three-methane system and simplified chain models. The prevalence of anti-cooperativity under ambient conditions suggests that driving forces other than hydrophobicity also play critical roles in protein thermodynamic cooperativity. We evaluate the effectiveness of several implicit-solvent potentials in mimicking explicit water simulated three-body PMFs. The favorability of the contact free energy minimum is found to be drastically overestimated by solvent accessible surface area (SASA). Both the SASA and a volume-based Gaussian solvent exclusion model fail to predict the desolvation barrier. However, this barrier is qualitatively captured by the molecular surface area model and a recent "hydrophobic force field." None of the implicit-solvent models tested are accurate for the entire range of three-methane configurations and several other thermodynamic signatures considered.

Interaction of processes on different length scales in a bioelastomer capable of performing energy conversion

This work concerns the aggregation properties of (Gly-Val-Gly-Val-Pro)(251) rec, a polypentapeptide reflecting a highly conserved repetitive unit of the bioelastomer, elastin. On raising the temperature of aqueous solutions above 25 degrees C, this polypeptide was already known to undergo concurrent conformational changes (hydrophobic folding), phase separation, and self-assembly with formation of aggregated three-stranded filaments composed of dynamic polypeptide helices, called beta-spirals. Aggregates obtained from the solution can be shaped into bands that acquire entropic elastic properties upon gamma-irradiation and can perform a variety of energy conversions. Previous studies have shown that aggregation is prompted by the (diverging) critical fluctuations of concentration occurring in the solution, in vicinity of its spinodal line. Here, we present combined circular dicroism (CD) and light scattering experiments, and independent fittings of experimental data to the theoretical spinodal and binodal (coexistence) lines. Results show the following logical and causal sequence of processes: (a) Smooth and progressive conformational changes promoted by concentration fluctuations occurring as temperature is raised "pull down" (in the temperature scale) the instability region of the solution. (b) This further promotes critical fluctuations. (c) The related locally high concentration prompts a further substantial conformational change ending in triple-helix formation and coacervation. (d) This intertwining of processes, covering different length scales (from that of individual peptides to the mesoscopic one of demixed regions), is related to the fact that solvent-induced interactions play a strong role over the entire scale span. These results concur with other recent ones in pointing out that process interactions over many length-scales probably reflect a frequent if not ubiquitous pattern in protein aggregation. This may be highly relevant to the desirable deep understanding of such phenomenon, whose interests cover many fields.

Effects of electric charges on hydrophobic forces. II

We study by molecular-dynamics simulations the effect of electric charges of either sign on hydrophobic interactions and on the dynamics of hydration water, using explicit water and very simplified solutes. Results show that the presence of a charged solute can disrupt the "hydrophobic contact bond" between two apolar solutes nearby, by forcing them towards a different configuration. As a consequence of different structural changes of the solvent caused by charges of opposite sign, the effect is markedly charge-sign-dependent. Analogous weaker effects appear to be induced by the presence of one additional apolar element. The dynamics of hydration water around each solute is also seen to be strongly influenced by the presence of other (charged or uncharged) nearby solutes. Comparison between our results on hydration water dynamics around charged solutes and available experimental data allows sorting out the effects of solute charge sign and size. Our results also offer a plain interpretation of the equivalence of the effects on water structure due to solute ions and to high pressures. These results reflect at a basic paradigmatic level the immensely more complex cases of well-known phenomena such as salting-in and salting-out, and of protein conformational changes caused, e.g., by the arrival of a charged or of an apolar group (phosphorilation or methylation). As it will be discussed, they help in the direction of Delbruck's desirable "progress towards a radical physical explanation" for this class of phenomena.

A complex clathrate hydrate structure showing bimodal guest hydration

Interactions between hydrophobic groups in water, as well as biomolecular hydration more generally, are intimately connected to the structure of liquid water around hydrophobic solutes. Such considerations have focused interest on clathrate hydrates: crystals in which a hydrogen-bonded network of water molecules encages hydrophobic guest molecules with which the water interacts only by non-directional van der Waals forces. Three structural families of clathrate hydrates have hitherto been recognized: cubic structure I (2M(S)-6M(L) x 46H2O), cubic structure II (16M(S) x 8M(L)-136H2O) and hexagonal structure H (M(L) x 3M(S) x 2M(S) x 34H2O) hydrates (here M(L) and M(S) are the hydrophobic guest sites associated with large and small cavities, respectively). Here we report a new hydrate structure: 1.67 choline hydroxide-tetra-n-propylammonium fluoride x 30.33H2O. This structure has a number of unusual features; in particular the choline guest exhibits both hydrophobic and hydrophilic modes of hydration. Formally the structure consists of alternating stacks of structure H and structure II hydrates, and might conceivably be found in those settings (such as seafloor deposits over natural-gas fields) in which clathrate hydrates form naturally.

The importance of hydration for the kinetics and thermodynamics of protein folding: simplified lattice models

BACKGROUND

Recent studies have proposed various sources for the origin of cooperativity in simplified protein folding models. Important contributions to cooperativity that have been discussed include backbone hydrogen bonding, sidechain packing and hydrophobic interactions. Related work has also focused on which interactions are responsible for making the free energy of the native structure a pronounced global minimum in the free energy landscape. In addition, two-flavor bead models have been found to exhibit poor folding cooperativity and often lack unique native structures. We propose a simple multibody description of hydration with expectations that it might modify the free energy surface in such a way as to increase the cooperativity of folding and improve the performance of two-flavor models.

RESULTS

We study the thermodynamics and kinetics of folding for designed 36-mer sequences on a cubic lattice using both our solvation model and the corresponding model without solvation terms. Degeneracies of the native states are studied by enumerating the maximally compact states. The histogram Monte Carlo method is used to obtain folding temperatures, densities of states and heat capacity curves. Folding kinetics are examined by accumulating mean first-passage times versus temperature. Sequences in the proposed solvation model are found to have more unique ground states, fold faster and fold with more cooperativity than sequences in the nonsolvation model.

CONCLUSIONS

We find that the addition of a multibody description of solvation can improve the poor performance of two-flavor lattice models and provide an additional source for more cooperative folding. Our results suggest that a better description of solvation will be important for future theoretical protein folding studies.