Hydration dependence of backbone and side chain polylysine dynamics: a 13C solid-state NMR and IR spectroscopy study

The molecular dynamics of solid poly-L-lysine has been studied by the following natural abundance (13)C-NMR relaxation methods: measurements of the relaxation times T(1) at two resonance frequencies, off-resonance T(1rho) at two spin-lock frequencies, and proton-decoupled T(1rho). Experiments were performed at different temperatures and hydration levels (up to 17% H(2)O by weight). The natural abundance (13)C-CPMAS spectrum of polylysine provides spectral resolution of all types of backbone and side chain carbons and thus, dynamic parameters could be determined separately for each of them. At the same time, the conformational properties of polylysine were investigated by Fourier transform infrared spectroscopy. The data obtained from the different NMR experiments were simultaneously analyzed using the correlation function formalism and model-free approach. The results indicate that in dry polylysine both backbone and side chains take part in two low amplitude motions with correlation times of the order of 10(-4) s and 10(-9) s. Upon hydration, the dynamic parameters of the backbone remain almost constant except for the amplitude of the slower process that increases moderately. The side chain dynamics reveals a much stronger hydration response: the amplitudes of both slow and fast motions increase significantly and the correlation time of the slow motion shortens by about five orders of magnitude, and at hydration levels of more than 10% H(2)O fast and slow side chain motions are experimentally indistinguishable. These changes in the molecular dynamics cannot be ascribed to any hydration-dependent conformational transitions of polylysine because IR spectra reveal almost no hydration dependence in either backbone or side chain absorption domains. The physical nature of the fast and slow motions, their correlation time distributions, and hydration dependence of microdynamic parameters are discussed.     

THz time scale structural rearrangements and binding modes in lysozyme-ligand interactions

Predicting the conformational changes in proteins that are relevant for substrate binding is an ongoing challenge in the aim of elucidating the functional states of proteins. The motions that are induced by protein-ligand interactions are governed by the protein global modes. Our measurements indicate that the detected changes in the global backbone motion of the enzyme upon binding reflect a shift from the large-scale collective dominant mode in the unbound state towards a functional twisting deformation that assists in closing the binding cleft. Correlated motion in lysozyme has been implicated in enzyme function in previous studies, but detailed characterization of the internal fluctuations that enable the protein to explore the ensemble of conformations that ultimately foster large-scale conformational change is yet unknown. For this reason, we use THz spectroscopy to investigate the picosecond time scale binding modes and collective structural rearrangements that take place in hen egg white lysozyme (HEWL) when bound by the inhibitor (NAG) ₃. These protein thermal motions correspond to fluctuations that have a role in both selecting and sampling from the available protein intrinsic conformations that communicate function. Hence, investigation of these fast, collective modes may provide knowledge about the mechanism leading to the preferred binding process in HEWL-(NAG) ₃. Specifically, in this work we find that the picosecond time scale hydrogen-bonding rearrangements taking place in the protein hydration shell with binding modify the packing density within the hydrophobic core on a local level. These localized, intramolecular contact variations within the protein core appear to facilitate the large cooperative movements within the interfacial region separating the α- and β- domain that mediate binding. The THz time-scale fluctuations identified in the protein-ligand system may also reveal a molecular mechanism for substrate recognition.     

Probing the collective vibrational dynamics of a protein in liquid water by terahertz absorption spectroscopy

Biological polymers are expected to exhibit functionally relevant, global, and subglobal collective modes in the terahertz (THz) frequency range (i.e., picosecond timescale). In an effort to monitor these collective motions, we have experimentally determined the absorption spectrum of solvated bovine serum albumin (BSA) from 0.3 to 3.72 THz (10-124 cm(-1)). We successfully extract the terahertz molar absorption of the solvated BSA from the much stronger attenuation of water and observe in the solvated protein a dense, overlapping spectrum of vibrational modes that increases monotonically with increasing frequency. We see no evidence of distinct, strong, spectral features, suggesting that no specific collective vibrations dominate the protein's spectrum of motions, consistent with the predictions of molecular dynamics simulations and normal mode analyses of a range of small proteins. The shape of the observed spectrum resembles the ideal quadratic spectral density expected for a disordered ionic solid, indicating that the terahertz normal mode density of the solvated BSA may be modeled, to first order, as that of a three-dimensional elastic nanoparticle with an aperiodic charge distribution. Nevertheless, there are important detailed departures from that of a disordered inorganic solid or the normal mode densities predicted for several smaller proteins. These departures are presumably the spectral features arising from the unique molecular details of the solvated BSA. The techniques used here and measurements have the potential to experimentally confront theoretical calculations on a frequency scale that is important for macromolecular motions in a biologically relevant water environment.     

Long-range DNA-water interactions

DNA functions only in aqueous environments and adopts different conformations depending on the hydration level. The dynamics of hydration water and hydrated DNA leads to rotating and oscillating dipoles that, in turn, give rise to a strong megahertz to terahertz absorption. Investigating the impact of hydration on DNA dynamics and the spectral features of water molecules influenced by DNA, however, is extremely challenging because of the strong absorption of water in the megahertz to terahertz frequency range. In response, we have employed a high-precision megahertz to terahertz dielectric spectrometer, assisted by molecular dynamics simulations, to investigate the dynamics of water molecules within the hydration shells of DNA as well as the collective vibrational motions of hydrated DNA, which are vital to DNA conformation and functionality. Our results reveal that the dynamics of water molecules in a DNA solution is heterogeneous, exhibiting a hierarchy of four distinct relaxation times ranging from ∼8 ps to 1 ns, and the hydration structure of a DNA chain can extend to as far as ∼18 Å from its surface. The low-frequency collective vibrational modes of hydrated DNA have been identified and found to be sensitive to environmental conditions including temperature and hydration level. The results reveal critical information on hydrated DNA dynamics and DNA-water interfaces, which impact the biochemical functions and reactivity of DNA.     

Excited state charge-transfer dynamics study of poplar plastocyanin by ultrafast pump-probe spectroscopy and molecular dynamics simulation

We have applied ultrafast pump-probe spectroscopy to investigate the excited state dynamics of the blue copper protein poplar plastocyanin, by exciting in the blue side of its 600-nm absorption band. The decay of the charge-transfer excited state occurs exponentially with a time constant of approximately 280 fs and is modulated by well visible oscillations. The Fourier transform of the oscillatory component, besides providing most of the vibrational modes found by conventional resonance Raman, presents additional bands in the low frequency region modes, which are reminiscent of collective motions of biological relevance. Notably, a high frequency mode at approximately 508 cm(-1), whose dynamics are consistent with that of the excited state and already observed for other blue copper proteins, is shown to be present also in poplar plastocyanin. This vibrational mode is reproduced by a molecular dynamics simulation involving the excited state of the copper site.     

Hydration dependent dynamics in RNA

The essential role played by local and collective motions in RNA function has led to a growing interest in the characterization of RNA dynamics. Recent investigations have revealed that even relatively simple RNAs experience complex motions over multiple time scales covering the entire ms–ps motional range. In this work, we use deuterium solid-state NMR to systematically investigate motions in HIV-1 TAR RNA as a function of hydration. We probe dynamics at three uridine residues in different structural environments ranging from helical to completely unrestrained. We observe distinct and substantial changes in ²H solid-state relaxation times and lineshapes at each site as hydration levels increase. By comparing solid-state and solution state ¹³C relaxation measurements, we establish that ns–μs motions that may be indicative of collective dynamics suddenly arise in the RNA as hydration reaches a critical point coincident with the onset of bulk hydration. Beyond that point, we observe smaller changes in relaxation rates and lineshapes in these highly hydrated solid samples, compared to the dramatic activation of motion occurring at moderate hydration.     

Analysis of Sub-τc and Supra-τc Motions in Protein Gβ1 Using Molecular Dynamics Simulations

The functions of proteins depend on the dynamical behavior of their native states on a wide range of timescales. To investigate these dynamics in the case of the small protein Gβ1, we analyzed molecular dynamics simulations with the model-free approach of nuclear magnetic relaxation. We found amplitudes of fast timescale motions (sub-τ_c, where τ_c is the rotational correlation time) consistent with S² obtained from spin relaxation measurements as well as amplitudes of slow timescale motions (supra-τ_c) in quantitative agreement with S² order parameters derived from residual dipolar coupling measurements. The slow timescale motions are associated with the large variations of the ³J couplings that follow transitions between different conformational substates. These results provide further characterization of the large structural fluctuations in the native states of proteins that occur on timescales longer than the rotational correlation time.     

Hydration dependence of conformational dielectric relaxation of lysozyme

Dielectric response of hen egg white lysozyme is measured in the far infrared (5-65 cm-1, 0.15-1.95 THz, 0.6-8.1 meV) as a function of hydration. The frequency range is associated with collective vibrational modes of protein tertiary structure. The observed frequency dependence of the absorbance is broad and glass-like. For the entire frequency range, there is a slight increase in both the absorbance and index of refraction with increasing hydration for <0.27 h (mass of H2O per unit mass protein). At 0.27 h, the absorbance and index begin to increase more rapidly. This transition corresponds to the point where the first hydration shell is filled. The abrupt increase in dielectric response cannot be fully accounted for by the additional contribution to the dielectric response due to bulk water, suggesting that the protein has not yet achieved its fully hydrated state. The broad, glass-like response suggests that at low hydrations, the low frequency conformational hen egg white lysozyme dynamics can be described by a dielectric relaxation model where the protein relaxes to different local minima in the conformational energy landscape. However, the low frequency complex permittivity does not allow for a pure relaxational mechanism. The data can best be modeled with a single low frequency resonance (nu approximately 120 GHz=4 cm-1) and a single Debye relaxation process (tau approximately .03-.04 ps). Terahertz dielectric response is currently being considered as a possible biosensing technique and the results demonstrate the required hydration control necessary for reliable biosensor applications.     

Equilibrium fluctuations in myoglobin and lysozyme

The angular dependencies of inelastic intensities of Rayleigh scattering of Moessbauer radiation were measured for myoglobin and lysozyme (in the hydration range h = 0.05-0.7). The data were fitted within the framework of model, when two types of intraglobular motions were taken into account: individual motions of small side-chain groups and cooperative motions of segments. The best agreement with the experiment at h > 0.05 was obtained when individual motions of small groups together with the cooperative motions of alpha-helices and beta-sheets for lysozyme, and alpha-helices for myoglobin were considered. At further hydration (h = 0.45), mean-square displacements (x2) of both types of motions strongly increase with the increase in hydration degree, while the motions with a large correlation radius (not less than macromolecule radius) remain nearly the same as for h = 0.05. The results of the study of the radial distribution function deduced by Fourier-transform from the diffuse x-ray measurements together with RSMR data allow one to conclude that the water during protein hydration competes with the intramolecular hydrogen bonds, loosens the protein and increases the internal dynamics. Concurrently, water arranges the ordering of macromolecule, which takes the native structure at h = 0.4-0.7. The analysis of auto and cross-correlation functions of bending fluctuations of alpha-helices in the large domain of lysozyme performed by molecular dynamics allows one to come to the final conclusion that it is the difference in the structural organization of myoglobin and lysozyme and not the presence of SS-bonds in lysozyme macromolecule that is responsible for different structural fluctuations in these proteins.     

Flexibility of ras lipid modifications studied by 2H solid-state NMR and molecular dynamics simulations

Human posttranslationally modified N-ras oncogenes are known to be implicated in numerous human cancers. Here, we applied a combination of experimental and computational techniques to determine structural and dynamical details of the lipid chain modifications of an N-ras heptapeptide in 1,2-dimyristoyl-sn-glycero-3-phosphocholine (DMPC) membranes. Experimentally, 2H NMR spectroscopy was used to study oriented membranes that incorporated ras heptapeptides with two covalently attached perdeuterated hexadecyl chains. Atomistic molecular dynamics simulations of the same system were carried out over 100 ns including 60 DMPC and 4 ras molecules. Several structural and dynamical experimental parameters could be directly compared to the simulation. Experimental and simulated 2H NMR order parameters for the methylene groups of the ras lipid chains exhibited a systematic difference attributable to the absence of collective motions in the simulation and to geometrical effects. In contrast, experimental 2H NMR spin-lattice relaxation rates for Zeeman order were well reproduced in the simulation. The lack of slower collective motions in the simulation did not appreciably influence the relaxation rates at a Larmor frequency of 115.1 MHz. The experimental angular dependence of the 2H NMR relaxation rates with respect to the external magnetic field was also relatively well simulated. These relaxation rates showed a weak angular dependence, suggesting that the lipid modifications of ras are very flexible and highly mobile in agreement with the low order parameters. To quantify these results, the angular dependence of the 2H relaxation rates was calculated by an analytical model considering both molecular and collective motions. Peptide dynamics in the membrane could be modeled by an anisotropic diffusion tensor with principal values of Dparallel=2.1x10(9) s(-1) and Dperpendicular=4.5x10(5) s(-1). A viscoelastic fitting parameter describing the membrane elasticity, viscosity, and temperature was found to be relatively similar for the ras peptide and the DMPC host matrix. Large motional amplitudes and relatively short correlation times facilitate mixing and dispersal with the lipid bilayer matrix, with implications for the role of the full-length ras protein in signal transduction and oncogenesis.     

Protein Collective Motions Coupled to Ligand Migration in Myoglobin

Ligand migration processes inside myoglobin and protein dynamics coupled to the migration were theoretically investigated with molecular dynamics simulations. Based on a linear response theory, we identified protein motions coupled to the transient migration of ligand, carbon monoxide (CO), through channels. The result indicates that the coupled protein motions involve collective motions extended over the entire protein correlated with local gating motions at the channels. Protein motions, coupled to opening of a channel from the distal pocket to a neighboring xenon site, were found to share the collective motion with experimentally observed protein motions coupled to a doming motion of the heme Fe atom upon photodissociation of the ligand. Analysis based on generalized Langevin dynamics elucidated slow and diffusive features of the protein response motions. Remarkably small transmission coefficients for rates of the CO migrations through myoglobin were found, suggesting that the CO migration dynamics are characterized as motions governed by the protein dynamics involving the collective motions, rather than as thermally activated transitions across energy barriers of well-structured channels.     

Resilience of the iron environment in heme proteins

Conformational flexibility is essential to the functional behavior of proteins. We use an effective force constant introduced by Zaccai, the resilience, to quantify this flexibility. Site-selective experimental and computational methods allow us to determine the resilience of heme protein active sites. The vibrational density of states of the heme Fe determined using nuclear resonance vibrational spectroscopy provides a direct experimental measure of the resilience of the Fe environment, which we compare quantitatively with values derived from the temperature dependence of atomic mean-squared displacements in molecular dynamics simulations. Vibrational normal modes in the THz frequency range dominate the resilience. Both experimental and computational methods find a higher resilience for cytochrome c than for myoglobin, which we attribute to the increased number of covalent links to the peptide in the former protein. For myoglobin, the resilience of the iron environment is larger than the average resilience previously determined for hydrogen sites using neutron scattering. Experimental results suggest a slightly reduced resilience for cytochrome c upon oxidation, although the change is smaller than reported in previous Mössbauer investigations on a bacterial cytochrome c, and is not reproduced by the simulations. Oxidation state also has no significant influence on the compressibility calculated for cyt c, although a slightly larger compressibility is predicted for myoglobin.     

Differences in internal dynamics of actin under different structural states detected by neutron scattering

F-actin, a helical polymer formed by polymerization of the monomers (G-actin), plays crucial roles in various aspects of cell motility. Flexibility of F-actin has been suggested to be important for such a variety of functions. Understanding the flexibility of F-actin requires characterization of a hierarchy of dynamical properties, from internal dynamics of the actin monomers through domain motions within the monomers and relative motions between the monomers within F-actin to large-scale motions of F-actin as a whole. As a first step toward this ultimate purpose, we carried out elastic incoherent neutron scattering experiments on powders of F-actin and G-actin hydrated with D₂O and characterized the internal dynamics of F-actin and G-actin. Well established techniques and analysis enabled the extraction of mean-square displacements and their temperature dependence in F-actin and in G-actin. An effective force constant analysis with a model consisting of three energy states showed that two dynamical transitions occur at ∼150 K and ∼245 K, the former of which corresponds to the onset of anharmonic motions and the latter of which couples with the transition of hydration water. It is shown that behavior of the mean-square displacements is different between G-actin and F-actin, such that G-actin is “softer” than F-actin. The differences in the internal dynamics are detected for the first time between the different structural states (the monomeric state and the polymerized state). The different behavior observed is ascribed to the differences in dynamical heterogeneity between F-actin and G-actin. Based on structural data, the assignment of the differences observed in the two samples to dynamics of specific loop regions involved in the polymerization of G-actin into F-actin is proposed.     

Do hydration dynamics follow the structural perturbation during thermal denaturation of a protein: a terahertz absorption study

We investigate the thermal denaturation of human serum albumin and the associated solvation using terahertz (THz) spectroscopy in aqueous buffer solution. Far- and near-ultraviolet circular dichroism spectroscopy reveal that the protein undergoes a native (N) to extended (E) state transition at temperature ≤55°C with a marginal change in the secondary and tertiary structure. At 70°C, the protein transforms into an unfolded (U) state with significant irreversible disruption of its structures. We measure the concentration- and temperature-dependent THz absorption coefficient (α) of the protein solution using a p-Ge THz difference spectrometer (2.1–2.8 THz frequency range), thereby probing the collective protein-water network dynamics. When the solvated protein is heated up to 55°C and cooled down again, a reversible change in THz absorption is observed. When increasing the temperature up to 70°C, we find a dramatic irreversible change of THz absorption. The increase in THz absorption compared to bulk water is attributed to a blue shift in the spectrum of the solvated protein compared to bulk water. This is supported by measurements of THz absorption coefficients using THz time-domain spectroscopy (0.1–1.2 THz frequency range). We also use picosecond-resolved fluorescence spectroscopy of the tryptophan 214 moiety of human serum albumin. All experimental observations can be explained by a change in the hydration dynamics of the solvated protein due to the additional exposure of hydrophobic residues upon unfolding.     

Classical molecular dynamics simulation of the photoinduced electron transfer dynamics of plastocyanin.

Classical molecular dynamics simulations are used to investigate the nuclear motions associated with photoinduced electron transfer in plastocyanin. The blue copper protein is modeled using a molecular mechanics potential; potential parameters for the copper-protein interactions are determined using an x-ray crystallographic structure and absorption and resonance Raman spectra. Molecular dynamics simulations yield a variety of information about the ground (oxidized) and optically excited (charge-transfer) states: 1) The probability distribution of the potential difference between the states, which is used to determine the coordinate and energy displacements, places the states well within the Marcus inverted region. 2) The two-time autocorrelation function of the difference potential in the ground state and the average of the difference potential after instantaneous excitation to the excited state are very similar (confirming linear response in this system); their decay indicates that vibrational relaxation occurs in about 1 ps in both states. 3) The spectral densities of various internal coordinates begin to identify the vibrations that affect the optical transition; the spectral density of the difference potential correlation function should also prove useful in quantum simulations of the back electron transfer. 4) Correlation functions of the protein atomic motions with the difference potential show that the nuclear motions are correlated over a distance of more than 20 A, especially along proposed electron transport paths.     

Perturbations at the chloride site during the photosynthetic oxygen-evolving cycle

Photosystem II (PSII) catalyzes the oxidation of water to O₂ at the manganese-containing, oxygen-evolving complex (OEC). Photoexcitation of PSII results in the oxidation of the OEC; four sequential oxidation reactions are required for the generation and release of molecular oxygen. Therefore, with flash illumination, the OEC cycles among five S _n states. Chloride depletion inhibits O₂ evolution. However, the binding site of chloride in the OEC is not known, and the role of chloride in oxygen evolution has not as yet been elucidated. We have employed reaction-induced FT-IR spectroscopy and selective flash excitation, which cycles PSII samples through the S state transitions. On the time scale employed, these FT-IR difference spectra reflect long-lived structural changes in the OEC. Bromide substitution supports oxygen evolution and was used to identify vibrational bands arising from structural changes at the chloride-binding site. Contributions to the vibrational spectrum from bromide-sensitive bands were observed on each flash. Sulfate treatment led to an elimination of oxygen evolution activity and of the FT-IR spectra assigned to the S₃ to S₀ (third flash) and S₀ to S₁ transitions (fourth flash). However, sulfate treatment changed, but did not eliminate, the FT-IR spectra obtained with the first and second flashes. Solvent isotope exchange in chloride-exchanged samples suggests flash-dependent structural changes, which alter protein dynamics during the S state cycle. Supported by NSF MCB 03-55421.     

Intrinsic Dynamics of Restriction Endonuclease EcoO109I Studied by Molecular Dynamics Simulations and X-Ray Scattering Data Analysis

EcoO109I is a type II restriction endonuclease that functions as a dimer in solution. Upon DNA binding to the enzyme, the two subunits rotate counterclockwise relative to each other, as the two catalytic domains undergo structural changes to capture the cognate DNA. Using a 150-ns molecular dynamics simulation, we investigated the intrinsic dynamics of the DNA-free enzyme in solution to elucidate the relationship between enzyme dynamics and structural changes. The simulation revealed that the enzyme is considerably flexible, and thus exhibits large fluctuations in the radius of gyration. The small-angle x-ray scattering profile calculated from the simulation, including scattering from explicit hydration water, was in agreement with the experimentally observed profile. Principal component analysis revealed that the major dynamics were represented by the open-close and counterclockwise motions: the former is required for the enzyme to access DNA, whereas the latter corresponds to structural changes upon DNA binding. Furthermore, the intrinsic dynamics in the catalytic domains were consistent with motions capturing the cognate DNA. These results indicate that the structure of EcoO109I is intrinsically flexible in the direction of its functional movement, to facilitate effective structural changes for sequence-specific DNA recognition and processing.     

A view of dynamics changes in the molten globule-native folding step by quasielastic neutron scattering

In order to understand the changes in protein dynamics that occur in the final stages of protein folding, we have used neutron scattering to probe the differences between a protein in its folded state and the molten globule states. The internal dynamics of bovine alpha-lactalbumin (BLA) and its molten globules (MBLA) have been examined using incoherent, quasielastic neutron scattering (IQNS). The IQNS results show length scale dependent, pico-second dynamics changes on length scales from 3.3 to 60 A studied. On shorter-length scales, the non-exchangeable protons undergo jump motions over potential barriers, as those involved in side-chain rotamer changes. The mean potential barrier to local jump motions is higher in BLA than in MBLA, as might be expected. On longer length scales, the protons undergo spatially restricted diffusive motions with the diffusive motions being more restricted in BLA than in MBLA. Both BLA and MBLA have similar mean square amplitudes of high frequency motions comparable to the chemical bond vibrational motions. Bond vibrational motions thus do not change significantly upon folding. Interestingly, the quasielastic scattering intensities show pronounced maxima for both BLA and MBLA, suggesting that "clusters" of atoms are moving collectively within the proteins on picosecond time scales. The correlation length, or "the cluster size", of such atom clusters moving collectively is dramatically reduced in the molten globules with the correlation length being 6.9 A in MBLA shorter than that of 18 A in BLA. Such collective motions may be important for the stability of the folded state, and may influence the protein folding pathways from the molten globules.     

Incoherent neutron scattering of copper azurin: a comparison with molecular dynamics simulation results

The low-frequency dynamics of copper azurin has been studied at different temperatures for a dry and deuterium hydrated sample by incoherent neutron scattering and the experimental results have been compared with molecular dynamics (MD) simulations carried out in the same temperature range. Experimental Debye-Waller factors are consistent with a dynamical transition at approximately 200 K which appears partially suppressed in the dry sample. Inelastic and quasielastic scattering indicate that hydration water modulates both vibrational and diffusive motions. The low-temperature experimental dynamical structure factor of the hydrated protein shows an excess of inelastic scattering peaking at about 3 meV and whose position is slightly shifted downwards in the dry sample. Such an excess is reminiscent of the “boson peak” observed in glass-like materials. This vibrational peak is quite well reproduced by MD simulations, although at a lower energy. The experimental quasielastic scattering of the two samples at 300 K shows a two-step relaxation behaviour with similar characteristic times, while the corresponding intensities differ only by a scale factor. Also, MD simulations confirm the two-step diffusive trend, but the slow process seems to be characterized by a decay faster than the experimental one. Comparison with incoherent neutron scattering studies carried out on proteins having different structure indicates that globular proteins display common elastic, quasielastic and inelastic features, with an almost similar hydration dependence, irrespective of their secondary and tertiary structure. Received: 12 October 1998 / Revised version: 19 February 1999 / Accepted: 1 March 1999     

Water's potential role: Insights from studies of the p53 core domain

Soluble proteins with amyloidogenic propensity such as the tumor suppressor protein p53 have high proportion of incompletely desolvated backbone H bonds (HB). Such bonds are vulnerable to water attack, thus potentially leading to the misfolding of these proteins. However, it is still not clear how the surrounding solvent influences the protein native states. To address this, systematic surveys by molecular dynamics simulations and entropy analysis were performed on the p53 core domain in this work. We examined seven wild/mutant X-ray structures and observed two types of water-network hydration in three "hot hydration centers" (DNA- or small molecule- binding surfaces of the p53 core domain). The "tight" water, resulting from the local collective hydrogen-bond interactions, is probably fundamental to the protein structural stability. The second type of water is highly "dynamical" and exchanges very fast within the bulk solution, which is unambiguously assisted by the local protein motions. An entropy mapping of the solvent around the protein and a temperature perturbation analysis further present the main features of the p53 hydration network. The particular environment created by different water molecules around the p53 core domain also partly explains the structural vulnerabilities of this protein.