Protein Hydration and Water Structure: X-ray Analysis of a Closely Packed Protein Crystal with Very Low Solvent Content

Low-humidity monoclinic lysozyme, resulting from a water-mediated transformation, has one of the lowest solvent contents (22% by volume) observed in a protein crystal. Its structure has been solved by the molecular replacement method and refined to an R value of 0.175 for 7684 observed reflections in the 10–1.75 Å resolution shell. 90% of the solvent in the well ordered crystals could be located. Favourable sites of hydration on the protein surface include side chains with multiple hydrogen-bonding centres, and regions between short hydrophilic side chains and the main-chain CO or NH groups of the same or nearby residues. Major secondary structural features are not disrupted by hydration. However, the free CO groups at the C terminii and, to a lesser extent, the NH groups at the N terminii of helices provide favourable sites for water interactions, as do reverse turns and regions which connect β-structure and helices. The hydration shell consists of discontinuous networks of water molecules, the maximum number of molecules in a network being ten. The substrate-binding cleft is heavily hydrated, as is the main loop region which is stabilized by water interactions. The protein molecules are close packed in the crystals with a molecular coordination number of 14. Arginyl residues are extensively involved in intermolecular hydrogen bonds and water bridges. The water molecules in the crystal are organized into discrete clusters. A distinctive feature of the clusters is the frequent occurrence of three-membered rings. The protein molecules undergo substantial rearrangement during the transformation from the native to the low-humidity form. The main-chain conformations in the two forms are nearly the same, but differences exist in the side-chain conformation. The differences are particularly pronounced in relation to Trp 62 and Trp 63. The shift in Trp 62 is especially interesting as it is also known to move during inhibitor binding.

Free Energy Barriers for Escape of Water Molecules from Protein Hydration Layer

Free energy barriers separating interfacial water molecules from the hydration layer at the surface of a protein to the bulk are obtained by using the umbrella sampling method of free energy calculation. We consider hydration layer of chicken villin head piece (HP-36) which has been studied extensively by molecular dynamics simulations. The free energy calculations reveal a strong sensitivity to the secondary structure. In particular, we find a region near the junction of first and second helix that contains a cluster of water molecules which are slow in motion, characterized by long residence times (of the order of 100 ps or more) and separated by a large free energy barrier from the bulk water. However, these ``slow'' water molecules constitute only about 5-10% of the total number of hydration layer water molecules. Nevertheless, they play an important role in stabilizing the protein conformation. Water molecules near the third helix (which is the important helix for biological function) are enthalpically least stable and exhibit the fastest dynamics. Interestingly, barrier height distributions of interfacial water are quite broad for water surrounding all the three helices (and the three coils), with the smallest barriers found for those near the helix-3. For the quasi-bound water molecules near the first and second helices, we use well-known Kramers' theory to estimate the residence time from the free energy surface, by estimating the friction along the reaction coordinate from the diffusion coefficient by using Einstein relation. The agreement found is satisfactory. We discuss the possible biological function of these slow, quasi-bound (but transient) water molecules on the surface.

Hydration dynamics of protein molecules in aqueous solution: Unity among diversity

Dielectric dispersion and NMRD experiments have revealed that a significant fraction of water molecules in the hydration shell of various proteins do not exhibit any slowing down of dynamics. This is usually attributed to the presence of the hydrophobic residues (HBR) on the surface, although HBRs alone cannot account for the large amplitude of the fast component. Solvation dynamics experiments and also computer simulation studies, on the other hand, repeatedly observed the presence of a non-negligible slow component. Here we show, by considering three well-known proteins (lysozyme, myoglobin and adelynate kinase), that the fast component arises partly from the response of those water molecules that are hydrogen bonded with the backbone oxygen (BBO) atoms. These are structurally and energetically less stable than those with the side chain oxygen (SCO) atoms. In addition, the electrostatic interaction energy distribution (EIED) of individual water molecules (hydrogen bonded to SCO) with side chain oxygen atoms shows a surprising two peak character with the lower energy peak almost coincident with the energy distribution of water hydrogen bonded to backbone oxygen atoms (BBO). This two peak contribution appears to be quite general as we find it for lysozyme, myoglobin and adenylate kinase (ADK). The sharp peak of EIED at small energy (at less than 2 k(B)T) for the BBO atoms, together with the first peak of EIED of SCO and the HBRs on the protein surface, explain why a large fraction (similar to 80%) of water in the protein hydration layer remains almost as mobile as bulk water Significant slowness arises only from the hydrogen bonds that populate the second peak of EIED at larger energy (at about 4 k(B)T). Thus, if we consider hydrogen bond interaction alone, only 15-20% of water molecules in the protein hydration layer can exhibit slow dynamics, resulting in an average relaxation time of about 5-10 ps. The latter estimate assumes a time constant of 20-100 ps for the slow component. Interestingly, relaxation of water molecules hydrogen bonded to back bone oxygen exhibit an initial component faster than the bulk, suggesting that hydrogen bonding of these water molecules remains frustrated. This explanation of the heterogeneous and non-exponential dynamics of water in the hydration layer is quantitatively consistent with all the available experimental results, and provides unification among diverse features.

Hydrophobic hydration driven self-assembly of curcumin in water: Similarities to nucleation and growth under large metastability, and an analysis of water dynamics at heterogeneous surfaces

As the beneficial effects of curcumin have often been reported to be limited to its small concentrations, we have undertaken a study to find the aggregation properties of curcumin in water by varying the number of monomers. Our molecular dynamics simulation results show that the equilibrated structure is always an aggregated state with remarkable structural rearrangements as we vary the number of curcumin monomers from 4 to 16 monomers. We find that the curcumin monomers form clusters in a very definite pattern where they tend to aggregate both in parallel and anti-parallel orientation of the phenyl rings, often seen in the formation of beta-sheet in proteins. A considerable enhancement in the population of parallel alignments is observed with increasing the system size from 12 to 16 curcumin monomers. Due to the prevalence of such parallel alignment for large system size, a more closely packed cluster is formed with maximum number of hydrophobic contacts. We also follow the pathway of cluster growth, in particular the transition from the initial segregated to the final aggregated state. We find the existence of a metastable structural intermediate involving a number of intermediate-sized clusters dispersed in the solution. We have constructed a free energy landscape of aggregation where the metatsable state has been identified. The course of aggregation bears similarity to nucleation and growth in highly metastable state. The final aggregated form remains stable with the total exclusion of water from its sequestered hydrophobic core. We also investigate water structure near the cluster surface along with their orientation. We find that water molecules form a distorted tetrahedral geometry in the 1st solvation layer of the cluster, interacting rather strongly with the hydrophilic groups at the surface of the curcumin. The dynamics of such quasi-bound water molecules near the surface of curcumin cluster is considerably slower than the bulk signifying a restricted motion as often found in protein hydration layer. (C) 2014 AIP Publishing LLC.

Characterization of lysozyme crystals with unusually low solvent content

Studies on the low-humidity (88%) forms of tetragonal and monoclinic lysozyme, resulting from water-mediated transformations, have provided a wealth of information on the variability in protein hydration, its structural consequences and the water structure associated with proteins, in addition to facilitating the delineation of the rigid and the flexible regions in the protein molecule and the invariant features in its hydration shell. Surprisingly, monoclinic lysozyme continues to diffract even when the environmental humidity is drastically reduced, thus permitting the structural study of the enzyme at different levels of hydration. As part of a study in this direction, three very low humidity forms, two of them occuring at a nominal relative humidity of 38% and the other at 5% relative humidity, have been characterized. These have unprecedented low solvent contents of 16.9, 17.6 and 9.4%, respectively, as determined by the Matthews method.

Sensitivity of Water Dynamics to Biologically Significant Surfaces of Monomeric Insulin: Role of Topology and Electrostatic Interactions

In addition to the biologically active monomer of the protein insulin circulating in human blood, the molecule also exists in dimeric and hexameric forms that are used as storage. The insulin monomer contains two distinct surfaces, namely, the dimer forming surface (DFS) and the hexamer forming surface (HFS), that are specifically designed to facilitate the formation of the dimer and the hexamer, respectively. In order to characterize the structural and dynamical behavior of interfacial water molecules near these two surfaces (DFS and HFS), we performed atomistic molecular dynamics simulations of insulin with explicit water. Dynamical characterization reveals that the structural relaxation of the hydrogen bonds formed between the residues of DFS and the interfacial water molecules is faster than those formed between water and that of the HFS. Furthermore, the residence times of water molecules in the protein hydration layer for both the DFS and HFS are found to be significantly higher than those for some of the other proteins studied so far, such as HP-36 and lysozyme. In particular, we find that more structured water molecules, with higher residence times (similar to 300-500 ps), are present near HFS than those near DFS. A significant slowing down is observed in the decay of associated rotational auto time correlation functions of O-H bond vector of water in the vicinity of HFS. The surface topography and the arrangement of amino acid residues work together to organize the water molecules in the hydration layer in order to provide them with a preferred orientation. HFS having a large polar solvent accessible surface area and a convex extensive nonpolar region, drives the surrounding water molecules to acquire predominantly an outward H-atoms directed, clathrate-like structure. In contrast, near the DFS, the surrounding water molecules acquire an inward H-atoms directed orientation owing to the flat curvature of hydrophobic surface and the interrupted hydrophilic residual alignment. We have followed escape trajectory of several such quasi-bound water molecules from both the surfaces that reveal the significant differences between the two hydration layers.

Modular design of synthetic protein mimics. Crystal structures, assembly, and hydration of two 15- and 16-residue apolar, leucyl-rich helical peptides

Two IS- and 16-residue peptides containing a-aminoisobutyric acid (Aib) have been synthesized, as part of a strategy to construct stereochemically rigid peptide helices, in a modular approach to design of protein mimics. The peptides Boc-(Val-Ala-Leu-Aib),-OMe ( I ) and Boc-Val-Ala-Leu-Aib-Val-Ala-Leu-(Val-Ala-Leu-Aib()11z)- OhaMvee been crystallized.Both crystals are stable only in the presence of mother liquor or water. The crystal data are as follows. I: C78H140N16019~2H20,P2,, a = 16.391 (3) A, b = 16.860 (3) A, c = 18.428 (3) A, p = 103.02 (I)O, Z = 2, R = 9.6% for 3445 data with lFol >30(F), resolution 0.93 A. 11: C7,Hl,,N,S018.7.5H,0, C2221, a = 18.348 ( 5 ) A, b = 47.382 (1 1) A, c = 24.157 ( 5 ) A, Z =8, R = l0,6%, for 3147 data with lFol > 3a(F), resolution 1.00 A. The 15-residue peptide (11) is entirely a helical, while the 16-residue peptide ( I ) has a short segment of 310 helix at the N terminus. The packing of the helices in the crystals is rather incfficicnt with no particular attractions between Leu-Leu side chains, or any other pair. Both crystals have fairly large voids, which are filled with water molecules in a disordered fashion. Water molecule sites near the polar head-to-tail regions are well detcrmined, those closer to the hydrophobic side chains less so and a number of possible water sites in the remaining "empty" space are not determined. No interdigitation of Leu side chains is observed in the crystal as is hypothesized in the "leucine zipper" class of DNA binding proteins.

Detection of the protein dimers, multiple monomeric states and hydrated forms of Plasmodium falciparum triosephosphate isomerase in the gas phase

Dimeric and monomeric forms of the enzyme triosephosphate isomerase (TIM) from Plasmodium falciparum (Pf) have been detected under conditions of nanoflow by electrospray mass spectrometry. The dimer (M = 55 663 Da) exhibits a narrow charge state distribution with intense peaks limited to values of 18(+) to 21(+), maximal intensity being observed for charge states 19(+) and 20(+). A monomeric species with a charge state distribution ranging from 11(+) to 16(+) is also observed, which may be assigned to folded dissociated subunits. Complete dimer dissociation results under normal electrospray condition. The effects of solution pH and source temperature have been investigated. The observation of four distinct charge state distributions which may be assigned to a dimer, folded monomer, partially folded monomer and unfolded monomer is reported. Circular dichromism and fluorescence studies of Pf TIM at low pH support the retention of substantial secondary and tertiary structures. Satellite peaks in mass spectra corresponding to hydrated species are also observed and isotope shift upon deuteration is demonstrated. The analysis of all available independent crystal structures of Pf TIM and TIMs from other organisms permits identification of structurally conserved water molecules. Hydration observed in the dimer and folded monomeric forms in the gas phase may correspond to these conserved sites.

Structural Transformations in Protein Crystals caused by controlled Dehydration

Recent experiments in this laboratory on structural transformations caused by controlled dehydration of protein crystals have been reviewed. X-ray diffraction patterns of the following crystals have been examined under varying conditions of environmental humidity in the relative humidity range of 100-75%: a new crystal form of bovine pancreatic ribonuclease A grown from acetone solution in tris buffer (I), the well-known monoclinic form of the protein grown from aqueous ethanol (II), the same form grown from a solution of 2-methyl pentan-2,4-diol in phosphate buffer (III), tetragonal (IV), orthorhombic (V), monoclinic (VI) and triclinic (VII) hen egg white lysozyme, porcine 2 Zn insulin (VIII), porcine 4 Zn insulin (IX) and the crystals of concanavalin A(X). I, II, IV, V and VI undergo one or more transformations as evidenced by discontinuous changes in the unit cell dimensions, the diffraction pattern and the solvent content. Such water-mediated transformations do not appear to occur in the remaining crystals in the relative humidity range explored. The relative humidity at which the transformation occurs is reduced when 2-methyl pentan-2,4-diol is present in the mother liquor. The transformations are affected by the crystal structure but not by the amount of solvent in the crystals. The X-ray investigations reviewed here and other related investigations emphasize the probable importance of water-mediated transformations in exploring hydration of proteins and conformational transitions in them.

Pi-turns in proteins and peptides: Classification, conformation, occurrence, hydration and sequence.

The i + 5-->i hydrogen bonded turn conformation (pi-turn) with the fifth residue adopting alpha L conformation is frequently found at the C-terminus of helices in proteins and hence is speculated to be a "helix termination signal." An analysis of the occurrence of i + 5-->i hydrogen bonded turn conformation at any general position in proteins (not specifically at the helix C-terminus), using coordinates of 228 protein crystal structures determined by X-ray crystallography to better than 2.5 A resolution is reported in this paper. Of 486 detected pi-turn conformations, 367 have the (i + 4)th residue in alpha L conformation, generally occurring at the C-terminus of alpha-helices, consistent with previous observations. However, a significant number (111) of pi-turn conformations occur with (i + 4)th residue in alpha R conformation also, generally occurring in alpha-helices as distortions either at the terminii or at the middle, a novel finding. These two sets of pi-turn conformations are referred to by the names pi alpha L and pi alpha R-turns, respectively, depending upon whether the (i + 4)th residue adopts alpha L or alpha R conformations. Four pi-turns, named pi alpha L'-turns, were noticed to be mirror images of pi alpha L-turns, and four more pi-turns, which have the (i + 4)th residue in beta conformation and denoted as pi beta-turns, occur as a part of hairpin bend connecting twisted beta-strands. Consecutive pi-turns occur, but only with pi alpha R-turns. The preference for amino acid residues is different in pi alpha L and pi alpha R-turns. However, both show a preference for Pro after the C-termini. Hydrophilic residues are preferred at positions i + 1, i + 2, and i + 3 of pi alpha L-turns, whereas positions i and i + 5 prefer hydrophobic residues. Residue i + 4 in pi alpha L-turns is mainly Gly and less often Asn. Although pi alpha R-turns generally occur as distortions in helices, their amino acid preference is different from that of helices. Poor helix formers, such as His, Tyr, and Asn, also were found to be preferred for pi alpha R-turns, whereas good helix former Ala is not preferred. pi-Turns in peptides provide a picture of the pi-turn at atomic resolution. Only nine peptide-based pi-turns are reported so far, and all of them belong to pi alpha L-turn type with an achiral residue in position i + 4. The results are of importance for structure prediction, modeling, and de novo design of proteins.

Hydration and distortion of peptide helices in crystals. α-Helical structure of a dodecapeptide, Boc-(Ala-Leu-Aib)4-OMe

The dodecapeptide Boc-(Ala-Leu-Aib)(4)-OMe crystallized with two independent helical molecules in a triclinic cell. The two molecules are very similar in conformation, with a 3(10)-helix turn at the N-terminus followed by an alpha-helix, except for an elongated N(7)...O(3) distance in both molecules. All the helices in the crystal pack in a parallel motif. Eleven water sites have been found in the head-to-tail region between the apolar helices that participate in peptide-water hydrogen bonds and a network of water-water hydrogen bonds. The crystal parameters are as follows: 2(C58H104N12O15)+ca. 10H(2)O, space group P1 with a = 12.946(2), b = 17.321(3), c = 20.465(4) Angstrom, alpha = 103.12(2), beta = 105.63(2), gamma = 107.50(2)degrees, Z = 2, R = 10.9% for 5152 data observed > 3 sigma(F), resolution 1.0 Angstrom. In contrast to the shorter sequences [Karle et al. (1988)Proc. Natl. Acad. Sci. USA 85, 299-303] and Boc-(Ala-Leu-Aib)(2)-OMe [Karle et al. (1989) Biopolymers 28, 773-781], no insertion of a water molecule into the helix is observed. However, the elongated N---O distance between Ala(7) NH and Aib(3) CO in both molecules (molecule A, 3.40 Angstrom; molecule B, 3.42 Angstrom) is indicative of an incipient break in the helices. (C) Munksgaard 1994.

Interaction of water molecules in non-identical protein structures

The interaction of the protein atoms with the surrounding water oxygen atoms has been computed for 392 protein chains from 369 protein structures belonging to 90% non-homologous high resolution (<= 1.5 angstrom) protein Structures with a crystallographic R-factor <= 20%. The percentage composition of the polar atoms is found to be 36.3%. An average of 82.55% of water oxygen atoms are found to be in the primary hydration shell and 15.12% in the secondary hydration shell. The average Percentage of interactions of water oxygen atoms with the polar atoms of the main chain and side chain are 54% and 46%. respectively. The interaction of the acidic residues, aspartate and glutamate, with the water oxygen atoms is more when compared to that of the other residues.

Rigid and flexible region in lysozyme and the invariant features in its hydration shell

Water-mediated transformations provide a useful handle for exploring the flexibility in protein molecules and the invariant features in their hydration shells. Low-humidity monoclinic hen egg white lysozyme, resulting from such a transformation, has perhaps the lowest solvent content observed in any protein crystal so far and has a well-ordered structure. A detailed comparison involving this structure, low-humidity tetragonal lysozyme, and the other available refined crystal structures of the enzyme permits the delineation of the relatively rigid, moderately flexible and highly flexible regions of the molecule. The relatively rigid region forms a contiguous structural unit close to the molecular centroid and encompasses parts of of the main beta-structure and three alpha-helices. The hydration shell of the protein contains 30 invariant water molecules. Many of them are involved in holding different parts of the molecule together or in stabilizing local structure. Five of the six invariant water molecules attached to the substrate-binding region form part of a water cluster contiguous with the side-chains of the catalytic residues Glu-35 and Asp-52.

Biological water: Femtosecond dynamics of macromolecular hydration

The unique features of a macromolecule and water as a solvent make the issue of solvation unconventional, with questions about the static versus dynamic nature of hydration and the, physics of orientational and translational diffusion at the boundary. For proteins, the hydration shell that covers the surface is critical to the stability of its structure and function. Dynamically speaking, the residence time of water at the surface is a signature of its mobility and binding. With femtosecond time resolution it is possible to unravel the shortest residence times which are key for the description of the hydration layer, static or dynamic. In this article we review these issues guided by experimental studies, from this laboratory, of polar hydration dynamics at the surfaces of two proteins (Subtilisin Carlsberg (SC) and Monellin). The natural probe tryptophan amino acid was used for the interrogation of the dynamics, and for direct comparison we also studied the behavior in bulk water - a complete hydration in 1 ps. We develop a theoretical description of solvation and relate the theory to the experimental observations. In this - theoretical approach, we consider the dynamical equilibrium in the hydration shell, defining the rate processes for breaking and making the transient hydrogen bonds, and the effective friction in the layer which is defined by the translational and orientational motions of water molecules. The relationship between the residence time of water molecules and the observed slow component in solvation dynamics is a direct one. For the two proteins studied, we observed a "bimodal decay" for the hydration correlation function, with two primary relaxation times: ultrafast, typically 1 ps or less, and longer, typically 15-40 ps, and both are related to the residence time at the protein surface, depending on the binding energies. We end by making extensions to studies of the denatured state of the protein, random coils, and the biomimetic micelles, and conclude with our thoughts on the relevance of the dynamics of native structures to their functions.