Crystal structures of the catalytic domain of human stromelysin-1 (MMP-3) and collagenase-3 (MMP-13) with a hydroxamic acid inhibitor SM-25453

Crystal structures of the catalytic domain of human stromelysin-1 (MMP-3) and collagenase-3 (MMP-13) with a hydroxamic acid inhibitor SM-25453 have been solved at 2.01 and 2.37A resolutions, respectively. The results revealed that the binding modes for this inhibitor to MMP-3 and -13 were quite similar. However, subtle comparative differences were observed at the bottom of S1' pockets, which were occupied with the guanidinomethyl moiety of the inhibitor. A remarkable feature of the inhibitor was the deep penetration of its long aliphatic chain into the S1' pocket and exposure of the guanidinomethyl moiety to the solvent.     

Crystal structures of MMP-1 and -13 reveal the structural basis for selectivity of collagenase inhibitors

The X-ray crystal structures of the catalytic domain of human collagenase-3 (MMP-13) and collagenase-1 (MMP-1) with bound inhibitors provides a basis for understanding the selectivity profile of a novel series of matrix metalloprotease (MMP) inhibitors. Differences in the relative size and shape of the MMP S1' pockets suggest that this pocket is a critical determinant of MMP inhibitor selectivity. The collagenase-3 S1' pocket is long and open, easily accommodating large P1' groups, such as diphenylether. In contrast, the collagenase-1 S1' pocket must undergo a conformational change to accommodate comparable P1' groups. The selectivity of the diphenylether series of inhibitors for collagenase-3 is largely determined by their affinity for the preformed S1' pocket of collagenase-3, as compared to the induced fit in collagenase-1.     

High-resolution solution structure of the catalytic fragment of human collagenase-3 (MMP-13) complexed with a hydroxamic acid inhibitor

The high-resolution solution structure of the catalytic fragment of human collagenase-3 (MMP-13) complexed with a sulfonamide derivative of a hydroxamic acid compound (WAY-151693) has been determined by multidimensional heteronuclear NMR. A total of 30 structures were calculated for residues 7-164 by means of hybrid distance geometry-simulated annealing using a total of 3280 experimental NMR restraints. The atomic rms distribution about the mean coordinate positions for the 30 structures is 0.43(+/-0.05) A for the backbone atoms, 0.80(+/-0.09) A for all atoms, and 0.47(+/-0.04) A for all atoms excluding disordered side-chains. The overall structure of MMP-13 is composed of a beta-sheet consisting of five beta-strands in a mixed parallel and anti-parallel arrangement and three alpha-helices where its overall fold is consistent with previously solved MMP structures. A comparison of the NMR structure of MMP-13 with the published 1.6 A resolution X-ray structure indicates that the major differences between the structures is associated with loop dynamics and crystal-packing interactions. The side-chains of some active-site residues for the NMR and X-ray structures of MMP-13 adopt distinct conformations. This is attributed to the presence of unique inhibitors in the two structures that encounter distinct interactions with MMP-13. The major structural difference observed between the MMP-13 and MMP-1 NMR structures is the relative size and shape of the S1' pocket where this pocket is significantly longer for MMP-13, nearly reaching the surface of the protein. Additionally, MMP-1 and MMP-13 exhibit different dynamic properties for the active-site loop and the structural Zn-binding region. The inhibitor WAY-151693 is well defined in the MMP-13 active-site based on a total of 52 distance restraints. The binding motif of WAY-151693 in the MMP-13 complex is consistent with our previously reported MMP-1:CGS-27023A NMR structure and is similar to the MMP-13: RS-130830 X-ray structure.     

Efficient conformational sampling of local side-chain flexibility

Side-chain flexibility of ligand-binding sites needs to be considered in the rational design of novel inhibitors. We have developed a method to generate conformational ensembles that efficiently sample local side-chain flexibility from a single crystal structure. The rotamer-based approach is tested here for the S1' pocket of human collagenase-1 (MMP-1), which is known to undergo conformational changes in multiple side-chains upon binding of certain inhibitors. First, a raw ensemble consisting of a large number of conformers of the S1' pocket was generated using an exhaustive search of rotamer combinations on a template crystal structure. A combination of principal component analysis and fuzzy clustering was then employed to successfully identify a core ensemble consisting of a low number of representatives from the raw ensemble. The core ensemble contained geometrically diverse conformers of stable nature, as indicated in several cases by a relative energy lower than that of the minimised template crystal structure. Through comparisons with X-ray crystallography and NMR structural data we show that the core ensemble occupied a conformational space similar to that observed under experimental conditions. The synthetic inhibitor RS-104966 is known to induce a conformational change in the side-chains of the S1' pocket of MMP-1 and could not be docked in the template crystal structure. However, the experimental binding mode was reproduced successfully using members of the core ensemble as the docking target, establishing the usefulness of the method in drug design.     

Flexibility and variability of TIMP binding: X-ray structure of the complex between collagenase-3/MMP-13 and TIMP-2

The excessive activity of matrix metalloproteinases (MMPs) contributes to pathological processes such as arthritis, tumor growth and metastasis if not balanced by the tissue inhibitors of metalloproteinases (TIMPs). In arthritis, the destruction of fibrillar (type II) collagen is one of the hallmarks, with MMP-1 (collagenase-1) and MMP-13 (collagenase-3) being identified as key players in arthritic cartilage. MMP-13, furthermore, has been found in highly metastatic tumors. We have solved the 2.0 A crystal structure of the complex between the catalytic domain of human MMP-13 (cdMMP-13) and bovine TIMP-2. The overall structure resembles our previously determined MT1-MMP/TIMP-2 complex, in that the wedge-shaped TIMP-2 inserts with its edge into the entire MMP-13 active site cleft. However, the inhibitor is, according to a relative rotation of approximately 20 degrees, oriented differently relative to the proteinase. Upon TIMP binding, the catalytic zinc, the zinc-ligating side chains, the enclosing MMP loop and the S1' wall-forming segment move significantly and in concert relative to the rest of the cognate MMP, and the active site cleft constricts slightly, probably allowing a more favourable interaction between the Cys1(TIMP) alpha-amino group of the inhibitor and the catalytic zinc ion of the enzyme. Thus, this structure supports the view that the central N-terminal TIMP segment essentially defines the relative positioning of the TIMP, while the flanking edge loops determine the relative orientation, depending on the individual target MMP.     

Binding affinities for sulfonamide inhibitors with matrix metalloproteinase-2 using a linear response method

Due to their involvement in many pathological conditions, matrix metalloproteinases (MMPs), are very attractive therapeutic targets. Our study focuses on one of them, MMP-2, which is involved in tumor progression and metastasis. Recently, the solution structure of the catalytic domain of MMP-2 complexed with a hydroxamic acid inhibitor (SC-74020) was published by Feng et al. Using the Hanessian group published binding affinity data and the structure published by Feng as a basis, we have built a binding affinity model by targeting the S(2)' pocket of the enzyme with a set of nine alpha-N-sulfonylamino hydroxamic acid derivatives. Two binding geometries of each ligand have been generated corresponding to two binding modes denoted A and B, respectively, of which the first one is targeting the S(2)' pocket and the second one the S(1) pocket. For the binding affinity model developed for mode A the computed activities show a rmsd of 0.583 kcal/mol as compared with the experimental data, and a correlation coefficient r(2) of 0.779, while in the case of the binding mode B a rmsd of 0.834 kcal/mol and correlation coefficient r(2) of 0.500, respectively, were obtained. In conclusion, our data suggest a higher probability for the Phe(76) gated S(2)' open form pocket to accommodate the substituent alpha versus the wide solvent exposed S(1) subsite, probability which some research groups could have overlooked due to extensive use in their calculations of non revealing S(2)' pocket open state crystallographic structures instead of NMR ones.     

Bioactive conformation of a potent stromelysin inhibitor determined by X-nucleus filtered and multidimensional NMR spectroscopy

The biologically active conformation of a novel, very potent, nonpeptidic stromelysin inhibitor was determined by X-nucleus filtered and multidimensional NMR spectroscopy. This bound conformer was subsequently docked into the stromelysin catalytic domain (SCD) using intermolecular distance constraints derived from NOE data. The complex showed the S1′ pocket of stromelysin to be the major site of enzyme-inhibitor interaction with other portions of the inhibitor spanning the S2′ and S1 binding sites. Theoretical predictions of SCD-inhibitor binding from molecular modeling studies were consistent with the NMR data. Comparison of modeled enzyme-inhibitor complexes for stromelysin and collagenase revealed an alternate binding mode for the inhibitor in collagenase, suggesting a similar binding interaction might also be possible for stromelysin. The NMR results, however, revealed a single SCD-inhibitor binding mode and provided a structural template for the design of more potent stromelysin inhibitors.     

Tetrahydroisoquinoline-3-carboxylate based matrix-metalloproteinase inhibitors: design,synthesis and structure-activity relationship

The design, synthesis and structure-activity relationship (SAR) of a series of nonpeptidic 2-arylsulfonyl-1,2,3,4-tetrahydro-isoquinoline-3-carboxylates and-hydroxamates as inhibitors of the matrix metalloproteinase human neutrophil collagenase (MMP-8) is described here. Based on available X-ray structures of MMP-8/inhibitor complexes, our structure-based design strategy was directed to complement major protein-ligand interaction regions mainly in the S1' hydrophobic specificity pocket close to the catalytic zinc ion. Here, the rigid 1,2,3,4-tetrahydroisoquinoline scaffold (Tic) provides ideal geometry to combine hydroxamates and carboxylates as typical zinc complexing functionalities, with a broad variety of S1' directed mono- and biaryl substituents consisting of aromatic rings perfectly accommodated within this more hydrophobic region of the MMP-8 inhibitor binding site. The effect of different S1' directed substituents, zinc-complexing groups, chirality and variations of the tetrahydroisoquinoline ring-system is investigated by systematic studies. X-ray structure analyses in combination with 3D-QSAR studies provided an additional understanding of key determinants for MMP-8 affinity in this series. The hypothetical binding mode for a typical molecule as basis for our inhibitor design was found in good agreement with a 1.7 A X-ray structure of this candidate in complex with the catalytic domain of human MMP-8. After analysis of all systematic variations, 3D-QSAR and X-ray structure analysis, novel S1' directed substituents were designed and synthesized and biologically evaluated. This finally results in inhibitors, which do not only show high biological affinity for MMP-8, but also exhibit good oral bioavailability in several animal species.     

Structure of recombinant mouse collagenase-3 (MMP-13).

The matrix metalloproteinases are crucial in the physiological and pathological degradation of the mammalian extracellular matrix, including breast tumours, and osteoarthritic cartilage. These enzymes are classified according to their matrix substrate specificity. Collagenase-3 (MMP-13) is a member of this family and preferentially cleaves type II collagen, cartilage, fibronectin and aggrecan. Collagenase-3 is normally expressed in hypertrophic chondrocytes, periosteal cells, and osteoblasts during bone development. The structure of the catalytic domain of recombinant mouse collagenase-3, complexed to the hydroxamate inhibitor (RS-113456), is reported at 2.0 A resolution. Molecular replacement and weak phasing information from a single derivative determined the structure. Neither molecular replacement nor derivative methods had a sufficient radius of convergence to yield a refinable structure. The structure illuminates the atomic zinc ion interactions with functional groups in the active site, emphasizing zinc ligation and the very voluminous hydrophobic P1' group for the inhibitor potency. The structure provides insight into the specificity of this enzyme, facilitating design of specific inhibitors to target various diseases.     

Encounter with unexpected collagenase-1 selective inhibitor: switchover of inhibitor binding pocket induced by fluorine atom

Phosphonamide-based inhibitors having trifluoromethyl moiety showed highly selective inhibition against MMP-1. A possible mechanism of the selectivity of MMP-1 inhibitors through the switchover of the binding pocket was speculated by computational calculations. As a consequence of the unexpected selectivity, the specific interaction of CF3 group of the inhibitor and Arg214 in the S1' pocket of MMP-1 conducted a low binding energy.     

Inhibition of MMP-1 and MMP-13 with phosphinic acids that exploit binding in the S2 pocket

Through the use of empirical and computational methods, phosphinate-based inhibitors of MMP-1 and MMP-13 that bind into the S2 pocket of these enzymes were designed. The synthesis and testing of 2 suggested that binding was occurring as hypothesized. Structure determination of a co-crystal of 2 bound to the catalytic domain of MMP-1 confirmed the binding mode. Substituents binding into S2, S1', S2' and S3', were optimized yielding compounds with low double-digit nM IC50's against these enzymes.     

Reduced-bond tight-binding inhibitors of HIV-1 protease. Fine tuning of the enzyme subsite specificity.

Truncation of a peptide substrate in the N-terminus and replacement of its scissile amide bond with a non-cleavable reduced bond results in a potent inhibitor of HIV-1 protease. A series of such inhibitors has been synthesized, and S2-S3' subsites of the protease binding cleft mapped. The S2 pocket requires bulky Boc or PIV groups, large aromatic Phe residues are preferred in P1 and P1' and Glu in P2'. The S3' pocket prefers Phe over small Ala or Val. Introduction of a Glu residue into the P2' position yields a tight-binding inhibitor of HIV-1 protease, Boc-Phe-[CH2-NH]-Phe-Glu-Phe-OMe, with a subnanomolar inhibition constant. The relevant peptide derived from the same amino acid sequence binds to the protease with a Ki of 110 nM, thus still demonstrating a good fit of the amino acid residues into the protease binding pockets and also the importance of the flexibility of P1-P1' linkage for proper binding. A new type of peptide bond mimetic, N-hydroxylamine -CH2-N(OH)-, has been synthesized. Binding of hydroxylamino inhibitor of HIV-1 protease is further improved with respect to reduced-bond inhibitor.     

Structural evidence that brain cyclic nucleotide phosphodiesterase is a member of the 2H phosphodiesterase superfamily

2',3'-Cyclic-nucleotide 3'-phosphodiesterase (CNP) is an enzyme abundantly present in the central nervous system of mammals and some vertebrates. In vitro, CNP specifically catalyzes the hydrolysis of 2',3'-cyclic nucleotides to produce 2'-nucleotides, but the physiologically relevant in vivo substrate remains obscure. Here, we report the medium resolution NMR structure of the catalytic domain of rat CNP with phosphate bound and describe its binding to CNP inhibitors. The structure has a bilobal arrangement of two modules, each consisting of a four-stranded beta-sheet and two alpha-helices. The beta-sheets form a large cavity containing a number of positively charged and aromatic residues. The structure is similar to those of the cyclic phosphodiesterase from Arabidopsis thaliana and the 2'-5' RNA ligase from Thermus thermophilus, placing CNP in the superfamily of 2H phosphodiesterases that contain two tetrapeptide HX(T/S)X motifs. NMR titrations of the CNP catalytic domain with inhibitors and kinetic studies of site-directed mutants reveal a protein conformational change that occurs upon binding.     

Heterocycle-based MMP inhibitors with P2' substituents

Potent and selective inhibition of matrix metalloproteinases was demonstrated for a series of sulfonamide-based hydroxamic acids. The design of the heterocyclic sulfonamides incorporates a six- or seven-member central ring with a P2' substituent that can be modified. Binding interactions of this substituent at the S2' site are believed to contribute to high inhibitory potency against stromelysin, collagenase-3 and gelatinases A and B, and to provide selectivity against collagenase-1 and matrilysin. An X-ray structure of a stromelysin inhibitor complex was obtained to provide insights into the SAR and selectivity trends observed for the series.     

The replacement of ATP by the competitive inhibitor emodin induces conformational modifications in the catalytic site of protein kinase CK2

The structure of a complex between the catalytic subunit of Zea mays CK2 and the nucleotide binding site-directed inhibitor emodin (3-methyl-1,6,8-trihydroxyanthraquinone) was solved at 2.6-A resolution. Emodin enters the nucleotide binding site of the enzyme, filling a hydrophobic pocket between the N-terminal and the C-terminal lobes, in the proximity of the site occupied by the base rings of the natural co-substrates. The interactions between the inhibitor and CK2 alpha are mainly hydrophobic. Although the C-terminal domain of the enzyme is essentially identical to the ATP-bound form, the beta-sheet in the N-terminal domain is altered by the presence of emodin. The structural data presented here highlight the flexibility of the kinase domain structure and provide information for the design of selective ATP competitive inhibitors of protein kinase CK2.     

The discovery of fluoropyridine-based inhibitors of the factor VIIa/TF complex--Part 2

The activated factor VII/tissue factor complex (FVIIa/TF) is known to play a key role in the formation of blood clots. Inhibition of this complex may lead to new antithrombotic drugs. A fluoropyridine-based series of FVIIa/TF inhibitors was discovered which utilized a diisopropylamino group for binding in the S2 and S3 binding pockets of the active site of the enzyme complex. In this series, an enhancement in binding affinity was observed by substitution at the 5-position of the hydroxybenzoic acid sidechain. An X-ray crystal structure indicates that amides at this position may increase inhibitor binding affinity through interactions with the S1'/S2' pocket.     

Structural characterization of the RNase E S1 domain and identification of its oligonucleotide-binding and dimerization interfaces

S1 domains occur in four of the major enzymes of mRNA decay in Escherichia coli: RNase E, PNPase, RNase II, and RNase G. Here, we report the structure of the S1 domain of RNase E, determined by both X-ray crystallography and NMR spectroscopy. The RNase E S1 domain adopts an OB-fold, very similar to that found with PNPase and the major cold shock proteins, in which flexible loops are appended to a well-ordered five-stranded beta-barrel core. Within the crystal lattice, the protein forms a dimer stabilized primarily by intermolecular hydrophobic packing. Consistent with this observation, light-scattering, chemical crosslinking, and NMR spectroscopic measurements confirm that the isolated RNase E S1 domain undergoes a specific monomer-dimer equilibrium in solution with a K(D) value in the millimolar range. The substitution of glycine 66 with serine dramatically destabilizes the folded structure of this domain, thereby providing an explanation for the temperature-sensitive phenotype associated with this mutation in full-length RNase E. Based on amide chemical shift perturbation mapping, the binding surface for a single-stranded DNA dodecamer (K(D)=160(+/-40)microM) was identified as a groove of positive electrostatic potential containing several exposed aromatic side-chains. This surface, which corresponds to the conserved ligand-binding cleft found in numerous OB-fold proteins, lies distal to the dimerization interface, such that two independent oligonucleotide-binding sites can exist in the dimeric form of the RNase E S1 domain. Based on these data, we propose that the S1 domain serves a dual role of dimerization to aid in the formation of the tetrameric quaternary structure of RNase E as described by Callaghan et al. in 2003 and of substrate binding to facilitate RNA hydrolysis by the adjacent catalytic domains within this multimeric enzyme.