Entropy increases from different sources support the high-affinity binding of the N-terminal inhibitory domains of tissue inhibitors of metalloproteinases to the catalytic domains of matrix metalloproteinases-1 and -3

The avid binding of tissue inhibitors of metalloproteinases (TIMPs) to matrix metalloproteinases (MMPs) is crucial for the regulation of pericellular and extracellular proteolysis. The interactions of the catalytic domain (cd) of MMP-1 with the inhibitory domains of TIMP-1 and TIMP-2 (N-TIMPs) and MMP-3cd with N-TIMP-2 have been characterized by isothermal titration calorimetry and compared with published data for the N-TIMP-1/MMP-3cd interaction. All interactions are largely driven by increases in entropy but there are significant differences in the profiles for the interactions of both N-TIMPs with MMP-1cd as compared with MMP-3cd; the enthalpy change ranges from small for MMP-1cd to highly unfavorable for MMP-3cd (-0.1 ± 0.7 versus 6.0 ± 0.5 kcal mol(-1)). The heat capacity change (ΔC(p)) of binding to MMP-1cd (temperature dependence of ΔH) is large and negative (-210 ± 20 cal K(-1) mol(-1)), indicating a large hydrophobic contribution, whereas the ΔC(p) values for the binding to MMP-3cd are much smaller (-53 ± 3 cal K(-1) mol(-1)), and some of the entropy increase may arise from increased conformational entropy. Apart from differences in ionization effects, it appears that the properties of the MMP may have a predominant influence in the thermodynamic profiles for these N-TIMP/MMP interactions.     

Temperature dependence of dynamics and thermodynamics of the regulatory domain of human cardiac troponin C

Binding of Ca(2+) to the regulatory domain of troponin C (TnC) in cardiac muscle initiates a series of protein conformational changes and modified protein-protein interactions that initiate contraction. Cardiac TnC contains two Ca(2+) binding sites, with one site being naturally defunct. Previously, binding of Ca(2+) to the functional site in the regulatory domain of TnC was shown to lead to a decrease in conformational entropy (TDeltaS) of 2 and 0.5 kcal mol(-1) for the functional and nonfunctional sites, respectively, using (15)N nuclear magnetic resonance (NMR) relaxation studies [Spyracopoulos, L., et al. (1998) Biochemistry 37, 18032-18044]. In this study, backbone dynamics of the Ca(2+)-free regulatory domain are investigated by backbone amide (15)N relaxation measurements at eight temperatures from 5 to 45 degrees C. Analysis of the relaxation measurements yields an order parameter (S(2)) indicating the degree of spatial restriction for a backbone amide H-N vector. The temperature dependence of S(2) allows estimation of the contribution to protein heat capacity from pico- to nanosecond time scale conformational fluctuations on a per residue basis. The average heat capacity contribution (C(p,j)) from backbone conformational fluctuations for regions of secondary structure for the regulatory domain of cardiac apo-TnC is 6 cal mol(-1) K(-1). The average heat capacity for Ca(2+) binding site 1 is larger than that for site 2 by 1.3 +/- 0.8 cal mol(-1) K(-1), and likely represents a mechanism where differences in affinity between Ca(2+) binding sites for EF hand proteins can be modulated.     

Global orientation of bound MMP-3 and N-TIMP-1 in solution via residual dipolar couplings

Crystal structures of catalytic domains of MMP-3 and MT1-MMP bound to TIMP-1 or TIMP-2, respectively, differ in the orientation of the TIMP in the MMP active site. The orientation in solution of N-TIMP-1 in the MMP-3 active site has been investigated using residual dipolar couplings (RDCs). Fitting of the RDCs to the X-ray structures of the complexes suggests general agreement with the orientation of crystalline MMP-3(DeltaC) and TIMP-1 and a large disparity from the orientation of crystalline MT1-MMP(DeltaC) and TIMP-2. Rigid body docking of MMP-3 and N-TIMP-1 X-ray coordinates using RDCs and intermolecular NOEs provided a time-averaged orientation in solution differing from the crystal structure by a 5 degrees rotation toward the MT1-MMP(DeltaC)/TIMP-2 orientation. The slight discrepancy in orientations in solution and crystal lies within the experimental uncertainties. Intermolecular NOEs used in the docking corroborated the accuracy of mapping the interface by a paramagnetic NMR footprinting assay, a potential alternative source of contacts for docking. Some uncertainty in the N-TIMP-1 orientation in the MMP-3 active site, coupled with microsecond to millisecond fluctuations of the MMP-binding ridge of N-TIMP-1 in the complex and flexibility in MMP-3(DeltaC) S(1)' to S(3)' subsites, leaves open the possibility that N-TIMP-1 might dynamically pivot a few degrees or more in the arc toward the MT1-MMP(DeltaC)/TIMP-2 orientation. Differing TIMP orientations in MMP active sites are more likely to result from structural differences in TIMP AB hairpin loops than from crystal packing artifacts.     

Constraining specificity in the N-domain of tissue inhibitor of metalloproteinases-1; gelatinase-selective inhibitors

The tissue inhibitors of metalloproteinases (TIMPs) are endogenous inhibitors of the matrix metalloproteinases (MMPs). Since unregulated MMP activities are linked to arthritis, cancer, and atherosclerosis, TIMP variants that are selective inhibitors of disease-related MMPs have potential therapeutic value. The structures of TIMP/MMP complexes reveal that most interactions with the MMP involve the N-terminal pentapeptide of TIMP and the C-D beta-strand connector which occupy the primed and unprimed regions of the active site. The loop between beta-strands A and B forms a secondary interaction site for some MMPs, ranging from multiple contacts in the TIMP-2/membrane type-1 (MT1)-MMP complex to none in the TIMP-1/MMP-1 complex. TIMP-1 and its inhibitory domain, N-TIMP-1, are weak inhibitors of MT1-MMP; inhibition is not improved by grafting the longer AB loop from TIMP-2 into N-TIMP-1, but this change impairs binding to MMP-3 and MMP-7. Mutational studies with N-TIMP-1 suggest that its weak inhibition of MT1-MMP, as compared to other N-TIMPs, arises from multiple (>3) sequence differences in the interaction site. Substitutions for Thr2 of N-TIMP-1 strongly influence MMP selectivity; Arg and Gly, that generally reduce MMP affinity, have less effect on binding to MMP-9. When the Arg mutation is added to the N-TIMP-1(AB2) mutant, it produces a gelatinase-specific inhibitor with Ki values of 2.8 and 0.4 nM for MMP-2 and -9, respectively. Interestingly, the Gly mutant has a Ki of 2.1 nM for MMP-9 and >40 muM for MMP-2, indicating that engineered TIMPs can discriminate between MMPs in the same subfamily.     

NMR structure of tissue inhibitor of metalloproteinases-1 implicates localized induced fit in recognition of matrix metalloproteinases

A high quality solution structure of the matrix metalloproteinase inhibitory N-terminal domain of recombinant human tissue inhibitor of metalloproteinases-1 (N-TIMP-1) has been determined. For the rigidly packed residues, the average RMSD to the mean structure is 0. 57 A for the backbone atoms and 1.00 A for all heavy atoms. Comparison of the solution structure of free N-TIMP-1 with the crystal structure of TIMP-1 bound to the catalytic domain of MMP-3 ( Gomis-R]uth et al., 1997 ) shows that the structural core of the beta barrel flanked by helices is nearly unchanged by the association with MMP-3, evident from a backbone RMSD of 1.15 A. However, clear differences in the conformation of the MMP-binding ridge of free and MMP-bound TIMP-1 suggest induced fit throughout the ridge. The MMP-dependent conformational changes in the ridge include a dramatic bending of AB loop residues Glu28 through Leu34, moderate hinge bending of the CD-loop about residues Ala65 and Cys70, and modest bending of the Cys1 through Pro6 segment. A large number of interresidue Nuclear Overhauser enhancements (NOEs) augmented by stereospecific assignments, torsion restraints, and dipolar couplings (an average of 18 non-trivial restraints per residue) engender confidence in these structural inferences. A tight cluster of three lysine residues and one arginine residue atop beta-strands A and B, and identical among TIMP sequences, form the heart of a highly conserved electropositive patch that may interact with anionic components of the extracellular matrix.     

The effect of matrix metalloproteinase complex formation on the conformational mobility of tissue inhibitor of metalloproteinases-2 (TIMP-2)

The backbone mobility of the N-terminal domain of tissue inhibitor of metalloproteinases-2 (N-TIMP-2) was determined both for the free protein and when bound to the catalytic domain of matrix metalloproteinase-3 (N-MMP-3). Regions of the protein with internal motion were identified by comparison of the T(1) and T(2) relaxation times and (1)H-(15)N nuclear Overhauser effect values for the backbone amide (15)N signals for each residue in the sequence. This analysis revealed rapid internal motion on the picosecond to nanosecond time scale for several regions of free N-TIMP-2, including the extended beta-hairpin between beta-strands A and B, which forms part of the MMP binding site. Evidence of relatively slow motion indicative of exchange between two or more local conformations on a microsecond to millisecond time scale was also found in the free protein, including two other regions of the MMP binding site (the CD and EF loops). On formation of a tight N-TIMP-2. N-MMP-3 complex, the rapid internal motion of the AB beta-hairpin was largely abolished, a change consistent with tight binding of this region to the MMP-3 catalytic domain. The extended AB beta-hairpin is not a feature of all members of the TIMP family; therefore, the binding of this highly mobile region to a site distant from the catalytic cleft of the MMPs suggests a key role in TIMP-2 binding specificity.     

Dynamic characterisation of the netrin-like domain of human type 1 procollagen C-proteinase enhancer and comparison to the N-terminal domain of tissue inhibitor of metalloproteinases (TIMP)

The backbone mobility of the C-terminal domain of procollagen C-proteinase enhancer (NTR PCOLCE1), part of a connective tissue glycoprotein, was determined using 15N NMR spectroscopy. NTR PCOLCE1 has been shown to be a netrin-like domain and adopts an OB-fold such as that found in the N-terminal domain of tissue inhibitors of metalloproteinases-1 (N-TIMP-1), N-TIMP-2, the laminin-binding domain of agrin and the C-terminal domain of complement protein C5. NMR relaxation dynamics of NTR PCOLCE1 highlight conformational flexibility in the N-terminus, strand A and the proximal CD loop. This region in N-TIMP is known to be essential for inhibitory activity against the matrix metalloproteinases and suggests that this region is of equal importance for NTR PCOLCE1, although the specific functional activity of the NTR PCOLCE1 domain is still unknown. Dynamics observed within the structural core of NTR PCOLCE1 that are not observed in N-TIMP molecules suggest that although the two domains have a similar architecture, the NTR PCOLCE1 domain will show different thermodynamic properties on binding and hence the target molecule could be somewhat different from that observed for the TIMPs. ModelFree order parameters show that NTR PCOLCE1 has more flexibility than both N-TIMP-1 and N-TIMP-2.     

Main chain and side chain dynamics of the ubiquitin conjugating enzyme variant human Mms2 in the free and ubiquitin-bound States

Protein ubiquitination involves a cascade of enzymatic steps where ubiquitin (Ub) is sequentially transferred as a thiolester intermediate from an E1 enzyme to an E2 enzyme and finally to the protein target with the help of a Ub-protein ligase. Protein ubiquitination brought about by the Ubc13-Mms2 (E2-E2) complex has a unique role in the cell, unrelated to protein degradation. The Mms2-Ubc13 heterodimer links Ub molecules to one another through an isopeptide bond between its own C-terminus and Lys-63 on another Ub. The role of Mms2 is to orient a target-bound Ub molecule such that its Lys-63 is proximal to the C-terminus of the Ub molecule that is covalently linked to the active site of Ubc13. To gain insight into the influence of protein dynamics on the affinity of Ub for Mms2, we have determined pico- to nanosecond time scale fluctuations of the main chain and methyl side chains of human Mms2 in the free and Ub-bound states using solution state (15)N and (2)H nuclear magnetic resonance relaxation measurements. Analysis of the relaxation data allows for a semiquantitative estimation of the conformational entropy change (TDeltaS(NMR)) for the main chain and side chain methyl groups of Mms2 upon binding Ub. The value of TDeltaS(NMR) for the main chain and side chain methyl groups of Mms2 is -8 +/- 2 and -2 +/- 2 kcal mol(-)(1), respectively. The experimental DeltaG(binding) for the Mms2.Ub complex is -6 kcal mol(-)(1). Estimation of DeltaG(binding) using an empirical structure-based approach that does not account for changes in main chain entropy yields a value of -17 +/- 2 kcal mol(-)(1). However, inclusion of TDeltaS(NMR) for the main chain of Mms2 increases the estimated DeltaG(binding) to -9 +/- 3 kcal mol(-)(1). Assuming that changes in Ub main chain dynamics contribute to TDeltaS(NMR) to the same extent as Mms2, the estimated DeltaG(binding) is further reduced to -1 +/- 4 kcal mol(-)(1), a value close to the experimental DeltaG(binding).     

Structure, dynamics, and thermodynamics of the structural domain of troponin C in complex with the regulatory peptide 1-40 of troponin I

The structure of the calcium-saturated C-domain of skeletal troponin C (CTnC) in complex with a regulatory peptide comprising residues 1-40 (Rp40) of troponin I (TnI) was determined using nuclear magnetic resonance (NMR) spectroscopy. The solution structure determined by NMR is similar to the structure of the C-domain from intact TnC in complex with TnI(1)(-)(47) determined by X-ray crystallography [Vassylyev, D. G., Takeda, S., Wakatsuki, S., Maeda, K., and Maeda, Y. (1998) Proc. Natl. Acad. Sci. U.S.A. 95, 4847-4852]. Changes in the dynamic properties of CTnC.2Ca2+ induced by Rp40 binding were investigated using backbone amide (15)N NMR relaxation measurements. Analysis of NMR relaxation data allows for extraction of motional order parameters on a per residue basis, from which the contribution of changes in picosecond to nanosecond time scale motions to the conformational entropy associated with complex formation can be estimated. The results indicate that binding of Rp40 decreases backbone flexibility in CTnC, particularly at the end of the C-terminal helix. The backbone conformational entropy change (-TDeltaS) associated with binding of Rp40 to CTnC.2Ca2+ determined from (15)N relaxation data is 9.6 +/- 0.7 kcal mol(-1) at 30 degrees C. However, estimation of thermodynamic quantities using a structural approach [Lavigne, P., Bagu, J. R., Boyko, R., Willard, L., Holmes, C. F., and Sykes, B. D. (2000) Protein Sci. 9, 252-264] reveals that the change in solvation entropy upon complex formation is dominant and overcomes the thermodynamic "cost" associated with "stiffening" of the protein backbone upon Rp40 binding. Additionally, backbone amide (15)N relaxation data measured at different concentrations of CTnC.2Ca2+.Rp40 reveal that the complex dimerizes in solution. Fitting of the apparent global rotational correlation time as a function of concentration to a monomer-dimer equilibrium yields a dimerization constant of approximately 8.3 mM.     

Conformational entropy changes upon lactose binding to the carbohydrate recognition domain of galectin-3

The conformational entropy of proteins can make significant contributions to the free energy of ligand binding. NMR spin relaxation enables site-specific investigation of conformational entropy, via order parameters that parameterize local reorientational fluctuations of rank-2 tensors. Here we have probed the conformational entropy of lactose binding to the carbohydrate recognition domain of galectin-3 (Gal3), a protein that plays an important role in cell growth, cell differentiation, cell cycle regulation, and apoptosis, making it a potential target for therapeutic intervention in inflammation and cancer. We used ¹⁵N spin relaxation experiments and molecular dynamics simulations to monitor the backbone amides and secondary amines of the tryptophan and arginine side chains in the ligand-free and lactose-bound states of Gal3. Overall, we observe good agreement between the experimental and computed order parameters of the ligand-free and lactose-bound states. Thus, the ¹⁵N spin relaxation data indicate that the molecular dynamics simulations provide reliable information on the conformational entropy of the binding process. The molecular dynamics simulations reveal a correlation between the simulated order parameters and residue-specific backbone entropy, re-emphasizing that order parameters provide useful estimates of local conformational entropy. The present results show that the protein backbone exhibits an increase in conformational entropy upon binding lactose, without any accompanying structural changes.     

Conformational dynamics and thermodynamics of protein-ligand binding studied by NMR relaxation

Protein conformational dynamics can be critical for ligand binding in two ways that relate to kinetics and thermodynamics respectively. First, conformational transitions between different substates can control access to the binding site (kinetics). Secondly, differences between free and ligand-bound states in their conformational fluctuations contribute to the entropy of ligand binding (thermodynamics). In the present paper, I focus on the second topic, summarizing our recent results on the role of conformational entropy in ligand binding to Gal3C (the carbohydrate-recognition domain of galectin-3). NMR relaxation experiments provide a unique probe of conformational entropy by characterizing bond-vector fluctuations at atomic resolution. By monitoring differences between the free and ligand-bound states in their backbone and side chain order parameters, we have estimated the contributions from conformational entropy to the free energy of binding. Overall, the conformational entropy of Gal3C increases upon ligand binding, thereby contributing favourably to the binding affinity. Comparisons with the results from isothermal titration calorimetry indicate that the conformational entropy is comparable in magnitude to the enthalpy of binding. Furthermore, there are significant differences in the dynamic response to binding of different ligands, despite the fact that the protein structure is virtually identical in the different protein-ligand complexes. Thus both affinity and specificity of ligand binding to Gal3C appear to depend in part on subtle differences in the conformational fluctuations that reflect the complex interplay between structure, dynamics and ligand interactions.     

Temperature dependence of the backbone dynamics of ribonuclease A in the ground state and bound to the inhibitor 5'-phosphothymidine (3'-5')pyrophosphate adenosine 3'-phosphate

The interaction of the dinucleotide inhibitor 5'-phosphothymidine(3',5')pyrophosphate adenosine 3'-phosphate (pTppAp) with bovine pancreatic ribonuclease A (RNase A) was characterized by calorimetry and solution NMR spectroscopy. Calorimetric data show that binding of pTppAp to RNase A is exothermic (DeltaH = -60.1 +/- 4.1 kJ/mol) with a dissociation constant of 16 nM at 298 K. At this temperature, the binding results in an entropy loss (TDeltaS = -16.8 +/- 7.3 kJ/mol) that is more favorable than that with the product analogue, 2'-CMP (TDeltaS = -31.3 +/- 0.9 kJ/mol). Temperature-dependent calorimetric experiments give a DeltaC(p) for ligand binding of -230 +/- 100 J/mol K. Binding of pTppAp results in noticeable effects on the backbone amide chemical shifts and dynamics. Amide backbone (15)N NMR spin-relaxation studies were performed on both apo RNase A and RNase A/pTppAp as a function of temperature. At each temperature, the model-free-determined order parameters, S(2), were significantly higher for RNase A/pTppAp than for the apo enzyme indicating a decrease in the conformational entropy of the protein upon ligand binding. Furthermore, the magnitude of this difference varies along the amino acid sequence specifically locating the entropic changes. The temperature dependence of S(2) at each residue enabled assessment of the local heat capacity changes (DeltaC(p)) from ligand binding. In an overall, average sense, DeltaC(p) for the protein backbone, determined from the NMR dynamics measurements, did not differ between apo RNase A and RNase A/pTppAp indicating that backbone dynamics contribute little to DeltaC(p) for protein-ligand interactions in this system. However, residue-by-residue comparison of the temperature-dependent change in entropy (DeltaS(B)) between free and bound forms reveals nonzero contributions to DeltaC(p) at individual sites. The balance of positive and negative changes reveals a redistribution of energetics upon binding. Furthermore, experiment and semiempirical estimates suggest that a large negative DeltaC(p) should accompany binding of pTppAp, and we conclude that this contribution must arise from factors other than amide backbone dynamics.     

High pressure NMR reveals a variety of fluctuating conformers in beta-lactoglobulin

High pressure 1H/15N two-dimensional NMR spectroscopy has been used to study conformational fluctuation in bovine beta-lactoglobulin at pH 2.0 and 36 degrees C. Pressure dependencies of 1H and 15N chemical shifts and cross-peak intensities were analyzed at more than 80 independent atom sites between 30 and 2000 bar. Unusually large and non-linear chemical shift pressure dependencies are found for residues centering in the hydrophobic core region, suggesting the existence of low-lying excited native states (N') of the protein. Measurement of 1H/15N cross-peak intensities at individual amide sites as a function of pressure suggests that unfolding events occur independently in two sides of the beta-barrel, i.e. the hydrophobic core side (betaF-H) (producing I2) and the non-core side (betaB-E) (producing I1). At 1 bar the stability is higher for the core region (DeltaG0 = 6.5(+/-2.0) kcal/mol) than for the non-core region (4.6(+/-1.3) kcal/mol), but at high pressure the stability is reversed due to a larger DeltaV value of unfolding for the core region (90.0(+/-35.2) ml/mol) than that for the non-core region (57.4(+/-14.4) ml/mol), possibly due to an uneven distribution of cavities. The DeltaG0 profile along the amino acid sequence obtained from the pressure experiment is found to coincide well with that estimated from hydrogen exchange experiments. Altogether, the high pressure NMR experiment has revealed a variety of fluctuating conformers of beta-lactoglobulin, notably N, N', I1, I2 and the totally unfolded conformer U. Fluctuation of N to I1 and I2 conformers with open barrel structures could be a common design of lipocalin family proteins which bind various hydrophobic compounds in its barrel structure.     

Solvation and the hidden thermodynamics of a zinc finger probed by nonstandard repair of a protein crevice

The classical Zn finger contains a phenylalanine at the crux of its three architectural elements: a beta-hairpin, an alpha-helix, and a Zn(2+)-binding site. Surprisingly, phenylalanine is not required for high-affinity Zn2+ binding, but instead contributes to the specification of a precise DNA-binding surface. Substitution of phenylalanine by leucine leads to a floppy but native-like structure whose Zn affinity is maintained by marked entropy-enthalpy compensation (DeltaDeltaH -8.3 kcal/mol and -TDeltaDeltaS 7.7 kcal/mol). Phenylalanine and leucine differ in shape, size, and aromaticity. To distinguish which features correlate with dynamic stability, we have investigated a nonstandard finger containing cyclohexanylalanine at this site. The structure of the nonstandard finger is similar to that of the native domain. The cyclohexanyl ring assumes a chair conformation, and conformational fluctuations characteristic of the leucine variant are damped. Although the nonstandard finger exhibits a lower affinity for Zn2+ than does the native domain (DeltaDeltaG -1.2 kcal/mol), leucine-associated perturbations in enthalpy and entropy are almost completely attenuated (DeltaDeltaH -0.7 kcal/mol and -TDeltaDeltaS -0.5 kcal/mol). Strikingly, global changes in entropy (as inferred from calorimetry) are in each case opposite in sign from changes in configurational entropy (as inferred from NMR). This seeming paradox suggests that enthalpy-entropy compensation is dominated by solvent reorganization rather than nominal molecular properties. Together, these results demonstrate that dynamic and thermodynamic perturbations correlate with formation or repair of a solvated packing defect rather than type of physical interaction (aromatic or aliphatic) within the core.     

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.     

Crystal structure of the catalytic domain of matrix metalloproteinase-1 in complex with the inhibitory domain of tissue inhibitor of metalloproteinase-1

The mammalian collagenases are a subgroup of the matrix metalloproteinases (MMPs) that are uniquely able to cleave triple helical fibrillar collagens. Collagen breakdown is an essential part of extracellular matrix turnover in key physiological processes including morphogenesis and wound healing; however, unregulated collagenolysis is linked to important diseases such as arthritis and cancer. The tissue inhibitors of metalloproteinases (TIMPs) function in controlling the activity of MMPs, including collagenases. We report here the structure of a complex of the catalytic domain of fibroblast collagenase (MMP-1) with the N-terminal inhibitory domain of human TIMP-1 (N-TIMP-1) at 2.54 A resolution. Comparison with the previously reported structure of the TIMP-1/stromelysin-1 (MMP-3) complex shows that the mechanisms of inhibition of both MMPs are generally similar, yet there are significant differences in the protein-protein interfaces in the two complexes. Specifically, the loop between beta-strands A and B of TIMP-1 makes contact with MMP-3 but not with MMP-1, and there are marked differences in the roles of individual residues in the C-D connector of TIMP-1 in binding to the two MMPs. Structural rearrangements in the bound MMPs are also strikingly different. This is the first crystallographic structure that contains the truncated N-terminal domain of a TIMP, which shows only minor differences from the corresponding region of the full-length protein. Differences in the interactions in the two TIMP-1 complexes provide a structural explanation for the results of previous mutational studies and a basis for designing new N-TIMP-1 variants with restricted specificity.     

Thermodynamics of folding and binding in an affibody:affibody complex

Affibody binding proteins are selected from phage-displayed libraries of variants of the 58 residue Z domain. Z(Taq) is an affibody originally selected as a binder to Taq DNA polymerase. The anti-Z(Taq) affibody was selected as a binder to Z(Taq) and the Z(Taq):anti-Z(Taq) complex is formed with a dissociation constant K(d)=0.1 microM. We have determined the structure of the Z(Taq):anti-Z(Taq) complex as well as the free state structures of Z(Taq) and anti-Z(Taq) using NMR. Here we complement the structural data with thermodynamic studies of Z(Taq) and anti-Z(Taq) folding and complex formation. Both affibody proteins show cooperative two-state thermal denaturation at melting temperatures T(M) approximately 56 degrees C. Z(Taq):anti-Z(Taq) complex formation at 25 degrees C in 50 mM NaCl and 20 mM phosphate buffer (pH 6.4) is enthalpy driven with DeltaH degrees (bind) = -9.0 (+/-0.1) kcal mol(-1)(.) The heat capacity change DeltaC(P) degrees (,bind)=-0.43 (+/-0.01) kcal mol(-1) K(-1) is in accordance with the predominantly non-polar character of the binding surface, as judged from calculations based on changes in accessible surface areas. A further dissection of the small binding entropy at 25 degrees C (-TDeltaS degrees (bind) = -0.6 (+/-0.1) kcal mol(-1)) suggests that a favourable desolvation of non-polar surface is almost completely balanced by unfavourable conformational entropy changes and loss of rotational and translational entropy. Such effects can therefore be limiting for strong binding also when interacting protein components are stable and homogeneously folded. The combined structure and thermodynamics data suggest that protein properties are not likely to be a serious limitation for the development of engineered binding proteins based on the Z domain.     

Tissue inhibitor of metalloproteinases-1 undergoes microsecond to millisecond motions at sites of matrix metalloproteinase-induced fit

The N-terminal, matrix metalloproteinase (MMP)-inhibitory fragment of recombinant, human tissue inhibitor of metalloproteinases (TIMP-1) exhibits varied backbone dynamics and rigidity. Most striking is the presence of chemical exchange in the MMP-binding ridge reported to undergo conformational change upon MMP binding. Conformational exchange fluctuations in microseconds to milliseconds map to the sites of MMP-induced fit at residues Val29 through Leu34 of the AB loop and to the Ala65 and Cys70 "hinges" of the CD loop of TIMP-1. Slow chemical exchange is also present at the type I turn of the EF loop at the base of the MMP-binding ridge. These functional slow motions and other fast internal motions are evident from backbone (15)N spin relaxation at 500 and 750 MHz, whether interpreted by the model-free formalism with axial diffusion anisotropy or by the reduced spectral density approach. The conformational exchange is confirmed by its deviation from the trend between R(2) and the cross-correlation rate eta. The magnetic field-dependence indicates that the chemical exchange broadening in the AB and CD loops is fast on the time-scale of chemical shift differences. The conformational exchange rates for most of these exchanging residues, which can closely approach MMP, appear to be a few thousand to several thousand per second. The slow dynamics of the TIMP-1 AB loop contrast the picosecond to nanosecond dynamics reported in the longer TIMP-2 AB loop.     

The role of backbone motions in ligand binding to the c-Src SH3 domain.

The Src homology 3 (SH3) domain of pp60(c-src) (Src) plays dual roles in signal transduction, through stabilizing the repressed form of the Src kinase and through mediating the formation of activated signaling complexes. Transition of the Src SH3 domain between a variety of binding partners during progression through the cell cycle requires adjustment of a delicate free energy balance. Although numerous structural and functional studies of SH3 have provided an in-depth understanding of structural determinants for binding, the origins of binding energy in SH3-ligand interactions are not fully understood. Considering only the protein-ligand interface, the observed favorable change in standard enthalpy (DeltaH=-9.1 kcal/mol) and unfavorable change in standard entropy (TDeltaS=-2.7 kcal/mol) upon binding the proline-rich ligand RLP2 (RALPPLPRY) are inconsistent with the predominantly hydrophobic interaction surface. To investigate possible origins of ligand binding energy, backbone dynamics of free and RLP2-bound SH3 were performed via (15)N NMR relaxation and hydrogen-deuterium (H/(2)H) exchange measurements. On the ps-ns time scale, assuming uncorrelated motions, ligand binding results in a significant reduction in backbone entropy (-1.5(+/-0.6) kcal/mol). Binding also suppresses motions on the micros-ms time scale, which may additionally contribute to an unfavorable change in entropy. A large increase in protection from H/(2)H exchange is observed upon ligand binding, providing evidence for entropy loss due to motions on longer time scales, and supporting the notion that stabilization of pre-existing conformations within a native state ensemble is a fundamental paradigm for ligand binding. Observed changes in motion on all three time scales occur at locations both near and remote from the protein-ligand interface. The propagation of ligand binding interactions across the SH3 domain has potential consequences in target selection through altering both free energy and geometry in intact Src, and suggests that looking beyond the protein-ligand interface is essential in understanding ligand binding energetics.