Backbone dynamics of the 269-residue protease Savinase determined from 15N-NMR relaxation measurements

Backbone dynamics of Savinase, a subtilisin of 269 residues secreted by Bacillus lentus, have been studied using 15N relaxation measurements derived from proton-detected dimensional 1H-15N-NMR spectroscopy. 15N spin-lattice rate constants (R1), spin-spin relaxation-rate constants(R2), and 1H-15N nuclear Overhauser effects (NOE) were determined for 84% of the backbone amide 15N nuclei. The model-free formalism [Lipari, G. & Szabo, A. (1982) J. Am. Chem. Soc. 104, 4546-4559] was used to derive values for a generalized order parameter, S2, interpretable as a measure of the amplitude of motion on the picosecond-nanosecond timescale, for each N-H bond vector. Additional terms used to fit the data include an effective correlation time for internal motions (taue) and an exchange term (Rex) to account for exchange contributions to R2. The overall rotational correlation time (taum) is 9.59 +/- 0.02 ns; the average order parameter (S2) is 0.90 +/- 0.07, indicative of a rigid structure consistent with Savinase's high degree of secondary structure and compact tertiary fold. Residues S125-S128, located in the substrate-binding region, represent the longest stretch of protein which exhibits disorder on the picosecond-nanosecond timescale. These residues also exhibit significant exchange terms, possibly indicative of motion on the microsecond-millisecond timescale, which could also be influenced by the proximity of the phenyl ring of the substituted aryl boronic acid inhibitor used in this study. S103 and G219 in the substrate-binding region, represent the longest stretch of protein which exhibits disorder on the picosecond-nanosecond timescale. These residues also exhibit significant exchange terms, possibly indicative of motion on the microsecond-millisecond timescale, which could also be influenced by the proximity of the phenyl ring of the substituted aryl boronic acid inhibitor used in this study. S103 and G219 in the substrate-binding region also show flexibility on the picosecond-nanosecond timescale. There is also significant motion in the turn, G258-T260, of a small solvent-exposed loop region which may make the protein vulnerable autolysis at that point. Some residues in both calcium-binding sites and nearby also show mobility.