Conformational variability of matrix metalloproteinases: beyond a single 3D structure.

The structures of the catalytic domain of matrix metalloproteinase 12 in the presence of acetohydroxamic acid and N-isobutyl-N-[4-methoxyphenylsulfonyl]glycyl hydroxamic acid have been solved by x-ray diffraction in the crystalline state at 1.0 and 1.3-A resolution, respectively, and compared with the previously published x-ray structure at 1.2-A resolution of the adduct with batimastat. The structure of the N-isobutyl-N-[4-methoxyphenylsulfonyl]glycyl hydroxamic acid adduct has been solved by NMR in solution. The three x-ray structures and the solution structure are similar but not identical to one another, the differences being sizably higher in the loops. We propose that many of the loops show a dynamical behavior in solution on a variety of time scales. Different conformations of some flexible regions of the protein can be observed as "frozen" in different crystalline environments. The mobility in solution studied by NMR reveals conformational equilibria in accessible time scales, i.e., from 10(-5) s to ms and more. Averaging of some residual dipolar couplings is consistent with further motions down to 10(-9) s. Finally, local thermal motions of each frozen conformation in the crystalline state at 100 K correlate well with local motions on the picosecond time scale. Flexibility/conformational heterogeneity in crucial parts of the catalytic domain is a rule rather than an exception in matrix metalloproteinases, and its extent may be underestimated by inspection of one x-ray structure. Backbone flexibility may play a role in the difficulties encountered in the design of selective inhibitors, whereas it may be a requisite for substrate binding and broad substrate specificity.

Cross correlation rates between Curie spin and dipole-dipole relaxation in paramagnetic proteins: the case of cerium substituted calbindin D9k.

Cross correlation rates between Curie spin relaxation and H-N dipole-dipole coupling (gamma(HM,HN)CS,DD) have been determined for a calcium binding protein, Calbindin D9k, in which one of the two calcium ions is substituted with cerium(III). Gamma(HM,HN)CS,DD values depend on both the metal-to-proton distances and the M-H-N angles and can be used as an additional constraint in order to refine the solution structure of paramagnetic metalloproteins. For this purpose, we have implemented a new module (CCR-DYANA) in a version of the program DYANA (PARAMAGNETIC-DYANA), which can be used together with other paramagnetism-based constraints such as pseudocontact shifts, residual dipolar couplings and hyperfine based Karplus relationships. This integrated structure calculation protocol has the advantage that different paramagnetic-based constraints are treated by the same algorithm in a way that the efficiency of each class of constraints can be analyzed and compared.