Crystal structure of hen egg-white lysozyme at a hydrostatic pressure of 1000 atmospheres

The crystal structure of tetragonal hen egg-white lysozyme at a hydrostatic pressure of 1000 atmospheres has been determined by X-ray diffraction to a nominal resolution of 2 A. The crystals, originally grown in 0.83 M-NaCl, had to be transferred to 1.4 M-NaCl to prevent crystal cracking at 300 to 400 atm. The a and b axes of the unit cell contracted by 0.6%, whilst the c axis increased by 0.1%. The unit cell volume contracted by 1.1%. Both the 1 atm and the 1000 atm structures were refined by restrained least-squares to yield final R factors of 14.9% in each case. Since the data were collected by an accurate difference protocol, the change in structure is considered to be more accurate than the absolute structure. The probable accuracy of the atomic shifts is shown to be +/- 0.06 A. The estimated volume decrease of the whole molecule corresponded to an isothermal compressibility of 4.7 X 10(-3) kbar-1. The contraction was non-uniformly distributed. Domain 2 (residues 40 to 88) was essentially incompressible, whilst domain 1 (residues 1 to 39, 89 to 129) had a compressibility of 5.7 X 10(-3) kbar-1. The interdomain region was also compressible. The average B factor decreased about 1 A2 at 1000 atm, but there was a wide range of decreases and increases in individual values. The pressure-induced deformation was analyzed with difference distance matrices. The beta-sheet (residues 42 to 60) and helix 2 (residues 24 to 36) were deformed the least under pressure. The other helices were more deformed and one loop region (residues 61 to 87) actually appeared to expand. The main-chain atoms of the beta-sheet and helix 2 were used to perform a least-squares superposition of the 1 atm and 1000 atm models. The root-mean-square pressure-induced shift for all atoms was 0.2 A, with a few atoms moving more than 1 A. There was no evidence for co-ordinated movement about the hinge axis defined by alpha carbon atoms 38 and 97. The 1 atm and 1000 atm refined structures included 151 and 163 ordered water molecules, respectively. The changes in these ordered water molecules and the mean compressibility of all of the solvent in the crystal will be described elsewhere.     

Conservation of water molecules in an antibody–antigen interaction

The solvation of the antibody–antigen Fv D1.3–lysozyme complex is investigated through a study of the conservation of water molecules in crystal structures of the wild-type Fv fragment of antibody D1.3, 5 free lysozyme, the wild-type Fv D1.3–lysozyme complex, 5 Fv D1.3 mutants complexed with lysozyme and the crystal structure of an idiotope (Fv D1.3)-abti-idiotope (Fv E5.2) complex. In all, there are 99 water molecules common to the wild-type and mutant antibody–lysozyme complexes. The antibody–lysozyme interface includes 25 well-ordered solvent molecules, conserved among the wild-type and mutant Fv D1.3–lysozyme complexes, which are bound directly or through other water molecules to both antibody and antigen. In addition to contributing hydrogen bonds to the antibody–antigen interaction the solvent molecules fill many interface cavities. Comparison with x-ray crystal structures of free Fv D1.3 and free lysozyme shows that 20 of these conserved interface waters in the complex were bound to one of the free proteins. Uo to 23 additional water molecules are also found in the antibody–antigen interface, however these waters do no bridge antibody and antigen and their temperature factors are much higher than those of the 25 well-ordered waters. Fifteen water molecules are displaced to form the complex, some of which are substituted by hydrophilic protein atoms, and 5 water molecules are added at the antibody–antigen interface with the formation of the complex. While the current crystal models of the D1.3–lysozyme complex do not demonstrate the increase in bound waters found in a physico-chemical study of the interaction at decreased water activities, the 25 well-ordered interface water contribute a net gain of 10 hydrogen bonds to complex stability.     

Display and interpretation of solvent electron density distributions in insulin crystals

In macromolecular crystallography, three-dimensional contour surfaces are useful for interactive computer graphics displays of the protein electron density but are less effective for presenting static images of large volumes of solvent density. A raster-based computer graphics program which displays depth-cued projections of continuous density distributions has been developed to analyze the distribution of solvent atoms in macromolecular crystals. Maps of the water distribution in the cubic insulin crystal show some well-ordered waters, which are bound to surrounding protein atoms by multiple hydrogen bonds, and an ill-defined solvent structure at a greater distance from the protein surface. Molecular dynamics calculations were used to assist in the interpretation of the time-varying solvent structure within two enclosed cavities in the crystal. Two water molecules that ligate a sodium ion were almost immobile during the simulation but the majority of water molecules were found to move rapidly between the density maxima identified from the crystallographic refinement.     

Solvating atomic level fine-grained proteins in supra-molecular level coarse-grained water for molecular dynamics simulations

Simulation of the dynamics of a protein in aqueous solution using an atomic model for both the protein and the many water molecules is still computationally extremely demanding considering the time scale of protein motions. The use of supra-atomic or supra-molecular coarse-grained (CG) models may enhance the computational efficiency, but inevitably at the cost of reduced accuracy. Coarse-graining solvent degrees of freedom is likely to yield a favourable balance between reduced accuracy and enhanced computational speed. Here, the use of a supra-molecular coarse-grained water model that largely preserves the thermodynamic and dielectric properties of atomic level fine-grained (FG) water in molecular dynamics simulations of an atomic model for four proteins is investigated. The results of using an FG, a CG, an implicit, or a vacuum solvent environment of the four proteins are compared, and for hen egg-white lysozyme a comparison to NMR data is made. The mixed-grained simulations do not show large differences compared to the FG atomic level simulations, apart from an increased tendency to form hydrogen bonds between long side chains, which is due to the reduced ability of the supra-molecular CG beads that represent five FG water molecules to make solvent-protein hydrogen bonds. But, the mixed-grained simulations are at least an order of magnitude faster than the atomic level ones.     

Structural and thermodynamic analysis of the binding of solvent at internal sites in T4 lysozyme

To investigate the structural and thermodynamic basis of the binding of solvent at internal sites within proteins a number of mutations were constructed in T4 lysozyme. Some of these were designed to introduce new solvent-binding sites. Others were intended to displace solvent from preexisting sites. In one case Val-149 was replaced with alanine, serine, cysteine, threonine, isoleucine, and glycine. Crystallographic analysis shows that, with the exception of isoleucine, each of these substitutions results in the binding of solvent at a polar site that is sterically blocked in the wild-type enzyme. Mutations designed to perturb or displace a solvent molecule present in the native enzyme included the replacement of Thr-152 with alanine, serine, cysteine, valine, and isoleucine. Although the solvent molecule was moved in some cases by up to 1.7 A, in no case was it completely removed from the folded protein. The results suggest that hydrogen bonds from the protein to bound solvent are energy neutral. The binding of solvent to internal sites within proteins also appears to be energy neutral except insofar as the bound solvent may prevent a loss of energy due to potential hydrogen bonding groups that would otherwise be unsatisfied. The introduction of a solvent-binding site appears to require not only a cavity to accommodate the water molecule but also the presence of polar groups to help satisfy its hydrogen-bonding potential. It may be easier to design a site to accommodate two or more water molecules rather than one as the solvent molecules can then hydrogen-bond to each other. For similar reasons it is often difficult to design a point mutation that will displace a single solvent molecule from the core of a protein.     

Effect of various humidity on local dynamic structure of lysozyme in a spin-labeled tetragonal crystal

The dynamics of the side groups of amino acid residues and local conformational changes in the lysozyme molecule upon dehydration and rehydration of lysozyme crystals were studied by the methods of spin label, X-ray diffraction, and molecular dynamics. The His15 residue of lysozyme from chicken egg white was modified by spin label, and spin-labeled tetragonal crystals of the protein were grown. The spatial structure of the covalently bound spin label and its immediate surroundings in the lysozyme tetragonal crystal was determined. The conformation of a fragment of the lysozyme molecule with the spin label on His15, optimized by the method of molecular dynamics, closely agreed with X-ray data. It was found by the X-ray diffraction analysis that a decrease in relative humidity to 40% is accompanied by both a decrease in the unit cell volume by 27% and a change in the diffraction field of roentgenograms from 0.23 to 0.60 HM. The dehydration of spin-labeled lysozyme crystals leads to an anomalous widening of EPR peaks without changes in their position. The dehydration in the humidity range studied has a two-stage character. The decrease in humidity to 75% is accompanied by a sharp change in the parameters measured, and on further decrease in humidity to 40% they change insignificantly. The first stage is caused by the removal of the greater part of molecules of bulk water, and the second stage is due to the removal of the remaining bulk water and possible changes in the dynamics of weakly bound water molecules and their position. The simulation of experimental EPR spectra showed that the anomalous broadening of the spectrum upon dehydration is related to an increase in the dispersion of spin label orientations induced by changes in the network of hydrogen bonds generated by water molecules in the vicinity of the spin label and a possible turn (by no more than 5 degrees) of the entire protein molecule. After rehydration, the physical state of the lysozyme crystal did not return to the starting point.     

Neutron diffraction study of carbonmonoxymyoglobin.

Neutron diffraction data from a crystal of carbonmonoxymyoglobin were refined by PROLSQ, a modern restrained least-squares procedure in reciprocal space, in conjunction with a solvent analysis technique, to a final R-factor of 11.3%. The ligand CO occupies two sites and its binding conformations are distorted from the linear conformation. The N epsilon atom of the distal histidine residue is deprotonated (not deuterated), and a water molecule is bound to the N delta atom of the distal histidine. The side-chain of Lys56 (D6) exists in two alternative charge-binding sites. His24 (B5) and His119 (GH1) share a hydrogen atom. His12 (A10) and His36 (C1) are deprotonated. The deprotonated imidazole ring of His12 (A10) may act as a hydrogen-bond acceptor. The heme group is planar within 0.09 A root-mean-square (r.m.s.) deviation from planarity. The solvent environments for the two propionic acid groups are different. The side-chain of Arg45 (CD3) forms hydrogen bonds with the side-chain of Asp60 (E3) and one of the two propionic acid groups. An average N-2H . . . O angle in helical regions is 147 (+/- 11) degrees. Eleven main-chain amide hydrogen atoms from hydrophobic residues do not exchange with deuterium. The overall atomic occupancy factors for the main-chain and side-chain atoms are quite uniform, at 0.97 (+/- 0.07) and 0.93 (+/- 0.10), respectively, as shown by an occupancy analysis made at the end of the refinement procedure.     

On the edge of the denaturation process: application of X-ray diffraction to barnase and lysozyme cross-linked crystals with denaturants in molar concentrations

Structural data about the early step of protein denaturation were obtained from cross-linked crystals for two small proteins: barnase and lysozyme. Several denaturant agents like urea, bromoethanol or thiourea were used at increasing concentrations up to a limit leading to crystal disruption (>or=2 to 6 M). Before the complete destruction of the crystal order started, specific binding sites were observed at the protein surfaces, an indication that the preliminary step of denaturation is the disproportion of intermolecular polar bonds to the benefit of the agent "parasiting" the surface. The analysis of the thermal factors first agree with a stabilization effect at low or moderate concentration of denaturants rapidly followed by a destabilization at specific weak points when the number of sites increase (overflooding effect).     

Structural studies of mutants of T4 lysozyme that alter hydrophobic stabilization

Multiple replacements at amino acid position 3 of bacteriophage T4 lysozyme have shown that the conformational stability of the protein is directly governed by the hydrophobicity of the residue substituted (Matsumura, M., Becktel, W. J., and Matthews, B. W. (1988) Nature 334, 406-410). Of the 13 mutant lysozymes made by site-directed mutagenesis, two variants, one with valine (I3V) and the other with tyrosine (I3Y), were crystallized and their structures solved. In this report we describe the crystal structures of these variants at 1.7 A resolution. While the structure of the I3V mutant is essentially the same as that of wild-type lysozyme, the I3Y mutant has substantial changes in its structure. The most significant of these are that the side chain of the tyrosine is not accommodated within the interior of the protein and the amino-terminal polypeptide (residues 1-9) moves 0.6-1.1 A relative to the wild-type structure. Using coordinates based on the wild-type and available mutant structures, solvent accessible surface area of residue 3 as well as the adjacent 9 residues in the folded form were calculated. The free energy of stabilization based on the transfer of these residues from a fully extended form to the interior to the folded protein was found to correlate well with the protein stability determined by thermodynamic analysis. The enhanced thermostability of the variant Ile-3----Leu, relative to wild-type lysozyme, can also be rationalized by surface-area calculations based on a model-built structure. Noncrystallization of most lysozyme variants at position 3 appears to be due to disruption of intermolecular contacts in the crystal. The Ile-3----Val variant is closely isomorphous with wild-type and maintains the same crystal contacts. In the Ile-3----Tyr variant, however, a new set of contacts is made in which direct protein-protein hydrogen bonds are replaced by protein-water-protein hydrogen bonds as well as a novel hydrogen bond involving the phenolic hydroxyl of the substituted tyrosine.     

A molecular dynamics study of solvent behavior around a protein

The solvent structure and behavior around a protein were examined by analyzing a trajectory of molecular dynamics simulation of thetrp-holorepressor in a periodic box of water. The calculated selfdiffusion coefficient indicated that the solvent within 10 Å of the protein had lower mobility. Examination of the solvent diffusion around different atoms of different kinds of residues showed no general tendency. Thisfact suggested that the solvent mobility is not influenced significantly bythe kind of the atom or residue they solvated. Distribution analysis aroundthe protein revealed two peaks of water oxygen: a sharp one at 2.8 Å around polar and charged atoms and a broad one at ～3.4 Å aroundapolar atoms. The former was stabilized by water–protein hydrogen bonds, and the latter was stabilized by water-lwater hydrogen bonds, suggesting the existence of a hydrophobic shell. An analysis of protein atom–water radial distribution functions confirmed these shell structures around polar or charged atoms and apolar ones. © 1993 Wiley-Liss, Inc.     

Hydrogen bonding between sugar and protein is responsible for inhibition of dehydration-induced protein unfolding.

The nature of the interaction responsible for the inhibition of protein unfolding and subsequent damage by sugars during dehydration is unclear. The relationship between sample moisture content measured by coulometric Karl Fischer titration and the apparent moisture content predicted by the area of the protein side chain carboxylate band at approximately 1580 cm-1 in infrared spectra of dried protein-sugar samples was examined. For samples in which a high level of native protein structure was retained in the dried solid, the apparent moisture content predicted by the carboxylate band area was greater than the actual moisture content, indicating that protection results from direct sugar-protein hydrogen bonding and not entrapment of water at the protein surface. Further, we show that the degree of structural protection conferred by sucrose and trehalose apparent in second derivative, amide I infrared spectra, correlates with the extent of hydrogen bonding between sugar and protein. The failure of dextran to inhibit dehydration-induced lysozyme unfolding is shown to result from the inability of the polymer to hydrogen bond adequately to the protein. Therefore, formation of an amorphous phase alone is not sufficient to maintain protein structure during dehydration. Glucose hydrogen bonds to a high degree with dried lysozyme, but is incapable of inhibiting lyophilization-induced protein unfolding in the absence of an effective cryoprotectant. However, the addition of polyethylene glycol, which is known to protect proteins during freezing, but not drying, to glucose protected lysozyme structure during lyophilization. Together, these results show that hydrogen bonding between carbohydrate and protein is necessary to prevent dehydration-induced protein damage. However, hydrogen bonding alone is not sufficient to protect proteins during lyophilization in the absence of adequate freezing protection.     

Contribution of hydrogen bonds to the conformational stability of human lysozyme: calorimetry and X-ray analysis of six Ser --> Ala mutants.

To further examine the contribution of hydrogen bonds to the conformational stability of the human lysozyme, six Ser to Ala mutants were constructed. The thermodynamic parameters for denaturation of these six Ser mutant proteins were investigated by differential scanning calorimetry (DSC), and the crystal structures were determined by X-ray analysis. The denaturation Gibbs energy (DeltaG) of the Ser mutant proteins was changed from 2.0 to -5.7 kJ/mol, compared to that of the wild-type protein. With an analysis in which some factors that affected the stability due to mutation were considered, the contribution of hydrogen bonds to the stability (Delta DeltaGHB) was extracted on the basis of the structures of the mutant proteins. The results showed that hydrogen bonds between protein atoms and between a protein atom and a water bound with the protein molecule favorably contribute to the protein stability. The net contribution of one intramolecular hydrogen bond to protein stability (DeltaGHB) was 8.9 +/- 2.6 kJ/mol on average. However, the contribution to the protein stability of hydrogen bonds between a protein atom and a bound water molecule was smaller than that for a bond between protein atoms.     

Continuum electrostatic analysis of preferred solvation sites around proteins in solution

To understand water-protein interactions in solution, the electrostatic field is calculated by solving the Poisson-Boltzmann equation, and the free energy surface of water is mapped by translating and rotating an explicit water molecule around the protein. The calculation is applied to T4 lysozyme with data available on the conservation of solvent binding sites in 18 crystallographically independent molecules. The free energy maps around the ordered water sites provide information on the relationship between water positions in crystal structure and in solution. Results show that almost all conserved sites and the majority of nonconserved sites are within 1.3 A of local free energy minima. This finding is in sharp contrast to the behavior of randomly placed water molecules in the boundary layer, which, on the average, must travel more than 3 A to the nearest free energy minimum. Thus, the solvation sites are at least partially determined by protein-water interactions rather than by crystal packing alone. The characteristic water residence times, obtained from the free energies at the local minima, are in good agreement with nuclear magnetic resonance experiments. Only about half of the potential sites show up as ordered water in the 1.7 A resolution X-ray structure. Crystal packing interactions can stabilize weak or mobile potential sites (in fact, some ordered water positions are not close to free energy minima) or can prevent water from occupying certain sites. Apart from a few buried water molecules that are strong binders, the free energies are not very different for conserved and nonconserved sites. We show that conservation of a water site between two crystals occurs if the positions of protein atoms, primarily contributing to the free energy at the local minimum, do not substantially change from one structure to the other. This requirement can be correlated with the nature of the side chain contacting the water molecule in the site.     

Conservation of solvent-binding sites in 10 crystal forms of T4 lysozyme.

Solvent-binding sites were compared in 10 different crystal forms of phage T4 lysozyme that were refined using data from 2.6 A to 1.7 A resolution. The sample included 18 crystallographically independent lysozyme molecules. Despite different crystallization conditions, variable crystal contacts, changes due to mutation, and varying attention to solvent during crystallographic refinement, 62% of the 20 most frequently occupied sites were conserved. Allowing for potential steric interference from neighboring molecules in the crystal lattice, this fraction increased to 79% of the sites. There was, however, no solvent-binding site that was occupied in all 18 lysozyme molecules. A buried double site was occupied in 17 instances and 2 other internal sites were occupied 15 times. Apart from these buried sites, the most frequently occupied sites were often at the amino-termini of alpha-helices. Solvent molecules at the most conserved sites tended to have crystallographic thermal factors lower than average, but atoms with low B-factors were not restricted to these sites. Although superficial inspection may suggest that only 50-60% (or less) of solvent-binding sites are conserved in different crystal forms of a protein, it appears that many sites appear to be empty either because of steric interference or because the apparent occupancy of a given site can vary from crystal to crystal. The X-ray method of identifying sites is somewhat subjective and tends to result in specification only of those solvent molecules that are well ordered and bound with high occupancy, even though there is clear evidence for solvent bound at many additional sites.(ABSTRACT TRUNCATED AT 250 WORDS)     

Downscaling cake filtration—Correlations between compressibility of lysozyme crystal filter cakes and filter area

Crystallization is a commonly used method for the purification and formulation of proteins. In order to harmonize crystallization with the mechanical separation process, it is necessary to estimate early the separation conditions of protein crystals. To fulfill this requirement, the feed material for filtration is minimized by reducing the filter area. The filtration behavior of lysozyme crystals in pressure filters with two different filter areas was compared. Crystal slurries with different mean crystal sizes and shapes were produced and the influence of the size and shape on the scalability of filtration data was examined. It was found that for different aggregated crystal structures and isometric crystals, the larger the compressibility of the cake was, the larger were the deviations for the two considered filter areas. For needle‐shaped crystals, the compressibility was not subject to deviations, but the absolute filter cake resistance was. Furthermore, the deviations of the cake resistances increased with higher pressures for all product systems. It was indicated that upscaling from the small filter area is possible for low compressible product systems and at low pressures. With high compressible products and for needle‐like crystals, the filtration time is underestimated with the small filter area.     

Vibrational frequency shifts as a probe of hydrogen bonds: thermal expansion and glass transition of myoglobin in mixed solvents

The contribution of hydrogen bonds to protein-solvent interactions and their impact on structural flexibility and dynamics of myoglobin are discussed. The shift of vibrational peak frequencies with the temperature of myoglobin in sucrose/water and glycerol/water solutions is used to probe the expansion of the hydrogen bond network. We observe a characteristic change in the temperature slope of the O–H stretching frequency at the glass transition which correlates with the discontinuity of the thermal expansion coefficient. The temperature-difference spectra of the amide bands show the same tendency, indicating that stronger hydrogen bonding in the bulk affects the main-chain solvent interactions in parallel. However, the hydrogen bond strength decreases relative to the bulk solvent with increasing cosolvent concentration near the protein surface, which suggests preferential hydration. Weaker and/or fewer hydrogen bonds are observed at low degrees of hydration. The central O–H stretching frequency of protein hydration water is red-shifted by 40 cm^–1 relative to the bulk. The shift increases towards lower temperatures, consistent with contraction and increasing strength of the protein-water bonds. The temperature slope shows a discontinuity near 180 K. The contraction of the network has reached a critical limit which leads to frozen-in structures. This effect may represent the molecular mechanism underlying the dynamic transition observed for the mean square displacements of the protein atoms and the heme iron of myoglobin. Received: 10 July 1996 / Accepted: 10 April 1997     

The heme redox center of chloroplast cytochrome f is linked to a buried five-water chain.

The crystal structure of the 252-residue lumen-side domain of reduced cytochrome f, a subunit of the proton-pumping integral cytochrome b6f complex of oxygenic photosynthetic membranes, was determined to a resolution of 1.96 A from crystals cooled to -35 degrees. The model was refined to an R-factor of 15.8% with a 0.013-A RMS deviation of bond lengths from ideality. Compared to the structure of cytochrome f at 20 degrees, the structure at -35 degrees has a small change in relative orientation of the two folding domains and significantly lower isotropic temperature factors for protein atoms. The structure revealed an L-shaped array of five buried water molecules that extend in two directions from the N delta 1 of the heme ligand His 25. The longer branch extends 11 A within the large domain, toward Lys 66 in the prominent basic patch at the top of the large domain, which has been implicated in the interaction with the electron acceptor, plastocyanin. The water sites are highly occupied, and their temperature factors are comparable to those of protein atoms. Virtually all residues that form hydrogen bonds with the water chain are invariant among 13 known cytochrome f sequences. The water chain has many features that optimize it as a proton wire, including insulation from the protein medium. It is suggested that this chain may function as the lumen-side exit port for proton translocation by the cytochrome b6f complex.     

Structural studies of mutants of the lysozyme of bacteriophage T4. The temperature-sensitive mutant protein Thr157----Ile

To understand the roles of individual amino acids in the folding and stability of globular proteins, a systematic structural analysis of mutants of the lysozyme of bacteriophage T4 has been undertaken. The isolation, characterization, crystallographic refinement and structural analysis of a temperature-sensitive lysozyme in which threonine 157 is replaced by isoleucine is reported here. This mutation reduces the temperature of the midpoint of the reversible thermal denaturation transition by 11 deg.C at pH 2.0. Electron density maps showing differences between the wild-type and mutant X-ray crystal structures have obvious features corresponding to the substitution of threonine 157 by isoleucine. There is little difference electron density in the remainder of the molecule, indicating that the structural changes are localized to the site of the mutation. High-resolution crystallographic refinement of the mutant lysozyme structure confirms that it is very similar to wild-type lysozyme. The largest conformational differences are in the gamma-carbon of residue 157 and in the side-chain of Asp159, which shift 1.0 A and 1.1 A, respectively. In the wild-type enzyme, the gamma-hydroxyl group of Thr157 participates in a network of hydrogen bonds. Substitution of Thr157 with an isoleucine disrupts this set of hydrogen bonds. A water molecule bound in the vicinity of Thr155 partially restores the hydrogen bond network in the mutant structure, but the buried main-chain amide of Asp159 is not near a hydrogen bond acceptor. This unsatisfied hydrogen-bonding potential is the most obvious reason for the reduction in stability of the temperature-sensitive mutant protein.     

Crystallographic studies of protein denaturation and renaturation. 1. Effects of denaturants on volume and X-ray pattern of cross-linked triclinic lysozyme crystals

Triclinic crystals of hen egg-white lysozyme cross-linked with glutaraldehyde have been treated with various denaturants and found to be susceptible to x-ray structure analysis even after major conformational changes in the protein. Cross-linked crystals were isomorphous with the native form, and electron density difference maps indicated the locations of intermolecular corss-links, but showed no appreciable differences in the protein conformation. Soaking of the cross-linked crystals in danaturant solutions of increasing concentrations caused corresponding increases in crystal volume and decreases in minimum observable x-ray spacings. These changes proved partly reversible on diluting the solutions, and measurements of crystal volume and minimums x-ray spacing were used to follow denaturation and renaturation as a function of concentration for several denaturants. Some of these, including bromoethanol and sodium dodecyl sulfate, had little effect on the crystals below critical concentrations at which there was a sharp volume increase and loss of x-ray pattern, which could, however, be regenerated to about 3.2-A resolution. Others, including KCNS and urea, caused more gradual changes, but with a smaller degree of recovery. It is suggested that at least two different denaturation mechanisms are involved with detergent-like reagents disrupting the hydrophobic interactions joining the two wings of the lysozyme molecule and hydrophilic denaturants interacting primarily with polar groups on the molecular surface.     

Electrophoresis in protein crystal: nonequilibrium molecular dynamics simulations

Electrophoresis of a mixture of NaCl and CaCl₂ in a lysozyme crystal is investigated using nonequilibrium molecular dynamics (MD) simulations. Upon exposure to an electric field, the stability of lysozyme is found to decrease slightly. This finding is demonstrated by increases in the root mean-square deviations of the heavy atoms of lysozyme, in the solvent-accessible surface area of hydrophobic residues, and in the number of hydrogen bonds between lysozyme and water. The solvent-accessible surface area of hydrophilic residues changes marginally, and the number of hydrogen bonds between lysozyme molecules decreases. Water molecules tend to align preferentially parallel to the electric field, and the dipole moment along the pore axis increases linearly with increasing field strength. Two pronounced layered structures are observed for Na⁺ and Ca²⁺ in the vicinity of protein surface, but only one enriched layer is observed for Cl⁻. The number distributions of all ions are nearly independent of the electric field. The water coordination numbers of all ions are smaller in the crystal than in aqueous bulk solution; however, the reverse is found for the Cl⁻ coordination numbers of cations. Both the water and the Cl⁻ coordination numbers are insensitive to the electric field. Ion diffusivities in the crystal are ∼2 orders of magnitude smaller than those in aqueous bulk solution. The drift velocities of ions increase proportionally to the electric field, particularly at high strengths, and depend on ionic charge and coordination with oppositely charged ions. Electrical current exhibits a linear relationship with the field strength. The zero-field electrical conductivity is estimated to be 0.56 S/m, which is very close to 0.61 S/m as predicted by the Nernst-Einstein equation.