Determination of a high-quality nuclear magnetic resonance solution structure of the bovine pancreatic trypsin inhibitor and comparison with three crystal structures.

A high-quality three-dimensional structure of the bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution was determined by 1H nuclear magnetic resonance (n.m.r.) spectroscopy and compared to the three available high-resolution X-ray crystal structures. A newly collected input of 642 distance constraints derived from nuclear Overhauser effects and 115 dihedral angle constraints was used for the structure calculations with the program DIANA, followed by restrained energy minimization with the program AMBER. The BPTI solution structure is represented by a group of 20 conformers with an average root-mean-square deviation (RMSD) relative to the mean solution structure of 0.43 A for backbone atoms and 0.92 A for all heavy atoms of residues 2 to 56. The pairwise RMSD values of the three crystal structures relative to the mean solution structure are 0.76 to 0.85 A for the backbone atoms and 1.24 to 1.33 A for all heavy atoms of residues 2 to 56. Small local differences in backbone atom positions between the solution structure and the X-ray structures near residues 9, 25 to 27, 46 to 48 and 52 to 58, and conformational differences for individual amino acid side-chains were analyzed for possible correlations with intermolecular protein-protein contacts in the crystal lattices, using the pairwise RMSD values among the three crystal structures as a reference.     

Packing interface energetics in different crystal forms of the λ Cro dimer

Variation among crystal structures of the λ Cro dimer highlights conformational flexibility. The structures range from a wild type closed to a mutant fully open conformation, but it is unclear if each represents a stable solution state or if one may be the result of crystal packing. Here we use molecular dynamics (MD) simulation to investigate the energetics of crystal packing interfaces and the influence of site‐directed mutagenesis on them in order to examine the effect of crystal packing on wild type and mutant Cro dimer conformation. Replica exchange MD of mutant Cro in solution shows that the observed conformational differences between the wild type and mutant protein are not the direct consequence of mutation. Instead, simulation of Cro in different crystal environments reveals that mutation affects the stability of crystal forms. Molecular Mechanics Poisson‐Boltzmann Surface Area binding energy calculations reveal the detailed energetics of packing interfaces. Packing interfaces can have diverse properties in strength, energetic components, and some are stronger than the biological dimer interface. Further analysis shows that mutation can strengthen packing interfaces by as much as ～5 kcal/mol in either crystal environment. Thus, in the case of Cro, mutation provides an additional energetic contribution during crystal formation that may stabilize a fully open higher energy state. Moreover, the effect of mutation in the lattice can extend to packing interfaces not involving mutation sites. Our results provide insight into possible models for the effect of crystallization on Cro conformational dynamics and emphasize careful consideration of protein crystal structures. Proteins 2014; 82:1128–1141. © 2013 Wiley Periodicals, Inc.     

Sampling of the native conformational ensemble of myoglobin via structures in different crystalline environments

Proteins sample multiple conformational substates in their native environment, but the process of crystallization selects the conformers that allow for close packing. The population of conformers can be shifted by varying the environment through a range of crystallization conditions, often resulting in different space groups and changes in the packing arrangements. Three high resolution structures of myoglobin (Mb) in different crystal space groups are presented, including one in a new space group P6(1)22 and two structures in space groups P2(1)2(1)2(1) and P6. We compare coordinates and anisotropic displacement parameters (ADPs) from these three structures plus an existing structure in space group P2(1). While the overall changes are small, there is substantial variation in several external regions with varying patterns of crystal contacts across the space group packing arrangements. The structural ensemble containing four different crystal forms displays greater conformational variance (Calpha rmsd of 0.54-0.79 A) in comparison to a collection of four Mb structures with different ligands and mutations in the same crystal form (Calpha rmsd values of 0.28-0.37 A). The high resolution of the data enables comparison of both the magnitudes and directions of ADPs, which are found to be suppressed by crystal contacts. A composite dynamic profile of Mb structural variation from the four structures was compared with an independent structural ensemble developed from NMR refinement. Despite the limitations and biases of each method, the ADPs of the crystallographic ensemble closely match the positional variance from the solution NMR ensemble with linear correlation of 0.8. This suggests that crystal packing selects conformers representative of the solution ensemble, and several different crystal forms give a more complete view of the plasticity of a protein structure.     

'Crystal lattice engineering,' an approach to engineer protein crystal contacts by creating intermolecular symmetry: crystallization and structure determination of a mutant human RNase 1 with a hydrophobic interface of leucines

A protein crystal lattice consists of surface contact regions, where the interactions of specific groups play a key role in stabilizing the regular arrangement of the protein molecules. In an attempt to control protein incorporation in a crystal lattice, a leucine zipper-like hydrophobic interface (comprising four leucine residues) was introduced into a helical region (helix 2) of the human pancreatic ribonuclease 1 (RNase 1) that was predicted to form a suitable crystallization interface. Although crystallization of wild-type RNase 1 has not yet been reported, the RNase 1 mutant having four leucines (4L-RNase 1) was successfully crystallized under several different conditions. The crystal structures were subsequently determined by X-ray crystallography by molecular replacement using the structure of bovine RNase A. The overall structure of 4L-RNase 1 is quite similar to that of the bovine RNase A, and the introduced leucine residues formed the designed crystal interface. To characterize the role of the introduced leucine residues in crystallization of RNase 1 further, the number of leucines was reduced to three or two (3L- and 2L-RNase 1, respectively). Both mutants crystallized and a similar hydrophobic interface as in 4L-RNase 1 was observed. A related approach to engineer crystal contacts at helix 3 of RNase 1 (N4L-RNase 1) was also evaluated. N4L-RNase 1 also successfully crystallized and formed the expected hydrophobic packing interface. These results suggest that appropriate introduction of a leucine zipper-like hydrophobic interface can promote intermolecular symmetry for more efficient protein crystallization in crystal lattice engineering efforts.     

Heterogeneity and inaccuracy in protein structures solved by X-ray crystallography

Proteins are dynamic molecules, exhibiting structural heterogeneity in the form of anisotropic motion and discrete conformational substates, often of functional importance. In protein structure determination by X-ray crystallography, the observed diffraction pattern results from the scattering of X-rays by an ensemble of heterogeneous molecules, ordered and oriented by packing in a crystal lattice. The majority of proteins diffract to resolutions where heterogeneity is difficult to identify and model, and are therefore approximated by a single, average conformation with isotropic variance. Here we show that disregarding structural heterogeneity introduces degeneracy into the structure determination process, as many single, isotropic models exist that explain the diffraction data equally well. The large differences among these models imply that the accuracy of crystallographic structures has been widely overestimated. Further, it suggests that analyses that depend on small differences in the relative positions of atoms may be flawed.     

Ligand-induced conformational changes in the crystal structures of Pneumocystis carinii dihydrofolate reductase complexes with folate and NADP+

Structural data from two independent crystal forms (P212121 and P21) of the folate (FA) binary complex and from the ternary complex with the oxidized coenzyme, NADP+, and recombinant Pneumocystis carinii dihydrofolate reductase (pcDHFR) refined to an average of 2.15 A resolution, show the first evidence of ligand-induced conformational changes in the structure of pcDHFR. These data are also compared with the crystal structure of the ternary complex of methotrexate (MTX) with NADPH and pcDHFR in the monoclinic lattice with data to 2.5 A resolution. Comparison of the data for the FA binary complex of pcDHFR with those for the ternary structures reveals significant differences, with a >7 A movement of the loop region near residue 23 that results in a new "flap-open" position for the binary complex, and a "closed" position in the ternary complexes, similar to that reported for Escherichia coli (ec) DHFR complexes. In the orthorhombic lattice for the binary FA pcDHFR complex, there is also an unwinding of a short helical region near residue 47 that places hydrophobic residues Phe-46 and Phe-49 toward the outer surface, a conformation that is stabilized by intermolecular packing contacts. The pyrophosphate moiety of NADP+ in the ternary folate pcDHFR complexes shows significant differences in conformation compared with that observed in the MTX-NADPH-pcDHFR ternary complex. Additionally, comparison of the conformations among these four pcDHFR structures reveals evidence for subdomain movement that correlates with cofactor binding states. The larger binding site access in the new "flap-open" loop 23 conformation of the binary FA complex is consistent with the rapid release of cofactor from the product complex during catalysis as well as the more rapid release of substrate product from the binary complex as a result of the weaker contacts of the closed loop 23 conformation, compared to ecDHFR.     

Conformational transitions of flanking purines in HIV-1 RNA dimerization initiation site kissing complexes studied by CHARMM explicit solvent molecular dynamics

Dimerization of HIV-1 genomic RNA is initiated by kissing loop interactions at the Dimerization Initiation Site (DIS). Dynamics of purines that flank the 5' ends of the loop-loop helix in HIV-1 DIS kissing complex were explored using explicit solvent molecular dynamics (MD) simulations with the CHARMM force field. Multiple MD simulations (200 ns in total) of X-ray structures for HIV-1 DIS Subtypes A, B, and F revealed conformational variability of flanking purines. In particular, the flanking purines, which in the starting X-ray structures are bulged-out and stack in pairs, formed a consecutive stack of four bulged-out adenines at the beginning of several simulations. This conformation is seen in the crystal structure of DIS Subtype F with no interference from crystal packing, and was frequently reported in our preceding MD studies performed with the AMBER force field. However, as CHARMM simulations progressed, the four continuously stacked adenines showed conformational transitions from the bulged-out into the bulged-in geometries. Although such an arrangement has not been seen in any X-ray structure, it has been suggested by a recent NMR investigation. In CHARMM simulations, in the longer time scale, the flanking purines display the tendency to move to bulged-in conformations. This is in contrast with the AMBER simulations, which indicate a modest prevalence for bulged-out flanking base positions in line with the X-ray data. The simulations also suggest that the intermolecular stacking between purines from the opposite hairpins can additionally stabilize the kissing complex.     

Packing forces in nine crystal forms of cutinase

During the characterization of mutants and covalently inhibited complexes of Fusarium solani cutinase, nine different crystal forms have been obtained so far. Protein mutants with a different surface charge distribution form new intermolecular salt bridges or long-range electrostatic interactions that are accompanied by a change in the crystal packing. The whole protein surface is involved in the packing contacts and the hydrophobicities of the protein surfaces in mutual contact turned out to be noncorrelated, which indicates that the packing interactions are nonspecific. In the case of the hydrophobic variants, the packing contacts showed some specificity, as the protein in the crystal tends to form either crystallographic or noncrystallographic dimers, which shield the hydrophobic surface from the solvent. The likelihood of surface atoms to be involved in a crystal contact is the same for both polar and nonpolar atoms. However, when taking areas in the 200–600 Ų range, instead of individual atoms, the either highly hydrophobic or highly polar surface regions were found to have an increased probability of establishing crystal lattice contacts. The protein surface surrounding the active-site crevice of cutinase constitutes a large hydrophobic area that is involved in packing contacts in all the various crystalline contexts. Proteins 31:320–333, 1998. © 1998 Wiley-Liss, Inc.     

A few low-frequency normal modes predominantly contribute to conformational responses of hen egg white lysozyme in the tetragonal crystal to variations of molecular packing controlled by environmental humidity

The structures of proteins in crystals are fixed by molecular interactions with neighboring molecules, except in non-contacting flexible regions. Thus, it is difficult to imagine what conformational changes occur in solution. However, if molecular interactions can be changed by manipulating molecular packing in crystals, it may be possible to visualize conformational responses of proteins at atomic resolution by diffraction experiments. For this purpose, it is suitable to control the molecular packing in protein crystals by changing the volume of solvent channels through variation of the environmental relative humidity. Here, we studied conformational responses of hen egg white lysozyme (HEWL) in the tetragonal crystal by X-ray diffraction experiments using a humidity-control apparatus, which provided air flow of 20-98%rh at 298 K. First, we monitored the lattice parameters and crystalline order during dehydration and rehydration of HEWL crystal between 61 and 94%rh at 300 K. Then two crystal structures at a resolution of 2.1 ? using diffraction data obtained at 84.2 and 71.9%rh were determined to discuss the conformational responses of HEWL against the external perturbation induced by changes in molecular packing. The structure at 71.9%rh displayed a closure movement that was likely induced by the molecular contacts formed during dehydration and could be approximated by ten low-frequency normal modes for the crystal structure obtained at 84.2%rh. In addition, we observed reorganization of hydration structures at the molecular interfaces between symmetry neighbors. These findings suggest that humidity-controlled X-ray crystallography is an effective tool to investigate the responses of inherent intramolecular motions of proteins to external perturbations.     

Structure of crambin in solution, crystal and in the trajectories of molecular dynamics simulations

The mechanisms of the three-dimensional crambin structure alterations in the crystalline environments and in the trajectories of the molecular dynamics simulations in the vacuum and crystal surroundings have been analyzed. In the crystalline state and in the solution the partial regrouping of remote intramolecular packing contacts, involved in the formation and stabilization of the tertiary structure of the crambin molecule, occurs in NMR structures. In the crystalline state it is initiated by the formation of the intermolecular contacts, the conformational influence of its appearance is distributed over the structure. The changes of the conformations and positions of the residues of the loop segments, where the intermolecular contacts of the crystal surroundings are preferably concentrated, are most observable. Under the influence of these contacts the principal change of the regular secondary structure of crambin is taking place: extension of the two-strand β structure to the three-strand structure with the participation of the single last residue N46 of the C-terminal loop. In comparison with the C-terminal loop the more profound changes are observed in the conformation and the atomic positions of the backbone atoms and in the solvent accessibility of the residues of the interhelical loop. In the solution of the ensemble of the 8 NMR structures relative accessibility to the solvent differs more noticeably also in the region of the loop segments and rather markedly in the interhelical loop. In the crambin cryogenic crystal structures the positions of the atoms of the backbone and/or side chain of 14–18 of 46 residues are discretely disordered. The disorganizations of at least 8 of 14 residues occur directly in the regions of the intermolecular contacts and another 5 residues are disordered indirectly through the intramolecular contacts with the residues of the intermolecular contacts. Upon the molecular dynamics simulation in the vacuum surrounding as in the solution of the crystalline structure of crambin the essential changes of the backbone conformation are caused by the intermolecular contacts absence, but partly masked by the structure changes owing to the nonpolar H atoms absence on the simulated structure. The intermolecular contact absence is partly manifested upon the molecular dynamics simulation of the crambin crystal with one protein molecule. Compared to the crystal structure the lengths of the interpeptide hydrogen bonds and other interresidue contacts in an average solution NMR structure are somewhat shorter and accordingly the energy of the interpeptide hydrogen bonds is better. This length shortening can occur at the stage of the refinement of the NMR structures of the crambin and other proteins by its energy minimizations in the vacuum surroundings and not exist in the solution protein structures.     

Crystal packing analysis of Rhodopsin crystals

Oligomerization has been proposed as one of several mechanisms to regulate the activity of G protein-coupled receptors (GPCRs), but little is known about the structure of GPCR oligomers. Crystallographic analyses of two new crystal forms of rhodopsin reveal an interaction surface which may be involved in the formation of functional dimers or oligomers. New crystallization conditions lead to the formation of two crystal forms with similar rhodopsin-rhodopsin interactions, but changes in the crystal lattice are induced by the addition of different surfactant additives. However, the intermolecular interactions between rhodopsin molecules in these crystal structures may reflect the contacts necessary for the maintenance of dimers or oligomers in rod outer segment membranes. Similar contacts may assist in the formation of dimers or oligomers in other GPCRs as well. These new dimers are compared with other models proposed by crystallography or EM and AFM studies. The inter-monomer surface contacts are different for each model, but several of these models coincide in implicating helix I, II, and H-8 as contributors to the main contact surface stabilizing the dimers.     

Reprint of "Crystal packing analysis of Rhodopsin crystals" [J. Struct. Biol. 158 (2007) 455-462

Oligomerization has been proposed as one of several mechanisms to regulate the activity of G protein-coupled receptors (GPCRs), but little is known about the structure of GPCR oligomers. Crystallographic analyses of two new crystal forms of rhodopsin reveal an interaction surface which may be involved in the formation of functional dimers or oligomers. New crystallization conditions lead to the formation of two crystal forms with similar rhodopsin-rhodopsin interactions, but changes in the crystal lattice are induced by the addition of different surfactant additives. However, the intermolecular interactions between rhodopsin molecules in these crystal structures may reflect the contacts necessary for the maintenance of dimers or oligomers in rod outer segment membranes. Similar contacts may assist in the formation of dimers or oligomers in other GPCRs as well. These new dimers are compared with other models proposed by crystallography or EM and AFM studies. The inter-monomer surface contacts are different for each model, but several of these models coincide in implicating helix I, II, and H-8 as contributors to the main contact surface stabilizing the dimers.     

Crystal structures of five bovine chymotrypsin complexes with P1 BPTI variants

The bovine chymotrypsin-bovine pancreatic trypsin inhibitor (BPTI) interaction belongs to extensively studied models of protein-protein recognition. The accommodation of the inhibitor P1 residue in the S1 binding site of the enzyme forms the hot spot of this interaction. Mutations introduced at the P1 position of BPTI result in a more than five orders of magnitude difference of the association constant values with the protease. To elucidate the structural aspects of the discrimination between different P1 residues, crystal structures of five bovine chymotrypsin-P1 BPTI variant complexes have been determined at pH 7.8 to a resolution below 2 A. The set includes polar (Thr), ionizable (Glu, His), medium-sized aliphatic (Met) and large aromatic (Trp) P1 residues and complements our earlier studies of the interaction of different P1 side-chains with the S1 pocket of chymotrypsin. The structures have been compared to the complexes of proteases with similar and dissimilar P1 preferences, including Streptomyces griseus proteases B and E, human neutrophil elastase, crab collagenase, bovine trypsin and human thrombin. The S1 sites of these enzymes share a common general shape of significant rigidity. Large and branched P1 residues adapt in their complexes similar conformations regardless of the polarity and size differences between their S1 pockets. Conversely, long and flexible residues such as P1 Met are present in the disordered form and display a conformational diversity despite similar inhibitory properties with respect to most enzymes studied. Thus, the S1 specificity profiles of the serine proteases appear to result from the precise complementarity of the P1-S1 interface and minor conformational adjustments occurring upon the inhibitor binding.     

Calculations on crystal packing of a flexible molecule, Leu-enkephalin

Crystal packing calculations have been carried out on a substantial number of conformations of Leu-enkephalin; namely, those obtained both from crystal structures and from energy minimizations on isolated molecules, and with and without waters of crystallization. The known crystal structures represent the most energetically stable packings found. The conformations of the enkephalin molecules in the crystal are not the most stable for an isolated molecule; i.e. intermolecular interactions force the isolated molecule to change conformation in order to achieve a small packing volume and an optimal packing energy in the crystal. It is found that the packing energy of an enkephalin molecule is a reasonably smooth function of its molecular volume in the unit cell, if structures with intermolecular hydrogen bonding are excluded, and is substantially independent of other details of the molecular conformation or of the crystal packing. Hydrogen bonding provides additional stabilization of the crystal structure, and would likely permit crystallization of the system if it is sufficiently dense. Solvent molecules further stabilize the structure when they can also provide intermolecular hydrogen bonds.     

Structure and internal dynamics of the bovine pancreatic trypsin inhibitor in aqueous solution from long-time molecular dynamics simulations

Structural and dynamic properties of bovine pancreatic trypsin inhibitor (BPTI) in aqueous solution are investigated using two molecular dynamics (MD) simulations: one of 1.4 ns length and one of 0.8 ns length in which atom-atom distance bounds derived from NMR spectroscopy are included in the potential energy function to make the trajectory satisfy these experimental data more closely. The simulated properties of BPTI are compared with crystal and solution structures of BPTI, and found to be in agreement with the available experimental data. The best agreement with experiment was obtained when atom-atom distance restraints were applied in a time-averaged manner in the simulation. The polypeptide segments found to be most flexible in the MD simulations coincide closely with those showing differences between the crystal and solution structures of BPTI. © 1995 Wiley-Liss, Inc.     

A large data set comparison of protein structures determined by crystallography and NMR: statistical test for structural differences and the effect of crystal packing

The existence of a large number of proteins for which both nuclear magnetic resonance (NMR) and X-ray crystallographic coordinates have been deposited into the Protein Data Bank (PDB) makes the statistical comparison of the corresponding crystal and NMR structural models over a large data set possible, and facilitates the study of the effect of the crystal environment and other factors on structure. We present an approach for detecting statistically significant structural differences between crystal and NMR structural models which is based on structural superposition and the analysis of the distributions of atomic positions relative to a mean structure. We apply this to a set of 148 protein structure pairs (crystal vs NMR), and analyze the results in terms of methodological and physical sources of structural difference. For every one of the 148 structure pairs, the backbone root-mean-square distance (RMSD) over core atoms of the crystal structure to the mean NMR structure is larger than the average RMSD of the members of the NMR ensemble to the mean, with 76% of the structure pairs having an RMSD of the crystal structure to the mean more than a factor of two larger than the average RMSD of the NMR ensemble. On average, the backbone RMSD over core atoms of crystal structure to the mean NMR is approximately 1 A. If non-core atoms are included, this increases to 1.4 A due to the presence of variability in loops and similar regions of the protein. The observed structural differences are only weakly correlated with the age and quality of the structural model and differences in conditions under which the models were determined. We examine steric clashes when a putative crystalline lattice is constructed using a representative NMR structure, and find that repulsive crystal packing plays a minor role in the observed differences between crystal and NMR structures. The observed structural differences likely have a combination of physical and methodological causes. Stabilizing attractive interactions arising from intermolecular crystal contacts which shift the equilibrium of the crystal structure relative to the NMR structure is a likely physical source which can account for some of the observed differences. Methodological sources of apparent structural difference include insufficient sampling or other issues which could give rise to errors in the estimates of the precision and/or accuracy.     

Exploring the dynamic information content of a protein NMR structure: Comparison of a molecular dynamics simulation with the NMR and X‐ray structures of Escherichia coli ribonuclease HI

The multiconformer nature of solution nuclear magnetic resonance (NMR) structures of proteins results from the effects of intramolecular dynamics, spin diffusion and an uneven distribution of structural restraints throughout the molecule. A delineation of the former from the latter two contributions is attempted in this work for an ensemble of 15 NMR structures of the protein Escherichia coli ribonuclease HI (RNase HI). Exploration of the dynamic information content of the NMR ensemble is carried out through correlation with data from two crystal structures and a 1.7‐ns molecular dynamics (MD) trajectory of RNase HI in explicit solvent. Assessment of the consistency of the crystal and mean MD structures with nuclear Overhauser effect (NOE) data showed that the NMR ensemble is overall more compatible with the high‐resolution (1.48 Å) crystal structure than with either the lower‐resolution (2.05 Å) crystal structure or the MD simulation. Furthermore, the NMR ensemble is found to span more conformational space than the MD simulation for both the backbone and the sidechains of RNase HI. Nonetheless, the backbone conformational variability of both the NMR ensemble and the simulation is especially consistent with NMR relaxation measurements of two loop regions that are putative sites of substrate recognition. Plausible side‐chain dynamic information is extracted from the NMR ensemble on the basis of (i) rotamericity and syn‐pentane character of variable torsion angles, (ii) comparison of the magnitude of atomic mean‐square fluctuations (msf) with those deduced from crystallographic thermal factors, and (iii) comparison of torsion angle conformational behavior in the NMR ensemble and the simulation. Several heterogeneous torsion angles, while adopting non‐rotameric/syn‐pentane conformations in the NMR ensemble, exist in a unique conformation in the simulation and display low X‐ray thermal factors. These torsions are identified as sites whose variability is likely to be an artifact of the NMR structure determination procedure. A number of other torsions show a close correspondence between the conformations sampled in the NMR and MD ensembles, as well as significant correlations among crystallographic thermal factors and atomic msf calculated from the NMR ensemble and the simulation. These results indicate that a significant amount of dynamic information is contained in the NMR ensemble. The relevance of the present findings for the biological function of RNase HI, protein recognition studies, and previous investigations of the motional content of protein NMR structures are discussed. Proteins 1999;36:87–110. © 1999 Wiley‐Liss, Inc.     

Reversible lattice repacking illustrates the temperature dependence of macromolecular interactions

Flash-freezing, which has become routine in macromolecular X-ray crystallography, causes the crystal to contract substantially. In the case of Escherichia coli beta-galactosidase the changes are reversible and are shown to be due to lattice repacking. On cooling, the area of the protein surface involved in lattice contacts increases by 50 %. There are substantial alterations in intermolecular contacts, these changes being dominated by the long, polar side-chains. For entropic reasons such side-chains, as well as surface solvent molecules, tend to be somewhat disordered at room temperature but can form extensive hydrogen-bonded networks on cooling. Low-temperature density measurements suggest that, at least in some cases, the beneficial effect of cryosolvents may be due to a density increase on vitrification which reduces the volume of bulk solvent that needs to be expelled from the crystal. Analysis of beta-galactosidase and several other proteins suggests that both intramolecular and intermolecular contact interfaces can be perturbed by cryocooling but that the changes tend to be more dramatic in the latter case. The temperature-dependence of the intermolecular interactions suggests that caution may be necessary in interpreting protein-protein and protein-nucleic acid interactions based on low-temperature crystal structures.