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

Pressure and temperature jump relaxation kinetics of the conformational change in Salmonella typhimurium tryptophan synthase L-serine complex: large activation compressibility and heat capacity changes demonstrate the contribution of solvation

Tryptophan synthase is an alpha2beta2 multienzyme complex that exhibits coupling of the alpha- and beta-subunit reactions by tightly controlled allosteric interactions. A wide range of parameters can affect the allosteric interactions, including monovalent cations, pH, alpha-site and beta-site ligands, temperature, and pressure. Rapid changes in hydrostatic pressure (P-jump) and temperature (T-jump) were used to examine the effects of pressure and temperature on the rates of the interconversion of external aldimine and aminoacrylate intermediates in the Tryptophan synthase-L-Ser complex. The intense fluorescence emission of the Tryptophan synthase L-Ser external aldimine complex at 495 nm, with 420 nm excitation, provides a probe of the conformational state of Trp synthase. P-jump measurements allowed the determination of rate constants for the reactions in the presence of Na(+), Na(+) with benzimidazole (BZI), and NH4(+). The data require a compressibility term, beta(o)(double dagger), to obtain good fits, especially for the NH4(+) and BZI/Na(+) data. The compressibility changes are consistent with changes in solvation in the transition state. The transition state for the relaxation is more similar in volume to the closed aminoacrylate complex in the presence of Na(+), while it is more similar to the open external aldimine in the presence of NH4(+). Differences between the relaxations for positive and negative P-jumps may arise from changing relative populations of microstates with pressure. T-jump experiments of the Na(+) form of the tryptophan synthase-L-Ser complex show large changes in rate and amplitude over the temperature range from 7 to 45 degrees C. The Arrhenius plots show strong curvature, and hence require a heat capacity term, DeltaC(p)(double dagger), to obtain good fits. The values of DeltaC(p)(double dagger) are very large and negative (-3.6 to -4.4 kJ mol(-1) K(-1)). These changes are also consistent with large changes in solvation in the transition state for interconversion of external aldimine and aminoacrylate intermediates in the Tryptophan synthase-L-Ser complex.

Hydration dynamics in a partially denatured ensemble of the globular protein human alpha-lactalbumin investigated with molecular dynamics simulations

Atomistic molecular dynamics simulations are used to probe changes in the nature and subnanosecond dynamical behavior of solvation waters that accompany partial denaturation of the globular protein, human alpha-lactalbumin. A simulated ensemble of subcompact conformers, similar to the molten globule state of human alpha-lactalbumin, demonstrates a marginal increase in the amount of surface solvation relative to the native state. This increase is accompanied by subtle but distinct enhancement in surface water dynamics, less favorable protein-water interactions, and a marginal decrease in the anomalous behavior of solvation water dynamics. The extent of solvent influx is not proportional to the increased surface area, and the partially denatured conformers are less uniformly solvated compared to their native counterpart. The observed solvation in partially denatured conformers is lesser in extent compared to earlier experimental estimates in molten globule states, and is consistent with more recent descriptions based on nuclear magnetic relaxation dispersion studies.

Thickness of the hydration layer of a protein from molecular dynamics simulation

Water molecules around a protein exhibit slow dynamics with respect to that of pure bulk water. One important issue in protein hydration is the thickness of the hydration layer (i.e., the distance from the protein surface up to which the water dynamics is influenced by the protein). Estimation of thickness is crucial to understand better the properties of "biological water" and the role that it plays in guiding the protein's function. We have performed an atomistic molecular dynamics simulation of an aqueous solution of the protein villin headpiece subdomain or HP-36 to estimate the thickness of its hydration water. In particular, several dynamical properties of water around different segments (three alpha-helices) of the protein have been calculated by varying the thickness of the hydration layers. It is found that in general the influence of the helices on water properties extends beyond the first hydration layer. However, the heterogeneous nature of water among the first hydration layers of the three helices diminishes as the thickness is increased. It indicates that, for a small protein such as HP-36, the thickness of "biological water" is uniform for different segments of the protein.

Dynamics of water in the hydration layer of a partially unfolded structure of the protein HP-36

Atomistic molecular dynamics simulations of the folded native structure and a partially unfolded molten globule structure of the protein villin headpiece subdomain or HP-36 have been carried out with explicit solvent to explore the effects of unfolding on the dynamical behavior of water present in the hydration layers of different segments (three alpha-helices) of the protein. The calculations revealed that the unfolding of helix-2 influences the translational and rotational motions of water present in the hydration layers of the three helices in a heterogeneous manner. It is observed that a correlation exists between the unfolding of helix-2 and the microscopic kinetics of protein-water hydrogen bonds formed by its residues. This in turn has an influence on the rigidity of the hydration layers of the helices in the unfolded structure versus that in the folded native structure. These results should provide a microscopic explanation to recent solvation dynamics experiments on folded native and unfolded structures of proteins.

Molecular dynamics simulations for water and ions in protein crystals

The spatial and temporal properties of water and ions in bionanoporous materials-protein crystals-have been investigated using molecular dynamics simulations. Three protein crystals are considered systematically with different morphologies and chemical topologies: tetragonal lysozyme, orthorhombic lysozyme, and tetragonal thermolysin. It is found that the thermal fluctuations of C(alpha) atoms in the secondary structures of protein molecules are relatively weak due to hydrogen bonding. The solvent-accessible surface area per residue is nearly identical in the three protein crystals; the hydrophobic and hydrophilic residues in each crystal possess approximately the same solvent-accessible surface area. Water distributes heterogeneously and has different local structures within the biological nanopores of the three protein crystals. The mobility of water and ions in the crystals is enhanced as the porosity increases and also by the fluctuations of protein atoms particularly in the two lysozyme crystals. Anisotropic diffusion is found preferentially along the pore axis, as experimentally observed. The anisotropy of the three crystals increases in the order: tetragonal thermolysin < tetragonal lysozyme < orthorhombic lysozyme.

X-ray solution scattering (SAXS) combined with crystallography and computation: defining accurate macromolecular structures, conformations and assemblies in solution

Crystallography supplies unparalleled detail on structural information critical for mechanistic analyses; however, it is restricted to describing low energy conformations of macromolecules within crystal lattices. Small angle X-ray scattering (SAXS) offers complementary information about macromolecular folding, unfolding, aggregation, extended conformations, flexibly linked domains, shape, conformation, and assembly state in solution, albeit at the lower resolution range of about 50 A to 10 A resolution, but without the size limitations inherent in NMR and electron microscopy studies. Together these techniques can allow multi-scale modeling to create complete and accurate images of macromolecules for modeling allosteric mechanisms, supramolecular complexes, and dynamic molecular machines acting in diverse processes ranging from eukaryotic DNA replication, recombination and repair to microbial membrane secretion and assembly systems. This review addresses both theoretical and practical concepts, concerns and considerations for using these techniques in conjunction with computational methods to productively combine solution scattering data with high-resolution structures. Detailed aspects of SAXS experimental results are considered with a focus on data interpretation tools suitable to model protein and nucleic acid macromolecular structures, including membrane protein, RNA, DNA, and protein-nucleic acid complexes. The methods discussed provide the basis to examine molecular interactions in solution and to study macromolecular flexibility and conformational changes that have become increasingly relevant for accurate understanding, simulation, and prediction of mechanisms in structural cell biology and nanotechnology.

Effects of surface water on protein dynamics studied by a novel coarse-grained normal mode approach

Normal mode analysis (NMA) has received much attention as a direct approach to extract the collective motions of macromolecules. However, the stringent requirement of computational resources by classical all-atom NMA limits the size of the macromolecules to which the method is normally applied. We implemented a novel coarse-grained normal mode approach based on partitioning the all-atom Hessian matrix into relevant and nonrelevant parts. It is interesting to note that, using classical all-atom NMA results as a reference, we found that this method generates more accurate results than do other coarse-grained approaches, including elastic network model and block normal mode approaches. Moreover, this new method is effective in incorporating the energetic contributions from the nonrelevant atoms, including surface water molecules, into the coarse-grained protein motions. The importance of such improvements is demonstrated by the effect of surface water to shift vibrational modes to higher frequencies and by an increase in overlap of the coarse-grained eigenvector space (the motion directions) with that obtained from molecular dynamics simulations of solvated protein in a water box. These results not only confirm the quality of our method but also point out the importance of incorporating surface structural water in studying protein dynamics.

An extended dynamical hydration shell around proteins

The focus in protein folding has been very much on the protein backbone and sidechains. However, hydration waters make comparable contributions to the structure and energy of proteins. The coupling between fast hydration dynamics and protein dynamics is considered to play an important role in protein folding. Fundamental questions of protein hydration include, how far out into the solvent does the influence of the biomolecule reach, how is the water affected, and how are the properties of the hydration water influenced by the separation between protein molecules in solution? We show here that Terahertz spectroscopy directly probes such solvation dynamics around proteins, and determines the width of the dynamical hydration layer. We also investigate the dependence of solvation dynamics on protein concentration. We observe an unexpected nonmonotonic trend in the measured terahertz absorbance of the five helix bundle protein lambda(6-85)* as a function of the protein: water molar ratio. The trend can be explained by overlapping solvation layers around the proteins. Molecular dynamics simulations indicate water dynamics in the solvation layer around one protein to be distinct from bulk water out to approximately 10 A. At higher protein concentrations such that solvation layers overlap, the calculated absorption spectrum varies nonmonotonically, qualitatively consistent with the experimental observations. The experimental data suggest an influence on the correlated water network motion beyond 20 A, greater than the pure structural correlation length usually observed.

Temperature dependence of protein-hydration hydrodynamics by molecular dynamics simulations

The dynamics of water molecules near the protein surface are different from those of bulk water and influence the structure and dynamics of the protein itself. To elucidate the temperature dependence hydration dynamics of water molecules, we present results from the molecular dynamic simulation of the water molecules surrounding two proteins (Carboxypeptidase inhibitor and Ovomucoid) at seven different temperatures (T=273 to 303 K, in increments of 5 K). Translational diffusion coefficients of the surface water and bulk water molecules were estimated from 2 ns molecular dynamics simulation trajectories. Temperature dependence of the estimated bulk water diffusion closely reflects the experimental values, while hydration water diffusion is retarded significantly due to the protein. Protein surface induced scaling of translational dynamics of the hydration waters is uniform over the temperature range studied, suggesting the importance protein-water interactions.

Dehydration-driven solvent exposure of hydrophobic surfaces as a driving force in peptide folding

Recent work has shown that the nature of hydration of pure hydrophobic surfaces changes with the length scale considered: water hydrogen-bonding networks adapt to small exposed hydrophobic species, hydrating or "wetting" them at relatively high densities, whereas larger hydrophobic areas are "dewetted" [Chandler D (2005), Nature 29:640-647]. Here we determine whether this effect is also present in peptides by examining the folding of a beta-hairpin (the 14-residue amyloidogenic prion protein H1 peptide), using microsecond time-scale molecular dynamics simulations. Two simulation models are compared, one explicitly including the water molecules, which may thus adapt locally to peptide configurations, and the other using a popular continuum approximation, the generalized Born/surface area implicit solvent model. The results obtained show that, in explicit solvent, peptide conformers with high solvent-accessible hydrophobic surface area indeed also have low hydration density around hydrophobic residues, whereas a concomitant higher hydration density around hydrophilic residues is observed. This dewetting effect stabilizes the fully folded beta-hairpin state found experimentally. In contrast, the implicit solvent model destabilizes the fully folded hairpin, tending to cluster hydrophobic residues regardless of the size of the exposed hydrophobic surface. Furthermore, the rate of the conformational transitions in the implicit solvent simulation is almost doubled with respect to that of the explicit solvent. The results suggest that dehydration-driven solvent exposure of hydrophobic surfaces may be a significant factor determining peptide conformational equilibria.

Simulation studies of the protein-water interface. I. Properties at the molecular resolution

We report molecular dynamics simulations of three globular proteins: ubiquitin, apo-calbindin D(9K), and the C-terminal SH2 domain of phospholipase C-gamma1 in explicit water. The proteins differ in their overall charge and fold type and were chosen to represent to some degree the structural variability found in medium-sized proteins. The length of each simulation was at least 15 ns, and larger than usual solvent boxes were used. We computed radial distribution functions, as well as orientational correlation functions about the surface residues. Two solvent shells could be clearly discerned about charged and polar amino acids. Near apolar amino acids the water density near such residues was almost devoid of structure. The mean residence time of water molecules was determined for water shells about the full protein, as well as for water layers about individual amino acids. In the dynamic properties, two solvent shells could be characterized as well. However, by comparison to simulations of pure water it could be shown that the influence of the protein reaches beyond 6 A, i.e., beyond the first two shells. In the first shell (r < or =3.5 A), the structural and dynamical properties of solvent waters varied considerably and depended primarily on the physicochemical properties of the closest amino acid side chain, with which the waters interact. By contrast, the solvent properties seem not to depend on the specifics of the protein studied (such as the net charge) or on the secondary structure element in which an amino acid is located. While differing considerably from the neat liquid, the properties of waters in the second solvation shell (3.5< r < or =6 A) are rather uniform; a direct influence from surface amino acids are already mostly shielded.

Coupling between hydration layer dynamics and unfolding kinetics of HP-36

We have performed atomistic molecular dynamics simulations of aqueous solutions of HP-36 at 300 K in its native state, as well as at high temperatures to explore the unfolding dynamics of the protein and its correlation with the motion of water around it. On increasing the temperature a partially unfolded molten globule state is formed where the smallest alpha helix (helix 2) unfolds into a coil. It is observed that the unfolding is initiated around the residue Phe-18 which shows a sharp displacement during unfolding. We have noticed that the unfolding of the protein affects the density of water near the protein surface. Besides, the dynamics of water in the protein hydration layer has been found to be strongly correlated with the time evolution of the unfolding process. We have introduced and calculated a displacement time correlation function to monitor the change in water motion relative to the protein backbone during unfolding. We find that the unfolding of helix 2 is associated with an increase in mobility of water around it as compared to water around the other two helices. We have also explored the microscopic aspects of secondary structure specific and site specific solvation dynamics of the protein. The calculations reveal that unfolding influences the solvation dynamics of the protein molecule in a heterogeneous manner depending on the location of the polar probe residues. This seems to be in agreement with recent experimental findings.

Coupling between lysozyme and trehalose dynamics: microscopic insights from molecular-dynamics simulations

We have carried out molecular-dynamics simulations on fully flexible all-atom models of the protein lysozyme immersed in trehalose, an effective biopreservative, with the purpose of exploring the nature and extent of the dynamical coupling between them. Our study shows a strong coupling over a wide range of temperatures. We found that the onset of anharmonic behavior was dictated by changes in the dynamics and relaxation processes in the trehalose glass. The physical origin of protein-trehalose coupling was traced to the hydrogen bonds formed at the interface between the protein and the solvent. Moreover, protein-solvent hydrogen bonding was found to control the structural relaxation of the protein. The dynamics of the protein was found to be heterogeneous; the motions of surface and core atoms had different dependencies on temperature and, in addition, the surface atoms were more sensitive to the dynamics of the solvent than the core atoms. From the solvent perspective we found that the dynamics near the protein surface showed an unexpected enhanced mobility compared to the bulk. These results shed some light on the microscopic origins of the dynamical coupling in protein-solvent systems.

Correlation between the Dynamics of Hydrogen Bonds and the Local Density Reorganization in the Protein Hydration Layer

An atomistic molecular dynamics simulation of the protein villin headpiece subdomain or HP-36 has been carried out with explicit water to explore the microscopic inhomogeneity of local density reorganization of the hydration layers of the three alpha-helical segments of the protein. The density reorganization of the hydration layer of helix-3 is found to occur faster than that for the hydration layers of the other two helices. It is noticed that such inhomogeneous density reorganization at the surface of different secondary structures exhibits excellent correlation with the microscopic dynamics of hydrogen bonds between the protein residues and the hydration water. Further, it is observed that the reorientation of water molecules involved in the formation and breaking of protein-water or water-water hydrogen bonds plays an important role in determining the dynamics of local density of the hydration layer. The faster density reorganization of the hydration layer of helix-3 is also consistent with the functionality of HP-36, as helix-3 contains several active site residues.

Single particle and collective hydration dynamics for hydrophobic and hydrophilic peptides

We have conducted extensive molecular dynamics simulations to study the single particle and collective dynamics of water in solutions of N-acetyl-glycine-methylamide, a model hydrophilic protein backbone, and N-acetyl-leucine-methylamide, a model (amphiphilic) hydrophobic peptide, as a function of peptide concentration. Various analytical models commonly used in the analysis of incoherent quasielastic neutron scattering (QENS), are tested against the translational and rotational intermediate scattering function, the mean square displacement of the water molecule center of mass, and fits to the second-order rotational correlation function of water evaluated directly from the simulation data. We find that while the agreement between the model-free analysis and analytical QENS models is quantitatively poor, the qualitative feature of dynamical heterogeneity due to caging is captured well by all approaches. The center of mass collective and single particle intermediate scattering functions of water calculated for these peptide solutions show that the crossover from collective to single particle-dominated motions occurs at a higher value of Q for high concentration solutions relative to low concentration because of the greater restriction in movement of water molecules due to confinement. Finally, we have shown that at the same level of confinement of the two peptides, the aqueous amphiphilic amino acid solution shows the strongest deviation between single particle and collective dynamics relative to the hydrophilic amino acid, indicating that chemical heterogeneity induces even greater spatial heterogeneity in the water dynamics.

Proton momentum distribution in a protein hydration shell

The momentum distribution of protons in the hydration shell of a globular protein has been measured through deep inelastic neutron scattering at 180 and 290 K, below and above the crossover temperature Tc=1.23Tg, where Tg=219 K is the glass transition temperature. It is found that the mean kinetic energy of the water hydrogens shows no temperature dependence, but the measurements are accurate enough to indicate a sensible change of momentum distribution and effective potential felt by protons, compatible with the transition from a single to a double potential well. This could support the presence of tunneling effects even at room temperature, playing an important role in biological function.

Grand canonical Monte Carlo simulations of water in protein environments

The grand canonical simulation algorithm is considered as a general methodology to sample the configuration of water molecules confined within protein environments. First, the probability distribution of the number of water molecules and their configuration in a region of interest for biochemical simulations, such as the active site of a protein, is derived by considering a finite subvolume in open equilibrium with a large system serving as a bulk reservoir. It is shown that the influence of the bulk reservoir can be represented as a many-body potential of mean force acting on the atoms located inside the subvolume. The grand canonical Monte Carlo (GCMC) algorithm, augmented by a number of technical advances to increase the acceptance of insertion attempts, is implemented, and tested for simple systems. In particular, the method is illustrated in the case of a pure water box with periodic boundary conditions. In addition, finite spherical systems of pure water and containing a dialanine peptide, are simulated with GCMC while the influence of the surrounding infinite bulk is incorporated using the generalized solvent boundary potential [W. Im, S. Berneche, and B. Roux, J. Chem. Phys. 114, 2924 (2001)]. As a last illustration of water confined in the interior of a protein, the hydration of the central cavity of the KcsA potassium channel is simulated.

Conformational heterogeneity and low-frequency vibrational modes of proteins

Molecular dynamics simulation and normal mode analysis are used to calculate the vibrational density of states of dihydrofolate reductase complexed with nicotinamide adenine dinucleotide phosphate at 120 K and the results are compared with the experimental spectrum derived from inelastic neutron scattering. The simulation results indicate that the experimental spectrum arises from an average over proteins trapped in different conformations with structural differences mainly in the loop regions, and that these conformations have significantly different low-frequency (<20 cm(-1)) spectra. Thus, the experimentally measured spectrum is an average over the vibrational modes of different protein conformations and is thus inhomogeneously broadened. The implications of this broadening for future neutron scattering experiments and ligand binding calculations are discussed.

Properties of Hydration Shells of Protein Molecules at their Pressure- and Temperature-Induced Native-Denatured Transition

Properties of water at the surface of biomolecules are important for their conformational stability. The behaviour of hydrating water at protein transition (t) pressures P(t) and temperatures T(t) , with the points (P(t),T(t) ) lying in the Native-Denatured (N-D) transition line, is studied. Hydration shells at the hydrophilic regions of protein molecules with surface charge density sigma are investigated with the help of the equation of state of water in an open system. The local values of sigma rather close to each other (sigma(D) approximately 0.3 C m(-2)) are found for six different experimental lines of the N-D transition found in the literature. The values sigma(D) correspond to the crossings of the total pressure (P(t)+Pi) vs sigma isotherms at different T(t) (Pi-electrostriction pressure). The pressures P(t) and temperatures T(t) appear to be related with some selected sites at the surfaces of the protein molecules.

PDB_Hydro: incorporating dipolar solvents with variable density in the Poisson-Boltzmann treatment of macromolecule electrostatics

We describe a new way to calculate the electrostatic properties of macromolecules which eliminates the assumption of a constant dielectric value in the solvent region, resulting in a Generalized Poisson–Boltzmann–Langevin equation (GPBLE). We have implemented a web server () that both numerically solves this equation and uses the resulting water density profiles to place water molecules at preferred sites of hydration. Surface atoms with high or low hydration preference can be easily displayed using a simple PyMol script, allowing for the tentative prediction of the dimerization interface in homodimeric proteins, or lipid binding regions in membrane proteins. The web site includes options that permit mutations in the sequence as well as reconstruction of missing side chain and/or main chain atoms. These tools are accessible independently from the electrostatics calculation, and can be used for other modeling purposes. We expect this web server to be useful to structural biologists, as the knowledge of solvent density should prove useful to get better fits at low resolution for X-ray diffraction data and to computational biologists, for whom these profiles could improve the calculation of interaction energies in water between ligands and receptors in docking simulations.

Protein stability and dynamics in the pressure–temperature plane

The pressure-temperature stability diagram of proteins and the underlying assumptions of the elliptical shape of the diagram are discussed. Possible extensions, such as aggregation and fibril formation, are considered. An important experimental observation is the extreme pressure stability of the mature fibrils. Molecular origins of the diagram in terms of models of the partial molar volume of a protein focus on cavities and hydration. Changes in thermal expansivity, compressibility and heat capacity in terms of fluctuations of the enthalpy and volume change of the unfolding should also focus on these parameters. It is argued that the study of water-soluble polymers might further our understanding of the stability diagram. Whereas the role of water in protein behaviour is unquestioned, the role of cavities is less clear.

A molecular dynamics simulation of SNase and its hydration shell at high temperature and high pressure

Temperature- and pressure-induced unfolding of staphylococcal nuclease (SNase) was studied by Royer, Winter et al. using a variety of experimental techniques (SAXS, FT-IR and fluorescence spectroscopy, DSC, PPC, densimetry). For a more detailed understanding of the underlying mechanistic processes of the different unfolding scenarios, we have carried out a series of molecular dynamics (MD) computer simulations on SNase. We investigated the initial changes of the structure of the protein upon application of pressure (up to 5 kbar) and discuss volumetric and structural differences between the native and pressure pre-denatured state. Additionally, we have obtained the compressibility of the protein and hydration water and compare these data with experimental results. As water plays a crucial role in determining the structure, dynamics and function of proteins, we undertook a detailed analysis of the structure of the interfacial water and the protein-solvent H-bond network as well. Moreover, we report here also MD results on the temperature-induced unfolding of SNase. The time evolution of the protein volume and solvent accessible surface area during thermal unfolding have been investigated, and we present a detailed discussion of the temperature-induced unfolding pathway of SNase in terms of secondary and tertiary structural changes.

Protein unfolding, amyloid fibril formation and configurational energy landscapes under high pressure conditions

High hydrostatic pressure induces conformational changes in proteins ranging from compression of the molecules to loss of native structure. In this tutorial review we describe how the interplay between the volume change and the compressibility leads to pressure-induced unfolding of proteins and dissociation of amyloid fibrils. We also discuss the effect of pressure on protein folding and free energy landscapes. From a molecular viewpoint, pressure effects can be rationalised in terms of packing and hydration of proteins.

New force field for calcium binding sites in annexin–membrane complexes

For accurate classical molecular dynamics (MD) simulations of the calcium mediated bound complexes of annexin and membrane we have developed new force-field parameters correctly describing the interaction of the Ca ion with its environment. We have used quantum chemical calculations to investigate the potential energy surface experienced by the Ca ion within the three different binding sites found in domain 1 of annexin V (ANX V/1). Based on these calculations we were able to quantify the charge polarization of atoms within the binding sites, and to determine the geometry and force constants of harmonic restraints between the Ca ion and its coordinating oxygen atoms. Harmonic restraints were introduced to compensate for the deviations between the quantum mechanical potential energy surface and that of the classical force field. Our analysis has shown that using the refined force field for the Ca binding sites enables long-time MD simulations that conserve the initial structure of ANX V/1 significantly better than MD simulations using the standard force field.

Pressure perturbation calorimetric studies of the solvation properties and the thermal unfolding of proteins in solution—experiments and theoretical interpretation

We used pressure perturbation calorimetry (PPC), a relatively new and efficient technique, to study the solvation and volumetric properties of amino acids and peptides as well as of proteins in their native and unfolded state. In PPC, the coefficient of thermal expansion of the partial volume of the protein is deduced from the heat consumed or produced after small isothermal pressure jumps, which strongly depends on the interaction of the protein with the solvent or cosolvent at the protein-solvent interface. Furthermore, the effects of various chaotropic and kosmotropic cosolvents on the volume and expansivity changes of proteins were measured over a wide concentration range with high precision. Depending on the type of cosolvent and its concentration, specific differences were found for the solvation properties and unfolding behaviour of the proteins, and the volume change upon unfolding may even change sign. To yield a molecular interpretation of the different terms contributing to the partial protein volume and its temperature dependence, and hence a better understanding of the PPC data, molecular dynamics computer simulations on SNase were also carried out and compared with the experimental data. The PPC studies introduced aim to obtain more insight into the basic thermodynamic properties of protein solvation and volume effects accompanying structural transformations of proteins in various cosolvents on one hand, as these form the basis for understanding their physiological functions and their use in drug designing and formulations, but also to initiate further valuable applications in studies of other biomolecular and chemical systems.

Pressure-dependent transition in protein dynamics at about revealed by molecular dynamics simulation

Molecular dynamics simulations of a crystalline protein, Staphylococcal nuclease, over the pressure range 1 bar to 15 kbar reveal a qualitative change in the internal protein motions at approximately 4 kbar. This change involves the existence of two linear regimes in the mean-square displacement for internal protein motion, <mu2>(P) with a twofold decrease in the slope for P>4 kbar. The major effect of pressure on the dynamics is a loss, with increasing pressure of large amplitude, collective protein modes below 2 THz effective frequency, accompanied by restriction of large-scale solvent translational motion.

A salt-bridge motif involved in ligand binding and large-scale domain motions of the maltose-binding protein

The uptake of nutrients is essential for the survival of bacterial cells. Many specialized systems have evolved, such as the maltose-dependent ABC transport system that transfers oligosaccharides through the cytoplasmic membrane. The maltose/maltodextrin-binding protein (MBP) serves as an initial high-affinity binding component in the periplasm that delivers the bound sugar into the cognate ABC transporter MalFGK(2). We have investigated the domain motions induced by the binding of the ligand maltotriose into the binding cleft using molecular dynamics simulations. We find that MBP is predominantly in the open state without ligand and in the closed state with ligand bound. Oligosaccharide binding induces a closure motion (30.0 degrees rotation), whereas ligand removal leads to domain opening (32.6 degrees rotation) around a well-defined hinge affecting key areas relevant for chemotaxis and transport. Our simulations suggest that a "hook-and-eye" motif is involved in the binding. A salt bridge between Glu-111 and Lys-15 forms that effectively locks the protein-ligand complex in a semiclosed conformation inhibiting any further opening and promoting complete closure. This previously unrecognized feature seems to secure the ligand in the binding site and keeps MBP in the closed conformation and suggests a role in the initial steps of substrate transport.

How thick is the layer of thermal volume surrounding the protein?

Investigation on the volume properties of protein hydration layers is reported. Presented results are based on combination of Monte Carlo modeling and available experimental data. Six globular proteins with known data are chosen for analysis. Analyzing the model and the experimental results we found that water molecules bound to proteins by hydrogen bond are preferentially located at the places with local depressions on the protein surface. Consequently, the hydration level is not strictly proportional to the area of charged and polar surfaces, but also depends on the shape of the molecular surface. The thickness of the thermal volume layer as calculated in the framework of the scaled particle theory is 0.6-0.65 A for chosen proteins. The obtained value is significantly lower than that presented for proteins in earlier papers (where proportionality between the hydration level and the area of charged and polar surfaces was assumed), but is close to the value published for small solute molecules. Discussion including the influence of protein size and the thermal motion of the surface is presented.

High-Density Hydration Layer of Lysozymes: Molecular Dynamics Decomposition of Solution Scattering Data

A characterization of the physical properties of protein hydration water is critical for understanding protein structure and function. Recent small-angle X-ray and neutron scattering data indicate that the density of water on the surface of lysozyme is significantly higher than in bulk water. Here, we provide an interpretation of the scattering results using a molecular dynamics simulation, which allows us to make quantitative predictions about density variations in the first hydration shell. The perturbation relative to bulk water involves statistically significant changes in the average water structure in the first hydration layer. The water density in the first hydration shell is increased by 5% with respect to the bulk. In regions of higher water density, the water dipoles align more parallel to each other and the number of hydrogen bonds per water molecule is higher. Increased water density is found for water molecules interacting with hydrogen and carbon atoms in the backbone or with nonpolar or negatively charged side-chain groups.

Molecular view of water dynamics near model peptides

Incoherent quasi-elastic neutron scattering (QENS) has been used to measure the dynamics of water molecules in solutions of a model protein backbone, N-acetyl-glycine-methylamide (NAGMA), as a function of concentration, for comparison with results for water dynamics in aqueous solutions of the N-acetyl-leucine-methylamide (NALMA) hydrophobic peptide at comparable concentrations. From the analysis of the elastic incoherent structure factor, we find significant fractions of elastic intensity at high and low concentrations for both solutes, which corresponds to a greater population of protons with rotational time scales outside the experimental resolution (>13 ps). The higher-concentration solutions show a component of the elastic fraction that we propose is due to water motions that are strongly coupled to the solute motions, while for low-concentration solutions an additional component is activated due to dynamic coupling between inner and outer hydration layers. An important difference between the solute types at the highest concentration studied is found from stretched exponential fits to their experimental intermediate scattering functions, showing more pronounced anomalous diffusion signatures for NALMA, including a smaller stretched exponent beta and a longer structural relaxation time tau than those found for NAGMA. The more normal water diffusion exhibited near the hydrophilic NAGMA provides experimental support for an explanation of the origin of the anomalous diffusion behavior of NALMA as arising from frustrated interactions between water molecules when a chemical interface is formed upon addition of a hydrophobic side chain, inducing spatial heterogeneity in the hydration dynamics in the two types of regions of the NALMA peptide. We place our QENS measurements on model biological solutes in the context of other spectroscopic techniques and provide both confirming as well as complementary dynamic information that attempts to give a unifying molecular view of hydration dynamics signatures near peptides and proteins.

Study on hydration layers near nanoscale silica dispersed in aqueous solutions through viscosity measurement

On the basis of the Einstein theory of viscosity of dispersion, a parameter, termed as solvation factor, is presented to evaluate the solvation degree of nanoscale particles dispersed in a liquid in this work. The value of the parameter is obtained through the measurements of relative viscosity of the dispersions as a function of the volume fraction of dry particles. The solvation factor has been used to study the hydration layers near nanoscale silica particles dispersed in water and aqueous electrolyte (NaCl and CaCl2) solutions in this work. The experimental results have shown that a strong hydration indeed applied to the silica surfaces in aqueous solutions, leaving a large volume of hydration layers on the surfaces. Also, it has been found that the hydration of the nanoscale silica particles could be greatly enhanced if they were dispersed in aqueous NaCl or CaCl2 solutions, which might be attributed to that the hydrated cations (Na+ or Ca2+) bind onto the silica/ water interface and thus increase the volume of the hydration layers.

Light-scattering studies of protein solutions: role of hydration in weak protein-protein interactions

We model the hydration contribution to short-range electrostatic/dispersion protein interactions embodied in the osmotic second virial coefficient, B(2), by adopting a quasi-chemical description in which water molecules associated with the protein are identified through explicit molecular dynamics simulations. These water molecules reduce the surface complementarity of highly favorable short-range interactions, and therefore can play an important role in mediating protein-protein interactions. Here we examine this quasi-chemical view of hydration by predicting the interaction part of B(2) and comparing our results with those derived from light-scattering measurements of B(2) for staphylococcal nuclease, lysozyme, and chymotrypsinogen at 25 degrees C as a function of solution pH and ionic strength. We find that short-range protein interactions are influenced by water molecules strongly associated with a relatively small fraction of the protein surface. However, the effect of these strongly associated water molecules on the surface complementarity of short-range protein interactions is significant, and must be taken into account for an accurate description of B(2). We also observe remarkably similar hydration behavior for these proteins despite substantial differences in their three-dimensional structures and spatial charge distributions, suggesting a general characterization of protein hydration.

Coupling between lysozyme and glycerol dynamics: Microscopic insights from molecular-dynamics simulations

We explore possible molecular mechanisms behind the coupling of protein and solvent dynamics using atomistic molecular-dynamics simulations. For this purpose, we analyze the model protein lysozyme in glycerol, a well-known protein-preserving agent. We find that the dynamics of the hydrogen bond network between the solvent molecules in the first shell and the surface residues of the protein controls the structural relaxation (dynamics) of the whole protein. Specifically, we find a power-law relationship between the relaxation time of the aforementioned hydrogen bond network and the structural relaxation time of the protein obtained from the incoherent intermediate scattering function. We demonstrate that the relationship between the dynamics of the hydrogen bonds and the dynamics of the protein appears also in the dynamic transition temperature of the protein. A study of the dynamics of glycerol as a function of the distance from the surface of the protein indicates that the viscosity seen by the protein is not the one of the bulk solvent. The presence of the protein suppresses the dynamics of the surrounding solvent. This implies that the protein sees an effective viscosity higher than the one of the bulk solvent. We also found significant differences in the dynamics of surface and core residues of the protein. The former is found to follow the dynamics of the solvent more closely than the latter. These results allowed us to propose a molecular mechanism for the coupling of the solvent-protein dynamics.

Stability of proteins: Temperature, pressure and the role of the solvent

We focus on the various aspects of the physics related to the stability of proteins. We review the pure thermodynamic aspects of the response of a protein to pressure and temperature variations and discuss the respective stability phase diagram. We relate the experimentally observed shape of this diagram to the low degree of correlation between the fluctuations of enthalpy and volume changes associated with the folding-denaturing transition and draw attention to the fact that one order parameter is not enough to characterize the transition. We discuss in detail microscopic aspects of the various contributions to the free energy gap of proteins and put emphasis on how a cosolvent may either enlarge or diminish this gap. We review briefly the various experimental approaches to measure changes in protein stability induced by cosolvents, denaturants, but also by pressure and temperature. Finally, we discuss in detail our own molecular dynamics simulations on cytochrome c and show what happens under high pressure, how glycerol influences structure and volume fluctuations, and how all this compares with experiments.

Retention of nativelike conformation by proteins embedded in high external electric fields

In this Communication, we show that proteins embedded in high external electric fields are capable of retaining a nativelike fold pattern. We have tested the metalloprotein azurin, immobilized onto SiO2 substrates in air with proper electrode configuration, by applying static fields up to 10(6)-10(7) Vm. The effects on the conformational properties of protein molecules have been determined by means of intrinsic fluorescence measurements. Experimental results indicate that no significant field-induced conformational alteration occurs. Such results are also discussed and supported by theoretical predictions of the inner protein fields.