Nuclear magnetic resonance observation and dynamics of specific amide protons in T4 lysozyme

We have produced T4 lysozyme using a bacterial expression system which allows efficient incorporation of isotopically labeled amino acids in lysozyme. By using conditions that repress the expression of various transaminases, we have incorporated 15N-labeled amino acid into the five phenylalanine residues of the protein. The relatively large spin--spin coupling (87 +/- 3 Hz) between the 15N nucleus and the phenylalanine amide protons may then be exploited in a variety of ways to selectively observe the five phenylalanine amide proton resonances. These include a simple "echo difference" technique which displays the amide proton resonances in one dimension and a "forbidden echo" technique [Bax, A., Griffey, R. H., & Hawkins, B.L. (1983) J. Magn. Reson. 55, 301-335] which gives two-dimensional information allowing the proton and 15N chemical shifts of each amide to be determined. With these approaches, all five phenylalanine amide protons give resolved resonances. Deuterium exchange experiments demonstrate that three of the five resonances are slow to exchange (half-times of about 1 week at pH 5.5 and 4 degrees C) while the other two are rapid with complete exchange in hours or less. These observations correlate well with the secondary structure of the protein which shows three residues in alpha-helical regions and two residues in surface-exposed environments. This approach of isotopic substitution on nitrogen or carbon atoms is of general utility and should allow virtually any proton on a protein of molecular weight 20 000 or thereabout to be selectively observed.     

Conformational analysis of Epac activation using amide hydrogen/deuterium exchange mass spectrometry

Exchange proteins directly activated by cAMP (Epac) play important roles in mediating the effects of cAMP through the activation of downstream small GTPases, Rap. To delineate the mechanism of Epac activation, we probed the conformation and structural dynamics of Epac using amide hydrogen/deuterium exchange and structural modeling. Our studies show that cAMP induces significant conformational changes that lead to a spatial rearrangement of the regulatory components of Epac and allows the exposure of the catalytic core for effector binding without imposing significant conformational change on the catalytic core. Homology modeling and comparative structural analyses of the cAMP binding domains of Epac and cAMP-dependent protein kinase (PKA) lead to a model of Epac activation, in which Epac and PKA activation by cAMP employs the same underlying principle, although the detailed structural and conformational changes associated with Epac and PKA activation are significantly different.     

Equilibrium amide hydrogen exchange and protein folding kinetics

The classical Linderstrøm-Lang hydrogen exchange (HX) model is extended to describe the relationship between the HX behaviors (EX1 and EX2) and protein folding kinetics for the amide protons that can only exchange by global unfolding in a three-state system including native (N), intermediate (I), and unfolded (U) states. For these slowly exchanging amide protons, it is shown that the existence of an intermediate (I) has no effect on the HX behavior in an off-pathway three-state system (IUN). On the other hand, in an on-pathway three-state system (UIN), the existence of a stable folding intermediate has profound effect on the HX behavior. It is shown that fast refolding from the unfolded state to the stable intermediate state alone does not guarantee EX2 behavior. The rate of refolding from the intermediate state to the native state also plays a crucial role in determining whether EX1 or EX2 behavior should occur. This is mainly due to the fact that only amide protons in the native state are observed in the hydrogen exchange experiment. These new concepts suggest that caution needs to be taken if one tries to derive the kinetic events of protein folding from equilibrium hydrogen exchange experiments.     

Progress in computation and amide hydrogen exchange for prediction of protein-protein complexes

The macromolecular docking problem that must be solved for experimental biologists is prediction of the structures of complexes for which the components are known or reliably modeled in the unbound state, but the structure of the complex is unknown. The current state of the art in macromolecular docking is such that solving this problem usually requires supplementary experimental chemical and/or biological information to evaluate computational predictions. Amide (1)H/(2)H exchange measured by mass spectroscopy is a promising approach for obtaining such information, because it can reveal interfacial regions of each member of the complex and identify regions of conformational flexibility in the structure. In a previous article (Anand et al., Proc Natl Acad Sci USA 2003;100:13264-13269), we used (1)H/(2)H exchange data to predict the structure of a complex between regulatory and catalytic subunits of protein kinase A. Comparison of the prediction with a recent crystal structure determination (Kim et al., Science 2005;307:690-696) showed large conformational change in the regulatory subunit on formation of the complex. Analysis of the prediction, previous CAPRI results, novel data processing methods for the (1)H/(2)H exchange data, and new fragment docking computations give grounds for cautious optimism that this method can be useful even in cases of substantial conformational change.     

Detection and characterization of an early folding intermediate of T4 lysozyme using pulsed hydrogen exchange and two-dimensional NMR.

Two-dimensional 1H-15N NMR techniques combined with pulsed hydrogen-deuterium exchange have been used to characterize the folding pathway of T4 lysozyme. In the unfolded state, there is little differential protection of the various amides from hydrogen exchange. In the native folded structure, 84 amides of the 164 residues are sufficiently spectrally resolved and protected from solvent exchange to serve as probes of the folding pathway. These probes are located in both the N-terminal and C-terminal domains of the native folded structure of the protein. The studies described here show that at least one intermediate is formed early during refolding at low denaturant concentrations. This intermediate (or intermediates) forms very rapidly (within the 10-ms temporal resolution of our mixing device) under the conditions used and is completed at least 10 times faster than the overall folding event. The intermediate(s) protect(s) from exchange a subset of amides in the N-terminal and C-terminal regions of the protein. In the final folded states these protected regions correspond to two alpha-helices and a beta-sheet region. These amides are protected from exchange by factors between 20 and 200 as compared to the fully unfolded protein. Protection of this magnitude is consistent with the formation of somewhat exposed secondary structure in these regions and could represent a "molten globule"-like or a "framework"-like structure for the intermediate(s) in which specific parts of the sequence form isolated secondary structures that are not stabilized by extensive tertiary interactions.     

The influence of glycerol on hydrogen isotope exchange in lysozyme

Hydrogen isotope exchange rates for lysozyme in glycerol cosolvent mixtures [D. G. Knox and A. Rosenberg (1980) Biopolymers 19 , 1049–1068] have been analyzed as functions of solvent viscosity and glycerol activity in an attempt to determine which solvent properties influence protein internal dynamics. The effect of glycerol on the fast- and slow-exchanging protons is different. Slow-exchanging protons [H(t) < 20] are slowed by ever-increasing amounts as H(t) decreases. However, comparison with data for the effect of glycerol on the thermal unfolding of lysozyme [K. Gekko (1982) J. Biochem. 19 , 1197–1204] indicates that the large decrease in exchange rates for the slow protons is not consistent with a local unfolding mechanism of exchange. These effects are also too large to be easily rationalized in terms of solvent viscosity. Instead, we suggest that the large effect of glycerol on exchange of the slow protons is due to a “compression” of the protein, as a result of thermodynamically unfavorable interactions of glycerol with the protein surface. This reduces the protein void volume, which in turn decreases the probability of conformational transitions required for exchange of the slowest protons. Present data do not allow a distinction to be made between thermodynamic (glycerol activity) and dynamic (solvent viscosity) influences on exchange rates for the fast-exchanging protons, although the effect of glycerol on these protons is also probably too large to be consistent with a local unfolding mechanism. In this case, glycerol decreases the rate of catalyst diffusion within the protein matrix, either by decreasing the probabilities or amplitudes of “gating” reactions that allow passage of the catalyst from the solvent to the exchange site, or by increasing the relaxation times for these conformational rearrangements.     

Kinetics of unfolding and folding from amide hydrogen exchange in native ubiquitin

Amide hydrogen (NH) exchange is one of the few experimental techniques with the potential for determining the thermodynamics and kinetics of conformational motions at nearly every residue in native proteins. Quantitative interpretation of NH exchange in terms of molecular motions relies on a simple two-state kinetic model: at any given slowly exchanging NH, a closed or exchange-incompetent conformation is in equilibrium with an open or exchange-competent conformation. Previous studies have demonstrated the accuracy of this model in measuring conformational equilibria by comparing exchange data with the thermodynamics of protein unfolding. We report here a test of the accuracy of the model in determining the kinetics of conformational changes in native proteins. The kinetics of folding and unfolding for ubiquitin have been measured by conventional methods and compared with those derived from a comprehensive analysis of the pH dependence of exchange in native ubiquitin. Rate constants for folding and unfolding from these two very different types of experiments show good agreement. The simple model for NH exchange thus appears to be a robust framework for obtaining quantitative information about molecular motions in native proteins.     

Analysis of conformational changes during activation of protein kinase Pak2 by amide hydrogen/deuterium exchange

During apoptotic stress, protein kinase Pak2 is cleaved by caspase 3 to form a heterotetramer that is constitutively activated following autophosphorylation. The active protein kinase migrates slightly slower than the inactive holoenzyme when analyzed by gel filtration, suggesting an expanded conformation. Activation of Pak2 comprises a series of structural changes resulting from caspase cleavage, ATP binding, and autophosphorylation of Pak2. Changes at each step were individually analyzed by amide hydrogen/deuterium exchange coupled with mass spectrometry and compared with inactive Pak2. The auto-inhibited form was shown to bind ATP in the active site, with minor changes in the glycine loop and the autoinhibitory domain (AID). Caspase cleavage produced significant changes in solvent accessibility in the AID and upper lobe of the catalytic domain. Cleavage of ATP-bound Pak2 relaxes the allosteric inhibition, as shown by increased solvent accessibility in the upper and lower lobes, including the G-helix, facilitating the autophosphorylation of two sites required for activation, Ser-141 in the regulatory domain and Thr-402 in the catalytic domain. Autophosphorylation increased the amide hydrogen/deuterium exchange solvent accessibility of the contact region between the AID and the G-helix, the E-F loop, and the N terminus. Thus, activation of Pak2 via caspase cleavage is associated with structural relaxation of Pak2 that allows for complete auto-phosphorylation, resulting in a more comprehensive solvent-exposed and conformationally dynamic enzyme.     

Atomic coordinates for T4 phage lysozyme.

Atomic coordinates are presented for the lysozyme from T4 bacteriophage. The coordinates were derived from a 2.4 Å resolution electron density map based on two isomorphous heavy-atom derivatives, interpreted in terms of the known amino acid sequence, and adjusted to have stereochemically acceptable bond lengths and angles.     

Flexibility and ligand exchange in a buried cavity mutant of T4 lysozyme studied by multinuclear NMR

The Leu99-->Ala mutant of T4 lysozyme contains a large internal cavity in the core of its C-terminal domain that is capable of reversibly binding small hydrophobic compounds. Although the cavity is completely buried, molecules such as benzene or xenon can exchange rapidly in and out. The dynamics of the unliganded protein have been compared to the wild-type protein by measuring the NMR spin relaxation rates of backbone amide and side chain methyl nuclei. Many residues surrounding the cavity were found to be affected by a chemical exchange process with a rate of 1500 +/- 200 s(-1), which is quenched upon addition of saturating amounts of the ligand xenon. The relationship between the structure, dynamics, and energetics of the T4 lysozyme mutant is discussed.     

Backbone dynamics of detergent-solubilized alamethicin from amide hydrogen exchange measurements.

Alamethicin is a 20 amino acid antibiotic peptide produced by the soil fungus Trichoderma viride. The peptide inserts into bacterial membranes and self-associates to form ion channels, but the details of this process are unknown. Residue-specific acid- and base-catalyzed exchange data were obtained for 16 of 18 backbone amides of alamethicin dissolved in sodium dodecyl sulfate micelles using high-resolution 2-dimensional heteronuclear nuclear magnetic resonance spectroscopy. To facilitate interpretation of the exchange data, we synthesized N-acetyl-alpha-aminoisobutyric acid-N'-methyl and N-acetyl-alanine-N'-methyl and measured the pD dependence of their hydrogen-deuterium exchange rates to determine the sequence-dependent inductive and steric effects of the alpha-aminoisobutyric acid residue. Intramolecular H-bonding in alamethicin was monitored through the exchange parameters kmin (minimum exchange rate) which indicate that the backbone is significantly more stable than the backbones of alanine-based helical peptides. Rapid exchange at Gly-11 suggests a highly local conformational flexibility in the middle of the peptide. Interactions with the detergent micelle were revealed by the exchange parameters pDmin (pD of minimum exchange) which suggest that the N-terminus of alamethicin interacts more strongly with the detergent micelle than does the C-terminus. A periodicity in pDmin difference data reveals that one surface of the helix interacts more strongly with the micelle. The surface consists of residues 1, 5, 9, 13, 16, and 20. The opposite face of the helix contains several polar residues (two glutamines and a glycine), suggesting that, on average, this face of the helix is directed toward the solvent. These results serve as a model for the interaction of the peptide with membranes containing anionic lipid. In combination with published molecular dynamics simulations [Gibbs et al. (1997) Biophys. J. 72, 2490-2495], the present results also offer insight into the mechanisms of hydrogen-deuterium exchange in helical peptides.     

Effects of co-operative ligand binding on protein amide NH hydrogen exchange

Amide protection factors have been determined from NMR measurements of hydrogen/deuterium amide NH exchange rates measured on assigned signals from Lactobacillus casei apo-DHFR and its binary and ternary complexes with trimethoprim (TMP), folinic acid and coenzymes (NADPH/NADP(+)). The substantial sizes of the residue-specific DeltaH and TDeltaS values for the opening/closing events in NH exchange for most of the measurable residues in apo-DHFR indicate that sub-global or global rather than local exchange mechanisms are usually involved. The amide groups of residues in helices and sheets are those most protected in apo-DHFR and its complexes, and the protection factors are generally related to the tightness of ligand binding. The effects of ligand binding that lead to changes in amide protection are not localised to specific binding sites but are spread throughout the structure via a network of intramolecular interactions. Although the increase in protein stability in the DHFR.TMP.NADPH complex involves increased ordering in the protein structure (requiring TDeltaS energy) this is recovered, to a large extent, by the stronger binding (enthalpic DeltaH) interactions made possible by the reduced motion in the protein. The ligand-induced protection effects in the ternary complexes DHFR.TMP.NADPH (large positive binding co-operativity) and DHFR.folinic acid.NADPH (large negative binding co-operativity) mirror the co-operative effects seen in the ligand binding. For the DHFR.TMP.NADPH complex, the ligand-induced protection factors result in DeltaDeltaG(o) values for many residues being larger than the DeltaDeltaG(o) values in the corresponding binary complexes. In contrast, for DHFR.folinic acid.NADPH, the DeltaDeltaG(o) values are generally smaller than many of those in the corresponding binary complexes. The results indicate that changes in protein conformational flexibility on formation of the ligand complex play an important role in determining the co-operativity in the ligand binding.     

Helix bending in alamethicin: molecular dynamics simulations and amide hydrogen exchange in methanol

Molecular dynamics simulations of alamethicin in methanol were carried out with either a regular alpha-helical conformation or the x-ray crystal structure as starting structures. The structures rapidly converged to a well-defined hydrogen-bonding pattern with mixed alpha-helical and 3(10)-helical hydrogen bonds, consistent with NMR structural characterization, and did not unfold throughout the 1-ns simulation, despite some sizable backbone fluctuations involving reversible breaking of helical hydrogen bonds. Bending of the helical structure around residues Aib10-Aib13 was associated with reversible flips of the peptide bonds involving G11 (Aib10-G11 or G11-L12 peptide bonds), yielding discrete structural states in which the Aib10 carbonyl or (rarely) the G11 carbonyl was oriented away from the peptide helix. These peptide bond reversals could be accommodated without greatly perturbing the adjacent helical structure, and intramolecular hydrogen bonding was generally maintained in bent states through the formation of new (non-alpha or 3[10]) hydrogen bonds with good geometries: G11 NH-V9 CO (inverse gamma turn), Aib13 NH-Aib8 CO (pi-helix) and, rarely, L12 NH- Q7 NH (pi-helix). These observations may reconcile potentially conflicting NMR structural information for alamethicin in methanol, in which evidence for conformational flexibility in the peptide sequence before P14 (G11-Aib13) contrasts with the stability of backbone amide NH groups to exchange with solvent. Similar reversible reorientation of the Thr11-Gly12 peptide bond of melittin is also observed in dynamics simulations in methanol (R. B. Sessions, N. Gibbs, and C. E. Dempsey, submitted). This phenomenon may have some role in the orientation of the peptide carbonyl in solvating the channel lumen in membrane ion channel states of these peptides.     

Measurement of amide hydrogen exchange rates with the use of radiation damping

A simple method for measuring amide hydrogen exchange rates is presented, which is based on the selective inversion of water magnetization with the use of radiation damping. Simulations show that accurate exchange rates can be measured despite the complications of radiation damping and cross relaxation to the exchange process between amide and water protons. This method cannot eliminate the contributions of the exchange-relayed NOE and direct NOE to the measured exchange rates, but minimize the direct NOE contribution. In addition, the amides with a significant amount of such indirect contributions are possible to be identified from the shape of the exchange peak intensity profiles or/and from the apparent relaxation rates of amide protons which are extracted from fitting the intensity profiles to an equation established here for our experiment. The method was tested on ubiquitin and also applied to an acyl carrier protein. The amide exchange rates for the acyl carrier protein at two pHs indicate that the entire protein is highly dynamic on the second timescale. Low protection factors for the residues in the regular secondary structural elements also suggest the presence of invisible unfolded species. The highly dynamic nature of the acyl carrier protein may be crucial for its interactions with its substrate and enzymes.     

Systematic mutation of bacteriophage T4 lysozyme

Amber mutations were introduced into every codon (except the initiating AUG) of the bacteriophage T4 lysozyme gene. The amber alleles were introduced into a bacteriophage P22 hybrid, called P22 e416, in which the normal P22 lysozyme gene is replaced by its T4 homologue, and which consequently depends upon T4 lysozyme for its ability to form a plaque. The resulting amber mutants were tested for plaque formation on amber suppressor strains of Salmonella typhimurium. Experiments with other hybrid phages engineered to produce different amounts of wild-type T4 lysozyme have shown that, to score as deleterious, a mutation must reduce lysozyme activity to less than 3% of that produced by wild-type P22 e416. Plating the collection of amber mutants covering 163 of the 164 codons of T4 lysozyme, on 13 suppressor strains that each insert a different amino acid substitutions at every position in the protein (except the first). Of the resulting 2015 single amino acid substitutions in T4 lysozyme, 328 were found to be sufficiently deleterious to inhibit plaque formation. More than half (55%) of the positions in the protein tolerated all substitutions examined. Among (N-terminal) amber fragments, only those of 161 or more residues are active. The effects of many of the deleterious substitutions are interpretable in light of the known structure of T4 lysozyme. Residues in the molecule that are refractory to replacements generally have solvent-inaccessible side-chains; the catalytic Glu11 and Asp20 residues are notable exceptions. Especially sensitive sites include residues involved in buried salt bridges near the catalytic site (Asp10, Arg145 and Arg148) and a few others that may have critical structural roles (Gly30, Trp138 and Tyr161).     

Thermal-induced unfolding domains in aldolase identified by amide hydrogen exchange and mass spectrometry.

Amide hydrogen exchange has been measured in short segments of intact rabbit muscle aldolase at temperatures of 14-50 degrees C by the protein fragmentation/mass spectrometry method (Zhang Z, Smith DL, 1993, Protein Sci 2:522-531). Deuterium levels in some segments did not change over the temperature range of the measurements, whereas deuterium levels in other segments increased rapidly with temperature. These results demonstrate that the equilibrium constant for local unfolding, Kunf, of some segments increases with temperature in the low temperature range (14-30 degrees C) of this study. Aldolase begins to lose activity at temperatures above 40 degrees C. In the 40-50 degrees C temperature range, Kunf is greater than 10(-4) in some regions and less than 10(-6) in other regions. This wide range of regional stability in the temperature range where aldolase begins to denature is interpreted in terms of cooperative unfolding/folding domains. Regions of highest stability were located along the hydrophobic subunit binding surface. It is proposed that hydrogen exchange might be used to identify unfolding domains in multidomain proteins whose thermodynamic properties have been determined by differential scanning calorimetry.