Multiple hydration layers in cubic insulin crystals

Dipolar Nanodomains in Protein Hydration Shells

The network of hydrogen bonds characteristic of bulk water is significantly disturbed at the protein-water interface, where local fields induce mutually frustrated dipolar domains with potentially novel structure and dynamics. Here the dipolar susceptibility of hydration shells of lysozyme is studied by molecular dynamics simulations in a broad range of temperatures, 140-300 K. The real part of the susceptibility passes through a broad maximum as a function of temperature. The maximum shifts to higher temperatures with increasing frequency of the dielectric experiment. This phenomenology is consistent with that reported for bulk relaxor ferroelectrics, where it is related to the formation of dipolar nanodomains. Nanodomains in the hydration shell extend 12-15 Å from the protein surface into the bulk. Their dynamics are significantly slower than the dynamics of bulk water. The domains dynamically freeze into a ferroelectric glass below 160 K, at which point the Arrhenius plot of the dipolar relaxation time becomes significantly steeper.

Decomposition of protein experimental compressibility into intrinsic and hydration shell contributions

The experimental determination of protein compressibility reflects both the protein intrinsic compressibility and the difference between the compressibility of water in the protein hydration shell and bulk water. We use molecular dynamics simulations to explore the dependence of the isothermal compressibility of the hydration shell surrounding globular proteins on differential contributions from charged, polar, and apolar protein-water interfaces. The compressibility of water in the protein hydration shell is accounted for by a linear combination of contributions from charged, polar, and apolar solvent-accessible surfaces. The results provide a formula for the deconvolution of experimental data into intrinsic and hydration contributions when a protein of known structure is investigated. The physical basis for the model is the variation in water density shown by the surface-specific radial distribution functions of water molecules around globular proteins. The compressibility of water hydrating charged atoms is lower than bulk water compressibility, the compressibility of water hydrating apolar atoms is somewhat larger than bulk water compressibility, and the compressibility of water around polar atoms is about the same as the compressibility of bulk water. We also assess whether hydration water compressibility determined from small compound data can be used to estimate the compressibility of hydration water surrounding proteins. The results, based on an analysis from four dipeptide solutions, indicate that small compound data cannot be used directly to estimate the compressibility of hydration water surrounding proteins.

Structure, dynamics and reactions of protein hydration water

The apparent simplicity of the water molecule belies the wide range of fascinating protein phenomena in which it participates. We review recent computer simulation work on buried, internal water molecules, discussing the thermodynamics of water molecule binding and the participation of water in proton transfer reactions. Surface water molecules are also considered, with emphasis on the modification of average solvent structure on a protein surface, the role of water in the protein dynamical 'glass' transition and a simplified description of the protein motions thereby activated.

Influence of pressure on the retention and separation of insulin variants under linear conditions

The effect of pressure on the retention behavior of insulin variants in RPLC on a YMC-ODS C18 column was investigated under linear conditions. The retention factors of these variants increase nearly 2-fold when the average column pressure is increased from 55 to 250 bar while their separation factors remain nearly unchanged. This effect is explained by a change of the partial molar volume of the insulin variants associated with their adsorption that decreases from -99 to -80 mL/mol for mobile-phase concentrations of acetonitrile increasing from 29 to 33% (v/v). This volume change is much larger than the one observed with low molecular weight compounds. For the same pressure variation, the average number Z of acetonitrile molecules displaced from the protein and the stationary phase upon adsorption increases from 22 to 23.3. The pressure-induced relative increase of the term b[S]/[D0]z (which corresponds to the initial slope of the adsorption isotherm) is approximately twice as large for Lispro than for porcine insulin. Because the binding constant of insulin decreases with increasing pressure, this suggests that the number of binding sites on the stationary phase increases even faster. Finally, it was observed that the column efficiency at flow rates higher than 0.6 mL/min increases slightly with increasing pressure. It is suggested that these observations are also valid for other proteins analyzed in RPLC.

Fast compaction of alpha-lactalbumin during folding studied by stopped-flow X-ray scattering

To monitor the fast compaction process during protein folding, we have used a stopped-flow small-angle X-ray scattering technique combined with a two-dimensional charge-coupled device-based X-ray detector that makes it possible to improve the signal-to-noise ratio of data dramatically, and measured the kinetic refolding reaction of alpha-lactalbumin. The results clearly show that the radius of gyration and the overall shape of the kinetic folding intermediate of alpha-lactalbumin are the same as those of the molten globule state observed at equilibrium. Thus, the identity between the kinetic folding intermediate and the equilibrium molten globule state is firmly established. The present results also suggest that the folding intermediate is more hydrated than the native state and that the hydrated water molecules are dehydrated when specific side-chain packing is formed during the change from the molten globule to the native state.

Is the first hydration shell of lysozyme of higher density than bulk water?

Characterization of the physical properties of protein surface hydration water is critical for understanding protein structure and folding. Here, using molecular dynamics simulation, we provide an explanation of recent x-ray and neutron solution scattering data that indicate that the density of water on the surface of lysozyme is significantly higher than that of bulk water. The simulation-derived scattering profiles are in excellent agreement with the experiment. In the simulation, the 3-A-thick first hydration layer is 15% denser than bulk water. About two-thirds of this increase is the result of a geometric contribution that would also be present if the water was unperturbed from the bulk. The remaining third arises from modification of the water structure and dynamics, involving approximately equal contributions from shortening of the average water-water O-O distance and an increase in the coordination number. Variation in the first hydration shell density is shown to be determined by topographical and electrostatic properties of the protein surface. On average, denser water is found in depressions on the surface in which the water dipoles tend to be aligned parallel to each other by the electrostatic field generated by the protein atoms.

Solvent-exposed tryptophans probe the dynamics at protein surfaces

The dynamics of single tryptophan (W) side chain of protease subtilisin Carlsberg (SC) and myelin basic protein (MBP) were used for probing the surface of these proteins. The W side chains are exposed to the solvent, as shown by the extent of quenching of their fluorescence by KI. Time-resolved fluorescence anisotropy measurements showed that the rotational motion of W is completely unhindered in the case of SC and partially hindered in the case of MBP. The rotational correlation time (phi) associated with the fast local motion of W did not scale linearly with the bulk solvent viscosity (eta) in glycerol-water mixtures. In contrast, phi values of either W side chains in the denatured proteins or the free W scaled almost linearly with eta, as expected by the Stokes-Einstein relationship. These results were interpreted as indicating specific partitioning of water at the surface of the proteins in glycerol-water mixtures.

Disordered water within a hydrophobic protein cavity visualized by x-ray crystallography

Water in the hydrophobic cavity of human interleukin 1beta, which was detected by NMR spectroscopy but was invisible by high resolution x-ray crystallography, has been mapped quantitatively by measurement and phasing of all of the low resolution x-ray diffraction data from a single crystal. Phases for the low resolution data were refined by iterative density modification of an initial flat solvent model outside the envelope of the atomic model. The refinement was restrained by the condition that the map of the difference between the electron density distribution in the full unit cell and that of the atomic model be flat within the envelope of the well ordered protein structure. Care was taken to avoid overfitting the diffraction data by maintaining phases for the high resolution data from the atomic model and by a resolution-dependent damping of the structure factor differences between data and model. The cavity region in the protein could accommodate up to four water molecules. The refined solvent difference map indicates that there are about two water molecules in the cavity region. This map is compatible with an atomic model of the water distribution refined by using XPLOR. About 70% of the time, there appears to be a water dimer in the central hydrophobic cavity, which is connected to the outside by two constricted channels occupied by single water molecules approximately 40% of the time on one side and approximately 10% on the other.

Cluster analysis of consensus water sites in thrombin and trypsin shows conservation between serine proteases and contributions to ligand specificity

Cluster analysis is presented as a technique for analyzing the conservation and chemistry of water sites from independent protein structures, and applied to thrombin, trypsin, and bovine pancreatic trypsin inhibitor (BPTI) to locate shared water sites, as well as those contributing to specificity. When several protein structures are superimposed, complete linkage cluster analysis provides an objective technique for resolving the continuum of overlaps between water sites into a set of maximally dense microclusters of overlapping water molecules, and also avoids reliance on any one structure as a reference. Water sites were clustered for ten superimposed thrombin structures, three trypsin structures, and four BPTI structures. For thrombin, 19% of the 708 microclusters, representing unique water sites, contained water molecules from at least half of the structures, and 4% contained waters from all 10. For trypsin, 77% of the 106 microclusters contained water sites from at least half of the structures, and 57% contained waters from all three. Water site conservation correlated with several environmental features: highly conserved microclusters generally had more protein atom neighbors, were in a more hydrophilic environment, made more hydrogen bonds to the protein, and were less mobile. There were significant overlaps between thrombin and trypsin conserved water sites, which did not localize to their similar active sites, but were concentrated in buried regions including the solvent channel surrounding the Na+ site in thrombin, which is associated with ligand selectivity. Cluster analysis also identified water sites conserved in thrombin but not trypsin, and vice versa, providing a list of water sites that may contribute to ligand discrimination. Thus, in addition to facilitating the analysis of water sites from multiple structures, cluster analysis provides a useful tool for distinguishing between conserved features within a protein family and those conferring specificity.

Protein hydration in solution: experimental observation by x-ray and neutron scattering

The structure of the protein-solvent interface is the subject of controversy in theoretical studies and requires direct experimental characterization. Three proteins with known atomic resolution crystal structure (lysozyme, Escherichia coli thioredoxin reductase, and protein R1 of E. coli ribonucleotide reductase) were investigated in parallel by x-ray and neutron scattering in H2O and D2O solutions. The analysis of the protein-solvent interface is based on the significantly different contrasts for the protein and for the hydration shell. The results point to the existence of a first hydration shell with an average density approximately 10% larger than that of the bulk solvent in the conditions studied. Comparisons with the results of other studies suggest that this may be a general property of aqueous interfaces.

Structure of cubic insulin crystals in glucose solutions

X-ray structures of cubic insulin crystals in high concentrations of glucose at different pH levels and temperatures have been refined to high resolution. We have identified one glucose-binding site near the N-terminus of the A-chain whose occupancy is pH dependent. The effects of reduced water activity on the ordered protein and solvent structures have been examined. Our analysis showed no notable conformational changes in the ordered protein structures or ordered solvent molecules near the protein surface, but the presence of glucose does have a significant effect on the overall density distribution of the bulk solvent in the solvent-accessible volume. We compared the structure of cubic insulin at room temperature and liquid-nitrogen temperature, under identical solvent conditions, using glucose as a cryoprotectant. In this case, we found that the average temperature factor of the protein is reduced and more water molecules can be identified, but there are no significant changes in the protein conformation.

Atomic and residue hydrophilicity in the context of folded protein structures

Water-protein interactions drive protein folding, stabilize the folded structure, and influence molecular recognition and catalysis. We analyzed the closest protein contacts of 10,837 water molecules in crystallographic structures to define a specific hydrophilicity scale reflecting specific rather than bulk solvent interactions. The tendencies of different atom and residue types to be the nearest protein neighbors of bound water molecules correlated with other hydrophobicity scales, verified the relevance of crystallographically determined water positions, and provided a direct experimental measure of water affinity in the context of the folded protein. This specific hydrophilicity was highly correlated with hydrogen-bonding capacity, and correlated better with experimental than computationally derived measures of partitioning between aqueous and organic phases. Atoms with related chemistry clustered with respect to the number of bound water molecules. Neutral and negatively charged oxygen atoms were the most hydrophilic, followed by positively-charged then neutral nitrogen atoms, followed by carbon and sulfur atoms. Agreement between observed side-chain specific hydrophilicity values and values derived from the atomic hydrophilicity scale showed that hydrophilicity values can be synthesized for different functional groups, such as unusual side or main chains, discontinuous epitopes, and drug molecules. Two methods of atomic hydrophilicity analysis provided a measure of complementarity in the interfaces of trypsin:pancreatic trypsin inhibitor and HIV protease:U-75875 inhibitor complexes.

Osmotic pressure method to measure salt induced folding/unfolding of bovine serum albumin

A new approach has been developed to monitor protein folding by utilizing osmotic pressure and a range of salt concentrations in a well characterized protein, bovine serum albumin (BSA). It is hypothesized that both the 'effective' osmotic molecular weight, Ae, and the solute/solvent interaction parameter, I, in the empirical relation Msolvent/Msolute = (RT rho/Ae)1/pi + I [1] can be used as measures of protein folding. I is a measure of solvent perturbed by the solute and is thought to depend directly upon the solvent accessible surface area (ASA). It is reasoned that larger solvent accessible surface area of an unfolded or denatured protein should perturb more water and produce larger I-values. Thus I-values allow calculation of a unfolded protein fraction, fua, due to changes in relative solvent accessible surface area. It has been observed that Ac decreases for filamentous, denatured proteins due to segmental motion of the molecule [2]. This allows calculation of unfolded protein fraction from the effective molecular weight, fum. Colloid osmotic pressure of BSA was measured in a range of salt concentrations at 25 degrees C, and pH = 7 (above the isoelectric point of BSA at pH = 5.4). Both S and I were used to monitor protein folding as the salt concentration was varied. In general, larger and variable I-values and smaller Ae were observed at salt concentrations less than 50 mmolal NaCl (Imax = 8.9), while constant I = 4.1 and Ae = 66,500 were observed above 50 mmolal NaCl. The two expressions for fractional unfolding (fua and fum) are in general agreement. Small differences in the parameters below 50 mmolal salt concentration are explained with well known shifts in the relative amounts of alpha-helix, beta-sheet and random coil in denatured BSA. The relative amounts of these shifts agree with predictions in the literature attributed to continuous BSA expansion rather than an 'all-or-none' conversion.

Structure and dynamics of the water around myoglobin

The interplay between simulations at various levels of hydration and experimental observables has led to a picture of the role of solvent in thermodynamics and dynamics of protein systems. One of the most studied protein-solvent systems is myoglobin, which serves as a paradigm for the development of structure-function relationships in many biophysical studies. We review here some aspects of the solvation of myoglobin and the resulting implications. In particular, recent theoretical and simulation studies unify much of the diverse set of experimental results on water near proteins.

Neutron diffraction analysis of the solvent accessible volume in cubic insulin crystals

The average contact distance between protein and solvent surface atoms in cubic insulin crystals has been determined from two sets of 15 A resolution neutron diffraction data. A contact distance between the water hydrogen sites and the protein surface that is significantly shorter than the average protein-water oxygen contact distance implies that many water molecules are oriented with hydrogen atoms pointed towards the protein surface. The shape of the protein/solvent interface is consistent with the protein envelope obtained from atomic co-ordinates.

A global model of the protein-solvent interface

The solvent structure and dynamics around myoglobin is investigated at the microscopic level of detail by computer simulation. We analyze a molecular dynamics trajectory in terms of solvent mobility and probability distribution. Local events, occurring in the protein-solvent interfacial region, which are often masked by other approaches are thus revealed. Specifically, the local solvent mobility is greatly enhanced for certain locations at the protein surface and in its interior. In addition, a strong correlation between the solvent mobility and density emerges on both global and local scales. We propose a simple model where the solvent distribution measured perpendicularly to the protein surface is utilized to reconstruct the simulated network of hydration within 6 A from the protein surface with a relative error of only 17%. The global precision of this solvation model matches results obtained with more complicated models usually used in refinement procedures in x-ray and neutron experiments but with far fewer parameters. The dramatically improved correspondence between observed and calculated x-ray intensities at low resolution relative to other methods both confirms the validity of the approach used in the MD (molecular dynamics) simulations and allows the results of this study to be implemented in solvent studies on real systems.