Comparison between the unfolding rate and structural fluctuations in native lysozyme--effects of denaturants, ligand binding, and intrachain cross-linking on hydrogen exchange and unfolding kinetics

The hydrogen-exchange reactions of peptide NH groups in lysozyme were studied by the change in the intensity of the amide II band in the ir spectrum. The slowest exchanging hydrogens, which are involved in intramolecular hydrogen bonding, are further divided into two groups at lower temperatures; half of them are exchanged through local unfolding and the other half through major cooperative unfolding. In order to study the correlation of the change in hydrogen-exchange rates with the change in the unfolding rate constant, we observed the effects of intrachain cross-linking, the addition of denaturant and ligand binding on the exchange rates through local unfolding. Although the exchange rate through major unfolding is greatly decreased by intrachain cross-linking between Glu 35 and Trp 108 (1/22000), the exchange rate through local unfolding is only slightly decreased (1/20). Even at higher temperatures, where most intact lysozyme molecules unfold, the folded conformation of cross-linked lysozyme remains compact, and no intermediate exists in which many side-chain atoms are packed loosely so that the hydrogen-exchange reaction occurs rapidly. Neither the addition of 2-PrOD molecules nor (NAG)₃ binding affects the exchange rates through local unfolding. Our experiments confirm that the change in the unfolding rate constant does not correlate with the change in fluctuations in the relatively flexible hydrogen-bonded structure through which the exchange of peptide hydrogens takes place.     

Hydrogen exchange of individual amide protons in the E. coli lac repressor DNA-binding domain: a nuclear magnetic resonance study

Proton exchange in lac repressor headpiece was studied by COSY and 2D NOE spectroscopy. The exchange rates of amide protons, stabilized by the hydrogen bonds of the three alpha-helices of the headpiece, could be determined quantitatively. The exchange rates in these helices showed repetitive patterns of about three to four residues. A correlation with the position of the amide proton in the interior or the exterior of the alpha-helix of the protein was found. The exchange data strongly support the validity of the three-dimensional structure, as determined recently (Kaptein, R. et al., J. Mol. Biol. 182, 179-182 (1985)).     

Unfolding and refolding of bovine beta-lactoglobulin monitored by hydrogen exchange measurements.

Bovine beta-lactoglobulin (beta-LG) is a widely studied protein belonging to the lipocalin family, whose structural characterisation has been reported by X-ray crystallography and NMR studies at physiological and acidic pH, respectively. Bovine beta-LG consists of nine antiparallel beta-sheets and a terminal alpha-helix segment. The beta-sheets form a calyx structure with a hydrophobic buried cluster conferring stability to the protein while a hydrophobic surface patch provides stabilising interactions between the barrel and the flanking terminal helix. Here, the stability and the folding properties of bovine beta-LG in the presence of a chemical denaturant is probed. The analysis of the NMR spectra recorded in aqueous solution with increasing amounts of urea revealed that the intensities of the backbone cross-peaks decrease upon increasing urea concentration, while their secondary shifts do not change significantly on going from 0 to 5 M urea, thus suggesting the presence of slow exchange between native and unfolded protein. Hydrogen exchange measurements at different urea concentrations were performed in order to evaluate the exchange rates of individual backbone amide protons. The opening reactions that determine protein exchange can be computed for the most slowly exchanging hydrogen atoms, and the measured exchange rates and the corresponding free energies can be expressed in terms of the equilibrium energetic for the global transition and the local fluctuations. Most of the residues converge to define a common isotherm identifying a unique cooperative folding unit, encompassing all the strands, except strand betaI, and the terminal region of the helix. The amides that do not join the same global unfolding isotherm are characterised by low DeltaGH20op and especially by low m values, indicating that they are already substantially exposed in the native state. A two-state unfolding model N <==> U is therefore proposed for this rather big protein, in agreement with CD data. Renaturation studies show that the unfolding mechanism is reversible up to 6 M urea and suggest a similar unfolding and refolding pathway. Present results are discussed in light of the hypothesis of an alpha-->beta transition proposed for bovine beta-LG refolding.     

Complex of lac repressor headpiece with a 14 base-pair lac operator fragment studied by two-dimensional nuclear magnetic resonance

Two-dimensional nuclear Overhauser enhancement spectra are presented of the complex of lac repressor headpiece with a 14 base-pair lac operator fragment. Analysis of nuclear Overhauser enhancements observed between protein and DNA shows that the second helix of the headpiece ("the recognition helix") binds in the major groove of DNA as has been suggested, but that the orientation of this helix is approximately 180 degrees different from the proposed models.     

H NMR study of a complex between the lac repressor headpiece and a 22 base pair symmetric lac operator

A complex between the lac repressor headpiece and a fully symmetric tight-binding 22 bp lac operator was studied by 2D NMR. Several 2D NOE spectra were recorded for the complex in both H2O and 2H2O. Many NOE cross-peaks between the headpiece and DNA could be identified, and changes in the chemical shift of the DNA protons upon complex formation were analyzed. Comparison of these data with those obtained for a complex between the headpiece and a 14 bp half-operator, studied previously [Boelens, R., Scheek, R. M., Lamerichs, R. M. J. N., de Vlieg, J., van Boom, J. H., & Kaptein, R. (1987) in DNA-ligand interactions (Guschlbauer, W., & Saenger, W., Eds.) pp 191-215, Plenum, New York], shows that two headpieces form a specific complex with the 22 bp lac operator in which each headpiece binds in the same way as found for the 14 bp complex. The orientation of the recognition helix in the major groove of DNA in these complexes is opposite with respect to the dyad axis to that found for other repressors.     

Insight into a Random Coil Conformation and an Isolated Helix: Structural and Dynamical Characterisation of the C-Helix Peptide from Hen Lysozyme

A 17 residue peptide corresponding to the C-helix of hen lysozyme (residues 86 to 102) has been investigated in detail to assess the factors that determine its conformation in both aqueous and trifluoroethanol (TFE) solutions. A thorough characterisation of the peptide by CD and NMR techniques under both conditions has been performed including the determination of complete NMR proton sequential assignments, and measurement of NOE effects,³J_HNαcoupling constants, temperature coefficients and residue-specific hydrogen-exchange rates. In water, the peptide adopts a largely unstructured conformation and NMR data, particularly coupling constants and chemical shift deviations, have been shown to agree closely with predictions from a model for a random coil based on the φ,ψ distributions in a protein database. This indicates that under these conditions the intrinsic conformational preferences of the individual amino acid residues are the dominating factors that determine the population of conformers adopted. With increasing concentrations of TFE a cooperative transition to an extensively helical conformation occurs and the resultant changes in CαH chemical shifts have been shown to correlate with the changes in φ,ψ populations. Using NOE and coupling constant data for this state, an ensemble of structures has been calculated and provides a model for a helix in the absence of tertiary interactions. In this model fluctuations, which increase in amplitude towards the termini, occur about the average helical φ,ψ angles and are responsible for increasing the values of³J_HNαcoupling constants above those anticipated for a static helix. The residue-specific rates of hydrogen exchange for the peptide in 50% TFE-d₃are consistent with such a model, the maximum protection from exchange being observed for residues in the centre of the helix.     

Investigation of the four cooperative unfolding units existing in the MICAL-1 CH domain

The solution structure of human MICAL-1 calpolnin homology (CH) domain is composed of six alpha helices and one 3(10) helix. To study the unfolding of this domain, we carry out native-state hydrogen exchange, intrinsic fluorescence and far-UV circular dichroism experiments. The free energy of unfolding, DeltaG(H2O), is calculated to be 7.11+/-0.58 kcal mol(-1) from GuHCl denaturation at pH 6.5. Four cooperative unfolding units are found using native-state hydrogen exchange experiment. Forty-seven slow-exchange residues can be studied by native-state hydrogen exchange experiments. From the concentration dependence of exchange rates, free energy of amide hydrogen with solvent, DeltaG(HX) and m-value (sensitivity of exposure to denaturant) are obtained, which reveal four cooperative unfolding units. The slowest exchanging protons are distributed throughout the whole hydrophobic core of the protein, which might be the folding core. These results will help us understand the structure of MICAL-1 CH domain more deeply.     

Hexafluoroacetone hydrate as a structure modifier in proteins: characterization of a molten globule state of hen egg-white lysozyme.

A molten globule-like state of hen egg-white lysozyme has been characterized in 25% aqueous hexafluoroacetone hydrate (HFA) by CD, fluorescence, NMR, and H/D exchange experiments. The far UV CD spectra of lysozyme in 25% HFA supports retention of native-like secondary structure while the loss of near UV CD bands are indicative of the overall collapse of the tertiary structure. The intermediate state in 25% HFA exhibits an enhanced affinity towards the hydrophobic dye, ANS, and a native-like tryptophan fluorescence quenching. 1-D NMR spectra indicates loss of native-like tertiary fold as evident from the absence of ring current-shifted 1H resonances. CD, fluorescence, and NMR suggest that the transition from the native state to a molten globule state in 25% HFA is a cooperative process. A second structural transition from this compact molten globule-like state to an "open" helical state is observed at higher concentrations of HFA (> or = 50%). This transition is characterized by a dramatic loss of ANS binding with a concomitant increase in far UV CD bands. The thermal unfolding of the molten globule state in 25% HFA is sharply cooperative, indicating a predominant role of side-chain-side-chain interactions in the stability of the partially folded state. H/D exchange experiments yield higher protection factors for many of the backbone amide protons from the four alpha-helices along with the C-terminal 3(10) helix, whereas little or no protection is observed for most of the amide protons from the triple-stranded antiparallel beta-sheet domain. This equilibrium molten globule-like state of lysozyme in 25% HFA is remarkably similar to the molten globule state observed for alpha-lactalbumin and also with the molten globule state transiently observed in the kinetic refolding experiments of hen lysozyme. These results suggest that HFA may prove generally useful as a structure modifier in proteins.     

Determinants of protein hydrogen exchange studied in equine cytochrome c.

The exchange of a large number of amide hydrogens in oxidized equine cytochrome c was measured by NMR and compared with structural parameters. Hydrogens known to exchange through local structural fluctuations and through larger unfolding reactions were separately considered. All hydrogens protected from exchange by factors greater than 10(3) are in defined H-bonds, and almost all H-bonded hydrogens including those at the protein surface were measured to exchange slowly. H-exchange rates do not correlate with H-bond strength (length) or crystallographic B factors. It appears that the transient structural fluctuation necessary to bring an exchangeable hydrogen into H-bonding contact with the H-exchange catalyst (OH(-)-ion) involves a fairly large separation of the H-bond donor and acceptor, several angstroms at least, and therefore depends on the relative resistance to distortion of immediately neighboring structure. Accordingly, H-exchange by way of local fluctuational pathways tends to be very slow for hydrogens that are neighbored by tightly anchored structure and for hydrogens that are well buried. The slowing of buried hydrogens may also reflect the need for additional motions that allow solvent access once the protecting H-bond is separated, although it is noteworthy that burial in a protein like cytochrome c does not exceed 4 angstroms. When local fluctuational pathways are very slow, exchange can become dominated by a different category of larger, cooperative, segmental unfolding reactions reaching up to global unfolding.     

Ligand Entry into Fatty Acid Binding Protein via Local Unfolding Instead of Gap Widening

Fatty acid binding proteins play an important role in the transportation of fatty acids. Despite intensive studies, how fatty acids enter the protein cavity for binding is still controversial. Here, a gap-closed variant of human intestinal fatty acid binding protein was generated by mutagenesis, in which the gap is locked by a disulfide bridge. According to its structure determined here by NMR, this variant has no obvious openings as the ligand entrance and the gap cannot be widened by internal dynamics. Nevertheless, it still takes up fatty acids and other ligands. NMR relaxation dispersion, chemical exchange saturation transfer, and hydrogen-deuterium exchange experiments show that the variant exists in a major native state, two minor native-like states, and two locally unfolded states in aqueous solution. Local unfolding of either βB–βD or helix 2 can generate an opening large enough for ligands to enter the protein cavity, but only the fast local unfolding of helix 2 is relevant to the ligand entry process.     

Exchange behavior of the H-bonded amide protons in the 3 to 13 helix of ribonuclease S

The preceding article shows that there are eight highly protected amide protons in the S-peptide moiety of RNAase S at pH 5, 0 degrees C. The residues with protected NH protons are 7 to 13, whose amide protons are H-bonded in the 3 to 13 alpha-helix, and Asp 14, whose NH proton is H-bonded to the CO group of Val47. We describe here the exchange behavior of these eight protected protons as a function of pH. Exchange rates of the individual NH protons are measured by 1H nuclear magnetic resonance in D2O. A procedure is used for specifically labeling with 1H only these eight NH protons. The resonance assignments of the eight protons are made chiefly by partial exchange, through correlating the resonance intensities in spectra taken when the peptide is bound and when it is dissociated from S-protein in 3.5 M-urea-d4, in D2O, pH 2.3, -4 degrees C. The two remaining assignments are made and some other assignments are checked by measurements of the nuclear Overhauser effect between adjacent NH protons of the alpha-helix. There is a transition in exchange behavior between pH 3, where the helix is weakly protected against exchange, and pH 5 where the helix is much more stable. At pH 3.1, 20 degrees C, exchange rates are uniform within the helix within a factor of two, after correction for different intrinsic exchange rates. The degree of protection within the helix is only 10 to 20-fold at this pH. At pH 5.1, 20 degrees C, the helix is more stable by two orders of magnitude and exchange occurs preferentially from the N-terminal end. At both pH values the NH proton of Asp 14, which is just outside the helix, is less protected by an order of magnitude than the adjacent NH protons inside the helix. Opening of the helix can be observed below pH 3.7 by changes in chemical shifts of the NH protons in the helix. At pH 2.4 the changes are 25% of those expected for complete opening. Helix opening is a fast reaction on the n.m.r. time scale (tau much less than 1 ms) unlike the generalized unfolding of RNAase S which is a slow reaction. Dissociation of S-peptide from S-protein in native RNAase S at pH 3.0 also is a slow reaction. Opening of the helix below pH 3.7 is a two-state reaction, as judged by comparing chemical shifts with exchange rates. The exchange rates at pH 3.1 are predicted correctly from the changes in chemical shift by assuming that helix opening is a two-state reaction. At pH values above 3.7, the nature of the helix opening reaction changes. These results indicate that at least one partially unfolded state of RNAase S is populated in the low pH unfolding transition.     

Orientation of the putative recognition helix in the DNA-binding domain of Hin recombinase complexed with the hix site

On the basis of sequence similarity with other known DNA-binding proteins, the DNA-binding domain of Hin recombinase, residues 139-190, is thought to bind DNA by a helix-turn-helix motif. Two models can be considered that differ in the orientation of the recognition helix in the major groove of DNA. One is based on the orientation of the recognition helix found in the 434 repressor (1-69) and lambda repressor-DNA cocrystals, and the other is based on the NMR studies of lac repressor headpiece. Cleavage by EDTA.Fe attached to a lysine side chain (Ser183----Lys183) near the COOH terminus of Hin(139-184) reveals that the putative recognition helix is oriented toward the center of the inverted repeats in a manner similar to that seen in the 434 and lambda repressor-DNA cocrystals.     

Pressure effect on the conformational fluctuation of apomyoglobin in the native state

We have investigated the effect of pressure on fluctuations of the native state of sperm whale apomyoglobin (apoMb) by H/D exchange, fluorescence, and limited proteolysis. The results from intrinsic fluorescence showed that a large fraction of apoMb molecules is in the native conformation in the pressure range from 0.1 to 150 MPa at 293 K and pH 6.0. The H/D exchange of protons of the individual backbone amino acids in this pressure range was monitored by NMR. The rate of H/D exchange was enhanced at high pressure, with the protection factors for some residues decreasing by factors of more than 100 compared to the values at 0.1 MPa. The amplitude of the decrease of the protection factor varied among the individual amino acids on the same secondary structure unit. This result suggests that H/D exchange in apoMb is explained best by the penetration model, in which solvent penetrates into the protein matrix via small motions. The result from limited proteolysis under high pressure showed that a pressure increase does not induce local unfolding of the secondary structure units of apoMb. Conformational fluctuations much smaller than local unfolding evidently provide pathways for water to diffuse into the protein interior, and are enhanced by an increase of pressure.     

Protein hydrogen exchange mechanism: local fluctuations

Experiments were done to study the dynamic structural motions that determine protein hydrogen exchange (HX) behavior. The replacement of a solvent-exposed lysine residue with glycine (Lys8Gly) in a helix of recombinant cytochrome c does not perturb the native structure, but it entropically potentiates main-chain flexibility and thus can promote local distortional motions and large-scale unfolding. The mutation accelerates amide hydrogen exchange of the mutated residue by about 50-fold, neighboring residues in the same helix by less, and residues elsewhere in the protein not at all, except for Leu98, which registers the change in global stability. The pattern of HX changes shows that the coupled structural distortions that dominate exchange can be several residues in extent, but they expose to exchange only one amide NH at a time. This "local fluctuation" mode of hydrogen exchange may be generally recognized by disparate near-neighbor rates and a low dependence on destabilants (denaturant, temperature, pressure). In contrast, concerted unfolding reactions expose multiple neighboring amide NHs with very similar computed protection factors, and they show marked destabilant sensitivity. In both modes, ionic hydrogen exchange catalysts attack from the bulk solvent without diffusing through the protein matrix.     

31P NMR spectra of an oligodeoxyribonucleotide duplex lac operator-repressor headpiece complex

The interaction of a symmetric lac operator duplex, d(TGTGAGCGCTCACA)2, with the N-terminal 56-residue headpiece fragment of the lac repressor protein was monitored by 31P NMR spectroscopy. The changes in the 31P chemical shifts upon addition of the headpiece demonstrated an end point of two headpiece fragments per symmetric 14-mer duplex with each headpiece binding to the T1pG2pT3pG4pA5 ends of the duplex. The specific phosphate 31P perturbations observed are consistent with those residues implicated in protein binding by previous NMR, molecular biological, and biochemical techniques. Upon complexation, the 31P signals of phosphates G2-A5 showed upfield or downfield shifts (less than 0.2 ppm) while most other residues were unperturbed. The interactions were dependent on ionic strength. The 31P NMR data provide direct evidence for predominant recognition of the 5' strand of the 5'-TGTGA/3'-ACACT binding site.     

NMR study of the interaction between the lac repressor and the lac operator

Binding of the lac repressor headpiece, the N-terminal region of the lac repressor, to the lac operator of Escherichia coli was studied by 1H-NMR spectroscopy. Two DNA fragments, of 51 base pairs and 62 base pairs, containing the lac operator region, were investigated. The signals of their hydrogen-bonded imino protons were well resolved in the 500-MHz NMR spectra. The spectra of the free lac operator DNA are similar to those obtained from ring-current-shift calculations for a B-DNA structure. Complex formation with the headpiece led to small but nevertheless characteristic changes in the spectra. The fact that very few imino resonances shifted upon addition of headpiece, as well as the variety in direction and size of these chemical shifts, indicate the formation of a specific complex between the lac repressor and the lac operator. The observed changes in the resonance positions exclude the intercalation of tyrosine residues of the headpiece between adjacent base pairs of the lac operator as well as the formation of a cruciform structure. They rather reflect a small conformational transition in the DNA itself, caused for example by an alteration in the tilt of a few base pairs or a shift of the keto-enol tautomeric equilibrium of the bases towards the enolic form.     

Hierarchy of local structural and dynamics perturbations due to subdenaturing urea in the native state ensemble of DLC8 dimer

Local structural and dynamic modulations due to small environmental perturbations reflect the adaptability of the protein to different interactors. We have investigated here the preferential local perturbations in Dynein light chain protein (DLC8), a cargo adapter, by sub-denaturing urea concentrations. Equilibrium unfolding experiments by optical spectroscopic methods indicated a two state like unfolding of DLC8 dimer, with the transition mid-point occurring around 8.6 M urea. NMR studies identified the β3 and β4 strands, N-, C- terminal regions, loops connecting β1 to α1, α1 to α2 and β3 to β4 as the soft targets of urea perturbation and thus indicated potential unfolding initiation sites. Native-state hydrogen exchange studies suggested the unfolding to traverse from the edges towards the centre of the secondary structural elements. At 6 M urea the whole protein chain acts like a cooperative unit. These observations are expected to have important implications for the protein's multiple functions.     

Assignment of the 1H-NMR spectrum of a lac repressor headpiece-operator complex in H2O and identification of NOEs. Consequences for protein-DNA interaction

A complex between the headpiece amino-terminal residues 1-56 of lac repressor (HP56) and an 11-bp lac operator fragment was studied by 1H NMR. The sequence specific assignment of the exchangeable and non-exchangeable protons has been accomplished. Several protons have favourable chemical shifts in the complex, therefore new intraprotein NOEs could be found that had not been unambigously identified in the free protein. By comparison, most of these intraprotein NOEs are also present in the spectra of the free headpiece but some are different. Furthermore, several new proteins DNA NOEs could be identified. The NOE between the side-chain amide protons of Gln18 and C5H of C7 confirms the specific contact between these residues which was proposed from genetic experiments [Ebright, R. M. (1985) J. Biomol. Struct. & Dyn. 3, 281-297]. The implications of the new data for the interaction between the lac repressor headpiece and its operator are discussed.     

Photo-CIDNP study of the interaction between lac repressor headpiece and lac operator DNA

Lac repressor headpiece (HP) and intact lac repressor have been studied using the photo-CIDNP method. At neutral pH histidine 29, tyrosines 7, 12 and 17 and methionine 1 are polarised. His-29 polarizations are weaker and broader in HP59 than in HP51 indicating that the C-terminal octapeptide in HP59 adopts a conformation that allows an interaction with His-29. The photo-CIDNP spectra of intact lac repressor and HP51 are very similar, showing that the same residues are accessible to the photo-excited flavin. An equimolar mixture of HP51 and a 14 base pair lac operator fragment strongly suppresses the photo-CIDNP effect of tyrosines 7 and 17 and abolishes the His-29 polarizations. The results are compared with earlier photo-CIDNP measurements on a complex of headpiece with poly[d(AT)] and with a model derived from a 2D NMR study on a lac headpiece-operator complex.     

Determination of ultrafast protein folding rates from loop formation dynamics

Quenching of the triplet state of tryptophan by contact with cysteine can be used to measure the kinetics of loop formation in unfolded proteins. Here we show that cysteine quenching dynamics also provide a novel method for measuring folding rates when the exchange between folded and unfolded states is faster than the unquenched triplet lifetime (approximately 100 micros). We use this technique to investigate folding/unfolding kinetics of the 35 residue headpiece subdomain of the protein villin, which contains a single tryptophan residue and was engineered to contain a cysteine residue at the N terminus. At intermediate concentrations of denaturant the time-course of the triplet decay consists of two relaxations, the rates and amplitudes of which reveal the fast kinetics for folding and unfolding of this protein. The folding rates extracted using a simple kinetic model are close to those reported previously from laser-induced temperature-jump experiments that employ the change in tryptophan fluorescence as a probe. However, the results differ significantly from those reported from dynamic NMR line shape analysis on a variant with methionine at the N terminus, an issue that remains to be resolved. The analysis of the triplet quenching kinetics also shows that the quenching rates in the unfolded state increase with decreasing denaturant concentration, indicating a compaction of the unfolded protein.