Does water drive protein folding?

Effect of Main and Side Chains on the Folding Mechanism of the Trp-Cage Miniprotein

Proteins that do not fold into their functional native state have been linked to diseases. In this study, the influence of the main and side chains of individual amino acids on the folding of the tryptophan cage (Trp-cage), a designed 20-residue miniprotein, was analyzed. For this purpose, we calculated the solvation free energy (SFE) contributions of individual atoms by using the 3D-reference interaction site model with the atomic decomposition method. The mechanism by which the Trp-cage is stabilized during the folding process was examined by calculating the total energy, which is the sum of the conformational energy and SFE. The folding process of the Trp-cage resulted in a stable native state, with a total energy that was 62.4 kcal/mol lower than that of the unfolded state. The solvation entropy, which is considered to be responsible for the hydrophobic effect, contributed 31.3 kcal/mol to structural stabilization. In other words, the contribution of the solvation entropy accounted for approximately half of the total contribution to Trp-cage folding. The hydrophobic core centered on Trp6 contributed 15.6 kcal/mol to the total energy, whereas the solvation entropy contribution was 6.3 kcal/mol. The salt bridge formed by the hydrophilic side chains of Asp9 and Arg16 contributed 10.9 and 5.0 kcal/mol, respectively. This indicates that not only the hydrophobic core but also the salt bridge of the hydrophilic side chains gain solvation entropy and contribute to stabilizing the native structure of the Trp-cage.

Analysis of Protein Folding Simulation with Moving Root Mean Square Deviation

We apply moving root-mean-square deviation (mRMSD), which does not require a reference structure, as a method for analyzing protein dynamics. This method can be used to calculate the root-mean-square deviation (RMSD) of structure between two specified time points and to analyze protein dynamics behavior through time series analysis. We applied this method to the Trp-cage trajectory calculated by the Anton supercomputer and found that it shows regions of stable states as well as the conventional RMSD. In addition, we extracted a characteristic structure in which the side chains of Asp1 and Arg16 form hydrogen bonds near the most stable structure of the Trp-cage. We also determined that ≥20 ns is an appropriate time interval to investigate protein dynamics using mRMSD. Applying this method to NuG2 protein, we found that mRMSD can be used to detect regions of metastable states in addition to the stable state. This method can be applied to molecular dynamics simulations of proteins whose stable structures are unknown.

Structural Stability Analysis of Proteins Using End-to-End Distance: A 3D-RISM Approach

The stability of a protein is determined from its properties and surrounding solvent. In our previous study, the total energy as a sum of the conformational and solvation free energies was demonstrated to be an appropriate energy function for evaluating the stability of a protein in a protein folding system. We plotted the various energies against the root mean square deviation, required as a reference structure. Herein, we replotted the various energies against the end-to-end distance between the N- and C-termini, which is not a required reference and is experimentally measurable. The solvation free energies for all proteins tend to be low as the end-to-end distance increases, whereas the conformational energies tend to be low as the end-to-end distance decreases. The end-to-end distance is one of interesting measures to study the behavior of proteins.

Effect of Humid Air Exposed to IR Radiation on Enzyme Activity

Water vapor absorbs well in the infra-red region of the electromagnetic spectrum. Absorption of radiant energy by water or water droplets leads to formation of exclusion zone water that possesses peculiar physico-chemical properties. In the course of this study, normally functioning and damaged alkaline phosphatase, horseradish peroxidase and catalase were treated with humid air irradiated with infrared light with a wavelength in the range of 1270 nm and referred to as coherent humidity (CoHu). One-minute long treatment with CoHu helped to partially protect enzymes from heat inactivation, mixed function oxidation, and loss of activity due to partial unfolding. Authors suggest that a possible mechanism underlying the observed effects involves altering the physicochemical properties of aqueous media while treatment of the objects with CoHu where CoHu acts as an intermediary.

Water-mediated interactions destabilize proteins

Proteins are folded to avoid exposure of the nonpolar groups to water because water-mediated interactions between nonpolar groups are a promising factor in the thermodynamic stabilities of proteins-which is a well-accepted view as one of the unique effects of hydrophobic interactions. This article poses a critical question for this classical view by conducting an accurate solvation free-energy calculation for a thermodynamic cycle of a protein folding using a liquid-state density functional theory. Here, the solvation-free energy for a leucine zipper formation was examined in the coiled-coil protein GCN4-p1, a typical model for hydrophobic interactions, which demonstrated that water-mediated interactions were unfavorable for the association of nonpolar groups in the native state, while the dispersion forces between them were, instead, responsible for the association. Furthermore, the present analysis well predicted the isolated helical state stabilized by pressure, which was previously observed in an experiment. We reviewed the problems in the classical concept and semiempirical presumption that the energetic cost of the hydration of nonpolar groups is a driving force of folding.

Highly precise characterization of the hydration state upon thermal denaturation of human serum albumin using a 65 GHz dielectric sensor

The biological functions of proteins depend on harmonization with hydration water surrounding them. Indeed, the dynamical transition of proteins, such as thermal denaturation, is dependent on the changes in the mobility of hydration water. However, the role of hydration water during dynamical transition is yet to be fully understood due to technical limitations in precisely characterizing the amount of hydration water. A state-of-the-art CMOS dielectric sensor consisting of 65 GHz LC resonators addressed this issue by utilizing the feature that oscillation frequency sensitively shifts in response to the complex dielectric constant at 65 GHz with extremely high precision. This study aimed to establish an analytical algorithm to derive the hydration number from the measured frequency shift and to demonstrate the transition of hydration number upon the thermal denaturation of human serum albumin. The determined hydration number in the native state drew a "global" hydration picture beyond the first solvation shell, with substantially reduced uncertainty of the hydration number (about ±1%). This allowed the detection of a rapid increase in the hydration number at about 55 °C during the heating process, which was in excellent phase with the irreversible rupture of the α-helical structure into solvent-exposed extended chains, whereas the hydration number did not trace the forward path in the subsequent cooling process. Our result indicates that the weakening of water hydrogen bonds trigger the unfolding of the protein structure first, followed by the changes in the number of hydration water as a consequence of thermal denaturation.

Mutation-induced change in chignolin stability from π-turn to α-turn

A mutation from threonine to proline at the eighth residue in chignolin changes π-turn to α-turn.

The Role of Hydrogen Bonding in the Folding/Unfolding Process of Hydrated Lysozyme: A Review of Recent NMR and FTIR Results

The biological activity of proteins depends on their three-dimensional structure, known as the native state. The main force driving the correct folding mechanism is the hydrophobic effect and when this folding kinetics is altered, aggregation phenomena intervene causing the occurrence of illnesses such as Alzheimer and Parkinson’s diseases. The other important effect is performed by water molecules and by their ability to form a complex network of hydrogen bonds whose dynamics influence the mobility of protein amino acids. In this work, we review the recent results obtained by means of spectroscopic techniques, such as Fourier Transform Infrared (FTIR) and Nuclear Magnetic Resonance (NMR) spectroscopies, on hydrated lysozyme. In particular, we explore the Energy Landscape from the thermal region of configurational stability up to that of the irreversible denaturation. The importance of the coupling between the solute and the solvent will be highlighted as well as the different behaviors of hydrophilic and hydrophobic moieties of protein amino acid residues.

Hydration Structure and Dynamics of Inhibitor-Bound HIV-1 Protease

Water is essential in many biological processes, and the hydration structure plays a critical role in facilitating protein folding, dynamics, and ligand binding. A variety of biophysical spectroscopic techniques have been used to probe the water solvating proteins, often complemented with molecular dynamics (MD) simulations to resolve the spatial and dynamic features of the hydration shell, but comparing relative water structure is challenging. In this study 1 μs MD simulations were performed to identify and characterize hydration sites around HIV-1 protease bound to an inhibitor, darunavir (DRV). The water density, hydration site occupancy, extent and anisotropy of fluctuations, coordinated water molecules, and hydrogen bonds were characterized and compared to the properties of bulk water. The water density of the principal hydration shell was found to be higher than bulk, dependent on the topology and physiochemical identity of the biomolecular surface. The dynamics of water molecules occupying principal hydration sites was highly dependent on the number of water-water interactions and inversely correlated with hydrogen bonds to the protein-inhibitor complex. While many waters were conserved following the symmetry of homodimeric HIV protease, the asymmetry induced by DRV resulted in asymmetric lower-occupancy hydration sites at the concave surface of the active site. Key interactions between water molecules and the protease, that stabilize the protein in the inhibited form, were altered in a drug resistant variant of the protease indicating that modulation of solvent-solute interactions might play a key role in conveying drug resistance. Our analysis provides insights into the interplay between an enzyme inhibitor complex and the hydration shell and has implications in elucidating water structure in a variety of biological processes and applications including ligand binding, inhibitor design, and resistance.

Analysis of Structural Stability of Chignolin

We discuss the stability of an entire protein and the influence of main chains and side chains of individual amino acids to investigate the protein-folding mechanism. For this purpose, we calculated the solvation free-energy contribution of individual atoms using the three-dimensional reference interaction site model with the atomic decomposition method. We generated structures of chignolin miniprotein by a molecular dynamics simulation and classified them into six types: native 1, native 2, misfolded 1, misfolded 2, intermediate, and unfolded states. The total energies of the native (-171.1 kcal/mol) and misfolded (-171.2 kcal/mol) states were almost the same and lower than those of the intermediate (-158.5 kcal/mol) and unfolded (-148.1 kcal/mol) states; however, their components were different. In the native state, the side-chain interaction between Thr6 and Thr8 is important for the formation of π-turn. On the other hand, the hydrogen bonds between the atoms of the main chains in the misfolded state become stronger than those in the intermediate state.

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.

Structure of the ordered hydration of amino acids in proteins: analysis of crystal structures

Crystallography provides unique information about the arrangement of water molecules near protein surfaces. Using a nonredundant set of 2818 protein crystal structures with a resolution of better than 1.8 Å, the extent and structure of the hydration shell of all 20 standard amino-acid residues were analyzed as function of the residue conformation, secondary structure and solvent accessibility. The results show how hydration depends on the amino-acid conformation and the environment in which it occurs. After conformational clustering of individual residues, the density distribution of water molecules was compiled and the preferred hydration sites were determined as maxima in the pseudo-electron-density representation of water distributions. Many hydration sites interact with both main-chain and side-chain amino-acid atoms, and several occurrences of hydration sites with less canonical contacts, such as carbon-donor hydrogen bonds, OH-π interactions and off-plane interactions with aromatic heteroatoms, are also reported. Information about the location and relative importance of the empirically determined preferred hydration sites in proteins has applications in improving the current methods of hydration-site prediction in molecular replacement, ab initio protein structure prediction and the set-up of molecular-dynamics simulations.

Investigation of pH-induced conformational change and hydration of poly(methacrylic acid) by analytical ultracentrifugation

Analytical ultracentrifugation was performed on poly(methacrylic acid) (PMAA) with a series of weight average molar masses (Mw) in aqueous solutions as a function of pH. The scales of the sedimentation coefficient (s) and the diffusion coefficient (D) to Mw at infinite dilutions were obtained at different pH values, indicating that PMAA chains adopt a collapsed structure at low pH values, and stretch at pH higher than 5.2. Our results show that the sedimentation coefficient exhibits a minimum at pH ∼ 6.0, presumably due to the effect of the conformational change and the hydration state of PMAA chains. When pH increases from 6.0 to 8.5, PMAA chains with high molar mass shrink a little bit, presumably because the sodium ions act as a bridging agent between nonadjacent carboxylate groups. Furthermore, the weight average molar mass of PMAA at pH 8.5 increases by one fold than that at pH 4.0, indicating the condensation of sodium ions and the increase in the number of hydration water molecules around carboxylate groups at high pH values.

Role of Internal Water on Protein Thermal Stability: The Case of Homologous G Domains

In this work, we address the question of whether the enhanced stability of thermophilic proteins has a direct connection with internal hydration. Our model systems are two homologous G domains of different stability: the mesophilic G domain of the elongation factor thermal unstable protein from E. coli and the hyperthermophilic G domain of the EF-1α protein from S. solfataricus. Using molecular dynamics simulation at the microsecond time scale, we show that both proteins host water molecules in internal cavities and that these molecules exchange with the external solution in the nanosecond time scale. The hydration free energy of these sites evaluated via extensive calculations is found to be favorable for both systems, with the hyperthermophilic protein offering a slightly more favorable environment to host water molecules. We estimate that, under ambient conditions, the free energy gain due to internal hydration is about 1.3 kcal/mol in favor of the hyperthermophilic variant. However, we also find that, at the high working temperature of the hyperthermophile, the cavities are rather dehydrated, meaning that under extreme conditions other molecular factors secure the stability of the protein. Interestingly, we detect a clear correlation between the hydration of internal cavities and the protein conformational landscape. The emerging picture is that internal hydration is an effective observable to probe the conformational landscape of proteins. In the specific context of our investigation, the analysis confirms that the hyperthermophilic G domain is characterized by multiple states and it has a more flexible structure than its mesophilic homologue.

Massively parallel implementation of 3D-RISM calculation with volumetric 3D-FFT

A new three-dimensional reference interaction site model (3D-RISM) program for massively parallel machines combined with the volumetric 3D fast Fourier transform (3D-FFT) was developed, and tested on the RIKEN K supercomputer. The ordinary parallel 3D-RISM program has a limitation on the number of parallelizations because of the limitations of the slab-type 3D-FFT. The volumetric 3D-FFT relieves this limitation drastically. We tested the 3D-RISM calculation on the large and fine calculation cell (2048(3) grid points) on 16,384 nodes, each having eight CPU cores. The new 3D-RISM program achieved excellent scalability to the parallelization, running on the RIKEN K supercomputer. As a benchmark application, we employed the program, combined with molecular dynamics simulation, to analyze the oligomerization process of chymotrypsin Inhibitor 2 mutant. The results demonstrate that the massive parallel 3D-RISM program is effective to analyze the hydration properties of the large biomolecular systems.