Amide proton hydrogen exchange rates for sperm whale myoglobin obtained from 15N-1H NMR spectra

The hydrogen exchange behavior of exchangeable protons in proteins can provide important information for understanding the principles of protein structure and function. The positions and exchange rates of the slowly-exchanging amide protons in sperm whale myoglobin have been mapped using 15N-1H NMR spectroscopy. The slowest-exchanging amide protons are those that are hydrogen bonded in the longest helices, including members of the B, E, and H helices. Significant protection factors were observed also in the A, C, and G helices, and for a few residues in the D and F helices. Knowledge of the identity of slowly-exchanging amide protons forms the basis for the extensive quench-flow kinetic folding experiments that have been performed for myoglobin, and gives insights into the tertiary interactions and dynamics in the protein.     

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.     

Folding subdomains of thioredoxin characterized by native-state hydrogen exchange

Native-state hydrogen exchange (HX) studies, used in conjunction with NMR spectroscopy, have been carried out on Escherichia coli thioredoxin (Trx) for characterizing two folding subdomains of the protein. The backbone amide protons of only the slowest-exchanging 24 amino acid residues, of a total of 108 amino acid residues, could be followed at pH 7. The free energy of the opening event that results in an amide hydrogen exchanging with solvent (DeltaG(op)) was determined at each of the 24 amide hydrogen sites. The values of DeltaG(op) for the amide hydrogens belonging to residues in the helices alpha(1), alpha(2), and alpha(4) are consistent with them exchanging with the solvent only when the fully unfolded state is sampled transiently under native conditions. The denaturant-dependences of the values of DeltaG(op) provide very little evidence that the protein samples partially unfolded forms, lower in energy than the unfolded state. The amide hydrogens belonging to the residues in the beta strands, which form the core of the protein, appear to have higher values of DeltaG(op) than amide hydrogens belonging to residues in the helices, suggesting that they might be more stable to exchange. This apparently higher stability to HX of the beta strands might be either because they exchange out their amide hydrogens in a high energy intermediate preceding the globally unfolded state, or, more likely, because they form residual structure in the globally unfolded state. In either case, the central beta strands-beta(3,) beta(2), and beta(4)-would appear to form a cooperatively folding subunit of the protein. The native-state HX methodology has made it possible to characterize the free energy landscape that Trx can sample under equilibrium native conditions.     

Nuclear magnetic resonance studies of the internal dynamics in Apo, (Cd2+)1 and (Ca2+)2 calbindin D9k. The rates of amide proton exchange with solvent.

The backbone dynamics of the EF-hand Ca(2+)-binding protein, calbindin D9k, has been investigated in the apo, (Cd2+)1 and (Ca2+)2 states by measuring the rate constants for amide proton exchange with solvent. 15N-1H correlation spectroscopy was utilized to follow direct 1H-->2H exchange of the slowly exchanging amide protons and to follow indirect proton exchange via saturation transfer from water to the rapidly exchanging amide protons. Plots of experimental rate constants versus intrinsic rate constants have been analyzed to give qualitative insight into the opening modes of the protein that lead to exchange. These results have been interpreted within the context of a progressive unfolding model, wherein hydrophobic interactions and metal chelation serve to anchor portions of the protein, thereby damping fluctuations and retarding amide proton exchange. The addition of Ca2+ or Cd2+ was found to retard the exchange of many amide protons observed to be in hydrogen-bonding environments in the crystal structure of the (Ca2+)2 state, but not of those amide protons that were not involved in hydrogen bonds. The largest changes in rate constant occur for residues in the ion-binding loops, with substantial effects also found for the adjacent residues in helices I, II and III, but not helix IV. The results are consistent with a reorganization of the hydrogen-bonding networks in the metal ion-binding loops, accompanied by a change in the conformation of helix IV, as metal ions are chelated. Further analysis of the results obtained for the three states of metal occupancy provides insight into the nature of the changes in conformational fluctuations induced by ion binding.     

Investigation of the structural stability of cardiotoxin analogue III from the Taiwan cobra by hydrogen-deuterium exchange kinetics.

The conformational stability of a small ( approximately 7 kDa), all beta-sheet protein, cardiotoxin analogue III (CTX III), from the venom of the Taiwan cobra has been investigated by hydrogen-deuterium (H/D) exchange using two-dimensional NMR spectroscopy. The H/D exchange kinetics of backbone amide protons in CTX III has been monitored at pD 3.6 and 6.6 (at 25 degrees C), for over 5000 h. Examination of H/D exchange kinetics in the protein showed that a number of slowly exchanging residues are in the hydrophobic core of the protein. The average protection factor of the amide protons of residues belonging to the triple-stranded beta-sheet domain is about 20 times greater than that of those in the double-stranded beta-sheet segment. The residues in the C-terminal tail of the molecule, though structureless, have been found to exhibit significant protection against H/D exchange. Comparison of the quenched-flow H/D exchange data on CTX III with those obtained in the present study reveals that the most slowly exchanging portion constitutes the folding core of the protein.     

Investigation of the structural stability of the human acidic fibroblast growth factor by hydrogen-deuterium exchange

The conformational stability of the human acidic fibroblast growth factor (hFGF-1) is investigated using amide proton exchange and temperature-dependent chemical shifts, monitored by two-dimensional NMR spectroscopy. The change in free energy of unfolding (DeltaG(u)) of hFGF-1 is estimated to be 5.00 +/- 0.09 kcal.mol(-)(1). Amide proton-exchange rates of 74 residues (in hFGF-1) have been unambiguously measured, and the exchange process occurs predominately according to the conditions of the EX2 limit. The exchange rates of the fast-exchanging amide protons exposed to the solvent have been measured using the clean SEA-HSQC technique. The amide proton protection factor and temperature coefficient estimates show reasonably good correlation. Residues in beta-strands II and VI appear to constitute the stability core of the protein. Among the 12 beta-strands constituting the beta-barrel architecture of hFGF-1, beta-strand XI, located in the heparin binding domain, exhibits the lowest average protection factor value. Amide protons involved in the putative folding nucleation site in hFGF-1, identified by quench-flow NMR studies, do not represent the slow-exchanging core. Residues in portions of hFGF-1 experiencing high conformational flexibility mostly correspond to those involved in receptor recognition and binding.     

The membrane-bound conformation of alpha-lactalbumin studied by NMR-monitored 1H exchange

The interaction of bovine alpha-lactalbumin (BLA) with negatively charged phospholipid bilayers was studied by NMR monitored 1H exchange to characterize the conformational transition that enables a water-soluble protein to associate with and partially insert into a membrane. BLA was allowed to exchange in deuterated buffer in the absence (reference) and the presence (membrane-bound) of acidic liposomes at pH 4.5, experimental conditions that allow efficient protein-membrane interaction. After adjusting the pH to 6.0, to dissociate the protein from the membrane, reference and membrane-released samples of BLA were analysed by (F1) band-selective homonuclear decoupled total correlation spectroscopy in the alphaH-NH region. The overall exchange behaviour of the membrane-bound state is molten globule-like, suggesting an overall destabilization of the polypeptide. Nevertheless, the backbone amide protons of residues R10, L12, C77, K94, K98, V99 and W104 show significant protection against solvent exchange in the membrane-bound protein. We propose a mechanism for the association of BLA with negatively charged membranes that includes initial protonation of acidic side-chains at the membrane interface, and formation of an interacting site with the membrane which involves helixes A and C. In the next step these helices would slide away from each other, adopting a parallel orientation to the membrane, and would rotate to maximize the interaction between their hydrophobic residues and the lipid bilayer.     

Periodicity of amide proton exchange rates in a coiled-coil leucine zipper peptide.

The two-stranded coiled-coil motif, which includes leucine zippers, is a simple protein structure that is well suited for studies of helix-helix interactions. The interaction between helices in a coiled coil involves packing of "knobs" into "holes", as predicted by Crick in 1953 and confirmed recently by X-ray crystallography for the GCN4 leucine zipper [O'Shea, E.K., Klemm, J.D., Kim, P.S., & Alber, T. (1991) Science 254, 539]. A striking periodicity, extending over six helical turns, is observed in the rates of hydrogen-deuterium exchange for amide protons in a peptide corresponding to the leucine zipper of GCN4. Protons at the hydrophobic interface show the most protection from exchange. The NMR chemical shifts of amide protons in the helices also show a pronounced periodicity which predicts a short H-bond followed by a long H-bond every seven residues. This variation was anticipated in 1953 by Pauling and is sufficient to give rise to a local left-handed superhelical twist characteristic of coiled coils. The amide protons that lie at the base of the "hole" in the "knobs-into-holes" packing show slow amide proton exchange rates and are predicted to have short H-bond lengths. These results suggest that tertiary interactions can lead to highly localized, but substantial, differences in stability and dynamics within a secondary structure element and emphasize the dominant nature of packing interactions in determining protein structure.     

Structure and dynamics of the potato carboxypeptidase inhibitor by 1H and 15N NMR

The solution structure and backbone dynamics of the recombinant potato carboxypeptidase inhibitor (PCI) have been characterized by NMR spectroscopy. The structure, determined on the basis of 497 NOE-derived distance constraints, is much better defined than the one reported in a previous NMR study, with an average pairwise backbone root-mean-square deviation of 0.5 A for the well-defined region of the protein, residues 7-37. Many of the side-chains show now well-defined conformations, both in the hydrophobic core and on the surface of the protein. Overall, the solution structure of free PCI is similar to the one that it shows in the crystal of the complex with carboxypeptidase A. However, some local differences are observed in regions 15-21 and 27-29. In solution, the six N-terminal and the two C-terminal residues are rather flexible, as shown by 15N backbone relaxation measurements. The flexibility of the latter segment may have implications in the binding of the inhibitor by the enzyme. All the remaining residues in the protein are essentially rigid (S2 > 0.8) with the exception of two of them at the end of a short 3/10 helix. Despite the small size of the protein, a number of amide protons are protected from exchange with solvent deuterons. The slowest exchanging protons are those in a small two-strand beta-sheet. The unfolding free energies, as calculated from the exchange rates of these protons, are around 5 kcal/mol. Other protected amide protons are located in the segment 7-12, adjacent to the beta-sheet. Although these residues are not in an extended conformation in PCI, the equivalent residues in structurally homologous proteins form a third strand of the central beta-sheet. The amide protons in the 3/10 helix are only marginally protected, indicating that they exchange by a local unfolding mechanism, which is consistent with the increase in flexibility shown by some of its residues. Backbone alignment-based programs for folding recognition, as opposite to disulfide-bond alignments, reveal new proteins of unrelated sequence and function with a similar structure.     

Deuterium/hydrogen exchange factors measured by solution nuclear magnetic resonance spectroscopy as indicators of the structure and topology of membrane proteins

Deuterium/hydrogen exchange factors (chi) were measured for the backbone amide sites of the membrane-bound forms of the 50-residue fd coat protein and the 23-residue magainin2 peptide in lipid micelles by solution nuclear magnetic resonance spectroscopy. By combining kinetic and thermodynamic effects, deuterium/hydrogen exchange factors overcome the principal limitations encountered in the measurements of kinetic protection factors and thermodynamic fractionation factors for membrane proteins. The magnitudes of the exchange factors can be correlated with the structure and topology of membrane-associated polypeptides. In fd coat protein, residues in the transmembrane helix have exchange factors that are substantially smaller than those in the amphipathic surface helix or the loop connecting the two helices. For the amphipathic helical peptide, magainin2, the exchange factors of residues exposed to the solvent are appreciably larger than those that face the hydrocarbon portion of membrane bilayers. These examples demonstrate that deuterium/hydrogen exchange factors can be measured by solution NMR spectroscopy and used to identify residues in transmembrane helices as well as to determine the polarity of amphipathic helices in membrane proteins.     

Microscopic stability of cold shock protein A examined by NMR native state hydrogen exchange as a function of urea and trimethylamine N-oxide

Native state hydrogen exchange of cold shock protein A (CspA) has been characterized as a function of the denaturant urea and of the stabilizing agent trimethylamine N-oxide (TMAO). The structure of CspA has five strands of beta-sheet. Strands beta1-beta4 have strongly protected amide protons that, based on experiments as a function of urea, exchange through a simple all-or-none global unfolding mechanism. By contrast, the protection of amide protons from strand beta5 is too weak to measure in water. Strand beta5 is hydrogen bonded to strands beta3 and beta4, both of which afford strong protection from solvent exchange. Gaussian network model (GNM) simulations, which assume that the degree of protection depends on tertiary contact density in the native structure, accurately predict the strong protection observed in strands beta1-beta4 but fail to account for the weak protection in strand beta5. The most conspicuous feature of strand beta5 is its low sequence hydrophobicity. In the presence of TMAO, there is an increase in the protection of strands beta1-beta4, and protection extends to amide protons in more hydrophilic segments of the protein, including strand beta5 and the loops connecting the beta-strands. TMAO stabilizes proteins by raising the free energy of the denatured state, due to highly unfavorable interactions between TMAO and the exposed peptide backbone. As such, the stabilizing effects of TMAO are expected to be relatively independent of sequence hydrophobicity. The present results suggest that the magnitude of solvent exchange protection depends more on solvent accessibility in the ensemble of exchange susceptible conformations than on the strength of hydrogen-bonding interactions in the native structure.     

Spectroscopic studies of bacteriorhodopsin fragments dissolved in organic solution.

Fourier transform infrared and UV fourth-derivative spectroscopies were used to study the secondary structure of bacteriorhodopsin and its chymotryptic and one of the sodium borohydride fragments dissolved in chloroform-methanol (1:1, v/v), 0.1 M LiClO4. The C1 fragment (helices C, D, E, F, and G) showed an alpha-helical content of about 53%, whereas C2 (helices A and B) had about 60%, and B2 (helices F and G) about 65% alpha-helix. The infrared main band indicated differences in alpha-helical properties between these fragments. These techniques were also used to obtain information on the interactions among helices. According to the results obtained from the hydrogen/deuterium exchange kinetics, about 40% of the amide protons of C2 are particularly protected against exchange, whereas for the C1 fragment this process is unexpectedly fast. UV fourth-derivative spectra of these samples were used to obtain information about the environment of Trp side chains. The results showed that the Trp residues of C2 are more shielded from the solvent than those of C1 or B2. The results of this work indicate that the specific interactions existing between the transmembrane segments induce different types of helical conformations in native bacteriorhodopsin.     

Solution structure and dynamics of the functional domain of Paracoccus denitrificans cytochrome c(552) in both redox states

A soluble and fully functional 10.5 kDa fragment of the 18.2 kDa membrane-bound cytochrome c(552) from Paracoccus denitrificans has been heterologously expressed and (13)C/(15)N-labeled to study the structural features of this protein in both redox states. Well-resolved solution structures of both the reduced and oxidized states have been determined using high-resolution heteronuclear NMR. The overall protein topology consists of two long terminal helices and three shorter helices surrounding the heme moiety. No significant redox-induced structural differences have been observed. (15)N relaxation rates and heteronuclear NOE values were determined at 500 and 600 MHz. Several residues located around the heme moiety display increased backbone mobility in both oxidation states, while helices I, III, and V as well as the two concatenated beta-turns between Leu30 and Arg36 apparently form a less flexible domain within the protein structure. Major redox-state-dependent differences of the internal backbone mobility on the picosecond-nanosecond time scale were not evident. Hydrogen exchange experiments demonstrated that the slow-exchanging amide proton resonances mainly belong to the helices and beta-turns, corresponding to the regions with high order parameters in the dynamics data. Despite this correlation, the backbone amide protons of the oxidized cytochrome c(552) exchange considerably faster with the solvent compared to the reduced protein. Using both differential scanning calorimetry as well as temperature-dependent NMR spectroscopy, a significant difference in the thermostabilities of the two redox states has been observed, with transition temperatures of 349.9 K (76.8 degrees C) for reduced and 307.5 K (34.4 degrees C) for oxidized cytochrome c(552). These results suggest a clearly distinct backbone stability between the two oxidation states.     

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)).     

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.     

Sequence-specific 1H and 15N resonance assignments and secondary structure of GDP-bound human c-Ha-Ras protein in solution

Summary All the backbone ¹H and ¹⁵N magnetic resonances (except for Pro residues) of the GDP-bound form of a truncated human c-Ha-ras proto-oncogene product (171 amino acid residues, the Ras protein) were assigned by ¹⁵N-edited two-dimensional NMR experiments on selectively ¹⁵N-labeled Ras proteins in combination with three-dimensional NMR experiments on the uniformly ¹⁵N-labeled protein. The sequence-specific assignments were made on the basis of the nuclear Overhauser effect (NOE) connectivities of amide protons with preceding amide and/or Cprotons. In addition to sequential NOEs, vicinal spin coupling constants for amide protons and C protons and deuterium exchange rates of amide protons were used to characterize the secondary structure of the GDP-bound Ras protein; six strands and five helices were identified and the topology of these elements was determined. The secondary structure of the Ras protein in solution was mainly consistent with that in crystal as determined by X-ray analyses. The deuterium exchange rates of amide protons were examined to elucidate the dynamic properties of the secondary structure elements of the Ras protein in solution. In solution, the -sheet structure in the Ras protein is rigid, while the second helix (A66-R73) is much more flexible, and the first and fifth helices (S17-124 and V152-L171) are more rigid than other helices. Secondary structure elements at or near the ends of the effector-region loop were found to be much more flexible in solution than in the crystalline state.     

The conformational characteristics in solution of the synthetic cyclic hexapeptide 5L-ala-D-ala

An attempt to elucidate the solution conformation(s) of the synthetic cyclic hexapeptide 5L -ala·D-ala is described. Nuclear magnetic resonance (nmr) spectra are recorded for the purpose of measuring the vicinal coupling constant between the amide and α-protons in each residue and to observe the deuterium exchange rate and temperature dependence of the chemical shift of each amide proton. Low-energy cyclic conformations, whose individual residues are in conformations consistent with the observed amide to α-proton coupling constant, are searched for in an approximate theoretical treatment. The two lowest energy, all trans peptide bond conformations generated are distinguishable by the presence or absence of a single intramolecular hydrogen bond. The observed temperature independence of the chemical shift of one of the amide protons is consistent with the presence of a single intramolecular hydrogen bond, while the observation of similar deuterium exchange rates for each of the amide protons indicates their comparable availability to solvent. Consequently, it is concluded that 5L -ala·D-ala is in rapid equilibrium between conformations with and without a single internal hydrogen bond and possesses considerable conformational flexibility in solution.     

Flexible and rigid structures in HIV-1 p17 matrix protein monitored by relaxation and amide proton exchange with NMR

The HIV-1 p17 matrix protein is a multifunctional protein that interacts with other molecules including proteins and membranes. The dynamic structure between its folded and partially unfolded states can be critical for the recognition of interacting molecules. One of the most important roles of the p17 matrix protein is its localization to the plasma membrane with the Gag polyprotein. The myristyl group attached to the N-terminus on the p17 matrix protein functions as an anchor for binding to the plasma membrane. Biochemical studies revealed that two regions are important for its function: D14–L31 and V84–V88. Here, the dynamic structures of the p17 matrix protein were studied using NMR for relaxation and amide proton exchange experiments at the physiological pH of 7.0. The results revealed that the α12-loop, which includes the 14–31 region, was relatively flexible, and that helix 4, including the 84–88 region, was the most protected helix in this protein. However, the residues in the α34-loop near helix 4 had a low order parameter and high exchange rate of amide protons, indicating high flexibility. This region is probably flexible because this loop functions as a hinge for optimizing the interactions between helices 3 and 4. The C-terminal long region of K113-Y132 adopted a disordered structure. Furthermore, the C-terminal helix 5 appeared to be slightly destabilized due to the flexible C-terminal tail based on the order parameters. Thus, the dynamic structure of the p17 matrix protein may be related to its multiple functions.     

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.     

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.