Domain flexibility in ligand-free and inhibitor-bound Escherichia coli adenylate kinase based on a mode-coupling analysis of 15N spin relaxation

Adenylate kinase from Escherichia coli (AKeco), consisting of a 23.6-kDa polypeptide chain folded into domains CORE, AMPbd, and LID catalyzes the reaction AMP + ATP <--> 2ADP. The domains AMPbd and LID execute large-amplitude movements during catalysis. Backbone dynamics of ligand-free and AP(5)A-inhibitor-bound AKeco is studied with slowly relaxing local structure (SRLS) (15)N relaxation, an approach particularly suited when the global (tau(m)) and the local (tau) motions are likely to be coupled. For AKeco tau(m) = 15.1 ns, whereas for AKeco*AP(5)A tau(m) = 11.6 ns. The CORE domain of AKeco features an average squared order parameter, , of 0.84 and correlation times tau(f) = 5-130 ps. Most of the AKeco*AP(5)A backbone features = 0.90 and tau(f) = 33-193 ps. These data are indicative of relative rigidity. Domains AMPbd and LID of AKeco, and loops beta(1)/alpha(1), alpha(2)/alpha(3), alpha(4)/beta(3), alpha(5)/beta(4), and beta(8)/alpha(7) of AKeco*AP(5)A, feature a novel type of protein flexibility consisting of nanosecond peptide plane reorientation about the C(i-1)(alpha)-C(i)(alpha) axis, with correlation time tau(perpendicular) = 5.6-11.3 ns. The other microdynamic parameters underlying this dynamic model include S(2) = 0.13-0.5, tau(parallel) on the ps time scale, and a diffusion tilt beta(MD) ranging from 12 to 21 degrees. For the ligand-free enzyme the tau(perpendicular) mode was shown to represent segmental domain motion, accompanied by conformational exchange contributions R(ex) < or = 4.4 s(-1). Loop alpha(4)/beta(3) and alpha(5)/beta(4) dynamics in AKeco*AP(5)A is related to the "energetic counter-balancing of substrate binding" effect apparently driving kinase catalysis. The other flexible AKeco*AP(5)A loops may relate to domain motion toward product release.     

Escherichia coli adenylate kinase dynamics: comparison of elastic network model modes with mode-coupling (15)N-NMR relaxation data

The dynamics of adenylate kinase of Escherichia coli (AKeco) and its complex with the inhibitor AP(5)A, are characterized by correlating the theoretical results obtained with the Gaussian Network Model (GNM) and the anisotropic network model (ANM) with the order parameters and correlation times obtained with Slowly Relaxing Local Structure (SRLS) analysis of (15)N-NMR relaxation data. The AMPbd and LID domains of AKeco execute in solution large amplitude motions associated with the catalytic reaction Mg(+2)*ATP + AMP --> Mg(+2)*ADP + ADP. Two sets of correlation times and order parameters were determined by NMR/SRLS for AKeco, attributed to slow (nanoseconds) motions with correlation time tau( perpendicular) and low order parameters, and fast (picoseconds) motions with correlation time tau( parallel) and high order parameters. The structural connotation of these patterns is examined herein by subjecting AKeco and AKeco*AP(5)A to GNM analysis, which yields the dynamic spectrum in terms of slow and fast modes. The low/high NMR order parameters correlate with the slow/fast modes of the backbone elucidated with GNM. Likewise, tau( parallel) and tau( perpendicular) are associated with fast and slow GNM modes, respectively. Catalysis-related domain motion of AMPbd and LID in AKeco, occurring per NMR with correlation time tau( perpendicular), is associated with the first and second collective slow (global) GNM modes. The ANM-predicted deformations of the unliganded enzyme conform to the functional reconfiguration induced by ligand-binding, indicating the structural disposition (or potential) of the enzyme to bind its substrates. It is shown that NMR/SRLS and GNM/ANM analyses can be advantageously synthesized to provide insights into the molecular mechanisms that control biological function.     

A novel view of domain flexibility in E. coli adenylate kinase based on structural mode-coupling (15)N NMR relaxation.

Adenylate kinase from Escherichia coli (AKeco), consisting of a single 23.6 kDa polypeptide chain folded into domains CORE, AMPbd and LID, catalyzes the reaction AMP+ATP-->2ADP. In the ligand-free enzyme the domains AMPbd and LID execute large-amplitude movements controlling substrate binding and product release during catalysis. Domain flexibility is investigated herein with the slowly relaxing local structure (SRLS) model for (15)N relaxation. SRLS accounts rigorously for coupling between the global and local N-H motions through a local ordering potential exerted by the protein structure at the N-H bond. The latter reorients with respect to its protein surroundings, which reorient on the slower time scale associated with the global protein tumbling. AKeco diffuses globally with correlation time tau(m)=15.1 ns, while locally two different dynamic cases prevail. The domain CORE features ordering about the equilibrium N-H bond orientation with order parameters, S(2), of 0.8-0.9 and local motional correlation times, tau, mainly between 5-130 ps. This represents a conventional rigid protein structure with rapid small-amplitude N-H fluctuations. The domains AMPbd and LID feature small parallel (Z(M)) ordering of S(2)=0.2-0.5 which can be reinterpreted as high perpendicular (Y(M)) ordering. M denotes the local ordering/local diffusion frame. Local motion about Z(M) is given by tau( parallel) approximately 5 ps and local motion of the effective Z(M) axis about Y(M) by tau( perpendicular)=6-11 ns. Z(M) is tilted at approximately 20 degrees from the N-H bond. The orientation of the Y(M) axis may be considered parallel to the C(alpha)(i-1)-C(alpha)(i) axis. The tau( perpendicular) mode reflects collective nanosecond peptide-plane motions, interpretable as domain motion. A powerful new model of protein flexibility/domain motion has been established. Conformational exchange (R(ex)) processes accompany the tau( perpendicular) mode. The SRLS analysis is compared with the conventional model-free analysis.     

The atomistic mechanism of conformational transition in adenylate kinase: a TEE-REX molecular dynamics study

We report on an atomistic molecular dynamics simulation of the complete conformational transition of Escherichia coli adenylate kinase (ADK) using the recently developed TEE-REX algorithm. Two phases characterize the transition pathway of ADK, which folds into the domains CORE and LID and the AMP binding domain AMPbd. Starting from the closed conformation, half-opening of the AMPbd precedes a partially correlated opening of the LID and AMPbd, defining the second phase. A highly stable salt bridge D118-K136 at the LID-CORE interface, contributing substantially to the total nonbonded LID-CORE interactions, was identified as a major factor that stabilizes the open conformation. Alternative transition pathways, such as AMPbd opening following LID opening, seem unlikely, as full transition events were not observed along this pathway. The simulation data indicate a high enthalpic penalty, possibly obstructing transitions along this route.     

Minimum free energy path of ligand-induced transition in adenylate kinase

Large-scale conformational changes in proteins involve barrier-crossing transitions on the complex free energy surfaces of high-dimensional space. Such rare events cannot be efficiently captured by conventional molecular dynamics simulations. Here we show that, by combining the on-the-fly string method and the multi-state Bennett acceptance ratio (MBAR) method, the free energy profile of a conformational transition pathway in Escherichia coli adenylate kinase can be characterized in a high-dimensional space. The minimum free energy paths of the conformational transitions in adenylate kinase were explored by the on-the-fly string method in 20-dimensional space spanned by the 20 largest-amplitude principal modes, and the free energy and various kinds of average physical quantities along the pathways were successfully evaluated by the MBAR method. The influence of ligand binding on the pathways was characterized in terms of rigid-body motions of the lid-shaped ATP-binding domain (LID) and the AMP-binding (AMPbd) domains. It was found that the LID domain was able to partially close without the ligand, while the closure of the AMPbd domain required the ligand binding. The transition state ensemble of the ligand bound form was identified as those structures characterized by highly specific binding of the ligand to the AMPbd domain, and was validated by unrestrained MD simulations. It was also found that complete closure of the LID domain required the dehydration of solvents around the P-loop. These findings suggest that the interplay of the two different types of domain motion is an essential feature in the conformational transition of the enzyme.     

How long does DNA keep the memory of its conformation? A time-dependent canonical correlation analysis of molecular dynamics simulation

The time dependence of the correlation between motions of different parts of DNA is analyzed from a 200 ps molecular dynamics simulation of the double-stranded self-complementary d(CTGATCAG) in the B form. Each nucleotide is decomposed into three subunits corresponding to the furanose ring (SU), the base (BA), and the backbone (SK). The motion of each subunit is considered as the superimposition of rigid body translation, rigid body rotation, and internal deformation. Canonical time-dependent correlation functions calculated with coordinates describing the different components of the subunits motion are defined and computed. This allows us to probe how long a particular type of motion of one subunit influences the other types of motions of other subunits (cross correlation functions) or how long a particular subunit keeps the memory of its own conformation or location (autocorrelation functions). From auto-correlation analysis it is found that deformation decorrelates within a few tenths of picoseconds, rotational correlation times are on the order of 8 ps, while translational motions are long-time correlated. The deformation of a subunit is not correlated to the deformation of another one (at the 200 ps time scale of our simulation), but influences slightly their translation and orientation as time increases. © 1996 John Wiley & Sons, Inc.     

Backbone dynamics of the oligomerization domain of p53 determined from 15N NMR relaxation measurements.

The backbone dynamics of the tetrameric p53 oligomerization domain (residues 319-360) have been investigated by two-dimensional inverse detected heteronuclear 1H-15N NMR spectroscopy at 500 and 600 MHz. 15N T1, T2, and heteronuclear NOEs were measured for 39 of 40 non-proline backbone NH vectors at both field strengths. The overall correlation time for the tetramer, calculated from the T1/T2 ratios, was found to be 14.8 ns at 35 degrees C. The correlation times and amplitudes of the internal motions were extracted from the relaxation data using the model-free formalism (Lipari G, Szabo A, 1982, J Am Chem Soc 104:4546-4559). The internal dynamics of the structural core of the p53 oligomerization domain are uniform and fairly rigid, with residues 327-354 exhibiting an average generalized order parameter (S2) of 0.88 +/- 0.08. The N- and C-termini exhibit substantial mobility and are unstructured in the solution structure of p53. Residues located at the N- and C-termini, in the beta-sheet, in the turn between the alpha-helix and beta-sheet, and at the C-terminal end of the alpha-helix display two distinct internal motions that are faster than the overall correlation time. Fast internal motions (< or = 20 ps) are within the extreme narrowing limit and are of uniform amplitude. The slower motions (0.6-2.2 ns) are outside the extreme narrowing limit and vary in amplitude.(ABSTRACT TRUNCATED AT 250 WORDS)     

Backbone dynamics of apocytochrome b5 in its native, partially folded state

The backbone dynamics in the native state of apocytochrome b5 were studied using 15N nuclear magnetic spin relaxation measurements. The field (11.7 and 14.1 T) and temperature (10-25 degrees C) dependence of the relaxation parameters (R1, R2, and R1rho) and the 1H-15N NOE established that the protein undergoes multiple time scale internal motions related to the secondary structure. The relaxation data were analyzed with the reduced spectral density mapping approach and within the extended model-free framework. The apoprotein was confirmed to contain a disordered heme-binding loop of approximately 30 residues with dynamics on the sub-nanosecond time scale (0.6 < S2 < 0.7, 100 ps < taue < 500 ps). This loop is attached to a structured hydrophobic core, rigid on the picosecond time scale (S2 > 0.75, taue < 50 ps). The inability to fit the data for several residues with the model-free protocol revealed the presence of correlated motion. An exchange contribution was detected in the transverse relaxation rate (R2) of all residues. The differential temperature response of R2 along the backbone supported slower exchange rates for residues in the loop (tauex > 300 micros) than for the folded polypeptide chain (tauex < 150 micros). The distribution of the reduced spectral densities at the 1H and 15N frequencies followed the dynamic trend and predicted the slowing of the internal motions at 10 degrees C. Comparison of the dynamics with those of the holoprotein [Dangi, B., Sarma, S., Yan, C., Banville, D. L., and Guiles, R. D. (1998) Biochemistry 37, 8289-8302] demonstrated that binding of the heme alters the time scale of motions both in the heme-binding loop and in the structured hydrophobic core.     

Hydration-coupled dynamics in proteins studied by neutron scattering and NMR: the case of the typical EF-hand calcium-binding parvalbumin

The influence of hydration on the internal dynamics of a typical EF-hand calciprotein, parvalbumin, was investigated by incoherent quasi-elastic neutron scattering (IQNS) and solid-state 13C-NMR spectroscopy using the powdered protein at different hydration levels. Both approaches establish an increase in protein dynamics upon progressive hydration above a threshold that only corresponds to partial coverage of the protein surface by the water molecules. Selective motions are apparent by NMR in the 10-ns time scale at the level of the polar lysyl side chains (externally located), as well as of more internally located side chains (from Ala and Ile), whereas IQNS monitors diffusive motions of hydrogen atoms in the protein at time scales up to 20 ps. Hydration-induced dynamics at the level of the abundant lysyl residues mainly involve the ammonium extremity of the side chain, as shown by NMR. The combined results suggest that peripheral water-protein interactions influence the protein dynamics in a global manner. There is a progressive induction of mobility at increasing hydration from the periphery toward the protein interior. This study gives a microscopic view of the structural and dynamic events following the hydration of a globular protein.     

Normal mode analysis as a method to derive protein dynamics information from the Protein Data Bank

Normal mode analysis (NMA) can facilitate quick and systematic investigation of protein dynamics using data from the Protein Data Bank (PDB). We developed an elastic network model-based NMA program using dihedral angles as independent variables. Compared to the NMA programs that use Cartesian coordinates as independent variables, key attributes of the proposed program are as follows: (1) chain connectivity related to the folding pattern of a polypeptide chain is naturally embedded in the model; (2) the full-atom system is acceptable, and owing to a considerably smaller number of independent variables, the PDB data can be used without further manipulation; (3) the number of variables can be easily reduced by some of the rotatable dihedral angles; (4) the PDB data for any molecule besides proteins can be considered without coarse-graining; and (5) individual motions of constituent subunits and ligand molecules can be easily decomposed into external and internal motions to examine their mutual and intrinsic motions. Its performance is illustrated with an example of a DNA-binding allosteric protein, a catabolite activator protein. In particular, the focus is on the conformational change upon cAMP and DNA binding, and on the communication between their binding sites remotely located from each other. In this illustration, NMA creates a vivid picture of the protein dynamics at various levels of the structures, i.e., atoms, residues, secondary structures, domains, subunits, and the complete system, including DNA and cAMP. Comparative studies of the specific protein in different states, e.g., apo- and holo-conformations, and free and complexed configurations, provide useful information for studying structurally and functionally important aspects of the protein.     

Spectral densities of nitrogen nuclei in Escherichia coli ribonuclease HI obtained by 15N NMR relaxation and molecular dynamics

Summary Spectral densities of the ¹⁵N amide in Escherichia coli ribonuclease HI, obtained from NMR relaxation experiments, were compared with those calculated using a molecular dynamics (MD) simulation. All calculations and comparisons assumed that the auto-correlation function describing the internal motions of the molecule was independent of the auto-correlation function associated with overall rotational diffusion. Comparisons were limited to those residues for which the auto-correlation function of internal motions rapidly relaxed and reached a steady state within 205 ps. The results show the importance of frequency components as well as amplitudes of internal motions in order to obtain a meaningful comparison of MD simulations with NMR data.     

15N NMR relaxation studies of backbone dynamics in free and steroid-bound Delta 5-3-ketosteroid isomerase from Pseudomonas testosteroni

The backbone dynamics of Delta(5)-3-ketosteroid isomerase (KSI) from Pseudomonas testosteroni has been studied in free enzyme and its complex with a steroid ligand, 19-nortestosterone hemisuccinate (19-NTHS), by (15)N relaxation measurements. The relaxation data were analyzed using the model-free formalism to extract the model-free parameters (S(2), tau(e), and R(ex)) and the overall rotational correlation time (tau(m)). The rotational correlation times were 19.23 +/- 0.08 and 17.08 +/- 0.07 ns with the diffusion anisotropies (D( parallel)/D( perpendicular)) of 1.26 +/- 0.03 and 1.25 +/- 0.03 for the free and steroid-bound KSI, respectively. The binding of 19-NTHS to free KSI causes a slight increase in the order parameters (S(2)) for a number of residues, which are located mainly in helix A1 and strand B4. However, the majority of the residues exhibit reduced order parameters upon ligand binding. In particular, strands B3, B5, and B6, which have most of the residues involved in the dimer interaction, have the reduced order parameters in the steroid-bound KSI, indicating the increased high-frequency (pico- to nanosecond) motions in the intersubunit region of this homodimeric enzyme. Our results differ from those of previous studies on the backbone dynamics of monomeric proteins, in which high-frequency internal motions are typically restricted upon ligand binding.     

Multinuclear magnetic resonance studies of Escherichia coli adenylate kinase in free and bound forms. Resonance assignment, secondary structure and ligand binding.

The crystal structure of Escherichia coli adenylate kinase (AKe) revealed three main components: a CORE domain, composed of a five-stranded parallel beta-sheet surrounded by alpha-helices, and two peripheral domains involved in covering the ATP in the active site (LID) and binding of the AMP (NMPbind). We initiated a long-term NMR study aiming to characterize the solution structure, binding mechanism and internal dynamics of the various domains. Using single (15N) and double-labeled (13C and 15N) samples and double- and triple-resonance NMR experiments we assigned 97% of the 1H, 13C and 15N backbone resonances, and proton and 13Cbeta resonances for more than 40% of the side chains in the free protein. Analysis of a 15N-labeled enzyme in complex with the bi-substrate analogue [P1,P5-bis(5'-adenosine)-pentaphosphate] (Ap5A) resulted in the assignment of 90% of the backbone 1H and 15N resonances and 42% of the side chain resonances. Based on short-range NOEs and 1H and 13C secondary chemical shifts, we identified the elements of secondary structure and the topology of the beta-strands in the unliganded form. The alpha-helices and the beta-strands of the parallel beta-sheet in solution have the same limits (+/- 1 residue) as those observed in the crystal. The first helix (alpha1) appears to have a frayed N-terminal side. Significant differences relative to the crystal were noticed in the LID domain, which in solution exhibits four antiparallel beta-strands. The secondary structure of the nucleoside-bound form, as deduced from intramolecular NOEs and the 1Halpha chemical shifts, is similar to that of the free enzyme. The largest chemical shift differences allowed us to map the regions of protein-ligand contacts. 1H/2H exchange experiments performed on free and Ap5A-bound enzymes showed a general decrease of the structural flexibility in the complex which is accompanied by a local increased flexibility on the N-side of the parallel beta-sheet.     

Crystal structure of human adenylate kinase 4 (L171P) suggests the role of hinge region in protein domain motion

It is well known that motion of LID and NMP-binding (NMP_bind) domains in adenylate kinase (AK) is important in ligand binding and catalysis. However, the nature of such domain motions is poorly characterized. One of the critical hinge regions is hinge IV, which connects the CORE and LID domains. In addition, the hinge IV contains a strictly conserved residue, L171, in the AK family. To investigate the role of hinge IV, crystal structure of human adenylate kinase 4 (AK4) L171P mutant was determined. This mutation dramatically changes the orientation of the LID domain, which could be described as a novel twisted-and-closed conformation in contrast to the open and closed conformations in other AKs. This mutant provides a new example of domain motions in AK family.     

Solution conformation and dynamics of a fungal cell wall polysaccharide isolated from Microsporum gypseum

The conformational and dynamical features of a branched mannan isolated from a fungal cell wall have been analysed by homo and heteronuclear NMR methods, employing different magnetic fields. 1HNMR cross relaxation times have been obtained for this polysaccharide and have been interpreted qualitatively using different motional models. 13C NMR relaxation parameters (T1, T2, NOE) have also been measured and interpreted using different approximations based on the Lipari and Szabo model free approach. The analysis of the data indicate the existence of important flexibility for the different linkages of the polysaccharide. Motions in the range of 4–6 ns contribute to the relaxation of the macromolecule, although faster internal motions in the 500 ps and 100 ps timescales are also present. These time scales indicate that segmental motions as well as internal motions around the glycosidic linkages are the major sources of relaxation for this molecule at 318 K. Molecular dynamics simulations have also been performed. The obtained results also indicate that the polysaccharide possess a substantial amount of conformational freedom.     

The effects of ligand binding on the backbone dynamics of the kringle 1 domain of human plasminogen

The internal motions of the backbone nitrogen atoms of the kringle 1 domain of human plasminogen (K1(Pg)) were examined in the absence and presence of the ligand, epsilon-aminocaproic acid. These dynamic properties were determined from (15)N NMR relaxation data in terms of the extended model-free parameters. The model of isotropic reorientation was found sufficient to account for overall molecular tumbling for both apo and EACA-bound K1(Pg). The global rotational correlation time (tau(m)) for apo-K1(Pg) was 5.87(+/-0.01) ns, while the tau(m) for ligand-bound K1(Pg) was 5.20(+/-0.01) ns, suggesting that perhaps some small degree of aggregation occurred in the apo form of the kringle module. Complexation of K1(Pg) with ligand mainly reduced those internal motions that occurred on a 100 ps to 5 ns time-scale. The magnitude of the chemical exchange was also attenuated upon ligand binding. These data are consistent with studies employing other approaches that suggest that the binding pocket is preformed in K1(Pg).     

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

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

Dynamics of the three methionyl side chains of Streptomyces subtilisin inhibitor. Deuterium NMR studies in solution and in the solid state.

Streptomyces subtilisin inhibitor (SSI) contains three methionine residues in a subunit: two (at positions 73 and 70) in the crucial enzyme-recognition sites P1 and P4, respectively, and one (Met 103) in the hydrophobic core. The motions of the side chains of these three Met residues and the changes in mobility on binding with subtilisin were studied by deuterium NMR spectroscopy in solution and in crystalline and powder solids. For this purpose, the wild-type SSI was deuterium-labeled at the methyl groups of all three Met residues, and three artificial mutant proteins were labeled at only one specific Met methyl group each. In solution, for methionines 73 and 70, the effective correlation times were only 0.8-1.0 x 10(-10)s indicating that the two side chains on the surface fluctuate almost freely. On formation of a complex with subtilisin, however, these high mobilities were quenched, giving a correlation time of 1.1 x 10(-8)s for the side chains of methionines 70 and 73. The correlation time of Met 103, located in the hydrophobic core, was at least 1.0 x 10(-8)s in free SSI, showing that its side chain motion is highly restricted. The nature of the internal motions of the three Met side chains was examined in more detail by deuterium NMR spectroscopy of powder and crystalline samples. The spectral patterns of the powder samples depended critically on hydration: immediately after lyophilization, the side-chain motions of the three Met residues were nearly quenched. With gradual hydration to 0.20 gram of water per gram protein-water, the orientational fluctuation of the methyl axes of methionines 70 and 73 was selectively enhanced in both amplitude and frequency (to about 1 MHz) and, at nearly saturating hydration (0.60 gram of water per gram protein-water), became extremely high in amplitude and frequency (> 10 MHz). In contrast, the polycrystalline wild-type SSI spectrum showed fine structures, reflecting characteristic motions of the Met side chains. The polycrystalline spectrum could be reproduced reasonably well by the same motion models and parameters used to simulate the powder spectrum at the final level of hydration, suggesting that the side-chain motions are similar in the fully hydrated powder and in crystals. Spin-lattice relaxation measurements gave evidence that, even in crystals, the methyl axes of all three Met residues undergo rapid motions with correlation times between 10(-8) and 10(-10)s, comparable to the correlation times in solution. Finally, in the hydrated stoichiometric complex of SSI with subtilisin BPN' in the solid state, large-amplitude motions are absent, but the side chains of methionines 70 and/or 73 are likely to have small-amplitude motions.