Theoretical analysis on changes in thermodynamic quantities upon protein folding: Essential role of hydration

Effect of Conformational Equilibrium on Solvation Properties of 1,2-DCE in Water: A Solvation Thermodynamics and 3D-RISM Study

The contributions of the enthalpy and entropy of solvation for the study of chemical and biological systems are important in the prediction, interpretation, and manipulation of these processes. The relation between solvation Gibbs energies, enthalpies, and entropies of solvation, and their rigorous relation with the conformational equilibrium, are derived for the first time and applied with a computational method, in accordance with the Solvation Thermodynamics previous results, to 1,2-dichloroethane solvation in water. The rigid conformer calculations in solution were performed by using PC+/3D-RISM approach, with the conformational averaged results for enthalpy and solvation Gibbs energy reproducing the experimental results quite successfully. A qualitative agreement in the entropy of solvation predictions was also observed.

Molten globule-like transition state of protein barnase measured with calorimetric force spectroscopy

SignificanceUnderstanding the molecular forces driving the unfolded polypeptide chain to self-assemble into a functional native structure remains an open question. However, identifying the states visited during protein folding (e.g., the transition state between the unfolded and native states) is tricky due to their transient nature. Here, we introduce calorimetric force spectroscopy in a temperature jump optical trap to determine the enthalpy, entropy, and heat capacity of the transition state of protein barnase. We find that the transition state has the properties of a dry molten globule, that is, high free energy and low configurational entropy, being structurally similar to the native state. This experimental single-molecule study characterizes the thermodynamic properties of the transition state in funneled energy landscapes.

Solvent Composition Effects on the Structural Properties of the Aβ42 Monomer from the 3D-RISM-KH Molecular Theory of Solvation

Structural characterization of amyloid (A)β peptides implicated in Alzheimer's disease is a challenging problem due to their intrinsically disordered nature and their high propensity for aggregation. Only limited information is currently available from experiments on conformational properties and aggregation pathways of the peptides in cellular environments. In silico modeling complements experimental information, providing atomistic insight into structure and dynamics of different Aβ species. All-atom explicit solvent molecular dynamics (MD) simulations with a properly selected force field can deliver reliable structural and dynamic information. In the case of intrinsically disordered Aβ peptides, enhanced sampling simulations beyond the nanosecond time scale are required to obtain statistically meaningful results even for simple solvent conditions. To overcome the challenges of conformational sampling in crowded cellular environments, alternative approaches have to be used, including postprocessing of MD data. In this study, we employ the statistical-mechanical, three-dimensional reference interaction site model with the Kovalenko-Hirata closure integral equation molecular theory of solvation to describe solvent composition effects on the conformational equilibrium in a structural ensemble of the Aβ42 (covering residues 1-42) monomer based on a statistical reweighting technique. The methodology enables a computationally efficient prediction on how different factors in the cellular environment, such as solvent composition, nonpolar solvation, and macromolecular crowding, affect the structural properties of the monomer. Similarities have been identified between changes in the structural ensemble caused by nonpolar solvation and crowded environments modeled by ionic solution with large negative ions. In particular, both solvent conditions reduce the random coil content and enhance the helical structure content of the monomer. In contrast to the previous studies, which reported increased α-helical content of peptides in crowded environments, this work attributes these structural features to the difference in solvent exposure of hydrophilic residues of the monomer for different secondary structure elements, rather than to (entropic) excluded volume effects.

Statistical thermodynamics for the unexpectedly large difference between disaccharide stereoisomers in terms of solubility in water

We unravel the physical origins of the large difference between cellobiose and maltose, which consist of two β-1,4 and α-1,4 linked d-glucose units, respectively, in terms of the solubility in water. We construct a thermodynamic theory where the chemical-potential difference between disaccharides in water and in vacuum is identified as the key free-energy function. Its energetic and entropic components are calculated for cellobiose and maltose by statistical-mechanical theories for solute hydration. The disaccharide structures are taken into account at the atomic level and a molecular model is adopted for water. Molecular dynamics simulations are used to account for the conformational fluctuation of a disaccharide molecule, which also enables us to estimate the conformational entropy. We show that the cellobiose/maltose solubility ratio calculated is in good agreement with the experimental value. The solubility becomes much lower for cellobiose due to conformational-entropy and water-entropy effects. The former effect is relevant to higher stability of the intramolecular hydrogen bond between oxygen atoms in the six-membered ring and in the neighboring hydroxyl group: the hydration alters the fluctuation of a molecular conformation to a larger or less regular one, but the degree of this alteration is smaller. The latter effect is attributed to more separation of two hydroxymethyl groups in a molecule, causing lower probability of the overlap of excluded volumes generated by the groups for water molecules. We suggest that physicochemical properties of disaccharides in water become variable depending on the stereoisomerism through hydration effects and the origins of the variety are entropic.

Structure and permeability of ionomers studied by atomistic molecular simulation combined with the theory of solutions in the energy representation

Ionomers play a key role in forming the catalyst layer of polymer electrolyte fuel cells. In the present work, we performed atomistic molecular dynamics simulations and free-energy calculations with the energy-representation method for sulfonated polyethersulfone (SPES) and its derivatives toward the rational design of ionomers for carbon alloy catalysts. It was observed that H₂O aggregates strongly in the branched SPES systems with fluorocarbons and is located homogeneously in the systems without fluorocarbons. The O₂ permeability was then examined within the framework of the solubility-diffusion mechanism. The permeability was seen to be large for the branched SPES with fluorocarbons, indicating that the performance of ionomers as a permeation medium for O₂ may be tuned by the flexibility and branching of the polymer chain.

Statistical Thermodynamics for Actin-Myosin Binding: The Crucial Importance of Hydration Effects

Actomyosin is an important molecular motor, and the binding of actin and myosin is an essential research target in biophysics. Nevertheless, the physical factors driving or opposing the binding are still unclear. Here, we investigate the role of water in actin-myosin binding using the most reliable statistical-mechanical method currently available for assessing biomolecules immersed in water. This method is characterized as follows: water is treated not as a dielectric continuum but as an ensemble of molecules; the polyatomic structures of proteins are taken into consideration; and the binding free energy is decomposed into physically insightful entropic and energetic components by accounting for the hydration effect to its full extent. We find that the actin-myosin binding brings large gains of electrostatic and Lennard-Jones attractive interactions. However, these gains are accompanied by even larger losses of actin-water and myosin-water electrostatic and LJ attractive interactions. Although roughly half of the energy increase due to the losses is cancelled out by the energy decrease arising from structural reorganization of the water released upon binding, the remaining energy increase is still larger than the energy decrease brought by the gains mentioned above. Hence, the net change in system energy is positive, which opposes binding. Importantly, the binding is driven by a large gain of configurational entropy of water, which surpasses the positive change in system energy and the conformational entropy loss occurring for actin and myosin. The principal physical origin of the large water-entropy gain is as follows: the actin-myosin interface is closely packed with the achievement of high shape complementarity on the atomic level, leading to a large increase in the total volume available to the translational displacement of water molecules in the system and a resultant reduction of water crowding (i.e., entropic correlations among water molecules).

Water based on a molecular model behaves like a hard-sphere solvent for a nonpolar solute when the reference interaction site model and related theories are employed

For neutral hard-sphere solutes, we compare the reduced density profile of water around a solute g(r), solvation free energy μ, energy U, and entropy S under the isochoric condition predicted by the two theories: dielectrically consistent reference interaction site model (DRISM) and angle-dependent integral equation (ADIE) theories. A molecular model for water pertinent to each theory is adopted. The hypernetted-chain (HNC) closure is employed in the ADIE theory, and the HNC and Kovalenko-Hirata (K-H) closures are tested in the DRISM theory. We also calculate g(r), U, S, and μ of the same solute in a hard-sphere solvent whose molecular diameter and number density are set at those of water, in which case the radial-symmetric integral equation (RSIE) theory is employed. The dependences of μ, U, and S on the excluded volume and solvent-accessible surface area are analyzed using the morphometric approach (MA). The results from the ADIE theory are in by far better agreement with those from computer simulations available for g(r), U, and μ. For the DRISM theory, g(r) in the vicinity of the solute is quite high and becomes progressively higher as the solute diameter d U increases. By contrast, for the ADIE theory, it is much lower and becomes further lower as d U increases. Due to unphysically positive U and significantly larger |S|, μ from the DRISM theory becomes too high. It is interesting that μ, U, and S from the K-H closure are worse than those from the HNC closure. Overall, the results from the DRISM theory with a molecular model for water are quite similar to those from the RSIE theory with the hard-sphere solvent. Based on the results of the MA analysis, we comparatively discuss the different theoretical methods for cases where they are applied to studies on the solvation of a protein.

SAMPL5: 3D-RISM partition coefficient calculations with partial molar volume corrections and solute conformational sampling

Implicit solvent methods for classical molecular modeling are frequently used to provide fast, physics-based hydration free energies of macromolecules. Less commonly considered is the transferability of these methods to other solvents. The Statistical Assessment of Modeling of Proteins and Ligands 5 (SAMPL5) distribution coefficient dataset and the accompanying explicit solvent partition coefficient reference calculations provide a direct test of solvent model transferability. Here we use the 3D reference interaction site model (3D-RISM) statistical-mechanical solvation theory, with a well tested water model and a new united atom cyclohexane model, to calculate partition coefficients for the SAMPL5 dataset. The cyclohexane model performed well in training and testing ([Formula: see text] for amino acid neutral side chain analogues) but only if a parameterized solvation free energy correction was used. In contrast, the same protocol, using single solute conformations, performed poorly on the SAMPL5 dataset, obtaining [Formula: see text] compared to the reference partition coefficients, likely due to the much larger solute sizes. Including solute conformational sampling through molecular dynamics coupled with 3D-RISM (MD/3D-RISM) improved agreement with the reference calculation to [Formula: see text]. Since our initial calculations only considered partition coefficients and not distribution coefficients, solute sampling provided little benefit comparing against experiment, where ionized and tautomer states are more important. Applying a simple [Formula: see text] correction improved agreement with experiment from [Formula: see text] to [Formula: see text], despite a small number of outliers. Better agreement is possible by accounting for tautomers and improving the ionization correction.

Effects of monohydric alcohols and polyols on the thermal stability of a protein

The thermal stability of a protein is lowered by the addition of a monohydric alcohol, and this effect becomes larger as the size of hydrophobic group in an alcohol molecule increases. By contrast, it is enhanced by the addition of a polyol possessing two or more hydroxyl groups per molecule, and this effect becomes larger as the number of hydroxyl groups increases. Here, we show that all of these experimental observations can be reproduced even in a quantitative sense by rigid-body models focused on the entropic effect originating from the translational displacement of solvent molecules. The solvent is either pure water or water-cosolvent solution. Three monohydric alcohols and five polyols are considered as cosolvents. In the rigid-body models, a protein is a fused hard spheres accounting for the polyatomic structure in the atomic detail, and the solvent is formed by hard spheres or a binary mixture of hard spheres with different diameters. The effective diameter of cosolvent molecules and the packing fractions of water and cosolvent, which are crucially important parameters, are carefully estimated using the experimental data of properties such as the density of solid crystal of cosolvent, parameters in the pertinent cosolvent-cosolvent interaction potential, and density of water-cosolvent solution. We employ the morphometric approach combined with the integral equation theory, which is best suited to the physical interpretation of the calculation result. It is argued that the degree of solvent crowding in the bulk is the key factor. When it is made more serious by the cosolvent addition, the solvent-entropy gain upon protein folding is magnified, leading to the enhanced thermal stability. When it is made less serious, the opposite is true. The mechanism of the effects of monohydric alcohols and polyols is physically the same as that of sugars. However, when the rigid-body models are employed for the effect of urea, its addition is predicted to enhance the thermal stability, which conflicts with the experimental fact. We then propose, as two essential factors, not only the solvent-entropy gain but also the loss of protein-solvent interaction energy upon protein folding. The competition of changes in these two factors induced by the cosolvent addition determines the thermal-stability change.

Statistical thermodynamics of aromatic–aromatic interactions in aqueous solution

To elucidate the interactions between aromatic rings, which are believed to play essential roles in a variety of biological processes, we analyze the water-mediated interactions between toluene molecules along face-to-face stacked (FF) and point-to-face T-shaped (TS) paths using a statistical-mechanical theory of liquids combined with a molecular model for water.

Mechanism of One-to-Many Molecular Recognition Accompanying Target-Dependent Structure Formation: For the Tumor Suppressor p53 Protein as an Example

The new type of molecular recognition, in which an intrinsically disordered region (IDR) of a protein binds to many different target proteins with target-dependent structure formation, is indispensable to the expression of life phenomena and also implicated in a number of diseases. According to the prevailing view, the physicochemical factors responsible for the binding are also target dependent. Here we consider an IDR of the tumor suppressor p53 protein, p53CTD, as an important example related to carcinogenesis and analyze its binding to four targets accompanying the formation of target-dependent structures (i.e., helix, sheet, and two different coils) using our statistical-mechanical method combined with molecular models for water. We find that all of the seemingly different binding processes are driven by a large gain of the translational, configurational entropy of water in the system. The gain originates from sufficiently high shape complementarity on the atomic level within the p53CTD-target interface. It is also required that the electrostatic complementarity be ensured as much as possible to compensate for the dehydration. Such complementarities are achieved in harmony with the portion of the target to which p53CTD binds, leading to a large diversity of structures of p53CTD formed upon binding: If they are not achievable, the binding does not occur. This finding is made possible only by calculating the changes in thermodynamic quantities upon binding and decomposing them into physically insightful components.

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.

Statistical thermodynamics for functionally rotating mechanism of the multidrug efflux transporter AcrB

AcrB, a homotrimer, is the pivotal part of a multidrug efflux pump. A "functionally rotating" picture has been proposed for the drug transport by AcrB, but its mechanism remains unresolved. Here, we investigate the energetics of the whole functional rotation cycle using our theoretical methods. We find that the packing efficiency of AcrB is ununiform, and this ununiformity plays imperative roles primarily through the solvent-entropy effect. When a proton binds to or dissociates from a protomer, the packing properties of this protomer and its two interfaces are perturbed overall in the direction that the solvent translational entropy is lowered. The packing properties of the other two protomers are then reorganized with the recovery or maintenance of closely packed interfaces, so that the solvent-entropy loss can be compensated. The functional structural change by an isolated protomer would cause a seriously large free-energy increase. By forming a trimer, any free-energy increase caused by a protomer is always canceled out by the free-energy decrease brought by the other two protomers via the mechanism mentioned above. The functional structural rotation is thus accomplished using the free-energy decrease arising from the transfer of only a single proton per cycle. The similarities to F1-ATPase are also discussed.

Role of water in ligand binding to maltose-binding protein: insight from a new docking protocol based on the 3D-RISM-KH molecular theory of solvation

Maltose-binding protein is a periplasmic binding protein responsible for transport of maltooligosaccarides through the periplasmic space of Gram-negative bacteria, as a part of the ABC transport system. The molecular mechanisms of the initial ligand binding and induced large scale motion of the protein's domains still remain elusive. In this study, we use a new docking protocol that combines a recently proposed explicit water placement algorithm based on the 3D-RISM-KH molecular theory of solvation and conventional docking software (AutoDock Vina) to explain the mechanisms of maltotriose binding to the apo-open state of a maltose-binding protein. We confirm the predictions of previous NMR spectroscopic experiments on binding modes of the ligand. We provide the molecular details on the binding mode that was not previously observed in the X-ray experiments. We show that this mode, which is defined by the fine balance between the protein-ligand direct interactions and solvation effects, can trigger the protein's domain motion resulting in the holo-closed structure of the maltose-binding protein with the maltotriose ligand in excellent agreement with the experimental data. We also discuss the role of water in blocking unfavorable binding sites and water-mediated interactions contributing to the stability of observable binding modes of maltotriose.

ATP-induced conformational changes of nucleotide-binding domains in an ABC transporter. Importance of the water-mediated entropic force

ATP binding cassette (ABC) proteins belong to a superfamily of active transporters. Recent experimental and computational studies have shown that binding of ATP to the nucleotide binding domains (NBDs) of ABC proteins drives the dimerization of NBDs, which, in turn, causes large conformational changes within the transmembrane domains (TMDs). To elucidate the active substrate transport mechanism of ABC proteins, it is first necessary to understand how the NBD dimerization is driven by ATP binding. In this study, we selected MalKs (NBDs of a maltose transporter) as a representative NBD and calculated the free-energy change upon dimerization using molecular mechanics calculations combined with a statistical thermodynamic theory of liquids, as well as a method to calculate the translational, rotational, and vibrational entropy change. This combined method is applied to a large number of snapshot structures obtained from molecular dynamics simulations containing explicit water molecules. The results suggest that the NBD dimerization proceeds with a large gain of water entropy when ATP molecules bind to the NBDs. The energetic gain arising from direct NBD-NBD interactions is canceled by the dehydration penalty and the configurational-entropy loss. ATP hydrolysis induces a loss of the shape complementarity between the NBDs, which leads to the dissociation of the dimer, due to a decrease in the water-entropy gain and an increase in the configurational-entropy loss. This interpretation of the NBD dimerization mechanism in concert with ATP, especially focused on the water-mediated entropy force, is potentially applicable to a wide variety of the ABC transporters.

Binding of an RNA aptamer and a partial peptide of a prion protein: crucial importance of water entropy in molecular recognition

It is a central issue to elucidate the new type of molecular recognition accompanied by a global structural change of a molecule upon binding to its targets. Here we investigate the driving force for the binding of R12 (a ribonucleic acid aptamer) and P16 (a partial peptide of a prion protein) during which P16 exhibits the global structural change. We calculate changes in thermodynamic quantities upon the R12–P16 binding using a statistical-mechanical approach combined with molecular models for water which is currently best suited to studies on hydration of biomolecules. The binding is driven by a water-entropy gain originating primarily from an increase in the total volume available to the translational displacement of water molecules in the system. The energy decrease due to the gain of R12–P16 attractive (van der Waals and electrostatic) interactions is almost canceled out by the energy increase related to the loss of R12–water and P16–water attractive interactions. We can explain the general experimental result that stacking of flat moieties, hydrogen bonding and molecular-shape and electrostatic complementarities are frequently observed in the complexes. It is argued that the water-entropy gain is largely influenced by the geometric characteristics (overall shapes, sizes and detailed polyatomic structures) of the biomolecules.

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.

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.

Density functional theory for protein transfer free energy

We cast the problem of protein transfer free energy within the formalism of density functional theory (DFT), treating the protein as a source of external potential that acts upon the solvent. Solvent excluded volume, solvent-accessible surface area, and temperature dependence of the transfer free energy all emerge naturally within this formalism, and may be compared with simplified "back of the envelope" models, which are also developed here. Depletion contributions to osmolyte induced stability range from 5 to 10 kBT for typical protein lengths. The general DFT transfer theory developed here may be simplified to reproduce a Langmuir isotherm condensation mechanism on the protein surface in the limits of short-ranged interactions, and dilute solute. Extending the equation of state to higher solute densities results in non-monotonic behavior of the free energy driving protein or polymer collapse. Effective interaction potentials between protein backbone or side chains and TMAO are obtained, assuming a simple backbone/side chain two-bead model for the protein with an effective 6-12 potential with the osmolyte. The transfer free energy δg shows significant entropy: d(δg)/dT ≈ 20 kB for a 100-residue protein. The application of DFT to effective solvent forces for use in implicit-solvent molecular dynamics is also developed. The simplest DFT expressions for implicit-solvent forces contain both depletion interactions and an "impeded-solvation" repulsive force at larger distances.

Extended solvent-contact model for protein solvation: test cases for dipeptides

Solvation effects are critically important in the structural stabilization and functional optimization of proteins. Here, we propose a new solvation free energy function for proteins, and test its applicability in predicting the solvation free energies of dipeptides. The present solvation model involves the improvement of the previous solvent-contact model assuming that the molecular solvation free energy could be given by the sum over the individual atomic contributions. In addition to the existing solvent-contact term, the modified solvation free energy function includes the self-solvation term that reflects the effects of intramolecular interactions in the solute molecule on solute-solvent interactions. Four kinds of atomic parameters should be determined in this solvation model: atomic fragmental volume, maximum atomic occupancy, atomic solvation, and atomic self-solvation parameters. All of these parameters for 16 atom types are optimized with a standard genetic algorithm in such a way to minimize the difference between the solvation free energies of dipeptides obtained from high-level quantum chemical calculations and those predicted by the solvation free energy function. The solvation free energies of dipeptides estimated from the new solvation model are in good agreement with the quantum chemical results. Therefore, the optimized solvation free energy function is expected to be useful for examining the structural and energetic features of proteins in aqueous solution.

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.

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.

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.

Atomic decomposition of the protein solvation free energy and its application to amyloid-beta protein in water

We report the development of an atomic decomposition method of the protein solvation free energy in water, which ascribes global change in the solvation free energy to local changes in protein conformation as well as in hydration structure. So far, empirical decomposition analyses based on simple continuum solvation models have prevailed in the study of protein-protein interactions, protein-ligand interactions, as well as in developing scoring functions for computer-aided drug design. However, the use of continuum solvation model suffers serious drawbacks since it yields the protein free energy landscape which is quite different from that of the explicit solvent model and since it does not properly account for the non-polar hydrophobic effects which play a crucial role in biological processes in water. Herein, we develop an exact and general decomposition method of the solvation free energy that overcomes these hindrances. We then apply this method to elucidate the molecular origin for the solvation free energy change upon the conformational transitions of 42-residue amyloid-beta protein (Aβ42) in water, whose aggregation has been implicated as a primary cause of Alzheimer's disease. We address why Aβ42 protein exhibits a great propensity to aggregate when transferred from organic phase to aqueous phase.

Integral equations in the study of polar and ionic interaction site fluids

In this review article we consider some of the current integral equation approaches and application to model polar liquid mixtures. We consider the use of multidimensional integral equations and in particular progress on the theory and applications of three dimensional integral equations. The IEs we consider may be derived from equilibrium statistical mechanical expressions incorporating a classical Hamiltonian description of the system. We give example including salt solutions, inhomogeneous solutions and systems including proteins and nucleic acids.

Structural and thermodynamic investigations on the aggregation and folding of acylphosphatase by molecular dynamics simulations and solvation free energy analysis

Protein engineering method to study the mutation effects on muscle acylphosphatase (AcP) has been actively applied to describe kinetics and thermodynamics associated with AcP aggregation as well as folding processes. Despite the extensive mutation experiments, the molecular origin and the structural motifs for aggregation and folding kinetics as well as thermodynamics of AcP have not been rationalized at the atomic resolution. To this end, we have investigated the mutation effects on the structures and thermodynamics for the aggregation and folding of AcP by using the combination of fully atomistic, explicit-water molecular dynamics simulations, and three-dimensional reference interaction site model theory. The results indicate that the A30G mutant with the fastest experimental aggregation rate displays considerably decreased α1-helical contents as well as disrupted hydrophobic core compared to the wild-type AcP. Increased solvation free energy as well as hydrophobicity upon A30G mutation is achieved due to the dehydration of hydrophilic side chains in the disrupted α1-helix region of A30G. In contrast, the Y91Q mutant with the slowest aggregation rate shows a non-native H-bonding network spanning the mutation site to hydrophobic core and α1-helix region, which rigidifies the native state protein conformation with the enhanced α1-helicity. Furthermore, Y91Q exhibits decreased solvation free energy and hydrophobicity compared to wild type due to more exposed and solvated hydrophilic side chains in the α1-region. On the other hand, the experimentally observed slower folding rates in both mutants are accompanied by decreased helicity in α2-helix upon mutation. We here provide the atomic-level structures and thermodynamic quantities of AcP mutants and rationalize the structural origin for the changes that occur upon introduction of those mutations along the AcP aggregation and folding processes.

Ion and solvent density distributions around canonical B-DNA from integral equations

We calculate the water and ion spatial distributions around charged oligonucleotides using a renormalized three-dimensional reference interaction site theory coupled with the HNC closure. Our goal is to understand the balance between inter-DNA strand forces and solvation forces as a function of oligonucleotide length in the short strand limit. The DNA is considered in aqueous electrolyte solutions of 1 M KCl, 0.1 M KCl, or 0.1 M NaCl. The current theoretical results are compared to molecular dynamics (MD) simulations and experiments. It is found that the integral equation (IE) theory replicates the MD and the experimental results for the base-specific hydration patterns in both the major and the minor grooves. We are also able to discern characteristic structural pattern differences between Na(+) and K(+) ions. When compared to Poisson-Boltzmann methods, the IE theory, like simulation, predicts a richly structured ion environment, which is better described as multilayer rather than double layer.

Calculation of local water densities in biological systems: a comparison of molecular dynamics simulations and the 3D-RISM-KH molecular theory of solvation

Water plays a unique role in all living organisms. Not only is it nature's ubiquitous solvent, but it also actively takes part in many cellular processes. In particular, the structure and properties of interfacial water near biomolecules such as proteins are often related to the function of the respective molecule. It can therefore be highly instructive to study the local water density around solutes in cellular systems, particularly when solvent-mediated forces such as the hydrophobic effect are relevant. Computational methods such as molecular dynamics (MD) simulations seem well suited to study these systems at the atomic level. However, due to sampling requirements, it is not clear that MD simulations are, indeed, the method of choice to obtain converged densities at a given level of precision. We here compare the calculation of local water densities with two different methods: MD simulations and the three-dimensional reference interaction site model with the Kovalenko-Hirata closure (3D-RISM-KH). In particular, we investigate the convergence of the local water density to assess the required simulation times for different levels of resolution. Moreover, we provide a quantitative comparison of the densities calculated with MD and with 3D-RISM-KH and investigate the effect of the choice of the water model for both methods. Our results show that 3D-RISM-KH yields density distributions that are very similar to those from MD up to a 0.5 Å resolution, but for significantly reduced computational cost. The combined use of MD and 3D-RISM-KH emerges as an auspicious perspective for efficient solvent sampling in dynamical systems.

Towards a universal method for calculating hydration free energies: a 3D reference interaction site model with partial molar volume correction

We report a simple universal method to systematically improve the accuracy of hydration free energies calculated using an integral equation theory of molecular liquids, the 3D reference interaction site model. A strong linear correlation is observed between the difference of the experimental and (uncorrected) calculated hydration free energies and the calculated partial molar volume for a data set of 185 neutral organic molecules from different chemical classes. By using the partial molar volume as a linear empirical correction to the calculated hydration free energy, we obtain predictions of hydration free energies in excellent agreement with experiment (R = 0.94, σ = 0.99 kcal mol (- 1) for a test set of 120 organic molecules).