Correlated solvent-induced forces on a protein at single residue resolution: relation to conformation, stability, dynamics and 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.

Liquid-structure forces and electrostatic modulation of biomolecular interactions in solution

Molecular interactions in solution are controlled by the bulk medium and by the forces originating in the structured region of the solvent close to the solutes. In this paper, a model of electrostatic and liquid-structure forces for dynamics simulations of biomolecules is presented. The model introduces information on the microscopic nature of the liquid in the vicinity of polar and charged groups and the associated non-pairwise character of the forces, thus improving upon conventional continuum representations. The solvent is treated as a polar and polarizable medium, with dielectric properties described by an inhomogeneous version of the Onsager theory. This treatment leads to an effective position-dependent dielectric permittivity that incorporates saturation effects of the electric field and the spatial variation of the liquid density. The non-pairwise additivity of the liquid-structure forces is represented by centers of force located at specific points in the liquid phase. These out-of-the-solute centers are positioned at the peaks of liquid density and exert local, external forces on the atoms of the solute. The density is calculated from a barometric law, using a Lennard-Jones-type solute-liquid effective interaction potential. The conceptual aspects of the model and its exact numerical solutions are discussed for single alkali and halide ions and for ion-pair interactions. The practical aspects of the model and the simplifications introduced for efficient computation of forces in molecular solutes are discussed in the context of polar and charged amino acid dimers. The model reproduces the contact and solvent-separated minima and the desolvation barriers of intermolecular potentials of mean force of amino acid dimers, as observed in atomistic dynamics simulations. Possible refinements based on an improved treatment of molecular correlations are discussed.

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.

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.

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.