Free-energy function for discriminating the native fold of a protein from misfolded decoys

Universal effects of solvent species on the stabilized structure of a protein

We investigate the effects of solvent specificities on the stability of the native structure (NS) of a protein on the basis of our free-energy function (FEF). We use CPB-bromodomain (CBP-BD) and apoplastocyanin (apoPC) as representatives of the protein universe and water, methanol, ethanol, and cyclohexane as solvents. The NSs of CBP-BD and apoPC consist of 66% α-helices and of 35% β-sheets and 4% α-helices, respectively. In order to assess the structural stability of a given protein immersed in each solvent, we contrast the FEF of its NS against that of a number of artificially created, misfolded decoys possessing the same amino-acid sequence but significantly different topology and α-helix and β-sheet contents. In the FEF, we compute the solvation entropy using the morphometric approach combined with the integral equation theories, and the change in electrostatic (ES) energy upon the folding is obtained by an explicit atomistic but simplified calculation. The ES energy change is represented by the break of protein-solvent hydrogen bonds (HBs), formation of protein intramolecular HBs, and recovery of solvent-solvent HBs. Protein-solvent and solvent-solvent HBs are absent in cyclohexane. We are thus able to separately evaluate the contributions to the structural stability from the entropic and energetic components. We find that for both CBP-BD and apoPC, the energetic component dominates in methanol, ethanol, and cyclohexane, with the most stable structures in these solvents sharing the same characteristics described as an association of α-helices. In particular, those in the two alcohols are identical. In water, the entropic component is as strong as or even stronger than the energetic one, with a large gain of translational, configurational entropy of water becoming crucially important so that the relative contents of α-helix and β-sheet and the content of total secondary structures are carefully selected to achieve sufficiently close packing of side chains. If the energetic component is excluded for a protein in water, the priority is given to closest side-chain packing, giving rise to the formation of a structure with very low α-helix and β-sheet contents. Our analysis, which requires minimal computational effort, can be applied to any protein immersed in any solvent and provides robust predictions that are quite consistent with the experimental observations for proteins in different solvent environments, thus paving the way toward a more detailed understanding of the folding process.

Strong Dependence of Hydration State of F-Actin on the Bound Mg(2+)/Ca(2+) Ions

Understanding of the hydration state is an important issue in the chemomechanical energetics of versatile biological functions of polymerized actin (F-actin). In this study, hydration-state differences of F-actin by the bound divalent cations are revealed through precision microwave dielectric relaxation (DR) spectroscopy. G- and F-actin in Ca- and Mg-containing buffer solutions exhibit dual hydration components comprising restrained water with DR frequency f2 (<fw: DR frequency of bulk solvent, 17 GHz at 20 °C) and hypermobile water (HMW) with DR frequency f1 (>fw). The hydration state of F-actin is strongly dependent on the ionic composition. In every buffer tested, the HMW signal Dhyme (≡ (f1 - fw)δ1/(fwδw)) of F-actin is stronger than that of G-actin, where δw is DR-amplitude of bulk solvent and δ1 is that of HMW in a fixed-volume ellipsoid containing an F-actin and surrounding water in solution. Dhyme value of F-actin in Ca2.0-buffer (containing 2 mM Ca(2+)) is markedly higher than in Mg2.0-buffer (containing 2 mM Mg(2+)). Moreover, in the presence of 2 mM Mg(2+), the hydration state of F-actin is changed by adding a small fraction of Ca(2+) (∼0.1 mM) and becomes closer to that of the Ca-bound form in Ca2.0-buffer. This is consistent with the results of the partial specific volume and the Cotton effect around 290 nm in the CD spectra, indicating a change in the tertiary structure and less apparent change in the secondary structure of actin. The number of restrained water molecules per actin (N2) is estimated to be 1600-2100 for Ca2.0- and F-buffer and ∼2500 for Mg2.0-buffer at 10-15 °C. These numbers are comparable to those estimated from the available F-actin atomic structures as in the first water layer. The number of HMW molecules is roughly explained by the volume between the equipotential surface of -kT/2e and the first water layer of the actin surface by solving the Poisson-Boltzmann equation using UCSF Chimera.

Enhanced sampling simulations to construct free-energy landscape of protein-partner substrate interaction

Molecular dynamics (MD) simulations using all-atom and explicit solvent models provide valuable information on the detailed behavior of protein-partner substrate binding at the atomic level. As the power of computational resources increase, MD simulations are being used more widely and easily. However, it is still difficult to investigate the thermodynamic properties of protein-partner substrate binding and protein folding with conventional MD simulations. Enhanced sampling methods have been developed to sample conformations that reflect equilibrium conditions in a more efficient manner than conventional MD simulations, thereby allowing the construction of accurate free-energy landscapes. In this review, we discuss these enhanced sampling methods using a series of case-by-case examples. In particular, we review enhanced sampling methods conforming to trivial trajectory parallelization, virtual-system coupled multicanonical MD, and adaptive lambda square dynamics. These methods have been recently developed based on the existing method of multicanonical MD simulation. Their applications are reviewed with an emphasis on describing their practical implementation. In our concluding remarks we explore extensions of the enhanced sampling methods that may allow for even more efficient sampling.

Correlation analysis for heat denaturation of Trp-cage miniprotein with explicit solvent

Energetics was analyzed for Trp-cage miniprotein in water to elucidate the solvation effect in heat denaturation. The solvation free energy was computed for a set of protein structures at room and high temperatures with all-atom molecular dynamics simulation combined with the solution theory in the energy representation, and its correlations were investigated against the intramolecular (structural) energy of the protein and the average interaction energy of the protein with the solvent water. It was observed both at room and high temperatures that the solvation free energy is anticorrelated to the structural energy and varies in parallel to the electrostatic component of the protein-water interaction energy without correlations to the van der Waals and excluded-volume components. When the set of folded structures sampled at room temperature was compared with the set of unfolded ones at high temperature, it was found that the preference order of the two sets is in correspondence to the van der Waals and excluded-volume components in the sum of the protein intramolecular and protein-water intermolecular interactions and is not distinguished by the electrostatic component.

Physicochemical origin of high correlation between thermal stability of a protein and its packing efficiency: a theoretical study for staphylococcal nuclease mutants

There is an empirical rule that the thermal stability of a protein is related to the packing efficiency or core volume of the folded state and the protein tends to exhibit higher stability as the backbone and side chains are more closely packed. Previously, the wild type and its nine mutants of staphylococcal nuclease were compared by examining their folded structures. The results obtained were as follows: The stability was not correlated with the number of intramolecular hydrogen bonds, intramolecular electrostatic interaction energy, or degree of burial of the hydrophobic surface; though the empirical rule mentioned above held, it was not the proximate cause of higher stability; and the number of van der Waals contacts N vdW, or equivalently, the intramolecular van der Waals interaction energy was an important factor governing the stability. Here we revisit the wild type and its nine mutants of staphylococcal nuclease using our statistical-mechanical theory for hydration of a protein. A molecular model is employed for water. We show that the pivotal factor is the magnitude of the water-entropy gain upon folding. The gain originates from an increase in the total volume available to the translational displacement of water molecules coexisting with the protein in the system. The magnitude is highly correlated with the denaturation temperature T m. Moreover, the apparent correlation between N vdW and T m as well as the empirical rule is interpretable (i.e., their physicochemical meanings can be clarified) on the basis of the water-entropy effect.

Effects of sugars on the thermal stability of a protein

It is experimentally known that the heat-denaturation temperature of a protein is raised (i.e., its thermal stability is enhanced) by sugar addition. In earlier works, we proposed a physical picture of thermal denaturation of proteins in which the measure of the thermal stability is defined as the solvent-entropy gain upon protein folding at 298 K normalized by the number of residues. A multipolar-model water was adopted as the solvent. The polyatomic structures of the folded and unfolded states of a protein were taken into account in the atomic detail. A larger value of the measure implies higher thermal stability. First, we show that the measure remains effective even when the model water is replaced by the hard-sphere solvent whose number density and molecular diameter are set at those of real water. The physical picture is then adapted to the elucidation of the effects of sugar addition on the thermal stability of a protein. The water-sugar solution is modeled as a binary mixture of hard spheres. The thermal stability is determined by a complex interplay of the diameter of sugar molecules dC and the total packing fraction of the solution η: dC is estimated from the volume per molecule in the sugar crystal and η is calculated using the experimental data of the solution density. We find that the protein is more stabilized as the sucrose or glucose concentration becomes higher and the stabilization effect is stronger for sucrose than for glucose. These results are in accord with the experimental observations. Using a radial-symmetric integral equation theory and the morphometric approach, we decompose the change in the measure upon sugar addition into two components originating from the protein-solvent pair and protein-solvent many-body correlations, respectively. Each component is further decomposed into the excluded-volume and solvent-accessible-surface terms. These decompositions give physical insights into the microscopic origin of the thermal-stability enhancement by sugar addition. As an example, the higher stability of the native state relative to that of the unfolded state is found to be attributable primarily to an increase in the solvent crowding caused by sugar addition. Due to the hydrophilicity of sugar molecules, the addition of sugar by a larger amount or that with a larger molecular size leads to an increase in η which is large enough to make the solvent crowding more serious.

Structural stability of proteins in aqueous and nonpolar environments

A protein folds into its native structure with the α-helix and∕or β-sheet in aqueous solution under the physiological condition. The relative content of these secondary structures largely varies from protein to protein. However, such structural variability is not exhibited in nonaqueous environment. For example, there is a strong trend that alcohol induces a protein to form α-helices, and many of the membrane proteins within the lipid bilayer consists of α-helices. Here we investigate the structural stability of proteins in aqueous and nonpolar environments using our recently developed free-energy function F = (Λ - TS)∕(k(B)T(0)) = Λ∕(k(B)T(0)) - S∕k(B) (T(0) = 298 K and the absolute temperature T is set at T(0)) which is based on statistical thermodynamics. Λ∕(k(B)T(0)) and S∕k(B) are the energetic and entropic components, respectively, and k(B) is Boltzmann's constant. A smaller value of the positive quantity, -S, represents higher efficiency of the backbone and side-chain packing promoted by the entropic effect arising from the translational displacement of solvent molecules or the CH(2), CH(3), and CH groups which constitute nonpolar chains of lipid molecules. As for Λ, in aqueous solution, a transition to a more compact structure of a protein accompanies the break of protein-solvent hydrogen bonds: As the number of donors and acceptors buried without protein intramolecular hydrogen bonding increases, Λ becomes higher. In nonpolar solvent, lower Λ simply implies more intramolecular hydrogen bonds formed. We find the following. The α-helix and β-sheet are advantageous with respect to -S as well as Λ and to be formed as much as possible. In aqueous solution, the solvent-entropy effect on the structural stability is so strong that the close packing of side chains is dominantly important, and the α-helix and β-sheet contents are judiciously adjusted to accomplish it. In nonpolar solvent, the solvent-entropy effect is substantially weaker than in aqueous solution. Λ is crucial and the α-helix is more stable than the β-sheet in terms of Λ, which develops a tendency that α-helices are exclusively chosen. For a membrane protein, α-helices are stabilized as fundamental structural units for the same reason, but their arrangement is performed through the entropic effect mentioned above.

A new theoretical approach to biological self-assembly

Upon biological self-assembly, the number of accessible translational configurations of water in the system increases considerably, leading to a large gain in water entropy. It is important to calculate the solvation entropy of a biomolecule with a prescribed structure by accounting for the change in water-water correlations caused by solute insertion. Modeling water as a dielectric continuum is not capable of capturing the physical essence of the water entropy effect. As a reliable tool, we propose a hybrid of the angle-dependent integral equation theory combined with a multipolar water model and a morphometric approach. Using our methods wherein the water entropy effect is treated as the key factor, we can elucidate a variety of processes such as protein folding, cold, pressure, and heat denaturating of a protein, molecular recognition, ordered association of proteins such as amyloid fibril formation, and functioning of ATP-driven proteins.

Structural characteristics of yeast F1-ATPase before and after 16-degree rotation of the γ subunit: theoretical analysis focused on the water-entropy effect

We have recently proposed a novel picture of the rotation mechanism for F(1)-ATPase [T. Yoshidome, Y. Ito, M. Ikeguchi, and M. Kinoshita, J. Am. Chem. Soc. 133, 4030 (2011)]. In the picture, the asymmetric packing in F(1)-ATPase, originating from the water-entropy effect, plays the key role in the rotation. Here, we analyze the differences between the experimentally determined structures of yeast F(1)-ATPase before and after 16° rotation of the γ subunit with the emphasis on the water-entropy effect. For each of these structures, we calculate the hydration entropies of three sub-complexes comprising the γ subunit, one of the β subunits, and two α subunits adjacent to them. The β(E), β(TP), and β(DP) subunits are involved in sub-complexes I, II, and III, respectively. The calculation is performed using a hybrid of the angle-dependent integral equation theory combined with the molecular model for water and the morphometric approach. The absolute value of the hydration entropy is in the following order: sub-complex I > sub-complex II > sub-complex III. The packing efficiency of the sub-complex follows the opposite order. The rotation gives rise to less efficient packing in sub-complex III and a corresponding water-entropy loss. However, the other two sub-complexes, accompanying water-entropy gains, become more efficiently packed. These results are consistent with our picture of the rotation mechanism, supporting its validity. The water-entropy analysis shows that the interfaces of α(DP)-β(DP) and α(E)-β(E) become more open after the rotation, which is in accord with the experimental observation.

Physical origin of hydrophobicity studied in terms of cold denaturation of proteins: comparison between water and simple fluids

A clue to the physical origin of the hydrophobicity is in the experimental observations, which show that it is weakened at low temperatures. By considering a solvophobic model protein immersed in water and three species of simple solvents, we analyze the temperature dependence of the changes in free energy, energy, and entropy of the solvent upon protein unfolding. The angle-dependent and radial-symmetric integral equation theories and the morphometric approach are employed in the analysis. Each of the changes is decomposed into two terms, which depend on the excluded volume and on the area and curvature of the solvent-accessible surface, respectively. The excluded-volume term of the entropy change is further decomposed into two components representing the protein-solvent pair correlation and the protein-solvent-solvent triplet and higher-order correlation, respectively. We show that water crowding in the system becomes more serious upon protein unfolding but this effect becomes weaker as the temperature is lowered. If the hydrophobicity originated from the water structuring near a nonpolar solute, it would be strengthened upon lowering of the temperature. Among the three species of simple solvents, considerable weakening of the solvophobicity at low temperatures is observed only for the solvent where the particles interact through a strong attractive potential and the particle size is as small as that of water. Even in the case of this solvent, however, cold denaturation of a protein cannot be reproduced. It would be reproducible if the attractive potential was substantially enhanced, but such enhancement causes the appearance of the metastability limit for a single liquid phase.

Characterization of experimentally determined native-structure models of a protein using energetic and entropic components of free-energy function

We show how to characterize the native-structure models of a protein using our free-energy function F which is based on hydration thermodynamics. Ubiquitin is adopted as an example protein. We consider models determined by the X-ray crystallography and two types of NMR model sets. A model set of type 1 comprises candidate models for a fixed native structure, and that of type 2 forms an ensemble of structures representing the structural variability of the native state. In general, the X-ray models give lower F than the NMR models. There is a trend that, as a model deviates more from the model with the lowest F among the X-ray models, its F becomes higher. Model sets of type 1 and those of type 2, respectively, exhibit two different characteristics with respect to the correlation between the deviation and F. It is argued that the total amount of constraints such as NOEs effectively taken into account in constructing the NMR models can be examined by analyzing the behavior of F. We investigate structural characteristics of the models in terms of the energetic and entropic components of F which are relevant to intramolecular hydrogen bonding and to backbone and side-chain packing, respectively.

A scoring function based on solvation thermodynamics for protein structure prediction

We predict protein structure using our recently developed free energy function for describing protein stability, which is focused on solvation thermodynamics. The function is combined with the current most reliable sampling methods, i.e., fragment assembly (FA) and comparative modeling (CM). The prediction is tested using 11 small proteins for which high-resolution crystal structures are available. For 8 of these proteins, sequence similarities are found in the database, and the prediction is performed with CM. Fairly accurate models with average Cα root mean square deviation (RMSD) ∼ 2.0 Å are successfully obtained for all cases. For the rest of the target proteins, we perform the prediction following FA protocols. For 2 cases, we obtain predicted models with an RMSD ∼ 3.0 Å as the best-scored structures. For the other case, the RMSD remains larger than 7 Å. For all the 11 target proteins, our scoring function identifies the experimentally determined native structure as the best structure. Starting from the predicted structure, replica exchange molecular dynamics is performed to further refine the structures. However, we are unable to improve its RMSD toward the experimental structure. The exhaustive sampling by coarse-grained normal mode analysis around the native structures reveals that our function has a linear correlation with RMSDs < 3.0 Å. These results suggest that the function is quite reliable for the protein structure prediction while the sampling method remains one of the major limiting factors in it. The aspects through which the methodology could further be improved are discussed.