|Home | About | Journals | Submit | Contact Us | Français|
MMPs and TACE (ADAM-17) assume independent, parallel or opposite pathological roles in cancer, arthritis, and several other diseases. For therapeutic purposes, selective inhibition of individual MMPs and TACE is required in most cases due to distinct roles in diseases and the need to preserve activities in normal states. Toward this goal, we compared force-field interaction energies of five ubiquitous inhibitor atoms with flexible binding sites of 24 known human MMPs and TACE. The results indicate that MMPs 1−3, 10, 11, 13, 16, and 17 have at least one subsite very similar to TACE. S3 subsite is the best target for development of specific TACE inhibitors. Specific binding to TACE compared to most MMPs is promoted by placing a negatively charged ligand part at the bottom of S2 subsite, at the entrance of S1’ subsite, or the part of S3’ subsite that is close to catalytic zinc. Numerous other clues, consistent with available experimental data, are provided for design of selective inhibitors.
The members of the metzincin superfamily of zinc-dependent endopeptidases exhibit similar structural topology of the hexa-peptide catalytic sites due to the need to accommodate the peptide backbone and the presence of conserved zinc-binding motif HEbxHxbGbxHz (b indicates a bulky hydrophobic amino acid, x is a variable amino acid, and z is a family specific amino acid). The motif comprises three His residues that bind zinc and Glu forming H-bond to the water molecule coordinated by zinc. Three families belonging to the metzincin superfamily are expressed in humans: MMPs, disintegrins-metalloproteinases (ADAMs and ADAM-TS) and astacins.1 Metzincins, synthesized as proenzymes, are mostly secreted and anchored to extracellular matrix (ECM) or to the cell surface. This compartmentalization limits their activity radius to specific substrates within the pericellular space.2 After proteolytic activation, metzincins cleave ECM components and other proteins, including cytokines, bioactive peptides, receptors, and ligands.3 Membrane-anchored metzincins exhibit sheddase activities, which are most notably manifested in activation of the membrane-bound proforms of bioactive peptides. Proteolytic activities of metzincins are controlled by endogenous tissue inhibitors and, in plasma, by α2-macroglobulin. Despite shape similarity in the binding sites,4 all three human metzincin families have specific roles under both normal and pathological conditions.
The current 24-member MMP family exhibits cell-type dependent expression levels that are generally low, yet prone to a rapid increase when tissues undergo remodeling, such as in inflammation, wound healing, development, and cancer. Due to the potential for both adhesion and proteolysis, the ADAM family, presently comprising 34 members, is implicated in processes such as sperm-egg fusion, cell-cell adhesion, ectodomain shedding, myoblast fusion, and development.5 Astacins include bone morphogenetic protein-1, as well as meprins α and β.1 Meprins are found in mammalian epithelial cells in the kidneys and the intestine.6 In contrast to other metzincins, the latent homooligomeric meprin complexes can move through extracellular spaces and only become activated at sites where relevant proteases are present.7
Several MMPs and TACE (ADAM-17) are best-characterized human metzincins. Pathological roles of MMPs and TACE, due to either (i) increased or (ii) decreased activities, can be classified as (a) independent, (b) parallel, and (c) opposite. Category i.a is represented by participation of TACE in diseases such as diabetes,8 HIV cachexia,9 and sepsis.10 Deficiency of MMP-14 causing craniofacial dysmorphism, arthritis, osteopenia, dwarfism, and fibrosis of soft tissues11 can be assigned to category ii.a. Category i.b includes inflammatory diseases, in which TACE is implicated with MMPs 1 and 3 in rheumatoid arthritis,12 and with MMPs 8 and 9 in bacterial meningitis.13,14 In category c, Alzheimer's disease can be mentioned, with the activation of TACE showing a positive effect15 and that of MMP-2 being implicated in the neuropathology of the disease.16
The roles of MMPs and TACE in cancer are most complex and span all categories. MMPs are synthesized by tumor cells and frequently also by surrounding stromal cells.17 The ECM-degrading abilities of MMPs 1, 2, 3, 7, 9, 11, and 14 are important for tumor invasion into surrounding connective tissues, crossing blood vessel walls, and metastasis to distant organs.18 Cleavage of a variety of nonmatrix substrates forms the basis for a wide range of other roles of MMPs in the processes such as cell death, cell proliferation, malignant conversion, tumor promotion, and tumor-associated angiogenesis.19 For instance, an increase in activities of MMPs 7 and 9 suppresses tumor progression by generation of Angiostatin from plasminogen in plasma and subsequent reduction in tumor vascularization.20 TACE causes tumorigenesis in nude mice by shedding transforming growth factor α (TGF-α) from its transmembrane form.21 Soluble TGF-α activates the receptor for the epidermal growth factor, a tyrosine kinase family implicated in the development of tumors of epithelial origin. Inhibition of sheddase activities of ADAMs is considered the primary cause of side effects of some cancerostatic MMP inhibitors.22
For therapeutic use, selective inhibition of individual MMPs and TACE is required in most cases, although parallel inhibition of certain MMPs and TACE might occasionally be beneficial. Selectivity is needed because of distinct roles of the enzymes in a disease and, also, to preserve their functions in normal states.23 This task is complicated by shape similarity of binding sites of TACE and MMPs even though: (i) the sequence identity of the catalytic domains of TACE and MMPs is very low,24 only 12−15%; (ii) the polypeptide chain of TACE catalytic domain is significantly longer; (iii) unlike MMPs, TACE is stable in the absence of calcium;25 and (iv) some differences between MMP-2 and TACE in binding a selective MMP inhibitor have been demonstrated.26 The usual targeting of one subsite (mostly S1’ subsite that is often considered the MMP specificity pocket) may lead to limited selectivity toward TACE vs a small group of MMPs that fades away when a broader set of MMPs is analyzed.27
This study aims at a detailed comparison of binding sites of TACE and MMPs using force-field interaction energies of atomic probes in binding sites, which are treated as partially flexible to account for induced fit. To assess the similarities in individual parts of the binding sites, linear regression analysis (LRA) was preferred to previously used Principal Component Analysis,28-30 since the latter requires discarding of the repulsive energies that are an important factor for selectivity. The most pronounced differences between the binding sites were analyzed to estimate the properties of functional groups conferring selectivity.
The hexapeptide binding site of TACE was compared with those of 24 human MMPs in terms of their force-field interaction energies with probes representing the atoms most frequently occurring in inhibitors and substrates.
Catalytic domains of 24 human MMPs were represented by experimental structures and comparative models as previously described.31 For TACE, the Protein Data Bank32 file 1BKC was used. The residues of five subsites (S1, S1’, S2, S2’, and S3’) in TACE were identified based on the positions of catalytic zinc and coordinating His residues (His405, His409, His415), the bound ligand, and the peptide backbone structure in this part of the binding site (Figure 1). S2 subsite33 was identified as the loop between His409 and His415 residues that are coordinating catalytic zinc (Figure 1A). The third coordinating His405 residue and adjacent conserved Glu406 residue are at the entrance into S1’ subsite, which extends along the α-helix from the Glu406 and His405 residues down to Leu401 residue forming the bottom of S1’ subsite. The bottom part of S1’ subsite was identified from the relative positions of Leu401 in TACE and Tyr214 forming the bottom33 of S1’ subsite in MMP-7 (PDB file 1MMQ). Bound ligand is positioned between the catalytic zinc coordinated by His residues and a characteristic series of β-sheets. The β-sheet on the surface including Leu348 is forming one side of S3’ subsite, and extends34 to S1 subsite (Figure 1C). The stretch of residues outlining the other side of S3’ subsite, including Ala439, lies on the opposite side of the bound ligand and continues to S2’ subsite34 (Figure 1D). The opposing side chains of Ala439 and Leu348 are separating S3’ subsite from S1’ subsite.25 The structure of entire catalytic domain with highlighted residues forming the five subsites is shown in Figure 1E.
The superposition of TACE to MMPs was based on the structural similarity of their binding sites. The catalytic domain of MMP-7 showed the highest sequence identity with TACE, therefore one of the X-ray structures of MMP-7 (PDB file 1MMQ) was selected as a reference structure for superposition with TACE, using the C, N, and O backbone atoms. Remaining MMPs were superimposed on MMP-7 as described previously.31
Superimposed TACE and MMP-7 structures exhibit an almost perfect overlap of the backbones in individual subsites, except S3 subsite (Figure 2). However, the residues forming individual subsites have considerable effect on the overall shapes of the binding sites. As indicated by the different backbone conformations (Figure 2), the major differences will be found in the shapes of S3 subsites. The side chains of some of the amino acids from adjacent S1 subsite also contribute to these differences. In MMP structures, one of the residues outlining S1 subsite is conserved His that coordinates the structural zinc atom and is well embedded in the structure of the protein. This conformation creates only a weak steric barrier between S1 and S3 subsites of MMPs, opening the possibility of a conformational change of bound inhibitor from the bottom of S1 to S3 subsite and vice versa. In TACE, the same position in S1 subsite is occupied by Leu350 residue whose side chain is not incorporated in the protein surface and points toward the solvent. This extended conformation, along with the backbone of S3 subsite, separates S1 and S3 subsites of the TACE binding site making it less likely, even for very flexible inhibitors, to switch the bound conformation from the S1/S3 region (Figure 2) to S3 subsite. The Tyr352 residue at the bottom of S1 subsite in TACE has more open conformation than the corresponding residues in MMP structures (e.g. Phe in MMP-8), creating a larger cavity that is extending from S1 subsite toward S3 subsite (S1/S3 in Figure 2). Significant difference can be found in the shape of S2 subsite, which is broader and deeper in MMPs as compared to TACE. The side chain of Val353 in TACE is pointing toward the S2 subsite, which is making it shallower than in MMP-8, where the side chain of Gln residue is oriented toward the solvent forming a deeper S2 subsite. Another interesting structural difference of TACE is that S1’ and S3’ subsites, separated on the surface by the side chains of residues Ala439 and Leu348, are merged below the surface.25 This feature will be discussed in more detail in section Specificity Determinants and Comparison with Inhibitor Data.
Binding was examined for five probes representing the atom types most frequently occurring in the peptidic substrates and nonpeptidic inhibitors: negatively charged carbonyl oxygen, neutral sp3 carbon and sp3 oxygen, and positively charged sp3 nitrogen and hydrogen. For the calculation of interaction energies, the superimposed binding sites were enclosed in an irregular lattice of grid points taking actual contours of the binding sites. The lattice originated from a regularly spaced (1.5 Å) grid with the dimensions 18.0, 19.5 and 19.5 Å in x-, y-, and z-directions by dropping the grid points that were not involved in ligand binding inside the binding sites. The probe was placed in each of the grid points, and the Tripos force-field interaction energy35,36 was calculated for individual partially flexible binding sites.31 The interaction energies were compared with two goals: (i) to identify the MMPs and their subsites that are most similar to TACE and (ii) to delineate the properties of the regions responsible for specific interactions of TACE with individual probes.
The relatedness between parts of the binding sites of TACE and MMPs was evaluated by linear regression analysis (LRA) of the interaction energies of probes in pertinent grid points with TACE and MMPs and characterized by the slope, intercept, the correlation coefficient r, and the explained variance r2. For the specific MMP/probe/subsite combinations with strong linear relationships, similarity of the parts of binding sites was examined through the squared correlation coefficient R2 for the identity line between the interaction energies in TACE and MMP: R2 = 1 – SSR/SYY, with SSR = Σ(ETACE – EMMP)2 and SYY = Σ(ETACE-ĒTACE)2; both summations go through pertinent grid points; ETACE and EMMP are interaction energies of the probe with TACE and MMP, respectively; and ĒTACE is the average of interaction energies between the probe in the used grid points and TACE.
For the second goal, to define the regions of the binding sites, which are important for design of specific TACE inhibitors, the differences between interaction energies for TACE and all MMPs for each grid point were calculated and the grid points where the differences exceeded defined levels were considered responsible for specific interactions with ligands.
For all TACE/MMP pairs, the relationships between the interaction energies were examined for each probe or each subsite using the averages of explained variance r2 from LRA. The negatively charged probe was found to be the least distinctive (on average, 48% of the MMP data could be explained by TACE energies), while the neutral and positive probes were showing more or less similar levels of discrimination between TACE and MMPs (the average r2 was 0.35 for positive probes and 0.32 for neutral probes). For individual subsites, significantly higher relatedness, calculated as average r2 for all TACE/MMP pairs with combined probes, was found in subsites S3’, S1’ and S2 (average r2 = 0.47, 0.46, and 0.42, respectively) than in subsites S2’ and S1 (average r2 = 0.34 and 0.28, respectively), which were followed by the most different subsite S3 (average r2 = 0.17).
To compare the binding sites of TACE and MMPs, pairwise correlations were performed in two phases. First, relatedness characterized by the strength of the linear relationships between the interaction energies of the probe with TACE and MMP was evaluated by LRA. The explained variances (r2) of these relationships (Figure 3) can serve as a guide for finding the most distinct parts of the binding sites, which have weak relationships between the two sets of interactions energies. Such different spots are to be targeted for design of inhibitors specific for TACE and, if necessary, for a selected group of MMPs involved in the same disease state, such as MMPs 1 and 3 in rheumatoid arthritis.12
For the related TACE/MMP pairs with correlated interaction energies (r2 ≥ 0.7), similarity of subsites was evaluated by calculating the correlation coefficient R2 for the identity line between TACE/MMP interaction energies. The results led to identification of a number of MMP/probe/subsite combinations showing significant similarity to TACE (R2 > 0.7): (i) for MMP-1, negative oxygen probe in S1’, neutral sp3 oxygen probe in S2’ and S3’; (ii) for MMP-3, negative oxygen probe in S1’ and neutral sp3 oxygen probe in S3’; (iii) for MMP-10, negative oxygen probe in S3’ and positive hydrogen in S3’; (iv) for MMP-16, positive probes in S3’; (v) for MMPs 2, 11, 13, and 17, negative probe in S1’; and (vi) for MMP-8 neutral sp3 oxygen in S3’. For MMPs with multiple similar probe/subsite selections (MMPs 1, 3, 10, and 16), possible probe/subsite combinations have also been examined and all of them were similar as well.
Many related MMP/probe/subsite combinations with strong linear relationships of interaction energies to those of TACE (explained variance r2 > 0.7, Figure 3) are not significantly similar to their TACE counterparts (R2 < 0.7). This discrepancy is caused by the LRA-fitted line having either (i) the slope close to unity and a significant nonzero intercept or (ii) the slope different from unity. For electrostatic interactions, the first case is observed when the entire subsite has a higher polarity in TACE than in MMP, and in the second case polarity is continually changing along the subsite, i.e., the region in TACE is, compared to MMP, more positive on one side and more negative on the other side or vice versa.
This behavior is most frequently observed in S3’ subsite, where a number of TACE/MMP pairs exhibit strong linear relationships between the interaction energies with positive probes (hydrogen and sp3 nitrogen). However, the LRA-fitted lines often deviate from the identity line and their R2 values are low. A closer look at the data shows that the points could be split into two subsets fitted by different lines, one maintaining the slope close to unity and the other having a significantly lower slope (Figure 4). The spatial positions of the grid points associated with these two lines are different: the points defining the line with unity slope and a significant intercept are located in the part of S3’ subsite that is close, within 4 Å, to catalytic zinc. This result is in agreement with the analyses (see below) showing that this part of S3’ subsite is more polar in TACE than in MMPs. As pointed out previously for several MMPs,26 the second shell residues surrounding the catalytic zinc in TACE comprise three charged residues Glu406, Asp416, and Glu414, while MMPs 1, 3, 7−13, 17−23, and 25−28 contain one conserved Glu residue, and MMPs 2, 14−16, and 24 contain two Glu residues in the same region. As the probe moves further away from the catalytic zinc, polarity of the TACE S3’ subsite changes compared to that of MMPs, which results in the change of the slope. Despite the two outliers from the 4 Å region lying on the line representing the grid points farther from catalytic zinc (Figure 4), the described trends are obvious.
The relatedness r2 (Figure 3) of individual subsites, in combination with similarity R2 (see above), can serve as a guide for selection of the group of MMPs that should be used to assess specificity of novel TACE inhibitors. If a selective inhibition of TACE versus more similar MMPs is achieved, a good selectivity with respect to the rest of the MMPs can also be expected.
TACE and MMP binding sites were analyzed for differences that could aid in design of specific TACE inhibitors and substrates. For individual grid points, the differences between the interaction energies of the probe with MMP and TACE were calculated. The parts of the binding site where the energy differences were higher than a certain level, given as the fraction of the maximum interaction energy for TACE, were considered the points of specific interactions. When is the difference in interaction energies significant and will lead to specific binding? Our previous analysis of binding sites of MMPs showed that in order to obtain a good correlation between the experimental inhibitor or substrate data on two MMPs, the relatedness between the interaction energies of two binding sites had to be at least 70%.31 Here we detected that relatedness r2 ≤ 0.3 corresponds to regions that are showing clear preference for either MMP or TACE. The results for the difference levels 70% and higher are summarized in Figure 5.
For this analysis, four subsites (S1’, S2, S3, S3’) were split into two halves differing in the properties in TACE and MMPs. The parts marked as T (from top) are the halves of the deep pockets S1’ and S3 lying closer to the protein surface or the halves of the shallow cavities S2 and S3’ that are closer to catalytic zinc. The parts marked as B (from bottom) are the bottom parts of the deep pockets S1’ and S3 protruding into the interior of the protein or the parts of the shallow cavities S2 and S3’ lying farther away from the catalytic zinc. The exact positions of the T and B parts of the subsites are outlined in Figure 2. These subdivisions are in accord with the results of the LRA analyses that showed that the differences in properties of the subsite between TACE and MMP might not apply to the entire subsite but rather to smaller portions of the subsite (Figure 4).
The best target for design of specific TACE inhibitors is S3 subsite, as also indicated by simple superposition of the structures (Figure 2) and by LRA. In S3T part, the positively charged hydrogen and nitrogen probes have much stronger attractive interactions with TACE than with each of the MMPs. The result is corroborated by the much stronger attractive interactions of the negative oxygen probe with MMPs than with TACE in the same parts of the binding site. How are these differences associated with structures of TACE and MMPs? In the part of S3 subsite that continues from S1 subsite (S1/S3, Figure 2), the environment in MMPs is affected by the presence of structural zinc and a coordinating His residue; while in TACE, the Glu319 residue is playing the dominant role. The differences in the S3T part (Figure 2) are caused by the Asp344 residue in TACE and Pro in MMPs. On the other hand, the S3B part exhibits opposite tendencies: preference for the negative probe by TACE caused by the Gln341 residue and for positive probes by a number of MMPs. A common feature in both T and B parts of S3 subsite is the preference for neutral probe by TACE over a number of MMPs. However, as shown in the superposition of the structures (Figure 2), the main differences determining specific interactions with TACE or MMP are the shapes of S3 subsite.
Besides S3 subsite, the analysis of differences in interaction energies uncovered other parts of the binding site that could be targeted in the design of specific TACE inhibitors. Placement of the negatively charged functional group in part S3’T leads to a much stronger binding of ligand to TACE than to the majority of MMPs. As the ligand extends further toward the S3’B part, away from catalytic zinc, the charge should change to neutral or positive to preserve more attractive interactions with TACE over MMPs. The differences in the part S3’B are not as pronounced as in the immediate vicinity of the catalytic zinc. For most MMPs, the difference in attraction of the negative probe to MMP over TACE was only at 60% or lower levels (data not shown). The results obtained on the S3’T part of subsite are in a close agreement with the recently published experimental observation26 that the active site of TACE is more polar than that of MMP-2 as well as with the analysis of the residues surrounding catalytic zinc where TACE contains more charged residues than MMPs (see above).
Significant differences between the interaction of probes with TACE and MMPs are found in S2B part. To obtain preferential binding to TACE, a negatively charged functional group is required for MMPs 2, 3, and 7−18; in contrast, a positively charged substructure is needed for MMPs 20 and 23−25. The S2T part shows only mild preferences for the neutral probe in TACE over most MMPs.
To improve selectivity for TACE compared to specific MMPs, the ligand parts interacting with the S1’B part should be: (i) negatively charged for MMPs 7, 10 and 16; (ii) neutral or negatively charged for MMPs 19, 23, and 25; and (iii) neutral or positively charged for MMPs 20, 26, and 27. The negatively charged functional group in the S1’T part is preferred by TACE over all MMPs. This result is in agreement with several studies. Significant differences in polarity of binding sites, especially in the vicinity of catalytic zinc, have been shown.26 In the series of benzothiadiazepine inhibitors,37 the placement of the NH group in the S1’T part caused significant decrease in inhibition activity against TACE. In the series of peptidic inhibitors,38 the highest overall selectivity toward TACE, as compared to MMPs 1, 3 and 9, has been reported for compounds with thiophene ring in the S1’T part. The carbon atoms in the thiophene ring that are pointing toward the S1’ subsite have more negative partial charges39 than the carbons of the other functional groups, cyclohexane and isobutane, interacting with the S1’T part.
A comparison of interaction energies in S1 subsite indicates that, to obtain better ligand binding to TACE, a positive ligand part is required for MMPs 19−28, while for MMPs 1−10, 13, 14, 17, and 18 the part should be neutral or negative. The study on benzothiadiazepine inhibitors37 supports these results. The inhibitors with the ester or sulfonamide group in S1 subsite showed increased selectivity against TACE and MMPs 1, 2, and 9, while those with the amino group exhibited comparable inhibitory potencies.
The comparison of crystal structures of TACE and MMPs (Figure 6) uncovers another unique structural feature of TACE that is expected to play a significant role in specificity: S3’ subsite is, below the surface, merged with the S1’ pocket.25 The two subsites are on the surface separated by the side chains of residues Ala439 and Leu348.
Owing to this structural difference and the false31 reputation of S1’ subsite as an MMP specificity pocket, a number of attempts to design of specific TACE inhibitors focused on S1’ subsite. However, some inhibitor data indicate that other parts of the binding site might be involved in the specific interactions. A set of sixteen α-sulfonamide hydroxamates with butynyl group residing in S1’ subsite and variable functional groups interacting with S1 and S2’ subsites show selectivity toward TACE over MMP-1 ranging from none to ~800 fold. Large scatter in selectivity was also observed for TACE over MMP-9 (none to 100 fold) and MMP-13 (none to ~ 30) for the same set of compounds.40 A similar result was obtained from a series of lactams with (2-methyl-4-quinolinyl)-methoxy group interacting with S1’ subsite and variable groups interacting with other parts of the binding site.41 Selectivity of inhibitors against TACE versus MMPs 1, 2, and 9 ranges from single digits to thousands for each of the enzymes. The significance of the parts of the binding site other than S1’ subsite is also demonstrated by the data on a set of piperazine analogues.42 In a subset of seven compounds with varying functional groups interacting with S1’ subsite, there was only about a 3-fold difference in specificity of TACE over MMP-1. In contrast, another subset of seven compounds with invariant S1’ functional group and modifications in different parts of the molecule exhibited about a 70-fold change in specificity. Several studies38,43,44 have shown that, for the interaction with S1’ subsite, the size of the functional group alone is insufficient for obtaining selectivity to TACE, even when compared with MMP-1 which has a small S1’ pocket. Among α-substituted hydroxamates,45 selectivity against TACE vs MMP-1 significantly increased by replacing isobutyl in S1’ subsite with larger alkoxy-benzyl substituents. According to our results (Figure 5), a neutral probe, in this case a benzene ring, is highly preferred in TACE over MMP-1 in the region where S1’ subsite connects to S3’ subsite.
Presented results are useful for design of specific inhibitors of TACE or MMPs that is complicated by apparent high shape similarity of the binding sites. When chemical nature and flexibility of the binding sites are brought into analysis, overall similarity of binding sites of TACE and MMPs decreases to the point that can lead to an impression that the majority of synthesized inhibitors will be selective either to TACE or individual MMPs. However, interactions of inhibitors with multiple subsites require that similarities of all participating subsites be taken into account. The contributions from several subsites may mutually weaken or even cancel the overall effect. As shown in Figure 5, the selectivity of an inhibitor for TACE could be increased by placing positive or neutral functional group in S1 subsite, while in S3’ subsite the same functional groups would have no effect on selectivity. The binding modes of inhibitors should be identical in the targets assessed for selectivity. The physical barrier between S1 and S3’ subsites in MMP-7 is so negligible that simple, less anchored inhibitors could switch the orientations between the subsites. The inhibitors would be selective for TACE if binding in S1 subsites of both TACE and MMP-7. However, if they prefer S3’ subsite in MMP-7, selectivity for TACE is eliminated. We believe that only a careful selection of targeted subsites and properties of the functional groups interacting with them will lead to highly specific binding.
The guidelines for design of ligands specific for TACE or for TACE and a small group of MMPs, depending on the targeted disease, cannot be generalized into few simple rules. The results summarized in Figures 3 and and55 can be used to guide synthesis of inhibitors and selection of the group of MMPs that are similar to TACE and, consequently, provide a critical evaluation of the ligand selectivity. All analyses were performed without allowing the catalytic domains to do any large-scale motions. This fact limits the applicability of the results to small inhibitors and substrates. The bound conformations of inhibitors are determined by the most attractive interactions in the binding sites and may differ in individual enzymes. A difference in binding modes in the assessed enzymes precludes the selectivity estimation.
For TACE and nine types of MMPs, X-ray and/or NMR structures have been obtained from PDB database32 as follows: TACE – 1BKC, MMP-1 - 1AYK, 1CGL, 1HFC, 2AYK, 2TCL, 3AYK, 4AYK, 966C, 1CGE, 1CGF; MMP-2 - 1HOV, 1QIB; MMP-3 - 1B3D, 1B8Y, 1BIW, 1BM6, 1BQO, 1C3I, 1CAQ, 1CIZ, 1CQR, 1D5J, 1D7X, 1D8F, 1D8M, 1G49, 1G4K, 1GO5, 1HFS, 1HY7, 1SLM, 1SLN, 1UEA, 1UMS, 1UMT, 1USN, 2SRT, 2USN, 1QIA, 1QIC; MMP-7 - 1MMP, 1MMQ, 1MMR; MMP-8 – 1A85, 1A86, 1BZS, 1I73, 1I76, 1JAN, 1JAO, 1JAP, 1JAQ, 1JHI, 1JJ9, 1KBC, 1MMB, 1MNC; MMP-9 – 1GKC; MMP-12 – 1JIZ, 1JK3; MMP-13 – 1EUB, 1FLS, 1FM1, 456C, 830C; MMP-14 – 1BQQ, 1BUV. In the NMR files (shown in italics) containing multiple structures, only the first structure was selected. Comparative models of MMPs 2 and 9−18 were obtained from Mobashery et al.46 For MMPs 19−21 and 23−28, comparative models were constructed in our group as described in detail previously.31 For protonation states of ionizable residues, the Gasteiger charges for default atom definitions as implemented in the Sybyl/Biopolymer dictionary35 were used for all structures. Carboxyl groups of Glu and Asp are ionized, with delocalized electrons leading to two equivalent oxygen atoms with charges −0.550. The Lys residue contains a terminal sp3 nitrogen forming NH3+ group with overall charge +0.816. Arg has two equivalent terminal nitrogens with slightly positive partial charges (+0.205). All other residues are treated as nonionized.
The comparative modeling procedure has been validated by construction of the comparative model for MMP-2, for which X-ray structure as well as the comparative model prepared in different group46 was available. The final root-mean-squared deviation (RMSD) values between the X-ray structure and our comparative model, RMSD=1.614, and between X-ray structure and the published comparative model, RMSD=1.564, as well as the close correlation between their interaction energies with all five probes warrants their use for MMPs where no experimental structures are available. For a number of MMPs, multiple experimental structures, X-ray and/or NMR, were available from PDB database. For these types of MMPs, a representative set of interaction energies was created by averaging the interaction energies for all similar experimental structures. The experimental structures were considered similar if the correlation of their interaction energies was R2 > 0.70. For MMPs with multiple groups of similar structures, the groups including structures originating from more different laboratories were used to create the representative set of energies. This rule was implemented in order to capture the real similarities or differences between individual types of MMPs and TACE and not the subtle differences in individual experimental structures caused by different technique used for their determination.
The alignment of catalytic domains of TACE and MMPs using sequence homology in Sybyl 6.9135 shows very low sequence identity. Superposition of TACE (PDB file 1BKC) and MMP-8 (PDB file 1A85, used previously31 as a template for superposition of MMP structures) guided by the alignment results was not satisfactory: the RMSD of the entire catalytic domain was about 12 Å, while RMSDs of the aligned catalytic domains of MMP structures were between 1 and 2 Å. For a direct comparison of the interaction energies, the TACE structure (PDB file 1BKC) had to be superimposed manually on one of the MMP structures. The catalytic domain of MMP-7 showed highest sequence identity with TACE, hence it was selected as a reference structure for superposition with TACE. The residues of five subsites (S1, S1’, S2, S2’, and S3’) in TACE were easily identified based on the positions of basic structural elements such as catalytic zinc and coordinating His residues (His405, His409, His415), the ligand bound in the catalytic domain, and the peptide backbone structure in this part of the binding site (Figure 1). Due to significant differences in residues forming the five subsites in MMP-7 and TACE, only the heavy backbone atoms were used for superposition. The RMSD for the atoms used for superposition was 1.024 Å, compared to RMSD=3.587 Å of the same set of atoms when the entire catalytic domain was superimposed onto MMP-8 as was done for all MMP structures.31 The superposition of the side chains was also satisfactory as was verified through RMSD for the side chains that are identical in catalytic sites of TACE and MMP-7. When only heavy backbone atoms of the five subsites were superimposed, RMSD = 0.906 compared to RMSD = 3.420 for the whole domain superposition. Superposition with other MMPs was satisfactory, as indicated by RMSD=1.065 for side chains identical to MMP-8 (PDB file 1A85), which was used previously31 as a template for superposition of MMP structures. It is worth noting that some of the side chains, which are identical in TACE and MMP-7 were different from the side chains identical in TACE and MMP-8.
Five probes imitating the most frequently occurring atoms in metzincin substrates and inhibitors were selected: carbonyl oxygen with a charge of −1.5, neutral sp3 oxygen and sp3 carbon, and sp3 nitrogen and hydrogen with charges of +1. The interaction energies between the TACE catalytic domain and the probes were calculated as described previously.31 Briefly, the probes were placed into each grid point of regularly spaced (1.5 Å) grid enclosing the binding site and the geometry of the complex was optimized. Overall size of the grid was 18.0, 19.5 and 19.5 Å in x-, y-, and z-directions. Due to the irregular shape of the binding site, a large portion of the grid points was omitted during calculation and the interaction energies of probes with the binding site were calculated only in the grid points pertinent to ligand binding. The coordinates of actual grid points used in calculation are listed in Supporting Information. During optimization, all residues including the backbone atoms with at least one atom within 2 Å distance from the probe were treated as flexible, and the rest of the protein and the position of the probe were fixed. For geometry optimization, the Tripos force field35,36 with Gasteiger charges47 on protein atoms and a charge of +2 on the zinc atom were used. After geometry optimization and energy calculation, the probe was removed and energy of the protein without the probe was calculated. The interaction energy of the probe with the catalytic domain was obtained as the difference of the force-field energies of the optimized protein with and without the probe.
This work was supported in part by the NIH NCRR grants 1 PP20 RR 15566 and 1 P20 RR 16471, as well as by the access to resources of the Computational Chemistry and Biology Network and the Center for High Performance Computing, both at the North Dakota State University.
Supporting Information Available: Coordinates and charges of the used structures of MMP and TACE catalytic domains, superimposed for the calculations, in the mol2 format, and coordinates of the used grid points. This material is available free of charge via the Internet at http://pubs.acs.org.