(1)H, (13)C and (15)N NMR backbone assignments of the 269-residue serine protease PB92 from Bacillus alcalophilus

Altered flexibility in the substrate-binding site of related native and engineered high-alkaline Bacillus subtilisins

High-alkaline serine proteases have been successfully applied as protein degrading components of detergent formulations and are subject to extensive protein engineering efforts to improve their stability and performance. Dynamics has been suggested to play an important role in determining enzyme activity and specificity and it is therefore of interest to establish how local changes in internal mobility affect protein stability, specificity and performance. Here we present the dynamic properties of the 269 residue serine proteases subtilisin PB92 (Maxacal(TM)) and subtilisin BLS (Savinase(TM)), secreted by Bacillus lentus, and an engineered quadruple variant, DSAI, that has improved washing performance. T1, T2 and heteronuclear NOE measurements of the 15N nuclei indicate that for all three proteins the majority of the backbone is very rigid, with only a limited number of residues being involved in local mobility. Many of the residues that constitute the S1 and S4 pockets, determining substrate specificity, are flexible in solution. In contrast, the backbone amides of the residues that constitute the catalytic triad do not exhibit any motion. Subtilisins PB92, BLS and DSAI demonstrate similar but not identical NMR relaxation rates. A detailed analysis of local flexibility indicates that the motion of residues Thr143 and Ala194 becomes more restricted in subtilisin BLS and DSAI. Noteworthy, the loop regions involved in substrate binding become more structured in the engineered variant as compared with the two native proteases, suggesting a relation between altered dynamics and performance. Similar conclusions have been established by X-ray crystallograpic methods, as shown in the accompanying paper.

The solution structure of serine protease PB92 from Bacillus alcalophilus presents a rigid fold with a flexible substrate-binding site

BACKGROUND

Research on high-alkaline proteases, such as serine protease PB92, has been largely inspired by their industrial application as protein-degrading components of washing powders. Serine protease PB92 is a member of the subtilase family of enzymes, which has been extensively studied. These studies have included exhaustive protein engineering investigations and X-ray crystallography, in order to provide insight into the mechanism and specificity of enzyme catalysis. Distortions have been observed in the substrate-binding region of subtilisin crystal structures, due to crystal contacts. In addition, the structural variability in the substrate-binding region of subtilisins is often attributed to flexibility. It was hoped that the solution structure of this enzyme would provide further details about the conformation of this key region and give new insights into the functional properties of these enzymes.

RESULTS

The three-dimensional solution structure of the 269-residue (27 kDa) serine protease PB92 has been determined using distance and dihedral angle constraints derived from triple-resonance NMR data. The solution structure is represented by a family of 18 conformers which overlay onto the average structure with backbone and all-heavy-atom root mean square deviations (for the main body of the molecule) of 0.88 and 1.21 A, respectively. The family of structures contains a number of regions of relatively high conformational heterogeneity, including various segments that are involved in the formation of the substrate-binding site. The presence of flexibility within these segments has been established from NMR relaxation parameters and measurements of amide proton exchange rates.

CONCLUSIONS

The solution structure of the serine protease PB92 presents a well defined global fold which is rigid with the exception of a restricted number of sites. Among the limited number of residues involved in significant internal mobility are those of two pockets, termed S1 and S4, within the substrate-binding site. The presence of flexibility within the binding site supports the proposed induced fit mechanism of substrate binding.

Characterizing the use of perdeuteration in NMR studies of large proteins: 13C, 15N and 1H assignments of human carbonic anhydrase II

Perdeuteration of all non-exchangeable proton sites can significantly increase the size of proteins and protein complexes for which NMR resonance assignments and structural studies are possible. Backbone 1H, 15N, 13CO, 13C alpha and 13C beta chemical shifts and aliphatic side-chain 13C and 1H(N)/15N chemical shifts for human carbonic anhydrase II (HCA II), a 259 residue 29 kDa metalloenzyme, have been determined using a strategy based on 2D, 3D and 4D heteronuclear NMR experiments, and on perdeuterated 13C/15N-labeled protein. To date, HCA II is one of the largest monomeric proteins studied in detail by high-resolution NMR. Of the backbone resonances, 85% have been assigned using fully protonated 15N and 3C/15N-labeled protein in conjunction with established procedures based on now standard 2D and 3D NMR experiments. HCA II has been perdeuterated both to complete the backbone resonance assignment and to assign the aliphatic side-chain 13C and 1H(N)/15N resonances. The incorporation of 2H into HCA II dramatically decreases the rate of 13C and 1H(N)T2 relaxation. This, in turn, increases the sensitivity of several key 1H/13C/15N triple-resonance correlation experiments. Many otherwise marginal heteronuclear 3D and 4D correlation experiments, which are important to the assignment strategy detailed herein, can now be executed successfully on HCA II. Further analysis suggests that, from the perspective of sensitivity, perdeuteration should allow other proteins with rotational correlation times significantly longer than HCA II (tau c = 11.4 ns) to be studied successfully with these experiments. Two different protocols have been used to characterize the secondary structure of HCA II from backbone chemical-shift data. Secondary structural elements determined in this manner compare favorably with those elements determined from a consensus analysis of the HCA II crystal structure. Finally, having outlined a general strategy for assigning backbone and side-chain resonances in a perdeuterated large protein, we propose a strategy whereby this information can be used to glean more detailed structural information from the partially or fully protonated protein equivalent.

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.

Improved methods for structural studies of proteins using nuclear magnetic resonance spectroscopy

The past few years have seen the development of three- and four-dimensional heteronuclear nuclear magnetic resonance methods. Increased sophistication in labelling strategies, use of pulse-field gradients and the application of these methods at higher magnetic fields has, in combination with improved software, allowed studies of the structure, interactions and dynamics of significantly larger proteins (now up to approximately 270 amino acid residues).

The blind watchmaker and rational protein engineering

In the present review some scientific areas of key importance for protein engineering are discussed, such as problems involved in deducting protein sequence from DNA sequence (due to posttranscriptional editing, splicing and posttranslational modifications), modelling of protein structures by homology, NMR of large proteins (including probing the molecular surface with relaxation agents), simulation of protein structures by molecular dynamics and simulation of electrostatic effects in proteins (including pH-dependent effects). It is argued that all of these areas could be of key importance in most protein engineering projects, because they give access to increased and often unique information. In the last part of the review some potential areas for future applications of protein engineering approaches are discussed, such as non-conventional media, de novo design and nanotechnology.