Core-directed protein design. I. An experimental method for selecting stable proteins from combinatorial libraries.

The design of proteins represents a significant challenge to modern-day structural biology. A major obstacle here is the specification of well-packed hydrophobic cores to drive the folding and stabilization of the target. Computational approaches have been used to alleviate this by testing alternate sequences prior to the production and characterization of a few proteins. Here we present the experimental counterpart of this approach. We selected stable variants from a library of ubiquitin hydrophobic-core mutants as follows. Hexahistidine-tagged proteins were displayed on the surface of phage. These protein-phage were immobilized onto Ni-coated surfaces. The bound fusion-phage were treated with protease to remove unstable or poorly folded proteins. Stable phage fusions were eluted and infected into Escherichia coli, which allowed amplification for further selection, sequencing, or protein expression. Two Ni-derivatized supports were tested: Ni-NTA chips for surface plasmon resonance (SPR) and Ni-NTA agarose beads. SPR had an advantage in that the selection process could be monitored directly. This allowed individual clones and experimental conditions to be tested rapidly prior to preparative panning of the library, which was carried out using Ni-NTA agarose beads. We demonstrate the method by selecting stable core mutants of ubiquitin, the characterization of which is described in the following paper [Finucane, M. D., and Woolfson, D. N. (1999) Biochemistry 38, XXXXX-XXXXX]. As our method selects only on the basis of structure and stability, it will be of use in improving the stabilities and structural specificities of proteins of de novo design, and in establishing rules that link sequence and structure.

Core-directed protein design. II. Rescue of a multiply mutated and destabilized variant of ubiquitin.

We have applied the method described in the preceding paper [Finucane, M. D., et al. (1999) Biochemistry 38, 11604-11612], namely, stability-based selection using phage display, to explore the sequence requirements for packing in the hydrophobic core of ubiquitin. In contrast to the parent protein, which was a structurally compromised mutant, the selected variants could be overexpressed and purified in yields for structural studies. In particular, CD and NMR measurements showed that the selectants folded correctly to stable native-like structures. These points demonstrate the utility of our core-directed method for stabilizing and redesigning proteins. In addition and in contrast to foregoing studies on other proteins, which suggest that hydrophobic cores permit substitutions provided that hydrophobicity and core volumes are generally conserved, we find that the core of ubiquitin is surprisingly intolerant of amino acid substitutions; variants that survived our selection showed a clear consensus for the wild-type sequence. It is probable that our results differed from those from other groups for two reasons. First, ubiquitin may be unusual in that it has strict sequence requirements for its structure and stability. We discuss this result in light of sequence conservation in the eukaryotic ubiquitins and proteins of the ubiquitin structural superfamily. Second, our mutants were selected solely on the basis of stability, in contrast to the other studies that rely on function-based selection. The latter may lead to proteins that are more plastic and tolerant of substitutions.

The pH dependence of hydrogen-deuterium exchange in trp repressor: the exchange rate of amide protons in proteins reflects tertiary interactions, not only secondary structure.

The pH dependence of amide proton exchange rates have been measured for trp-repressor. One class of protons exchanges too fast to be measured in these experiments. Among the protons that have measurable hydrogen-deuterium exchange rates, two additional classes may be distinguished. The second class of protons are in elements of secondary structure that are mostly on the surface of the protein, and exchange linearly with increasing base concentration (log kex versus pH). The third class of amide protons is characterized by much higher protection against exchange at higher pH. These protons are located in the core of the protein, in helices B and C. The exchange rate in the core region does not increase linearly with pH, but rather goes through a minimum around pH 6. The mechanism of exchange for the slowly exchanging core protons is interpreted in terms of the two-process model of Hilton and Woodward (1979, Biochemistry 18:5834-5841), i.e., exchange through both a local mechanism that does not require unfolding of the protein, and a mechanism involving global unfolding of the protein. The increase in exchange rates at low pH is attributed to a partial unfolding of the repressor. It is concluded that the formation of secondary structure alone is insufficient to account for the high protection factors seen in the core of native proteins at higher pH, and that tertiary interactions are essential to stabilize the structure.

Mechanism of hydrogen-deuterium exchange in trp repressor studied by 1H-15N NMR.

Amide proton exchange rates have been measured for fast-exchanging amides in trp aporepressor, and compared with the rates measured in the holorepressor. The results indicate that the presence of the ligand stabilizes all of the amide protons in the molecule against exchange, not just those whose access to solvent it directly hinders. This global hindering of the exchange process by tryptophan implies that there is a non-random element in the transmission mechanism, so that damping of the exchange in one part of the molecule also damps exchange in another region. This damping at a distance is not associated with any measurable changes in the intervening average secondary structure. This suggests the existence of a concerted dynamic process in the protein backbone that is modulated by ligand binding and in turn affects the observed backbone proton exchange.

Investigation of protein amide-proton exchange by 1H longitudinal spin relaxation.

The general theory of the Linderstrøm-Lang model, which, in simplified form, is widely used for the interpretation of NMR data on backbone-proton exchange in proteins, is systematically discussed. An experimental protocol for testing the applicability of the customary simplifications is described and experimental data which require the application of the general, rather than of the simplified, theory are presented. The appearance of pH-dependent biexponential magnetization recovery for any amide group in a protein is an unequivocal indication that the conventional simplified versions of the model do not apply.

Solution dynamics of the trp repressor: a study of amide proton exchange by T1 relaxation.

The amide proton exchange rates of Escherichia coli trp repressor have been measured through their effects on the longitudinal relaxation rates of the amide protons. Three types of exchange regimes have been observed: (1) slow exchange (on a minute/hour time-scale), measurable by isotope exchange, but not by relaxation techniques in the core of the molecule; (2) relatively rapid exchange, with the rates on a T1 relaxation time-scale (seconds) in the DNA-binding region and (3) very fast exchange at the N and C termini. The results have been analyzed in terms of the two-site exchange model originally proposed by Linderstrøm-Lang, and of a three-site extension of the model. The values of the intrinsic exchange rates calculated using the two-state model agree with the values expected from the studies of Englander and co-workers for the very fast case of the chain terminals, but disagree with the literature values by two orders of magnitude in the intermediate case found in the DNA-binding region. The implication of these findings is that the "open" state of the two-state model in the DNA-binding region is not completely open and has an intrinsic exchange rate different from that of a random coil peptide. Alternatively, if the literature values of the intrinsic exchange rates are assumed to apply to the open states in all parts of the repressor molecule, two "closed" helical states have to be postulated, in slow exchange with each other, with only one of them in rapid exchange with the open state and hence with the solvent. Kinetically, the two models are indistinguishable.

A study of the stabilization of tetrahedral adducts by trypsin and delta-chymotrypsin.

delta-Chymotrypsin has been alkylated by 1-13C- and 2-13C-enriched tosylphenylalanylchloromethane. In the intact inhibitor derivative, signals due to the 1-13C- and 2-13C-enriched carbon atoms have chemical shifts which titrate from 55.10 to 59.50 p.p.m. and from 99.10 to 103.66 p.p.m. respectively with similar pKa values of 8.99 and 8.85 respectively. These signals are assigned to a tetrahedral adduct formed between the hydroxy group of serine-195 and the inhibitor. An additional signal at 58.09 p.p.m. and at 204.85 p.p.m. in the 1-13C- and 2-13C-enzyme-inhibitor derivatives respectively does not titrate when the pH is changed and it is assigned to alkylated methionine-192. On denaturation/autolysis of the 1-13C-enriched enzyme-inhibitor derivative these signals associated with the intact inhibitor derivative are no longer detected, and a new signal, which titrates from 56.28 to 54.84 p.p.m. with a pKa of 5.26, is detected. The titration shift of this signal is assigned to the deprotonation of the imidazolium cation of alkylated histidine-57 in the denatured/autolysed enzyme-inhibitor derivative. Model compounds which form stable hydrates and hemiketals in aqueous solutions have been synthesized. By comparing the 13C titration shifts of these model compounds with those of the 13C enriched trypsin- and delta-chymotrypsin-inhibitor derivatives, we deduce that, in both of the intact enzyme-inhibitor derivatives, the zwitterionic tetrahedral adduct containing the imidazolium cation of histidine-57 and the hemiketal oxyanion predominates at alkaline pH values. It is estimated that in both the trypsin and delta-chymotrypsin-inhibitor derivatives the concentration of this zwitterionic tetrahedral adduct is 10,000-fold greater than it would be in water. We conclude that the pKa of the oxyanion of the hemiketal in the presence of the imidazolium cation of histidine-57 is 7.9 and 8.9 in the trypsin and delta-chymotrypsin-inhibitor derivatives respectively and that the pKa of the imidazolium cation of histidine-57 is greater than 7.9 and greater than 8.9 when the oxyanion is present as its conjugate acid, whereas, when the oxyanion is present, the pKa of the imidazolium cation is greater than 11 in both enzyme-inhibitor derivatives. We discuss how these enzymes preferentially stabilize zwitterionic tetrahedral adducts in the intact enzyme-inhibitor derivatives and how they could stabilize similar tetrahedral intermediates during catalysis. It is suggested that substrate binding could raise the pKa of the imidazolium cation of histidine-57 before tetrahedral-intermediate formation which would explain the enhanced nucleophilicity of the hydroxy group of serine-195.(ABSTRACT TRUNCATED AT 400 WORDS)

A study of the relaxation parameters of a 13C-enriched methylene carbon and a 13C-enriched perdeuteromethylene carbon attached to chymotrypsin.

L-1-Chloro-4-phenyl-3-tosylamido[1-13C]butan-2-one (Tos-[1-13C]Phe-CH2Cl) and Tos-[1-13C,2H2]Phe-CH2Cl were prepared and used to alkylate delta-chymotrypsin. The relaxation parameters of the 13C-n.m.r. signal resulting from the alkylation of histidine-57 in both enzyme-inhibitor complexes were determined at 1.88 T and 6.34 T as well as the spin-lattice relaxation times of the backbone alpha-carbon atoms of the unenriched Tos-Phe-CH2-delta-chymotrypsin complex. It is concluded that the species examined do not have significant internal librational motions and that the rotational correlation time of the monomeric enzyme-inhibitor complex is 16.0 +/- 3.2 ns. The signal from the 13C-enriched atom of Tos-[1-13C,2H2]Phe-CH2Cl is split into a quintet (JCD = 23 Hz) whereas in the Tos-[1-13C,2H2]Phe-CH2-delta-chymotrypsin complex the signal from the 13C-enriched inhibitor carbon atom is decoupled. This decoupled signal had linewidths of 16 +/- 3 Hz and 52 +/- 2 Hz at 1.88 T and 6.34 T respectively, whereas linewidths at 40 +/- 2 Hz and 53 +/- 4 Hz were obtained for the same signal in the Tos-[1-13C]Phe-CH2-delta-chymotrypsin complex at 1.88 T and 6.34 T respectively. Therefore whereas deuteration produces a 2.5-fold reduction in linewidth at 1.88 T there is no significant decrease in the linewidth at 6.34 T. This result is explained by using the rigid rotor model, which predicts that the quadrupolar spin-lattice relaxation rate will be faster at low field strengths, resulting in more efficient deuterium decoupling by scalar relaxation of the second kind at lower field strengths. It is also predicted that deuterium decoupling by scalar relaxation will become less efficient as rotational correlation times increase. The consequences of these predictions for the detection of 13C-enriched atomic probes of proteins are discussed. It is also shown that a spin-echo pulse sequence can be used to remove signals due to protonated carbon atoms without attenuating the signal due to deuterated carbon atoms.

A 13C-n.m.r. investigation of the ionizations within an inhibitor--alpha-chymotrypsin complex. Evidence that both alpha-chymotrypsin and trypsin stabilize a hemiketal oxyanion by similar mechanisms.

13C-n.m.r. was used to investigate the structure of the inhibitor enzyme complex formed when alpha-chymotrypsin is alkylated by L-1-chloro-4-phenyl-3-tosylamido-[2-13C]butan-2-one. Two signals are detected. The one at 204.82 +/- 0.11 p.p.m. does not titrate from pH 3 to 9 and is assigned to alkylated methionine-192. The second signal titrates from 99.08 p.p.m. to 103.44 p.p.m. with pKa 8.67. This signal is assigned to a tetrahedral adduct formed between the hydroxy group of serine-195 and the inhibitor. The titration shift of the tetrahedral adduct is ascribed to the ionization of the hemiketal hydroxy group. It is proposed that the resulting oxyanion is stabilized by interaction with the imidazolium ion of histidine-57. It is argued that this interaction must raise the pKa of at least 70% of histidine-57 to greater than 11. On denaturation/autolysis of the inhibitor-enzyme complex neither of the signals associated with the intact complex is detected, but a new signal is observed that titrates from 203.52 p.p.m. to 206.08 p.p.m. with pKa = 5.27. This titration shift is assigned to the ionization of the imidazolium ion of alkylated histidine, confirming that the inhibitor has alkylated histidine-57. The significance of these results for the catalytic mechanism of the serine proteinases is discussed.