Design of potent and selective human cathepsin K inhibitors that span the active site

Potent and selective active-site-spanning inhibitors have been designed for cathepsin K, a cysteine protease unique to osteoclasts. They act by mechanisms that involve tight binding intermediates, potentially on a hydrolytic pathway. X-ray crystallographic, MS, NMR spectroscopic, and kinetic studies of the mechanisms of inhibition indicate that different intermediates or transition states are being represented that are dependent on the conditions of measurement and the specific groups flanking the carbonyl in the inhibitor. The species observed crystallographically are most consistent with tetrahedral intermediates that may be close approximations of those that occur during substrate hydrolysis. Initial kinetic studies suggest the possibility of irreversible and reversible active-site modification. Representative inhibitors have demonstrated antiresorptive activity both in vitro and in vivo and therefore are promising leads for therapeutic agents for the treatment of osteoporosis. Expansion of these inhibitor concepts can be envisioned for the many other cysteine proteases implicated for therapeutic intervention.

Structure-assisted design of mechanism-based irreversible inhibitors of human rhinovirus 3C protease with potent antiviral activity against multiple rhinovirus serotypes

Human rhinoviruses, the most important etiologic agents of the common cold, are messenger-active single-stranded monocistronic RNA viruses that have evolved a highly complex cascade of proteolytic processing events to control viral gene expression and replication. Most maturation cleavages within the precursor polyprotein are mediated by rhinovirus 3C protease (or its immediate precursor, 3CD), a cysteine protease with a trypsin-like polypeptide fold. High-resolution crystal structures of the enzyme from three viral serotypes have been used for the design and elaboration of 3C protease inhibitors representing different structural and chemical classes. Inhibitors having α,β-unsaturated carbonyl groups combined with peptidyl-binding elements specific for 3C protease undergo a Michael reaction mediated by nucleophilic addition of the enzyme’s catalytic Cys-147, resulting in covalent-bond formation and irreversible inactivation of the viral protease. Direct inhibition of 3C proteolytic activity in virally infected cells treated with these compounds can be inferred from dose-dependent accumulations of viral precursor polyproteins as determined by SDS/PAGE analysis of radiolabeled proteins. Cocrystal-structure-assisted optimization of 3C-protease-directed Michael acceptors has yielded molecules having extremely rapid in vitro inactivation of the viral protease, potent antiviral activity against multiple rhinovirus serotypes and low cellular toxicity. Recently, one compound in this series, AG7088, has entered clinical trials.

Effects of retroviral vector design on expression of human adenosine deaminase in murine bone marrow transplant recipients engrafted with genetically modified cells.

To determine which features of retroviral vector design most critically affect gene expression in hematopoietic cells in vivo, we have constructed a variety of different retroviral vectors which encode the same gene product, human adenosine deaminase (EC 3.5.4.4), and possess the same vector backbone yet differ specifically in transcriptional control sequences suggested by others to be important for gene expression in vivo. Murine bone marrow cells were transduced by each of the recombinant viruses and subsequently used to reconstitute the hematopoietic system of lethally irradiated recipients. Five to seven months after transplantation, analysis of the peripheral blood of animals transplanted with cells transduced by vectors which employ viral long terminal repeats (LTRs) for gene expression indicated that in 83% (77/93) of these animals, the level of human enzyme was equal to or greater than the level of endogenous murine enzyme. Even in bone marrow transplant recipients reconstituted for over 1 year, significant levels of gene expression were observed for each of the vectors in their bone marrow, spleen, macrophages, and B and T lymphocytes. However, derivatives of the parental MFG-ADA vector which possess either a single base mutation (termed B2 mutation) or myeloproliferative sarcoma virus LTRs rather than the Moloney murine leukemia virus LTRs led to significantly improved gene expression in all lineages. These studies indicate that retroviral vectors which employ viral LTRs for the expression of inserted sequences make it possible to obtain high levels of a desired gene product in most hematopoietic cell lineages for close to the lifetime of bone marrow transplant recipients.

Identification of receptor-binding domains on human interleukin 5 and design of an interleukin 5-derived receptor antagonist.

A detailed structure-function analysis of human interleukin 5 (hIL5) has been performed. The hIL5 receptor is composed of two different polypeptide chains, the alpha and beta subunits. The alpha subunit alone is sufficient for ligand binding, but association with the beta subunit leads to a 2- to 3-fold increase in binding affinity. The beta chain is shared with the receptors for IL3 and granulocyte/macrophage-colony-stimulating factor--hence the descriptor beta C (C for common). All hIL5 mutants were analyzed in a solid-phase binding assay for hIL5R alpha interaction and in a proliferation assay using IL5-dependent cell lines for receptor-complex activation. Most residues affecting binding to the receptor alpha subunit were clustered in a loop connecting beta-strand 1 and helix B (mutants H38A, K39A, and H41A), in beta-strand 2 (E89A and R91A; weaker effect for E90A) and close to the C terminus (T109A, E110A, W111S, and I112A). Mutations at one position, E13 (Glu13), caused a reduced activation of the hIL5 receptor complex. In the case of E13Q, only 0.05% bioactivity was detected on a hIL5-responsive subclone of the mouse promyelocytic cell line FDC-P1. Moreover, on hIL5-responsive TF1 cells, the same mutant was completely inactive and proved to have antagonistic properties. Interactions of this mutant with both receptor subunits were nevertheless indistinguishable from those of nonmutated hIL5 by crosslinking and Scatchard plot analysis of transfected COS-1 cells.

Protein structure-based design of potent orally bioavailable, nonpeptide inhibitors of human immunodeficiency virus protease.

A class of potent nonpeptidic inhibitors of human immunodeficiency virus protease has been designed by using the three-dimensional structure of the enzyme as a guide. By employing iterative protein cocrystal structure analysis, design, and synthesis the binding affinity of the lead compound was incrementally improved by over four orders of magnitude. An inversion in inhibitor binding mode was observed crystallographically, providing information critical for subsequent design and highlighting the utility of structural feedback in inhibitor optimization. These inhibitors are selective for the viral protease enzyme, possess good antiviral activity, and are orally available in three species.

The design, computer modeling, solution structure, and biological evaluation of synthetic analogs of bryostatin 1

The bryostatins are a unique family of emerging cancer chemotherapeutic candidates isolated from marine bryozoa. Although the biochemical basis for their therapeutic activity is not known, these macrolactones exhibit high affinities for protein kinase C (PKC) isozymes, compete for the phorbol ester binding site on PKC, and stimulate kinase activity in vitro and in vivo. Unlike the phorbol esters, they are not first-stage tumor promoters. The design, computer modeling, NMR solution structure, PKC binding, and functional assays of a unique class of synthetic bryostatin analogs are described. These analogs (7b, 7c, and 8) retain the putative recognition domain of the bryostatins but are simplified through deletions and modifications in the C4-C14 spacer domain. Computer modeling of an analog prototype (7a) indicates that it exists preferentially in two distinct conformational classes, one in close agreement with the crystal structure of bryostatin 1. The solution structure of synthetic analog 7c was determined by NMR spectroscopy and found to be very similar to the previously reported structures of bryostatins 1 and 10. Analogs 7b, 7c, and 8 bound strongly to PKC isozymes with Ki = 297, 3.4, and 8.3 nM, respectively. Control 7d, like the corresponding bryostatin derivative, exhibited weak PKC affinity, as did the derivative, 9, lacking the spacer domain. Like bryostatin, acetal 7c exhibited significant levels of in vitro growth inhibitory activity (1.8–170 ng/ml) against several human cancer cell lines, providing an important step toward the development of simplified, synthetically accessible analogs of the bryostatins.

Peptide design by artificial neural networks and computer-based evolutionary search

A technique for systematic peptide variation by a combination of rational and evolutionary approaches is presented. The design scheme consists of five consecutive steps: (i) identification of a “seed peptide” with a desired activity, (ii) generation of variants selected from a physicochemical space around the seed peptide, (iii) synthesis and testing of this biased library, (iv) modeling of a quantitative sequence-activity relationship by an artificial neural network, and (v) de novo design by a computer-based evolutionary search in sequence space using the trained neural network as the fitness function. This strategy was successfully applied to the identification of novel peptides that fully prevent the positive chronotropic effect of anti-β1-adrenoreceptor autoantibodies from the serum of patients with dilated cardiomyopathy. The seed peptide, comprising 10 residues, was derived by epitope mapping from an extracellular loop of human β1-adrenoreceptor. A set of 90 peptides was synthesized and tested to provide training data for neural network development. De novo design revealed peptides with desired activities that do not match the seed peptide sequence. These results demonstrate that computer-based evolutionary searches can generate novel peptides with substantial biological activity.

Design, synthesis, and biological characterization of a peptide-mimetic antagonist for a tethered-ligand receptor

Protease-activated receptors (PARs) represent a unique family of seven-transmembrane G protein-coupled receptors, which are enzymatically cleaved to expose a truncated extracellular N terminus that acts as a tethered activating ligand. PAR-1 is cleaved and activated by the serine protease α-thrombin, is expressed in various tissues (e.g., platelets and vascular cells), and is involved in cellular responses associated with hemostasis, proliferation, and tissue injury. We have discovered a series of potent peptide-mimetic antagonists of PAR-1, exemplified by RWJ-56110. Spatial relationships between important functional groups of the PAR-1 agonist peptide epitope SFLLRN were employed to design and synthesize candidate ligands with appropriate groups attached to a rigid molecular scaffold. Prototype RWJ-53052 was identified and optimized via solid-phase parallel synthesis of chemical libraries. RWJ-56110 emerged as a potent, selective PAR-1 antagonist, devoid of PAR-1 agonist and thrombin inhibitory activity. It binds to PAR-1, interferes with PAR-1 calcium mobilization and cellular function (platelet aggregation; cell proliferation), and has no effect on PAR-2, PAR-3, or PAR-4. By flow cytometry, RWJ-56110 was confirmed as a direct inhibitor of PAR-1 activation and internalization, without affecting N-terminal cleavage. At high concentrations of α-thrombin, RWJ-56110 fully blocked activation responses in human vascular cells, albeit not in human platelets; whereas, at high concentrations of SFLLRN-NH2, RWJ-56110 blocked activation responses in both cell types. Thus, thrombin activates human platelets independently of PAR-1, i.e., through PAR-4, which we confirmed by PCR analysis. Selective PAR-1 antagonists, such as RWJ-56110, should serve as useful tools to study PARs and may have therapeutic potential for treating thrombosis and restenosis.

Semirational design of active tumor suppressor p53 DNA binding domain with enhanced stability

We have designed a p53 DNA binding domain that has virtually the same binding affinity for the gadd45 promoter as does wild-type protein but is considerably more stable. The design strategy was based on molecular evolution of the protein domain. Naturally occurring amino acid substitutions were identified by comparing the sequences of p53 homologues from 23 species, introducing them into wild-type human p53, and measuring the changes in stability. The most stable substitutions were combined in a multiple mutant. The advantage of this strategy is that, by substituting with naturally occurring residues, the function is likely to be unimpaired. All point mutants bind the consensus DNA sequence. The changes in stability ranged from +1.27 (less stable Q165K) to −1.49 (more stable N239Y) kcal mol−1, respectively. The changes in free energy of unfolding on mutation are additive. Of interest, the two most stable mutants (N239Y and N268D) have been known to act as suppressors and restored the activity of two of the most common tumorigenic mutants. Of the 20 single mutants, 10 are cancer-associated, though their frequency of occurrence is extremely low: A129D, Q165K, Q167E, and D148E are less stable and M133L, V203A and N239Y are more stable whereas the rest are neutral. The quadruple mutant (M133LV203AN239YN268D), which is stabilized by 2.65 kcal mol−1 and Tm raised by 5.6°C is of potential interest for trials in vivo.

(4-Aminomethyl)phenylguanidine derivatives as nonpeptidic highly selective inhibitors of human urokinase

Increased expression of the serine protease urokinase-type plasminogen activator (uPA) in tumor tissues is highly correlated with tumor cell migration, invasion, proliferation, progression, and metastasis. Thus inhibition of uPA activity represents a promising target for antimetastatic therapy. So far, only the x-ray crystal structure of uPA inactivated by H-Glu-Gly-Arg-chloromethylketone has been reported, thus limited data are available for a rational structure-based design of uPA inhibitors. Taking into account the trypsin-like arginine specificity of uPA, (4-aminomethyl)phenylguanidine was selected as a potential P1 residue and iterative derivatization of its amino group with various hydrophobic residues, and structure–activity relationship-based optimization of the spacer in terms of hydrogen bond acceptor/donor properties led to N-(1-adamantyl)-N′-(4-guanidinobenzyl)urea as a highly selective nonpeptidic uPA inhibitor. The x-ray crystal structure of the uPA B-chain complexed with this inhibitor revealed a surprising binding mode consisting of the expected insertion of the phenylguanidine moiety into the S1 pocket, but with the adamantyl residue protruding toward the hydrophobic S1′ enzyme subsite, thus exposing the ureido group to hydrogen-bonding interactions. Although in this enzyme-bound state the inhibitor is crossing the active site, interactions with the catalytic residues Ser-195 and His-57 are not observed, but their side chains are spatially displaced for steric reasons. Compared with other trypsin-like serine proteases, the S2 and S3/S4 pockets of uPA are reduced in size because of the 99-insertion loop. Therefore, the peculiar binding mode of the new type of uPA inhibitors offers the possibility of exploiting optimized interactions at the S1′/S2′ subsites to further enhance selectivity and potency. Because crystals of the uPA/benzamidine complex allow inhibitor exchange by soaking procedures, the structure-based design of new generations of uPA inhibitors can rely on the assistance of x-ray analysis.

Identification of a second aryl phosphate-binding site in protein-tyrosine phosphatase 1B: A paradigm for inhibitor design

The structure of the catalytically inactive mutant (C215S) of the human protein-tyrosine phosphatase 1B (PTP1B) has been solved to high resolution in two complexes. In the first, crystals were grown in the presence of bis-(para-phosphophenyl) methane (BPPM), a synthetic high-affinity low-molecular weight nonpeptidic substrate (Km = 16 μM), and the structure was refined to an R-factor of 18.2% at 1.9 Å resolution. In the second, crystals were grown in a saturating concentration of phosphotyrosine (pTyr), and the structure was refined to an R-factor of 18.1% at 1.85 Å. Difference Fourier maps showed that BPPM binds PTP1B in two mutually exclusive modes, one in which it occupies the canonical pTyr-binding site (the active site), and another in which a phosphophenyl moiety interacts with a set of residues not previously observed to bind aryl phosphates. The identification of a second pTyr molecule at the same site in the PTP1B/C215S–pTyr complex confirms that these residues constitute a low-affinity noncatalytic aryl phosphate-binding site. Identification of a second aryl phosphate binding site adjacent to the active site provides a paradigm for the design of tight-binding, highly specific PTP1B inhibitors that can span both the active site and the adjacent noncatalytic site. This design can be achieved by tethering together two small ligands that are individually targeted to the active site and the proximal noncatalytic site.

Design of polyzinc finger peptides with structured linkers

Zinc finger domains are perhaps the most versatile of all known DNA binding domains. By fusing up to six zinc finger modules, which normally recognize up to 18 bp of DNA, designer transcription factors can be produced to target unique sequences within large genomes. However, not all continuous DNA sequences make good zinc finger binding sites. To avoid having to target unfavorable DNA sequences, we designed multizinc finger peptides with linkers capable of spanning long stretches of nonbound DNA. Two three-finger domains were fused by using either transcription factor IIIA for the Xenopus 5S RNA gene (TFIIIA) finger 4 or a non-sequence-specific zinc finger as a “structured” linker. Our gel-shift results demonstrate that these peptides are able to bind with picomolar affinities to target sequences containing 0–10 bp of nonbound DNA. Furthermore, these peptides display greater sequence selectivity and bind with higher affinity than similar six-finger peptides containing long, flexible linkers. These peptides are likely to be of use in understanding the behavior of polydactyl proteins in nature and in the targeting of human, animal, or plant genomes for numerous applications. We also suggest that in certain polydactyl peptides an individual finger can “flip” out of the major groove to allow its neighbors to bind shorter, nontarget DNA sequences.

Immucillin H, a powerful transition-state analog inhibitor of purine nucleoside phosphorylase, selectively inhibits human T lymphocytes

Transition-state theory has led to the design of Immucillin-H (Imm-H), a picomolar inhibitor of purine nucleoside phosphorylase (PNP). In humans, PNP is the only route for degradation of deoxyguanosine, and genetic deficiency of this enzyme leads to profound T cell-mediated immunosuppression. This study reports the biological effects and mechanism of action of Imm-H on malignant T cell lines and on normal activated human peripheral T cells. Imm-H inhibits the growth of malignant T cell leukemia lines with the induction of apoptosis. Imm-H also inhibits activated normal human T cells after antigenic stimulation in vitro. However, Imm-H did not inhibit malignant B cells, colon cancer cell lines, or normal human nonstimulated T cells, demonstrating the selective activity of Imm-H. The effects on leukemia cells were mediated by the cellular phosphorylation of deoxyguanosine and the accumulation of dGTP, an inhibitor of ribonucleotide diphosphate reductase. Cells were protected from the toxic effects of Imm-H when deoxyguanosine was absent or when deoxycytidine was present. Guanosine incorporation into nucleic acids was selectively blocked by Imm-H with no effect on guanine, adenine, adenosine, or deoxycytidine incorporation. Imm-H may have clinical potential for treatment of human T cell leukemia and lymphoma and for other diseases characterized by abnormal activation of T lymphocytes. The design of Imm-H from an enzymatic transition-state analysis exemplifies a powerful approach for developing high-affinity enzyme inhibitors with pharmacologic activity.

Structure-based design of selective and potent G quadruplex-mediated telomerase inhibitors

The telomerase enzyme is a potential therapeutic target in many human cancers. A series of potent inhibitors has been designed by computer modeling, which exploit the unique structural features of quadruplex DNA. These 3,6,9-trisubstituted acridine inhibitors are predicted to interact selectively with the human DNA quadruplex structure, as a means of specifically inhibiting the action of human telomerase in extending the length of single-stranded telomeric DNA. The anilino substituent at the 9-position of the acridine chromophore is predicted to lie in a third groove of the quadruplex. Calculated relative binding energies predict enhanced selectivity compared with earlier 3,6-disubstituted compounds, as a result of this substituent. The ranking order of energies is in accord with equilibrium binding constants for quadruplex measured by surface plasmon resonance techniques, which also show reduced duplex binding compared with the disubstituted compounds. The 3,6,9-trisubstututed acridines have potent in vitro inhibitory activity against human telomerase, with EC50 values of up to 60 nM.

The structure of the human βII-tryptase tetramer: Fo(u)r better or worse

Tryptases, the predominant serine proteinases of human mast cells, have recently been implicated as mediators in the pathogenesis of allergic and inflammatory conditions, most notably asthma. Their distinguishing features, their activity as a heparin-stabilized tetramer and resistance to most proteinaceous inhibitors, are perfectly explained by the 3-Å crystal structure of human βII-tryptase in complex with 4-amidinophenylpyruvic acid. The tetramer consists of four quasiequivalent monomers arranged in a flat frame-like structure. The active centers are directed toward a central pore whose narrow openings of approximately 40 Å × 15 Å govern the interaction with macromolecular substrates and inhibitors. The tryptase monomer exhibits the overall fold of trypsin-like serine proteinases but differs considerably in the conformation of six surface loops arranged around the active site. These loops border and shape the active site cleft to a large extent and form all contacts with neighboring monomers via two distinct interfaces. The smaller of these interfaces, which is exclusively hydrophobic, can be stabilized by the binding of heparin chains to elongated patches of positively charged residues on adjacent monomers or, alternatively, by high salt concentrations in vitro. On tetramer dissociation, the monomers are likely to undergo transformation into a zymogen-like conformation that is favored and stabilized by intramonomer interactions. The structure thus provides an improved understanding of the unique properties of the biologically active tryptase tetramer in solution and will be an incentive for the rational design of mono- and multifunctional tryptase inhibitors.