Understanding the role of hydrogen bonds in water dynamics and protein stability

Thermal tuning of protein hydration in a hyperthermophilic enzyme

Water at the protein surface is an active biological molecule that plays a critical role in many functional processes. Using NMR-restrained MD simulations, we here addressed how protein hydration is tuned at high biological temperatures by analysing homologous acylphosphatase enzymes (AcP) possessing similar structure and dynamics under very different thermal conditions. We found that the hyperthermophilic Sso AcP at 80°C interacts with a lower number of structured waters in the first hydration shell than its human homologous mt AcP at 37°C. Overall, the structural and dynamical properties of waters at the surface of the two enzymes resulted similar in the first hydration shell, including solvent molecules residing in the active site. By contrast the dynamical content of water molecules in the second hydration shell was found to diverge, with higher mobility observed in Sso AcP at 80°C. Taken together the results delineate the subtle differences in the hydration properties of mt AcP and Sso AcP, and indicate that the concept of corresponding states with equivalent dynamics in homologous mesophilic and hyperthermophylic proteins should be extended to the first hydration shell.

Structure and activity of a thermally stable mutant of Acanthamoeba actophorin

Actophorin, which was recently tested for crystallization under microgravity on the International Space Station, was subjected to mutagenesis to identify a construct with improved biophysical properties that were expected to improve the extent of diffraction. First, 20 mutations, including one C-terminal deletion of three residues, were introduced individually into actophorin, resulting in modest increases in thermal stability of between +0.5°C and +2.2°C. All but two of the stabilizing mutants increased both the rates of severing F-actin filaments and of spontaneous polymerization of pyrenyl G-actin in vitro. When the individual mutations were combined into a single actophorin variant, Acto-2, the overall thermal stability was 22°C higher than that of wild-type actophorin. When an inactivating S2P mutation in Acto-2 was restored, Acto-2/P2S was more stable by 20°C but was notably more active than the wild-type protein. The inactivating S2P mutation reaffirms the importance that Ser2 plays in the F-actin-severing reaction. The crystal structure of Acto-2 was solved to 1.7 Å resolution in a monoclinic space group, a first for actophorin. Surprisingly, despite the increase in thermal stability, the extended β-turn region, which is intimately involved in interactions with F-actin, is disordered in one copy of Acto-2 in the asymmetric unit. These observations emphasize the complex interplay among protein thermal stability, function and dynamics.

Advances in the study of supercooled water

In this review, we report recent progress in the field of supercooled water. Due to its uniqueness, water presents numerous anomalies with respect to most simple liquids, showing polyamorphism both in the liquid and in the glassy state. We first describe the thermodynamic scenarios hypothesized for the supercooled region and in particular among them the liquid-liquid critical point scenario that has so far received more experimental evidence. We then review the most recent structural indicators, the two-state model picture of water, and the importance of cooperative effects related to the fact that water is a hydrogen-bonded network liquid. We show throughout the review that water's peculiar properties come into play also when water is in solution, confined, and close to biological molecules. Concerning dynamics, upon mild supercooling water behaves as a fragile glass former following the mode coupling theory, and it turns into a strong glass former upon further cooling. Connections between the slow dynamics and the thermodynamics are discussed. The translational relaxation times of density fluctuations show in fact the fragile-to-strong crossover connected to the thermodynamics arising from the existence of two liquids. When considering also rotations, additional crossovers come to play. Mobility-viscosity decoupling is also discussed in supercooled water and aqueous solutions. Finally, the polyamorphism of glassy water is considered through experimental and simulation results both in bulk and in salty aqueous solutions. Grains and grain boundaries are also discussed.

Effect of trehalose on protein cryoprotection: Insights into the mechanism of slowing down of hydration water

We study, with molecular dynamics simulations, a lysozyme protein immersed in a water-trehalose solution upon cooling. The aim is to understand the cryoprotectant role played by this disaccharide through the modifications that it induces on the slow dynamics of protein hydration water with its presence. The α-relaxation shows a fragile to strong crossover about 20° higher than that in the bulk water phase and 15° higher than that in lysozyme hydration water without trehalose. The protein hydration water without trehalose was found to show a second slower relaxation exhibiting a strong to strong crossover coupled with the protein dynamical transition. This slower relaxation time importantly appears enormously slowed down in our cryoprotectant solution. On the other hand, this long-relaxation in the presence of trehalose is also connected with a stronger damping of the protein structural fluctuations than that found when the protein is in contact with the pure hydration water. Therefore, this appears to be the mechanism through which trehalose manifests its cryoprotecting function.

Protein folding and surface interaction phase diagrams in vitro and in cells

Protein stability is subject to environmental perturbations such as pressure and crowding, as well as sticking to other macromolecules and quinary structure. Thus, the environment inside and outside the cell plays a key role in how proteins fold, interact, and function on the scale from a few molecules to macroscopic ensembles. This review discusses three aspects of protein phase diagrams: first, the relevance of phase diagrams to protein folding and function in vitro and in cells; next, how the evolution of protein surfaces impacts on interaction phase diagrams; and finally, how phase separation plays a role on much larger length-scales than individual proteins or oligomers, when liquid phase-separated regions form to assist protein function and cell homeostasis.

Protein Unfolding and Aggregation near a Hydrophobic Interface

The behavior of proteins near interfaces is relevant for biological and medical purposes. Previous results in bulk show that, when the protein concentration increases, the proteins unfold and, at higher concentrations, aggregate. Here, we study how the presence of a hydrophobic surface affects this course of events. To this goal, we use a coarse-grained model of proteins and study by simulations their folding and aggregation near an ideal hydrophobic surface in an aqueous environment by changing parameters such as temperature and hydrophobic strength, related, e.g., to ions concentration. We show that the hydrophobic surface, as well as the other parameters, affect both the protein unfolding and aggregation. We discuss the interpretation of these results and define future lines for further analysis, with their possible implications in neurodegenerative diseases.

An assessment of water placement algorithms in quantum mechanics/molecular mechanics modeling: the case of rhodopsins' first spectral absorption band maxima

Quantum mechanics/molecular mechanics (QM/MM) models are a widely used tool to obtain detailed insight into the properties and functioning of proteins. The outcome of QM/MM studies heavily depends on the quality of the applied QM/MM model. Prediction and right placement of internal water molecules in protein cavities is one of the critical parts of any QM/MM model construction. Herein, we performed a systematic study of four protein hydration algorithms. We tested these algorithms for their ability to predict X-ray-resolved water molecules for a set of membrane photosensitive rhodopsin proteins, as well as the influence of the applied water placement algorithms on the QM/MM calculated absorption maxima (λ_max) of these proteins. We used 49 rhodopsins and their intermediates with available X-ray structures as the test set. We found that a proper choice of hydration algorithms and setups is needed to predict functionally important water molecules in the chromophore-binding cavity of rhodopsins, such as the water cluster in the N-H region of bacteriorhodopsin or two water molecules in the binding pocket of bovine visual rhodopsin. The QM/MM calculated λ_max of rhodopsins is also quite sensitive to the applied protein hydration protocols. The best methodology allows obtaining an 18.0 nm average value for the absolute deviation of the calculated λ_max from the experimental λ_max. Although the major effect of water molecules on λ_max originates from the water molecules located in the binding pocket, the water molecules outside the binding pocket also affect the calculated λ_max mainly by causing a reorganization of the protein structure. The results reported in this study can be used for the evaluation and further development of hydration methodologies, in general, and rhodopsin QM/MM models, in particular.

Proteins are Solitary! Pathways of Protein Folding and Aggregation in Protein Mixtures

We present a computational and experimental study on the folding and aggregation in solutions of multiple protein mixtures at different concentrations. We show how in protein mixtures each component is capable of maintaining its folded state at densities greater than the one at which they would precipitate in single-species solutions. We demonstrate the generality of our observation over many different proteins using computer simulations capable of fully characterizing the cross-aggregation phase diagram of all the mixtures. Dynamic light scattering experiments were performed to evaluate the aggregation of two proteins, bovine serum albumin (BSA) and consensus tetratricopeptide repeat (CTPR), in solutions of one or both proteins. The experiments confirm our hypothesis and the simulations. These findings elucidate critical aspects of the cross-regulation of expression and aggregation of proteins exerted by the cell and on the evolutionary selection of folding and non-aggregating protein sequences, paving the way for new experimental tests.

How Heterogeneous Are Trehalose/Glycerol Cryoprotectant Mixtures? A Combined Time-Resolved Fluorescence and Computer Simulation Investigation

Heterogeneity and molecular motions in representative cryoprotectant mixtures made of trehalose and glycerol are investigated in the temperature range 298 ≤ T (K) ≤ 353, via time-resolved fluorescence Stokes shift and anisotropy measurements, and molecular dynamics simulations of four-point density-time correlations and H-bond relaxations. Mixtures containing 5 and 20 wt % of trehalose along with neat glycerol are studied. Viscosity coefficients for these systems lie in the range 0.30 < η (P) < 23. Measured solute (Coumarin 153) rotation and solvation times reveal a substantial departure from the hydrodynamic viscosity dependence, suggesting the strong microheterogeneous nature of these systems. Fluorescence anisotropy decays are highly nonexponential, reflecting a non-Markovian character of the medium friction. A complete missing of the Stokes shift dynamics in these systems at 298 K but partial detection of it at other higher temperatures (shift magnitude being ∼400-600 cm^-1) indicates rigid solute environments. An amorphous solid-like feature emerges in the simulated radial distribution functions at these temperatures. Analyses of mean squared displacements reveal rattling-in-a-cage motion, non-Gaussian displacement distributions, and strong dynamic heterogeneity features. Simulated dynamic structure factors and four-point correlations hint, respectively, at very long α-relaxation and correlated time scales at 298 K. This explains the long solute rotation times (∼80-200 ns) measured at 298 K. Stretched exponential decay of the simulated H-bond relaxations with long time scales further highlights the strong temporal heterogeneity and slow dynamics inherent to these systems. In summary, this work provides the first insight into the molecular motions and interspecies interaction in a representative cryoprotectant mixture, and stimulates further study to investigate the interconnection between cryoprotection and dynamic heterogeneity.

Hydration of proteins and nucleic acids: Advances in experiment and theory. A review

BACKGROUND

Most biological processes involve water, and the interactions of biomolecules with water affect their structure, function and dynamics.

SCOPE OF REVIEW

This review summarizes the current knowledge of protein and nucleic acid interactions with water, with a special focus on the biomolecular hydration layer. Recent developments in both experimental and computational methods that can be applied to the study of hydration structure and dynamics are reviewed, including software tools for the prediction and characterization of hydration layer properties.

MAJOR CONCLUSIONS

In the last decade, important advances have been made in our understanding of the factors that determine how biomolecules and their aqueous environment influence each other. Both experimental and computational methods contributed to the gradually emerging consensus picture of biomolecular hydration.

GENERAL SIGNIFICANCE

An improved knowledge of the structural and thermodynamic properties of the hydration layer will enable a detailed understanding of the various biological processes in which it is involved, with implications for a wide range of applications, including protein-structure prediction and structure-based drug design.

Contribution of Water to Pressure and Cold Denaturation of Proteins

The mechanisms of cold and pressure denaturation of proteins are matter of debate and are commonly understood as due to water-mediated interactions. Here, we study several cases of proteins, with or without a unique native state, with or without hydrophilic residues, by means of a coarse-grain protein model in explicit solvent. We show, using Monte Carlo simulations, that taking into account how water at the protein interface changes its hydrogen bond properties and its density fluctuations is enough to predict protein stability regions with elliptic shapes in the temperature-pressure plane, consistent with previous theories. Our results clearly identify the different mechanisms with which water participates to denaturation and open the perspective to develop advanced computational design tools for protein engineering.

Role of cavities and hydration in the pressure unfolding of T4 lysozyme

It is well known that high hydrostatic pressures can induce the unfolding of proteins. The physical underpinnings of this phenomenon have been investigated extensively but remain controversial. Changes in solvation energetics have been commonly proposed as a driving force for pressure-induced unfolding. Recently, the elimination of void volumes in the native folded state has been argued to be the principal determinant. Here we use the cavity-containing L99A mutant of T4 lysozyme to examine the pressure-induced destabilization of this multidomain protein by using solution NMR spectroscopy. The cavity-containing C-terminal domain completely unfolds at moderate pressures, whereas the N-terminal domain remains largely structured to pressures as high as 2.5 kbar. The sensitivity to pressure is suppressed by the binding of benzene to the hydrophobic cavity. These results contrast to the pseudo-WT protein, which has a residual cavity volume very similar to that of the L99A-benzene complex but shows extensive subglobal reorganizations with pressure. Encapsulation of the L99A mutant in the aqueous nanoscale core of a reverse micelle is used to examine the hydration of the hydrophobic cavity. The confined space effect of encapsulation suppresses the pressure-induced unfolding transition and allows observation of the filling of the cavity with water at elevated pressures. This indicates that hydration of the hydrophobic cavity is more energetically unfavorable than global unfolding. Overall, these observations point to a range of cooperativity and energetics within the T4 lysozyme molecule and illuminate the fact that small changes in physical parameters can significantly alter the pressure sensitivity of proteins.

Critical behavior of a water monolayer under hydrophobic confinement

The properties of water can have a strong dependence on the confinement. Here, we consider a water monolayer nanoconfined between hydrophobic parallel walls under conditions that prevent its crystallization. We investigate, by simulations of a many-body coarse-grained water model, how the properties of the liquid are affected by the confinement. We show, by studying the response functions and the correlation length and by performing finite-size scaling of the appropriate order parameter, that at low temperature the monolayer undergoes a liquid-liquid phase transition ending in a critical point in the universality class of the two-dimensional (2D) Ising model. Surprisingly, by reducing the linear size L of the walls, keeping the walls separation h constant, we find a 2D-3D crossover for the universality class of the liquid-liquid critical point for L=h ≃ 50, i.e. for a monolayer thickness that is small compared to its extension. This result is drastically different from what is reported for simple liquids, where the crossover occurs for L=h ≃ 5, and is consistent with experimental results and atomistic simulations. We shed light on these findings showing that they are a consequence of the strong cooperativity and the low coordination number of the hydrogen bond network that characterizes water.

A facile process for soak-and-peel delamination of CVD graphene from substrates using water

We demonstrate a simple technique to transfer chemical vapour deposited (CVD) graphene from copper and platinum substrates using a soak-and-peel delamination technique utilizing only hot deionized water. The lack of chemical etchants results in cleaner CVD graphene films minimizing unintentional doping, as confirmed by Raman and electrical measurements. The process allows the reuse of substrates and hence can enable the use of oriented substrates for growth of higher quality graphene, and is an inherently inexpensive and scalable process for large-area production.

Carbohydrate recognition by RpfB from Mycobacterium tuberculosis unveiled by crystallographic and molecular dynamics analyses

Resuscitation of Mtb is crucial to the etiology of Tuberculosis, because latent tuberculosis is estimated to affect one-third of the world population. The resuscitation-promoting factor RpfB is mainly responsible for Mtb resuscitation from dormancy. Given the impact of latent Tuberculosis, RpfB represents an interesting target for tuberculosis drug discovery. However, no molecular models of substrate binding and catalysis are hitherto available for this enzyme. Here, we identified key interactions involved in substrate binding to RpfB by combining x-ray diffraction studies and computational approaches. The crystal structure of RpfB catalytic domain in complex with N,N',N"-triacetyl-chitotriose, as described here, provides the first, to our knowledge, atomic representation of ligand recognition by RpfB and demonstrates that the strongest interactions are established by the N-acetylglucosamine moiety in the central region of the enzyme binding cleft. Molecular dynamics analyses provided information on the dynamic behavior of protein-substrate interactions and on the role played by the solvent in RpfB function. These data combined with sequence conservation analysis suggest that Glu-292 is the sole residue crucial for catalysis, implying that RpfB acts via the formation of an oxocarbenium ion rather than a covalent intermediate. Present data represent a solid base for the design of effective drug inhibitors of RpfB. Moreover, homology models were generated for the catalytic domains of all members of the Mtb Rpf family (RpfA-E). The analysis of these models unveiled analogies and differences among the different members of the Rpf protein family.

Unifying microscopic mechanism for pressure and cold denaturations of proteins

We study the stability of globular proteins as a function of temperature and pressure through NPT simulations of a coarse-grained model. We reproduce the elliptical stability of proteins and highlight a unifying microscopic mechanism for pressure and cold denaturations. The mechanism involves the solvation of nonpolar residues with a thin layer of water. These solvated states have lower volume and lower hydrogen-bond energy compared to other conformations of nonpolar solutes. Hence, these solvated states are favorable at high pressure and low temperature, and they facilitate protein unfolding under these thermodynamical conditions.