Structure-based approaches for drug design generally do not incorporate solvent effects and dynamic information to predict inhibitor-binding affinity because of practical limitations. The matrix metalloproteinases (MMPs) have previously been demonstrated to exhibit significant mobility in their active sites. This dynamic characteristic significantly complicates the drug design process based on static structures, which was clearly observed for a class of hydroxamic acids containing a butynyl moiety. Compound 1 was expected to be selective against MMP-1 based on predicted steric clashes between the butynyl P1' group and the S1' pocket, but the observation of complex inhibitor dynamics in the NMR structure of MMP-1:1 provides an explanation for the low nanomolar binding to MMP-1.
A protocol is described for rapidly screening small organic molecules for their ability to bind a target protein while obtaining structure-related information as part of a structure-based drug discovery and design program. The methodology takes advantage of and combines the inherent strengths of size exclusion gel chromatography, mass spectrometry, and NMR to identify bound complexes in a relatively universal high-throughput screening approach. Size exclusion gel chromatography in the spin column format provides the high-speed separation of a protein-ligand complex from free ligands. The spin column eluent is then analyzed under denaturing conditions by electrospray ionization mass spectrometry (MS) for the presence of small molecular weight compounds formerly bound to the protein. Hits identified by MS are then individually assayed by chemical shift perturbations in a 2D 1H-15N HSQC NMR spectrum to verify specific interactions of the compound with the protein and identification of the binding site on the protein. The utility of the MS/NMR assay is demonstrated with the use of the catalytic fragment of human fibroblast collagenase (MMP-1) as a target protein and the screening of a library consisting of approximately 32 000 compounds for the identification of molecules that exhibit specific binding to the RGS4 protein.
Interleukin-13 has been implicated as a key factor in asthma, allergy, atopy and inflammatory response, establishing the protein as a valuable therapeutic target. The high-resolution solution structure of human IL-13 has been determined by multidimensional NMR. The resulting structure is consistent with previous short-chain left-handed four-helix bundles, where a significant similarity in the folding topology between IL-13 and IL-4 was observed. IL-13 shares a significant overlap in biological function with IL-4, a result of the common α chain component (IL-4Rα) in their respective receptors. Based on the available structural and mutational data, an IL-13/IL-4Rα model and a sequential mechanism for forming the signaling heterodimer is proposed for IL-13.
A novel series of anthranilic acid-based inhibitors of MMP-1, MMP-9, and MMP-13 was prepared and evaluated both in vitro and in vivo. The most potent compound, 6e, has in vivo activity in a rat sponge-wrapped cartilage model.
ZipA, an essential component of cell division in Escherichia coli, interacts with the FtsZ protein at the midcell in one of the initial steps of septum formation. The high-resolution solution structure of the 144-residue C-terminal domain of E. coli ZipA (ZipA(185)(-)(328)) has been determined by multidimensional heteronuclear NMR. A total of 30 structures were calculated by means of hybrid distance geometry-simulated annealing using a total of 2758 experimental NMR restraints. The atomic root means square distribution about the mean coordinate positions for residues 6-142 for the 30 structures is 0.37 +/- 0.04 A for the backbone atoms, 0. 78 +/- 0.05 A for all atoms, and 0.45 +/- 0.04 A for all atoms excluding disordered side chains. The NMR solution structure of ZipA(185)(-)(328) is composed of three alpha-helices and a beta-sheet consisting of six antiparallel beta-strands where the alpha-helices and the beta-sheet form surfaces directly opposite each other. A C-terminal peptide from FtsZ has been shown to bind ZipA(185)(-)(328) in a hydrophobic channel formed by the beta-sheet providing insight into the ZipA-FtsZ interaction. An unexpected similarity between the ZipA(185)(-)(328) fold and the split beta-alpha-beta fold observed in many RNA binding proteins may further our understanding of the critical ZipA-FtsZ interaction.
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.
Heterotrimeric guanine nucleotide-binding proteins (G-proteins) are transducers in many cellular transmembrane signaling systems where regulators of G-protein signaling (RGS) act as attenuators of the G-protein signal cascade by binding to the Galpha subunit of G-proteins (G(i)(alpha)(1)) and increasing the rate of GTP hydrolysis. The high-resolution solution structure of free RGS4 has been determined using two-dimensional and three-dimensional heteronuclear NMR spectroscopy. A total of 30 structures were calculated by means of hybrid distance geometry-simulated annealing using a total of 2871 experimental NMR restraints. The atomic rms distribution about the mean coordinate positions for residues 5-134 for the 30 structures is 0.47 +/- 0.05 A for the backbone atoms, 0. 86 +/- 0.05 A for all atoms, and 0.56 +/- 0.04 A for all atoms excluding disordered side chains. The NMR solution structure of free RGS4 suggests a significant conformational change upon binding G(i)(alpha)(1) as evident by the backbone atomic rms difference of 1. 94 A between the free and bound forms of RGS4. The underlying cause of this structural change is a perturbation in the secondary structure elements in the vicinity of the G(i)(alpha)(1) binding site. A kink in the helix between residues K116-Y119 is more pronounced in the RGS4-G(i)(alpha)(1) X-ray structure relative to the free RGS4 NMR structure, resulting in a reorganization of the packing of the N-terminal and C-terminal helices. The presence of the helical disruption in the RGS4-G(i)(alpha)(1) X-ray structure allows for the formation of a hydrogen-bonding network within the binding pocket for G(i)(alpha)(1) on RGS4, where RGS4 residues D117, S118, and R121 interact with residue T182 from G(i)(alpha)(1). The binding pocket for G(i)(alpha)(1) on RGS4 is larger and more accessible in the free RGS4 NMR structure and does not present the preformed binding site observed in the RGS4-G(i)(alpha)(1) X-ray structure. This observation implies that the successful complex formation between RGS4 and G(i)(alpha)(1) is dependent on both the formation of the bound RGS4 conformation and the proper orientation of T182 from G(i)(alpha)(1). The observed changes for the free RGS4 NMR structure suggest a mechanism for its selectivity for the Galpha-GTP-Mg(2+) complex and a means to facilitate the GTPase cycle.
The high-resolution NMR solution structure of the catalytic fragment of human collagenase-3 (MMP-13) was used as a starting point for structure-based design of selective inhibitors for MMP-13. The major structural difference observed between the MMP 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. On the basis of the extended nature of the MMP-13 S1' pocket an inhibitor potent and selective for MMP-13 was designed from an initial high throughput screening (HTS) lead. CL-82198 was identified as a weak (10 muM) inhibitor against MMP-13 while demonstrating no activity against MMP-1, MMP-9, or the related enzyme TACE. The drug-Like properties of CL-82198 made it an ideal candidate for optimization of enzyme potency and selectivity. On the basis of NMR binding studies, it was shown that inhibitor CL-82198 bound within the entire S1' pocket of MMP-13 which is the basis of its selectivity against MMP-1, MMP-9, and TACE. A strategy utilizing this information was devised for designing new inhibitors that showed enhanced selectivity toward MMP-13. Our design strategy combined the critical selectivity features of CL-82198 with the known potency features of a nonspecific MMP inhibitor (WAY-152177) to generate a potent and selective MMP-13 inhibitor (WAY-170523). WAY-170523 has an IC50 of 17 nM for MMP-13 and showed > 5800-, 56-, and >500-fold selectivity against MMP-1 MMP-9, and TACE, respectively.
The solution structure of the catalytic fragment of human fibroblast collagenase (MMP-1) complexed with a sulfonamide derivative of a hydroxamic acid compound (CGS-27023A) has been determined using two-dimensional and three-dimensional heteronuclear NMR spectroscopy. The solution structure of the complex was calculated by means of hybrid distance geometry-simulated annealing using a combination of experimental NMR restraints obtained from the previous refinement of the inhibitor-free MMP-1 (1) and recent restraints for the MMP-1:CGS-27023A complex. The hydroxamic acid moiety of CGS-27023A was found to chelate to the "right" of the catalytic zinc where the p-methoxyphenyl sits in the S1' active-site pocket, the isopropyl group is in contact with H83 and N80, and the pyridine ring is solvent exposed. The sulfonyl oxygens are in hydrogen-bonding distance to the backbone NHs of L81 and A82. This is similar to the conformation determined by NMR of the inhibitor bound to stromelysin (2, 3). A total of 48 distance restraints were observed between MMP-1 and CGS-27023A from 3D 13C-edited/12C-filtered NOESY and 3D 15N-edited NOESY experiments. An additional 18 intramolecular restraints were observed for CGS-27023A from a 2D 12C-filtered NOESY experiment. A minimal set of NMR experiments in combination with the free MMP-1 assignments were used to assign the MMP-1 (1)H, 13C, and 15N resonances in the MMP-1:CGS-27023A complex. The assignments of CGS-27023A in the complex were obtained from 2D 12C-filtered NOESY and 2D 12C-filtered TOCSY experiments.
The high-resolution solution structure of the inhibitor-free catalytic fragment of human fibroblast collagenase (MMP-1), a protein of 18.7 kDa, which is a member of the matrix metalloproteinase family, has been determined using three-dimensional heteronuclear NMR spectroscopy. A total of 30 structures were calculated by means of hybrid distance geometry-simulated annealing using a total of 3333 experimental NMR restraints, consisting of 2409 approximate interproton distance restraints, 84 distance restraints for 42 backbone hydrogen bonds, 426 torsion angle restraints, 125 3JNH alpha restraints, 153 C alpha restraints, and 136 C beta restraints. The atomic rms distribution about the mean coordinate positions for the 30 structures for residues 7-137 and 145-163 is 0.42 +/- 0.04 A for the backbone atoms, 0.80 +/- 0.04 A for all atoms, and 0.50 +/- 0.03 A for all atoms excluding disordered side chains. The overall structure of MMP-1 is composed of a beta-sheet consisting of five beta-strands in a mixed parallel and anti-parallel arrangement and three alpha-helices. A best-fit superposition of the NMR structure of inhibitor-free MMP-1 with the 1.56 A resolution X-ray structure by Spurlino et al. [Spurlino, J. C., Smallwood, A. M., Carlton, D. D., Banks, T. M., Vavra, K. J., Johnson, J. S., Cook, E. R., Falvo, J., and Wahl, R. C., et al. (1994) Proteins: Struct., Funct., Genet. 19, 98-109] complexed with a hydroxamate inhibitor yields a backbone atomic rms difference of 1.22 A. The majority of differences between the NMR and X-ray structure occur in the vicinity of the active site for MMP-1. This includes an increase in mobility for residues 138-144 and a displacement for the Ca(2+)-loop (residues 74-80). Distinct differences were observed for side-chain torsion angles, in particular, the chi 1 for N80 is -60 degrees in the NMR structure compared to 180 degrees in the X-ray. This results in the side chain of N80 occupying and partially blocking access to the active site of MMP-1.
Oncostatin M (OM) is a member of the cytokine family which regulates the proliferation and differentiation of a variety of cell types and includes interleukin-6 (IL-6), leukemia inhibitory factor (LIF), and granulocyte-colony stimulating factor (G-CSF). This family of proteins adopts a four-helix bundle fold with up-up-down-down topology and contains intramolecular disulfide bonds. Since an X-ray or NMR structure for OM is not currently available, a homology model for OM was determined from the X-ray structures of human growth hormone (hGH), LIF, and G-CSF where the alignment was based on secondary structure instead of sequence. The OM secondary structure was determined from NMR structural data, and the secondary structures for hGH, LIF, and G-CSF were obtained from the reported X-ray structures. The resulting homology model was refined using sequential NOE distance 13C restraints, chemical shift information, and a conformational database.
Interaction of basic fibroblast growth factor (FGF-2) with heparin or heparan sulfate proteoglycans (HSPGs) is required for receptor activation and initiation of biological responses. To gain insight into the mechanism of activation of the FGF receptor by FGF-2 and heparin, we have used NMR, dynamic light scattering, and HSPG-deficient cells and cell-free systems. The first 28 N-terminal residues in FGF-2 were found to be highly mobile and flexible, consistent with the disorder found in both the NMR and X-ray structures. The structure of an FGF-2-heparin-decasaccharide complex that binds to and activates the FGF receptor was compared to a heparin-tetrasaccharide-induced complex that does not promote an interaction with the receptor. The major change observed upon the addition of the tetrasaccharide to FGF-2 was an increase in the correlation time consistent with the formation of an FGF-2 dimer. The NMR line widths of FGF-2 in the presence of the decasaccharide are severely broadened relative to the tetrasaccharide, consistent with dynamic light scattering results which indicate FGF-2 is a tetramer. The interaction of these heparin species with FGF-2 does not induce a significant conformational change in the overall structure of FGF-2, but small chemical shift changes are observed in both heparin and receptor binding sites. A trans-oriented symmetric dimer of FGF-2 is formed in the presence of the tetrasaccharide whereas two cis-oriented dimers in a symmetric tetramer are formed in the presence of the decasaccharide. This suggests that the cis-oriented FGF-2 dimer is the minimal biologically active structural unit of FGF-2. These data allow us to propose a novel mechanism to explain the functional interaction of FGF-2 with heparin and its transmembrane receptor.