Partial molar volume of proteins studied by the three-dimensional reference interaction site model theory

Accurately Computing the Interacted Volume of Molecules over Their 3D Mesh Models

For quickly predicting the rational arrangement of catalysts and substrates, we previously proposed a method to calculate the interacted volumes of molecules over their 3D point cloud models. However, the nonuniform density in molecular point clouds may lead to incomplete contours in some slices, reducing the accuracy of the previous method. In this paper, we propose a two-step method for more accurately computing molecular interacted volumes. First, by employing a prematched mesh slicing method, we layer the 3D triangular mesh models of the electrostatic potential isosurfaces of two molecules globally, transforming the volume calculation into finding the intersecting areas in each layer. Next, by subdividing polygonal edges, we accurately identify intersecting parts within each layer, ensuring precise calculation of interacted volumes. In addition, we present a concise overview for computing intersecting areas in cases of multiple contour intersections and for improving computational efficiency by incorporating bounding boxes at three stages. Experimental results demonstrate that our method maintains high accuracy in different experimental data sets, with an average relative error of 0.16%. On the same experimental setup, our average relative error is 0.07%, which is lower than the previous algorithm's 1.73%, improving the accuracy and stability in calculating interacted volumes.

Theoretical Analysis of the Role of Water in Ligand Binding to Cucurbit[n]uril of Different Sizes

The role of water in host-ligand binding was investigated using a combination of molecular dynamics simulation and three-dimensional reference interaction site model theory. Three different hosts were selected (CB6, CB7, and CB8). Six organic molecules were used as representative ligands: dimethyl sulfoxide (DMSO), N,N-dimethylformamide (DMF), acetone, 2,3-diazabicyclo[2.2.2]oct-2-ene (DBO), cyclopentanone (CPN), and pyrrole. From the binding free energy and its components, we divided the ligands into two groups: those with relatively small molecular size (DMSO, DMF, acetone, and pyrrole) and those with relatively large molecular size (DBO and CPN). We established that the solvent water in the CB6 cavity can be completely displaced by small ligands, resulting in a greater binding affinity compared with larger CBs, except in the case of the small pyrrole ligand, due to outstanding intrinsic properties such as the relatively high hydrophobicity and low dipole moment. In the case of the large ligands, the solvent water can be displaced by DBO and CPN in both CB6 and CB7; there were similar tendencies in their binding affinities, with the greatest affinity in the CB7 complexes. However, the tendencies of the binding affinity components are completely different due to the difference between the complex structure and the solvation structure when a ligand binds with a CB structure. The binding affinities suggest that the size fit between the ligand and CB cannot guarantee the greatest binding affinity gain because the binding structure and intrinsic properties of CB and ligand equally play a crucial role.

gr Predictor: A Deep Learning Model for Predicting the Hydration Structures around Proteins

Among the factors affecting biological processes such as protein folding and ligand binding, hydration, which is represented by a three-dimensional water site distribution function around the protein, is crucial. The typical methods for computing the distribution functions, including molecular dynamics simulations and the three-dimensional reference interaction site model (3D-RISM) theory, require a long computation time ranging from hours to tens of hours. Here, we propose a deep learning (DL) model that rapidly estimates the distribution functions around proteins obtained using the 3D-RISM theory from the protein 3D structure. The distribution functions predicted using our DL model are in good agreement with those obtained using the 3D-RISM theory. Particularly, the coefficient of determination between the distribution function obtained by the DL model and that obtained using the 3D-RISM theory is approximately 0.98. Furthermore, using a graphics processing unit, the prediction by the DL model is completed in less than 1 min, more than 2 orders of magnitude faster than the calculation time of the 3D-RISM theory. The position of water molecules around the protein was estimated based on the distribution function obtained by our DL model, and the position of waters estimated by our DL model was in good agreement with that of water molecules estimated using the 3D-RISM theory and of crystallographic waters. Therefore, our DL model provides a practical and efficient way to calculate the three-dimensional water site distribution functions and to estimate the position of water molecules around the protein. The program called "gr Predictor" is available under the GNU General Public License from https://github.com/YoshidomeGroup-Hydration/gr-predictor.

Hydration Thermodynamics of Familial Parkinson's Disease-Linked Mutants of α-Synuclein

The hydration thermodynamics of different mutants of α-synuclein (α-syn) related to familial Parkinson's disease (PD) is explored using a computational approach that combines both molecular dynamics simulations in water and integral equation theory of molecular liquids. This analysis focuses on the change in conformational entropy, hydration free energy (HFE), and partial molar volume of α-syn upon mutation. The results show that A53T, A30P, E46K, and H50Q mutants aggregate more readily and display increased HFE and less negative interaction volume than the wild-type α-syn. In contrast, an opposite trend is observed for the G51D mutant with a lower experimental aggregation rate. The residuewise decomposition analysis of the HFE highlights that the dehydration/hydration of the hydrophilic residue-rich N- and C-termini of α-syn majorly contributes to the change upon mutation. The hydration shell contributions of different residues to the interaction volume are consistent with its increase/decrease upon mutation. This work shows that both HFE and interaction volume determine the aggregation kinetics of α-syn upon mutation and may serve as an appropriate benchmark for the treatment of PD.

Molecular Recognition and Self-Organization in Life Phenomena Studied by a Statistical Mechanics of Molecular Liquids, the RISM/3D-RISM Theory

There are two molecular processes that are essential for living bodies to maintain their life: the molecular recognition, and the self-organization or self-assembly. Binding of a substrate by an enzyme is an example of the molecular recognition, while the protein folding is a good example of the self-organization process. The two processes are further governed by the other two physicochemical processes: solvation and the structural fluctuation. In the present article, the studies concerning the two molecular processes carried out by Hirata and his coworkers, based on the statistical mechanics of molecular liquids or the RISM/3D-RISM theory, are reviewed.

Interaction Volume Is a Measure of the Aggregation Propensity of Amyloid-β

This study highlights the significance of the partial molar volume of amino acids in predicting the aggregation propensity of an intrinsically disordered protein, amyloid-β (Aβ), and its mutants in aqueous solution. The change in the interaction volume of the protein or mutant is quantitatively correlated with its calculated experimental aggregation propensity. This method also reveals how the interaction volume may be tuned by changing the charge and hydrophobicity of Aβ. While a positive change in the interaction volume and a higher aggregation propensity are observed for mutants with a decrease in the overall charge and/or an increase in hydrophobicity, a reverse trend is observed for the mutants with a decrease in the hydrophobicity and/or an increase in its charge. Hence, the interaction volume may be considered as a key parameter for monitoring protein aggregation that bridges the gap between the experimental aggregation kinetics and solvation thermodynamics.

The geometry of protein hydration

Based on molecular dynamics simulations of four globular proteins in dilute aqueous solution, with three different water models, we examine several, essentially geometrical, aspects of the protein-water interface that remain controversial or incompletely understood. First, we compare different hydration shell definitions, based on spatial or topological proximity criteria. We find that the best method for constructing monolayer shells with nearly complete coverage is to use a 5 Šwater-carbon cutoff and a 4 Šwater-water cutoff. Using this method, we determine a mean interfacial water area of 11.1 Ų which appears to be a universal property of the protein-water interface. We then analyze the local coordination and packing density of water molecules in the hydration shells and in subsets of the first shell. The mean polar water coordination number in the first shell remains within 1% of the bulk-water value, and it is 5% lower in the nonpolar part of the first shell. The local packing density is obtained from additively weighted Voronoi tessellation, arguably the most physically realistic method for allocating space between protein and water. We find that water in all parts of the first hydration shell, including the nonpolar part, is more densely packed than in the bulk, with a shell-averaged density excess of 6% for all four proteins. We suggest reasons why this value differs from previous experimental and computational results, emphasizing the importance of a realistic placement of the protein-water dividing surface and the distinction between spatial correlation and packing density. The protein-induced perturbation of water coordination and packing density is found to be short-ranged, with an exponential decay "length" of 0.6 shells. We also compute the protein partial volume, analyze its decomposition, and argue against the relevance of electrostriction.

Predicting the binding free energy of the inclusion process of 2-hydroxypropyl-β-cyclodextrin and small molecules by means of the MM/3D-RISM method

A protocol to calculate the binding free energy of a host-guest system is proposed based on the MM/3D-RISM method, taking cyclodextrin derivatives and their ligands as model systems. The protocol involves the procedure to identify the most probable binding mode (MPBM) of receptors and ligands by means of the umbrella sampling method. The binding free energies calculated by the MM/3D-RISM method for the complexes of the seven ligands with the MPBM of the cyclodextrin, and with the fluctuated structures around it, are in agreement with the corresponding experimental data in a semi-quantitative manner. It suggests that the protocol proposed here is promising for predicting the binding affinity of a small ligand to a relatively rigid receptor such as cyclodextrin.

Volume and Compressibility of Proteins

The partial specific (or molar) volume, expansibility, and compressibility of a protein are fundamental thermodynamic quantities for characterizing its structure in solution. We review the definitions, measurements, and implications of these volumetric quantities in relation to protein structural biology. The partial specific volumes under constant molality (isomolal) and chemical potential (isopotential) conditions of the cosolvent (multicomponent systems) are explained in terms of preferential solvent interactions relevant to the solubility and stability of proteins. The partial expansibility is briefly discussed in terms of the effects of temperature on protein-solvent interactions (hydration) and internal packing defects (cavities). We discuss the compressibility-structure-function relationships of proteins based on analyses of the correlations between the partial adiabatic compressibilities and the structures or functions of various globular proteins (including mutants), focusing on the roles of the internal cavities in structural fluctuations. The volume and compressibility changes associated with various conformational transitions are also discussed in terms of the changes in hydration and cavities in order to elucidate the nonnative structures and the transition mechanisms, especially those associated with pressure denaturation.

Probing “ambivalent” snug-fit sites in the KcsA potassium channel using three-dimensional reference interaction site model (3D-RISM) theory

The potassium channel is highly selective for K+ over Na+, and the mechanism underlying this selectivity remains unclear. We show the three-dimensional distribution functions (3D-DFs) of small cations (Li+, Na+, and K+) and the free energy profile of ions inside the open selectivity filter (SF) of the KcsA channel. Our previous results [S. Phongphanphanee, N. Yoshida, S. Oiki, F. Hirata. Abstract Book of 5th International Symposium on Molecular Science of Fluctuations toward Biological Functions, P062 (2012)] indicate that the 3D-DF for K+ exhibits distinct peaks at the sites formed by the eight carbonyl oxygen atoms belonging to the surrounding peptide-backbone and residues (the cage site). Li+ has sharp distributions in the 3D-DF at the center of a quadruplex composed of four carbonyl oxygen atoms (the plane site). Na+ has a rather diffuse distribution throughout the SF region with peaks both in the plane and in cage sites. The results provide microscopic evidence of the phenomenological findings that Li+ and Na+ are not excluded from the SF region and that the binding affinity alone does not cause the ion selectivity of KcsA. In the present study, with an ion placed explicitly along the pore axis, the free energy profiles of the ions inside the SF were calculated; from these profiles we suggest a new mechanism for selective K+ permeation. According to the model, a K+ ion must overcome a free energy barrier that is approximately half that of Na+ to exit from either of the SF mouths due to the existence of an intermediate local minimum along the route for climbing the barriers.

Analysis of biomolecular solvation sites by 3D-RISM theory

We derive, implement, and apply equilibrium solvation site analysis for biomolecules. Our method utilizes 3D-RISM calculations to quickly obtain equilibrium solvent distributions without either necessity of simulation or limits of solvent sampling. Our analysis of these distributions extracts highest likelihood poses of solvent as well as localized entropies, enthalpies, and solvation free energies. We demonstrate our method on a structure of HIV-1 protease where excellent structural and thermodynamic data are available for comparison. Our results, obtained within minutes, show systematic agreement with available experimental data. Further, our results are in good agreement with established simulation-based solvent analysis methods. This method can be used not only for visual analysis of active site solvation but also for virtual screening methods and experimental refinement.

Placevent: an algorithm for prediction of explicit solvent atom distribution-application to HIV-1 protease and F-ATP synthase

We have created a simple algorithm for automatically predicting the explicit solvent atom distribution of biomolecules. The explicit distribution is coerced from the three-dimensional (3D) continuous distribution resulting from a 3D reference interaction site model (3D-RISM) calculation. This procedure predicts optimal location of solvent molecules and ions given a rigid biomolecular structure and the solvent composition. We show examples of predicting water molecules near the KNI-272 bound form of HIV-1 protease and predicting both sodium ions and water molecules near the rotor ring of F-adenosine triphosphate (ATP) synthase. Our results give excellent agreement with experimental structure with an average prediction error of 0.39-0.65 Å. Further, unlike experimental methods, this method does not suffer from the partial occupancy limit. Our method can be performed directly on 3D-RISM output within minutes. It is extremely useful for examining multiple specific solvent-solute interactions, as a convenient method for generating initial solvent structures for molecular dynamics calculations, and may assist in refinement of experimental structures. © 2012 Wiley Periodicals, Inc.

Correlation of second virial coefficient with solubility for proteins in salt solutions

In this work, osmotic second virial coefficients (B(22)) were determined and correlated with the measured solubilities for the proteins, α-amylase, ovalbumin, and lysozyme. The B(22) values and solubilities were determined in similar solution conditions using two salts, sodium chloride and ammonium sulfate in an acidic pH range. An overall decrease in the solubility of the proteins (salting out) was observed at high concentrations of ammonium sulfate and sodium chloride solutions. However, for α-amylase, salting-in behavior was also observed in low concentration sodium chloride solutions. In ammonium sulfate solutions, the B(22) are small and close to zero below 2.4 M. As the ammonium sulfate concentrations were further increased, B(22) values decreased for all systems studied. The effect of sodium chloride on B(22) varies with concentration, solution pH, and the type of protein studied. Theoretical models show a reasonable fit to the experimental derived data of B(22) and solubility. B(22) is also directly proportional to the logarithm of the solubility values for individual proteins in salt solutions, so the log-linear empirical models developed in this work can also be used to rapidly predict solubility and B(22) values for given protein-salt systems.

Non-bulk-like solvent behavior in the ribosome exit tunnel

As nascent proteins are synthesized by the ribosome, they depart via an exit tunnel running through the center of the large subunit. The exit tunnel likely plays an important part in various aspects of translation. Although water plays a key role in many bio-molecular processes, the nature of water confined to the exit tunnel has remained unknown. Furthermore, solvent in biological cavities has traditionally been characterized as either a continuous dielectric fluid, or a discrete tightly bound molecule. Using atomistic molecular dynamics simulations, we predict that the thermodynamic and kinetic properties of water confined within the ribosome exit tunnel are quite different from this simple two-state model. We find that the tunnel creates a complex microenvironment for the solvent resulting in perturbed rotational dynamics and heterogenous dielectric behavior. This gives rise to a very rugged solvation landscape and significantly retarded solvent diffusion. We discuss how this non-bulk-like solvent is likely to affect important biophysical processes such as sequence dependent stalling, co-translational folding, and antibiotic binding. We conclude with a discussion of the general applicability of these results to other biological cavities.

Intrinsically disordered regions may lower the hydration free energy in proteins: a case study of nudix hydrolase in the bacterium Deinococcus radiodurans

The proteome of the radiation- and desiccation-resistant bacterium D. radiodurans features a group of proteins that contain significant intrinsically disordered regions that are not present in non-extremophile homologues. Interestingly, this group includes a number of housekeeping and repair proteins such as DNA polymerase III, nudix hydrolase and rotamase. Here, we focus on a member of the nudix hydrolase family from D. radiodurans possessing low-complexity N- and C-terminal tails, which exhibit sequence signatures of intrinsic disorder and have unknown function. The enzyme catalyzes the hydrolysis of oxidatively damaged and mutagenic nucleotides, and it is thought to play an important role in D. radiodurans during the recovery phase after exposure to ionizing radiation or desiccation. We use molecular dynamics simulations to study the dynamics of the protein, and study its hydration free energy using the GB/SA formalism. We show that the presence of disordered tails significantly decreases the hydration free energy of the whole protein. We hypothesize that the tails increase the chances of the protein to be located in the remaining water patches in the desiccated cell, where it is protected from the desiccation effects and can function normally. We extrapolate this to other intrinsically disordered regions in proteins, and propose a novel function for them: intrinsically disordered regions increase the "surface-properties" of the folded domains they are attached to, making them on the whole more hydrophilic and potentially influencing, in this way, their localization and cellular activity.

Molecular recognition in biomolecules studied by statistical-mechanical integral-equation theory of liquids

Recent progress in the theory of molecular recognition in biomolecules is reviewed, which has been made based on the statistical mechanics of liquids or the RISM/3D-RISM theory during the last five years in the authors' group. The method requires just the structure of protein and the potential energy parameters for the biomolecules and solutions as inputs. The calculation is carried out in two steps. The first step is to obtain the pair correlation functions for solutions consisting of water and ligands based on the RISM theory. Then, given the pair correlation functions prepared in the first step, we calculate the 3D-distribution functions of water and ligands around and inside protein based on the 3D-RISM theory. The molecular recognition of a ligand by the protein is realized by the 3D-distribution functions: if one finds some conspicuous peaks in the distribution of a ligand inside protein, then the ligand is regarded as "recognized" by the protein. Some biochemical processes are investigated, which are intimately related to the molecular recognition of small ligands including water, noble gases, and ions by a protein.

Microscopic Understanding of Preferential Exclusion of Compatible Solute Ectoine: Direct Interaction and Hydration Alteration

Ectoine, a zwitterionic compatible solute (CS), acts as an effective stabilizer of protein function. Using molecular dynamics simulation, solvent spatial distributions around both met-enkephalin (M-Enk) and chymotrypsin inhibitor 2 (CI2) were investigated at the molecular level in ectoine aqueous solution. An unexpected finding was that ectoine exhibits preferential binding, as an overall tendency, around both peptides. However, with the aid of the surficial Kirkwood-Buff parameter, it was clearly shown that the preferential exclusion of ectoine from the peptide surface was weaker in the smaller M-Enk than in the larger CI2. It is concluded that a denser and more structured hydration layer, such as that developed on the surface of CI2, is an important factor in the exclusion of ectoine.

Theoretical study of the partial molar volume change associated with the pressure-induced structural transition of ubiquitin

The partial molar volume (PMV) change associated with the pressure-induced structural transition of ubiquitin is analyzed by the three-dimensional reference interaction site model (3D-RISM) theory of molecular solvation. The theory predicts that the PMV decreases upon the structural transition, which is consistent with the experimental observation. The volume decomposition analysis demonstrates that the PMV reduction is primarily caused by the decrease in the volume of structural voids in the protein, which is partially canceled by the volume expansion due to the hydration effects. It is found from further analysis that the PMV reduction is ascribed substantially to the penetration of water molecules into a specific part of the protein. Based on the thermodynamic relation, this result implies that the water penetration causes the pressure-induced structural transition. It supports the water penetration model of pressure denaturation of proteins proposed earlier.

Theoretical study of volume changes accompanying xenon-lysozyme binding: implications for the molecular mechanism of pressure reversal of anesthesia

The change in partial molar volume (PMV) accompanying the xenon-lysozyme binding was investigated for elucidating the molecular mechanism of the pressure reversal of general anesthesia, using the three-dimensional reference interaction site model theory of molecular solvation. An increase of the PMV from xenon binding to the substrate binding site of lysozyme was found, and the binding is suppressed by pressure, while the internal site binding did not change the PMV. The PMV change was analyzed by decomposing it into several contributions from geometry and hydration. We also analyzed the hydration change due to the binding. From the results, we draw a molecular picture of the PMV change accompanying xenon-lysozyme binding, which gives a possible mechanism of pressure reversal of anesthesia.

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

The free energy change associated with the coil-to-native structural transition of protein G in aqueous solution is calculated by using the molecular theory of solvation, also known as the three-dimensional reference interaction site model theory, to uncover the molecular mechanism of protein folding. The free energy is decomposed into the protein intramolecular energy, the hydration energy, and the hydration entropy. The folding is accompanied with a large gain in the protein intramolecular energy. However, it is almost canceled by the correspondingly large loss in the hydration energy due to the dehydration, resulting in the total energy gain about an order of magnitude smaller than might occur in vacuum. The hydration entropy gain is found to be a substantial driving force in protein folding. It is comparable with or even larger than the total energy gain. The total energy gain coupled with the hydration entropy gain is capable of suppressing the conformational entropy loss in the folding. Based on careful analysis of the theoretical results, the authors present a challenging physical picture of protein folding where the overall folding process is driven by the water entropy effect.

Locating missing water molecules in protein cavities by the three-dimensional reference interaction site model theory of molecular solvation

Water molecules confined in protein cavities are of great importance in understanding the protein structure and functions. However, it is a nontrivial task to locate such water molecules in protein by the ordinary molecular simulation and modeling techniques as well as experimental methods. The present study proves that the three-dimensional reference interaction site model (3D-RISM) theory, a recently developed statistical-mechanical theory of molecular solvation, has an outstanding advantage in locating such water molecules. In this paper, we demonstrate that the 3D-RISM theory is able to reproduce the structure and the number of water molecules in cavities of hen egg-white lysozyme observed commonly in the X-ray structures of different resolutions and conditions. Furthermore, we show that the theory successfully identified a water molecule in a cavity, the existence of which has been ambiguous even from the X-ray results. In contrast, we confirmed that molecular dynamics simulation is helpless at present to find such water molecules because the results substantially depend on the initial coordinates of water molecules. Possible applications of the theory to problems in the fields of biochemistry and biophysics are also discussed.

Theoretical study of the cosolvent effect on the partial molar volume change of staphylococcal nuclease associated with pressure denaturation

We explain the molecular mechanism of the effect of urea and glycerol cosolvents on the partial molar volume (PMV) change associated with the pressure denaturation of staphylococcal nuclease (SNase) protein recently observed in experiments. Native and denatured conformations of SNase are produced by using molecular dynamics simulations in water, and the PMV is obtained from the integral equation theory of molecular liquids called 3D-RISM, which is based on statistical mechanics. The PMV of the native SNase in water predicted by 3D-RISM theory is in good agreement with experiment. The PMV changes associated with pressure denaturation in water and in water-urea and water-glycerol mixtures are qualitatively reproduced. By analyzing the results obtained, we found two interesting cosolvent effects on the PMV: (1) both urea and glycerol cosolvents increase the PMVs of both native and denatured SNase compared to those in water and (2) both urea and glycerol cosolvents increase the PMV of denatured SNase more than that of native SNase. We also showed that these two observations can be explained in terms of the thermal volume, which is related to the packing effect of solvent molecules.

Partial molar volume of anionic polyelectrolytes in aqueous solution

In this work the partial molar volumes (V) of different anionic polyelectrolytes and hydrophobically modified polyelectrolytes (PHM) were measured. Polymers like polymaleic acid-co-styrene, polymaleic acid-co-1-olefin, polymaleic acid-co-vinyl-2-pyrrolidone, and polyacrylic acid (abbreviated as MAS-n, PA-n-K2, AMVP, and PAA, respectively) were employed. These materials were investigated by density measurements in highly dilute aqueous solutions. The molar volume results allow us to discuss the effect of the carboxylic groups and the contributions from the comonomeric principal chain. The PAA presents the smaller V, while the largest V value was for AMVP. The V of PHM shows a linear relationship with the number of methylene groups in the lateral chain. It is found that the magnitude of the contribution per methylene group decreases as the hydrophobic character of the environment increases.

A theoretical analysis on hydration thermodynamics of proteins

The hydration free energy (HFE) of several proteins modeled using the all-atom force field is calculated by employing the three-dimensional reference interaction site model theory, a recently developed integral equation theory of molecular solvation. The HFE is decomposed into the energetic and entropic components under the isochoric condition. The former comprises the protein-water interaction energy and the water reorganization energy arising from the structural changes induced in water. Each component is further decomposed into the nonelectrostatic and electrostatic contributions. It is found that the HFE is governed by the nonelectrostatic hydration entropy and the electrostatic hydration energy. The nonelectrostatic hydration entropy is almost exclusively ascribed to the translational entropy loss of water upon the protein insertion. It asymptotically becomes proportional to the excluded volume (EV) for water molecules as the protein size increases. The hydration energy is determined by the protein-water interaction energy which is half compensated by the water reorganization energy. These energy terms are approximately proportional to the water-accessible surface area (ASA). The energetic and entropic contributions are balanced with each other and the HFE has no apparent linear relation with the EV and ASA.

Understanding high pressure stability of helical conformation of oligopeptides and helix bundle protein

The pressure effect on conformational equilibria of simple organic compounds and the pressure denaturation of proteins have been well investigated by using vibrational spectroscopy. However, there was no systematic investigation of the pressure effect on conformational equilibria of oligopeptides, which are located between the simple organic compounds and proteins. Here, we review the recent vibrational spectroscopic and theoretical studies of the pressure effect on conformational equilibria of model oligopeptides and helix bundle protein in aqueous solution.

Water molecules in a protein cavity detected by a statistical-mechanical theory

Four water molecules confined in a small cavity of hen egg-white lysozyme were detected by means of the three-dimensional (3D) RISM theory, a statistical-mechanical theory of molecular solutions. This is the first theoretical realization of confined molecules in a protein without making nonsense tricks, such as placing the molecules in the space a priori. Possible impacts which the result may have on biochemistry and biophysics, including the molecular recognition, enzymatic reactions, etc., are discussed.

Theoretical study of volume changes associated with the helix-coil transition of an alanine-rich peptide in aqueous solution

The changes in the partial molar volume (PMV) associated with the conformational transition of an alanine-rich peptide AK16 from the alpha-helix structure to various random coil structures are calculated by the three-dimensional interaction site model (3D-RISM) theory coupled with the Kirkwood-Buff theory. The volume change is analyzed by decomposing it into contributions from geometry and hydration: the changes in the van der Waals, void, thermal, and interaction volume. The total change in the PMV is positive. This is primarily due to the growth of void space within the peptide, which is canceled in part by the volume reduction resulting from the increase in the electrostatic interaction between the peptide and water molecules. The changes in the void and thermal volume of the coil structures are widely distributed and tend to compensate each other. Additionally, the relations between the hydration volume components and the surface properties are investigated. We categorize coil structures into extended coils with the PMV smaller than helix and general coils with the PMV larger than helix. The pressure therefore can both stabilize and destabilize the coil structures. The latter seems to be a more proper model of random coil structures of the peptide.