High-pressure protein crystallography of hen egg-white lysozyme

The crystal structure of hen egg-white lysozyme (HEWL) was analyzed under pressures of up to 950 MPa. The high pressure modified the conformation of the molecule and induced a novel phase transition in the tetragonal crystal of HEWL.


Introduction
High pressure induces structural changes in proteins, such as the dissociation of oligomers, conformational changes and denaturation. For this reason, high-pressure structural analysis has been widely used to study important characteristics of proteins such as their stability, folding and aggregation (see, for example, Balny, 2006). Because the effect of pressure on protein structures helps us to understand their dynamics, structural studies based on high-pressure nuclear magnetic resonance (NMR) and X-ray crystallography have recently evolved (Cioni & Gabellieri, 2011;Fourme et al., 2012). In the development of high-pressure X-ray protein crystallography (HPPX), the introduction of the diamond anvil cell (DAC) was a key step (Katrusiak & Dauter, 1996;Fourme et al., 2001;Girard et al., 2007). Using a DAC, the useful pressure range was extended by almost an order of magnitude relative to beryllium cells, which had previously been used to study protein crystals at pressures of less than 100 MPa. For example, the use of DACs permitted structural studies to be made of bovine copper zinc superoxide dismutase (SOD) at 1 GPa (Fourme et al., 2001) and of Cowpea mosaic virus at 330 MPa (Girard et al., 2005). Our group also developed a HPPX environment at the Photon Factory, Japan (Chavas et al., 2013), with which we studied the pressure-induced structural changes of 3-isopropylmalate dehydrogenase (IPMDH) up to 650 MPa (Nagae et al., 2012). In particular, the results for SOD and IPMDH show that the pressure-induced structural changes observed by HPPX are directly related to protein activity or function.
Compared with NMR, HPPX is a superior method because it is capable of directly detecting pressure-induced changes in hydration structure. Collins et al. (2005) reported water filling ISSN 1399-0047 of the nonpolar interior cavity of L99A mutant T4 lysozyme. As the initial steps of pressure denaturation, Nagae et al. (2012) also observed water penetration into the hydrophobic cavity at the dimer interface and a cleft on the surface of IPMDH. To further study high-pressure effects in proteins, including pressure-induced changes in hydration structure, we used HPPX to study the well known tetragonal crystal of HEWL. The first high-pressure studies of the crystal structure of HEWL were made at 100 MPa using a beryllium cell (Kundrot & Richards, 1987). A DAC was then used to measure the compressibility of orthorhombic crystals up to 1.0 GPa (Katrusiak & Dauter, 1996), and a higher pressure structure of HEWL at 820 MPa was determined using a DAC (Fourme et al., 2001(Fourme et al., , 2003. To our knowledge, however, no details of the high-pressure structure of HEWL have yet been reported. HEWL is one of the best characterized enzymes. Its reaction mechanism was discussed for many years, and was settled by evidence for a glycosyl-enzyme covalent intermediate at Asp52 (Vocadlo et al., 2001). However, basic questions remain for the other catalytic residue, Glu35. The protonated state of the carboxyl group of Glu35 or its abnormally elevated pK a is known to be important to its catalytic reaction (Phillips, 1967). The pK a of Glu35 has been experimentally confirmed to range from 6.0 to 6.8, in contrast to the standard value of 4.1 for glutamic acid (Kuramitsu & Hamaguchi, 1980;Webb et al., 2011). This abnormally high pK a for Glu35 is believed to result from a hydrophobic environment mainly caused by Trp108 (Inoue et al., 1992). It also provides a good benchmark for computer simulations of pK a (Wallace & Shen, 2011;Goh et al., 2014). Furthermore, Niimura et al. (1997) and Bon et al. (1999) studied the structure of lysozyme using neutron crystallography, which can observe H atoms. However, no structural arguments have yet appeared to explain the high pK a of Glu35 consistently. In contrast, the present study of the highpressure structure of lysozyme detects higher energy conformational states around Glu35 produced by pressure; these states cannot be detected in ambient crystallographic environments.
We have also observed a novel pressure-induced phase transition of the tetragonal crystal of HEWL. A disorderorder transition arising from pressure has previously been reported for Cowpea mosaic virus by Fourme et al. (2002). To the best of our knowledge, however, this is the first report of a pressure-induced crystal-to-crystal phase transition observed in a protein crystal.

Crystallization
For this study, we used tetragonal crystals (space group P4 3 2 1 2) because their higher symmetry accommodates the limitations imposed on the oscillation angle by the aperture of the DAC, which was approximately 70 for the entrance and exit apertures for the X-ray beams. Tetragonal crystals of HEWL (Wako Pure Chemical Industries, Osaka, Japan) were obtained by the batch crystallization method from a solution 7.8 (8.0) 6.0 (6.0) 6.2 (6.4) 7.6 (8.0) 5.0 (5.1) 6.3 (6.5) 6.7 (7.0) 6.1 (6.0) 5.9 (6. consisting of 40 mg ml À1 lysozyme, 50 mM sodium acetate buffer pH 4.5, 0.8 M sodium chloride at 293 K. In order to discuss the effects of substrate binding on the behaviour of Glu35, a complex of HEWL with tetra-N-acetylchitotetraose [(GlcNAc) 4 ; Dextra Laboratories, Berkshire, England] was also crystallized from 25 mg ml À1 lysozyme with a protein: substrate ratio of 1:1.2 and the same buffer conditions as the free enzyme. Because the solubility of tetragonal crystals of lysozyme increases with pressure (Suzuki et al., 2002), the salt concentration of the solution was increased to 1.5 M before pressurization.

High-pressure protein crystallography experiment
HPPX experiments were performed on beamline AR-NW12A at the Photon Factory. To collect high-completeness data sets at each pressure, we inserted two or three crystals into the DAC at different orientations. To prevent the crystals from moving in the DAC sample chamber, the crystals were fixed with tied fibres from cigarette filters (Nagae et al., 2012). The pressure in the DAC was gradually increased to the target pressure in about 10-15 min, and the DAC was left for about 10 min. The pressure was then measured before and after each X-ray measurement using the wavelength shift of ruby fluorescence. Diffraction data sets were collected from HEWL crystals at room temperature at pressures of up to 950 MPa. The crystals diffracted reasonably well at pressures as high as 950 MPa, but no diffraction was detected at 1 GPa. The threedimensional crystal structures were determined at a resolution of 1.5-1.9 Å . Data collection at atmospheric pressure was performed with a Rigaku FR-E Cu K X-ray source and a R-AXIS VII detector using a glass capillary. Table 1 summarizes the data-collection and scaling statistics for all diffraction data sets.

Structure analysis
The diffraction patterns were indexed, integrated and scaled using HKL-2000 (Otwinowski & Minor, 1997). During the integration, the mosaicity was used to check for radiation damage to the crystals. Below 890 MPa, the mosaicities before exposure were less than approximately 0.05 . At 950 MPa, however, the mosaicity increased to approximately 0.6 because a crystal-to-crystal phase transition occurs between 890 and 950 MPa, as discussed below. The structural analysis was based only on frames acquired at the lower mosaicity. The initial crystal structure of HEWL at 0.1 MPa (atmospheric pressure) was solved using MOLREP (Vagin & Teplyakov, 2010) with PDB entry 2lyz (Diamond, 1974) as the search model. For the complex with (GlcNAc) 4 , PDB entry 1lzc (Maenaka et al., 1995) was used. The structure was refined using REFMAC5 (Murshudov et al., 2011) implemented in the CCP4 suite  and then by manual fitting with Coot (Emsley & Cowtan, 2004). The crystal structures of HEWL at different high pressures were solved using the structure at atmospheric pressure as the starting model, and were then improved by several rounds of refinement with REFMAC5 and Coot. At 950 MPa (i.e. after the phase transition), the space group changes and the crystals were merohedrally twinned. The structure solution and refinement were therefore performed on an untwinned data set using DETWIN in CCP4 (the twin fraction was 0.36). The four structures were also searched by MOLREP and refined with REFMAC5 and Coot. Water molecules were added using the solvent mode of ARP/wARP (Lamzin & Wilson, 1993). We removed water with B factors greater than 60 Å 2 , map root-mean-square deviation (r.m.s.d.) levels less than 1.00 e Å À3 and closest contacts less than 2.30 Å or greater than 3.50 Å . The refinement statistics are compiled in Table 2.  jF obs j À jF calc j = P hkl jF obs j. ‡ R free is calculated using 5% of the data that were randomly chosen and excluded from the refinement.
The solvent-excluded volumes of HEWL were determined using VOIDOO with a probe radius of 1.4 Å (Kleywegt & Jones, 1994). The volumes of the internal cavities were calculated using CASTp with a probe radius of 0.9 Å (Dundas et al., 2006). The figures were generated using the PyMOL visualization software (DeLano, 2002). The displacement vectors shown in Figs. 2 and 4 were generated using the modevector.py script from the PyMOLWiki website (http:// www.pymolwiki.org) with a 0.5 Å cutoff, and the differencedistance maps were calculated using CMView (Vehlow et al., 2011). The internal cavities in Figs. 11 and 12 were generated using the HOLLOW tool (Ho & Gruswitz, 2008). Finally, the secondary-structure assignments of HEWL used in this paper are based on UniProt (http://www.uniprot.org/uniprot/ P00698).    Fig. 1(a) shows the relative variations in the unit-cell parameters and the unit-cell volume as a function of pressure. The unit-cell volume decreases from 2.38 Â 10 6 Å 3 at atmo-spheric pressure to 2.19 Â 10 6 Å 3 at 890 MPa. The compressibility of the unit cell is 9.4 Â 10 À2 GPa À1 , which is almost identical to the values of 9.8 Â 10 À2 GPa À1 reported by Fourme et al. (2001) and 9.4 Â 10 À2 GPa À1 reported by Ascone, Savino et al. (2010). As previously reported by Fourme et al. (2001)    anisotropic. The a axis decreases in size much more than the c axis. The compressibility of the a axis is 4.6 Â 10 À2 GPa À1 , whereas the length of the c axis remains almost constant.

Compressibility of the unit cell and the HEWL molecule
The molecular and cavity volumes of HEWL are 2.49 Â 10 5 and 7.48 Â 10 3 Å 3 at atmospheric pressure and 2.40 Â 10 5 and 3.67 Â 10 3 Å 3 at 890 MPa, respectively (Fig. 1b). Both the molecular and cavity volumes decrease approximately linearly with pressure up to 710 MPa. The compressibility of the molecular volume is 4.8 Â 10 À2 GPa À1 , which is almost equal to the published value of 4.7 Â 10 À2 GPa À1 (Kundrot & Richards, 1987), and the compressibility of the cavity volume is 0.80 Â 10 À2 GPa À1 . However, from 710 to 890 MPa, which is just below the pressure of the phase transition, the molecular and cavity volumes remain almost constant. Fig. 2 shows the structural changes of the molecule at 890 and 950 MPa from the structure at ambient pressure. Each structure was superimposed based on the C atoms of all residues 1-129, and the r.m.s.d. was analyzed with the CCP4 program SUPERPOSE (Krissinel & Henrick, 2004). Under high pressures, the structures become distorted with respect to the structure at ambient pressure. At 890 MPa, the -strand domain (residues 42-82) and the -helix domain (residues 1-41 and 83-129) move towards each other to close the active site (Fig. 3a).

Pressure-induced structural changes in molecules
Beyond this pressure range, a crystal-to-crystal phase transition occurs (Fig. 4). At 950 MPa, the crystal symmetry decreases with respect to that at lower pressure, and the space group changes from P4 3 2 1 2 to P4 3 . The a axis of the new P4 3 crystal corresponds to the diagonal (a + b) of the original P4 3 2 1 2 crystal lattice. Four molecules (A, B, C and D) are included in the asymmetric unit of the new P4 3 crystal, whereas these four molecules are symmetry mates in the original (low-pressure) P4 3 2 1 2 crystal, which has one molecule in the asymmetric unit. Figs. 4(b) and 4(c) show the general displacement of molecules on going from 0.1 to 950 MPa. The results also indicate that when the unit-cell volume is calculated in terms of the original P4 3 2 1 2 unit cell, its volume dramatically decreases after the phase transition (Fig. 1a). In particular, the c axis is significantly shortened. We estimate the average r.m.s.d. of the C -atom positions between 0.1 and 890 MPa to be 0.42 Å . However, the corresponding r.m.s.d.s for molecules A, B, C and D after the phase transition (at 950 MPa) are estimated to be 0.57, 0.70, 0.72 and 0.63 Å , respectively. One of the four molecules in the asymmetric unit, molecule C, is dramatically more distorted, and the active-site groove becomes significantly more closed (Figs. 2b and 3b).   The Wilson B factor as a function of pressure.
We analyzed the domain movement of molecule C between 0.1 and 950 MPa using DynDom (Hayward & Berendsen, 1998). The hinge regions of this movement are residues 41 and 42 and residues 82 and 83 (coloured green in Fig. 2c). The former occurs between 1 and 2 and the latter occurs in the middle of 4. The degree of closure motion is 59.2%. The rotational angle and conformational translation between the -domain and the -domain of the high-pressure (950 MPa) molecule C versus the ambient (0.1 MPa) structure are approximately 6.5 and 0.5 Å , respectively.

Reduction in temperature factor
As shown in Fig. 5, the Wilson B factors are reduced upon increasing the pressure from 0.1 to 890 MPa. In particular, the reduction between 0.1 to 190 MPa is significant. This initial decrease is also observed in the molecular and cavity volume reduction shown in Fig. 2. Fig. 6 shows the average temperature factors for all atoms in each residue at 0.1, 890 and 950 MPa. At 890 MPa, the decrease in the average B factor is uniform for the overall structure. At 950 MPa, despite the lower resolution than at 0.1-890 MPa, the average B factor also decreases across the molecules. However, the B factors around 2-3 and 6 do not decrease with respect to their values at 0.1 MPa.

Changes in hydration structure
3.4.1. Molecule surface. The number of water molecules observed gradually increases with pressure. For example, we assign 171 water molecules at 890 MPa compared with 100 at 0.1 MPa (Fig. 7). Among these water structures, 81 water molecules are conserved at both 0.1 and 890 MPa, and 90 new water molecules are found at 890 MPa, mostly located at the back surface of the active site. In addition, hydration develops on a hydrophobic surface of the protein, supported by hydrogen-bond networks. Examples of such water structures appear beside the hydrophobic residues Ile78 and Pro79 and Phe34, Phe38 and Trp123 (Fig. 8). The newly incorporated water molecules hydrogen-bond to protein residues or to the original water molecules that were present at 0.1 MPa. Sidechain conformational changes also help to anchor the network of new water molecules. For example, at 890 MPa, the sidechain conformations of Asn103 are modified to allow hydrogen-bonding to two new water molecules (Fig. 9). However, the number of observed water molecules decreases again at 950 MPa (Fig. 7), which may reflect an increase in disorder caused by the large movements incurred during the pressure-induced phase transition. This interpretation is also supported by the lower resolution of the diffraction data at 950 MPa.
3.4.2. Internal cavity. As the pressure increases, most of the internal cavities are compressed (Figs. 1b and 10). However, the volume of the cavity above Trp108 does not decrease monotonically with pressure. Instead, the volume of this cavity initially decreases from 0. Number of water molecules as a function of pressure. At 950 MPa, the number of water molecules is the average over the four molecules in the asymmetric unit.

Figure 8
Water molecules observed on the hydrophobic surface of HEWL at 890 MPa. (a) Hydration structure around Ile78 and Pro79. (b) Hydration structure around Phe34, Phe38 and Trp123. Upon increasing the pressure, several water molecules appear and form a hydrogen-bond network with each other over the hydrophobic surface. The surface in red is negatively charged, that in blue is positively charged and that in white is neutral. Water molecules found only at 0.1 MPa are coloured green, those found only at 890 MPa are coloured purple and those found at both 0.1 and 890 MPa are coloured cyan. 49.6, 48.2, 46.8 and 39.5 Å 3 at pressures of 0.1, 190, 380, 500 and 600 MPa, respectively). The cavity then swells to 44.6 Å 3 at 710 MPa and is compressed again to 42.4 and 43.0 Å 3 at 800 and 890 MPa, respectively. This increase in cavity volume comes from water penetrating into the cavity and from the movement of the surrounding residues. In fact, at pressures greater than 710 MPa a water molecule penetrates into the cavity (Fig. 11a). This water molecule appears to interact with Trp108 via a lone pair-interaction (Egli & Sarkhel, 2007;Jain et al., 2009). Upon increasing the pressure to 800 MPa, Glu35 becomes a dual conformer. The side chain of the new conformer is inserted into the cavity so that it hydrogen-bonds to the penetrating water molecule.
Above the phase-transition pressure, Glu35 of the four molecules adopts different conformations (Fig. 11b). In molecule A Glu35 takes on dual conformations, as in the 890 MPa structure. In molecules B, C and D, however, the side chains have a single conformation and protrude into the cavity to hydrogen-bond to the newly present water molecules. In molecules C and D additional water molecules penetrate above the previously present water molecule that interact with Internal cavities of HEWL from 0.1 to 890 MPa and at 920 MPa for the (GlcNAc) 4 complex. Most of the cavities are compressed as the pressure increases. The cavity for which the volume increases upon penetration by water is indicated by a black arrow in the 710 MPa structure. The bound (GlcNAc) 4 is coloured purple.

Figure 9
Change in the structure of HEWL upon hydration on going from 0.1 to 890 MPa: (a) 0.1 MPa and (b) 890 MPa. The conformation of Asn103 changes to hydrogen-bond to a newly arrived water molecule at 890 MPa. The protein surface is coloured as in Fig. 11. Trp108 and make a hydrogen-bond network between Trp108 and the main-chain carbonyl O atom of Trp28.
However, in the case of the complex of HEWL with (GlcNAc) 4 , the side chain of Glu35 is only in the outward conformation at 920 MPa and no water penetration into the cavity above Trp108 was observed, as shown in Fig. 11(a). A phase transition did not take place even when the crystal was pressurized to 920 MPa.

Pressure-induced phase transition
In the solution state, the HEWL molecule is expected to be denatured at pressures of several hundred MPa (Heremans & Wong, 1985;Jonas, 1990). In the crystal form, however, the HEWL molecule can tolerate pressures of up to 950 MPa. As previously mentioned by Katrusiak & Dauter (1996), the crystal packing might mitigate the effects of pressure. Upon increasing the pressure, a pressure-induced phase transition is observed between 890 and 950 MPa. This result is consistent with a previous HPPX study of HEWL, which reported a radical reduction of the unit-cell volume, diffuse scattering and the absence of diffraction after 15 min at 915 MPa (Fourme et al., 2001). This result of Fourme and coworkers may thus be the onset of the phase transition. At the phase transition, the four molecules in the asymmetric unit move towards each other (Fig. 4). In particular, molecules A and C approach each other along the c axis. Consequently, the c axis may become shorter, as shown in Fig. 1(a). After the phase transition, each structure in the asymmetric unit may embody slightly different high-energy conformations of HEWL, which are difficult to observe using the usual crystallographic methods because of their relatively low abundance at atmospheric pressure. As pointed out by Collins et al. (2011) andFourme et al. (2012), such structures are often essential for protein function. We discuss these aspects below.

Change in the volume and the structure of HEWL
In accordance with Le Chatelier's principle, an increase in pressure favours the reduction of the volume of a system. As shown in Fig. 1(a), the unit-cell volume of tetragonal crystals of HEWL decreases monotonically upon increasing the pressure. Note, however, that the molecular volume and the internal-cavity volume are almost constant from 710 to 890 MPa (Figs. 1b and 11a). At 710 MPa, the unit-cell and molecular volumes are compressed to 93.4 and 96.4% with respect to the corresponding volumes at atmospheric pressure. This is consistent with previous reports of the unit cell exhibiting a compressibility larger than that of the molecules Figure 11 Penetration of water and change in conformation of Glu35 in the internal cavity near Trp108: (a) from 0.1 to 890 MPa and at 920 MPa for the (GlcNAc) 4 complex and (b) for the four molecules at 950 MPa. The internal cavity is represented as a pink surface. The green mesh shows the difference electrondensity map calculated by omitting the water that penetrates and a Glu35 side chain. Maps are contoured at (a) 3.0 and (b) 2.5. Refaee et al., 2003). The partial molar volume (PMV) of a protein is expressed by the sum of the volumes of all the constituent atoms, the volumes of void cavities and the change in volume because of hydration. As viewed from the perspective of the PMV, the system shrinks in going from atmospheric pressure to 710 MPa by compressing the volume of the internal cavities. At pressures greater than 710 MPa, however, the system reduces the PMV mainly by filling the cavity above Trp108 with a water molecule from the bulk region (Fig. 11). This seems to explain why the molecular and internal-cavity volumes are outwardly constant from 710 to 890 MPa, whereas the unit-cell volume decreases mono-tonically. The same trend occurs in 3-isopropylmalate dehydrogenase from Shewanella oneidensis MR-1 (SoIPMDH; Nagae et al., 2012). This phenomenon might depend on the original size of the cavity in question. At 0.1 MPa, the volume of the cavity above Trp108 is 62.5 Å 3 and that at the dimer interface of SoIPMDH is 55.6 Å 3 . While, on the other hand, the hydrophobic cavity of the L99A mutant of T4 lysozyme is large as 160 Å 3 (Collins et al., 2005). In this case, the cavity is sufficiently large to accommodate water molecules, so its volume need not increase.
Despite the disadvantage of water localization in the hydrophobic environment, this penetration by water may be the main way to reduce the PMV of the molecule when the compression of the cavities reaches a physical limit. In studies involving a three-dimensional reference interaction-site model (RISM) theory of molecular solvation, Imai et al. (2007) reported that reduction of the PMV is associated with a pressure-induced structural change in ubiquitin. They attribute the reduction of the PMV primarily to the penetration of water into the hydrophobic core of ubiquitin. Imai & Sugita (2010) subsequently used a molecular-dynamics (MD) simulation to confirm the pressure-induced conformational changes linked to penetration by water. In lysozyme, to cancel out the drawback of water localization in a hydrophobic environment, the water that penetrates the cavity forms a lone pair-interaction with the indole ring of the Trp108 side chain (Figs. 11 and 12). At pressures greater than 800 MPa, the side chain of the active residue Glu35 becomes a dual conformer. This emerging conformer may also reduce the PMV of the molecule and be stabilized by hydrogen bonds to the water molecules that have penetrated the cavity.
Regarding the pressured-induced structural changes of HEWL, Kundrot & Richards (1987) previously reported a difference-distance map between the structure at ambient pressure and at 100 MPa, which indicates that with increasing pressure the distance between the -domains and the -domains decreases. These results provide evidence for a general tendency of the protein to become concentrated closer to its centre of mass, although some differences in directions persist. This tendency to bring the -domains and -domains into proximity corresponds to our results at 890 MPa (Figs. 2a and 3a). However, a large deformation of the -domain is detected only for molecule C at 950 MPa. In particular, molecule C undergoes a relatively large twist and closure (Fig. 2c), as detected by DynDom. This twist of the -domain is consistent with the results of high-pressure NMR studies (Refaee et al., 2003).  an increase in the number of water molecules upon pressurization (Fig. 7). Most new water molecules are localized on the surface of HEWL, which suggests that pressure suppresses atomic fluctuations, thereby reducing the PMV. Actually, the overall B factor decreases with pressure (by 8.7 Å 2 on going from 0.1 to 890 MPa; Fig. 5). A pressure-induced decrease in the B factor was also revealed by structural studies of Cowpea mosaic virus (Girard et al., 2005). MD simulations also indicate a pressure-induced decrease in atomic positional fluctuations (Brunne & van Gunsteren, 1993) and enhanced interactions between water molecules and protein (Marchi & Akasaka, 2001).
The present HPPX study reveals several water molecules on the hydrophobic surface of the protein forming a hydrogenbond network between themselves (Figs. 8 and 9). Such localization may be representative of the conformational change of water molecules driven by the pressure-induced necessity to reduce the PMV. Some residues change their side-chain conformation to interact with structured water molecules and also to secure hydrogen-bond networks at hydrophobic surfaces.

Pressure-induced conformational change of Glu35
As the pressure increases, Glu35, which is one of the two residues with a catalytic carboxyl group, changes its conformation. Between 0.1 and 600 MPa, Glu35 thrusts its side chain towards the active site. At 800 MPa, however, one of the dual conformers of Glu35 assumes a novel conformation in which its side chain faces an inner side of the hydrophobic cavity above Trp108. At higher pressure, for example in the structures of molecules B, C and D at 950 MPa, only the inward conformation is observed (Fig. 11b). The pK a of this glutamic acid is known to be 6.0-6.8, which is abnormally high (Kuramitsu et al., 1977;Webb et al., 2011) compared with the standard pK a of glutamic acid (4.1). This high pK a of Glu35 is critical for it to serve as a general acid catalyst in the glycosidase activity of lysozyme under physiological conditions. The hydrophobicity of Trp108 is considered to contribute to the high pK a of Glu35. Single mutations of Trp108 to the less hydrophobic residues Tyr and Gln decrease the activity of HEWL, and the pK a of Glu35 is reduced by 0.2 and 0.6, respectively (Inoue et al., 1992). To the best of our knowledge, these inward conformations have not yet been observed in crystalline lysozyme, including in neutron diffraction experiments (Niimura et al., 1997;Bon et al., 1999).
We used the empirical methods DEPTH (Tan et al., 2013) and PROPKA  to estimate the pK a of Glu35 in the structures determined by HPPX. As expected, the pK a values of the inward conformation are greater than those of the outward conformation (Table 3). Using DEPTH, the pK a values of the inward Glu35 are estimated to range from 6.3 to 6.5, which is consistent with the experimental results, whereas those of the outer conformation are less than 4.5. PROPKA also assigns greater pK a for the inward conformation than for the outward conformation, but the pK a of the glutamine side chain seems to be overestimated.
In an interesting twist, the refined NMR structure of HEWL shows that the inward or buried conformation of Glu35 is popular in solution (PDB entry 1e8l; Schwalbe et al., 2001). Surprisingly, 45 of the 50 low-energy structures include the inward conformer of Glu35. Thus, the inward conformation of Glu35 appears to be common in the solution state. Why conventional crystallography does not detect the inward conformation of the NMR structure at ambient pressure is not clear. Between NMR and X-ray structures, the root-meansquare deviation around Glu35 is not high; it is less than 1.0 Å (Schwalbe et al., 2001). One likely explanation invokes differences in salt concentrations and molecular packing in the crystal. Whether rightly or wrongly, the structure as determined by NMR does not explain why the side chain of Glu35 prefers an inward conformation. The HPPX results of this study show that water molecules penetrate into the hydrophobic cavity, thereby stabilizing the inward conformation of Glu35 by hydrogen-bonding to it. These results show that Trp108 is important not only for its hydrophobicity, but also for its ability to stabilize the water molecule in the hydrophobic cavity via a lone pair-interaction.
However, the inward conformation of Glu35 and water penetration above Trp108 were never observed in the crystal complexed with (GlcNAc) 4 even at 920 MPa. Pressureinduced phase transition with a large displacement of the molecule is not also observed for the complex crystal. These results show that the flexibility of the HEWL molecule is restricted when substrate is bound to the enzyme molecule. The inward conformation of Glu35 might occur only when the HEWL molecule is free from substrate. Thus, the present study shows that applying high pressure to crystallized HEWL changes the populations of the conformational substates, thereby revealing structures that are less prevalent at atmospheric pressure in the crystal (Collins et al., 2011;Fourme et al., 2012). Studies of these conformational substates can yield valuable information about the activity of proteins, for example, in the catalytic mechanism of enzymes.  Table 3 Estimated pK a value of Glu35 in each conformation using DEPTH and PROPKA.