Login| Sign Up| Help| Contact|

Patent Searching and Data


Title:
RICIN INHIBITORS AND METHODS FOR USE THEREOF
Document Type and Number:
WIPO Patent Application WO/2001/073438
Kind Code:
A1
Abstract:
Ricin A-chain is an N-glycosidase that attacks ribosomal RNA at a highly conserved adenine residue. Crystallographic studies show that not only adenine and formycin, but also pterin-based rings can bind in the ricin active site. For a better understanding of the recognition mode between ricin, and adenine-like rings, the interaction energies and geometries were calculated for a number of complexes. Shiga toxin, a compound essentially identical to the protein originally isolated from Shigella dysenteriae, has an active protein chain that is a homologue of the ricin active chain, and catalyzes the same depurination reaction. The present invention is drawn to identifying inhibitors of ricin and Shiga toxin, using methods molecular mechanics and ab initio methods and using the identified inhibitors as antidotes to ricin or Shiga toxin, or to facilitate immunotoxin treatment by controlling non-specific cytotoxicity.

Inventors:
ROBERTUS JON (US)
KERWIN SEAN MICHAEL (US)
YAN XINJIAN (US)
Application Number:
PCT/US2001/009400
Publication Date:
October 04, 2001
Filing Date:
March 23, 2001
Export Citation:
Click for automatic bibliography generation   Help
Assignee:
RES DEV FOUNDATION (US)
ROBERTUS JON (US)
KERWIN SEAN MICHAEL (US)
YAN XINJIAN (US)
International Classes:
C07D239/46; C07D239/48; C07D239/545; C07D285/22; C07D285/24; C07D473/18; C07D475/04; C07D487/04; C07D498/04; C07D513/04; C12Q1/34; G01N33/53; G01N33/569; C07D239/54; (IPC1-7): G01N33/566; G01N33/53; C07D471/02; C07D239/22; C07D473/24; C07D513/02
Foreign References:
US5219926A1993-06-15
US4371514A1983-02-01
Other References:
STECKER P.G. ET AL.: "The merck index of chemicals and drugs", 1960, MERCK AND CO., INC., RAHWAY, NJ, USA, XP002943916
YAN X. ET AL.: "Structure-based identification of a ricin inhibitor", J. MOL. BIOL., vol. 266, 1997, pages 1043 - 1049, XP002943917
OLSON M.A.: "Ricin A-chain structural determinant for binding substrate analogues: A molecular dynamics simulation analysis", PROTEINS: STRUCTURE, FUNCTION AND GENETICS, vol. 27, 1997, pages 80 - 95, XP002943918
SZEWCZAK A.A. ET AL.: "The conformation of the sarcin/ricin loop from 28S ribosomal RNA", PROC. NATL. ACAD. SCI. USA, vol. 90, October 1993 (1993-10-01), pages 9581 - 9585, XP002943919
RUTENBER E. ET AL.: "Crystallographic refinement of ricin to 2.5 A", PROTEINS: STRUCTURE, FUNCTION AND GENETICS, vol. 10, 1991, pages 240 - 250, XP002943920
Attorney, Agent or Firm:
Weiler, James F. (Suite 1560 Houston, TX, US)
Download PDF:
Claims:
WHAT IS CLAIMED IS :
1. A method of screening for a ricin inhibitor, comprising the steps of : identifying a test compound, wherein said test compound binds to the active site of the ricin A chain, has an aromatic heterocyclic molecular core resembling an adenine moiety in size and shape and is derivatized with polar substituents ; and measuring the ability of said test compound to displace binding of a second compound known to bind to the active site of the ricin A chain, wherein said second compound is selected from the group of pteroic acid, pteroic acid analogs, neopterin, pterin (3) tautomer, guanine (4) tautomer, 2-amino-5- (4\'- carboxyphenylaminomethyl)-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3- d] pyrimidine and 2-Amino-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine, wherein a competitive displacement of said second compound by the test compound indicates that the test compound is a ricin inhibitor.
2. The method of claim 1, wherein said test compound acts within the active site of ricin A chain, has nonpolar interactions with a side chain of an amino acid in said active site selected from the group Tyr 80, Ile 172 or Tyr 123, and has polar interactions with a side chain of an amino acid in the active site selected from the group of carbonyl oxygens of Gly 121 or Val 81, backbone amides of Val 81 or Tyr 123 and atoms on side chains of Arg 180, Tyr 80, Tyr 123 or Asn 78.
3. The method of claim 2, wherein said test compound further comprises at least one nonpolar moiety which interacts with side chains from Trp 211, Leu 45, Val 256, Tyr 257 or Thr 77 of said RTA chain.
4. The method of claim 2, wherein said test compound further comprises at least one polar moiety which interacts with the carbonyl oxygens of Thr 77 or Tyr 257, or side chains from Asn 47 and Arg 258 of said RTA chain.
5. The method of claim 2, wherein said test compound further comprises at least one polar moiety which interacts with an amino acid from a second pocket of said RTA chain.
6. The method of claim 1, wherein said identifying occurs by a technique selected from the group consisting of molecular modeling, crystallography, energy minimization protein synthesis assay and phage display.
7. A method of screening for a Shiga toxin inhibitor, comprising the steps of : identifying a test compound, wherein said test compound acts within an active site of Shiga toxin A chain, has an aromatic heterocyclic molecular core resembling an adenine moiety in size and shape and is derivatized with polar substituents such that said polar substituents interact in said active site of Shiga toxin A chain ; and measuring the ability of said test compound to displace binding of a second compound known to bind to the active site of the Shiga toxin, wherein said second compound is selected from the group of pteroic acid, pteroic acid analogs, neopterin, pterin (3) tautomer, guanine (4) tautomer, 2-amino-5- (4\'-carboxyphenylaminomethyl)-3, 4- dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine and 2-Amino-3, 4- dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine, wherein a competitive displacement of said second compound by the test compound indicates that the test compound is a Shiga toxin inhibitor.
8. The method of claim 7, wherein said identifying occurs by a technique selected from the group consisting of molecular modeling, crystallography, energy minimization protein synthesis assay and phage display.
9. A compound effective for inhibiting ricin, said compound able to act within an active site of RTA and having an aromatic heterocyclic molecular core, wherein said aromatic heterocyclic molecular core resembles an adenine moiety in size and shape and is derivatized with polar substituents such that said polar substituents interact in said active site of RTA.
10. The compound of claim 9, wherein said polar substituent is an amine group able to donate hydrogen bonds to a carbonyl oxygen of Val 81 or Gly 121.
11. The compound of claim 9, wherein said inhibitor further comprises at least one pendant group that binds an amino acid adjacent to said active site.
12. The compound of claim 9, wherein said inhibitor further comprises at least one moiety that reacts with a shallow channel in said RTA chain.
13. A compound effective for inhibiting ricin, said compound being able to act within an active site of RTA, having nonpolar interactions with a side chain of an amino acid in said active site selected from the group Tyr 80, Ile 172 or Tyr 123, and having polar interactions with a side chain of an amino acid in said active site selected from the group of carbonyl oxygens of Gly 121 or Val 81, backbone amides of Val 81 or Tyr 123, or atoms on side chains of Arg 180, Tyr 80, Tyr 123 or Asn 78.
14. The compound of claim 13, wherein said inhibitor further comprises at least one nonpolar moiety which interacts with side chains from Trp 211, Leu 45, Val 256, Tyr 257 or Thr 77 of said RTA chain.
15. The compound of claim 13, wherein said inhibitor further comprises at least one polar moiety which interacts with the carbonyl oxygens of Thr 77 or Tyr 257, or side chains from Asn 4 7 and Arg 258 of said RTA chain.
16. The compound of claim 13, wherein said inhibitor further comprises at least one polar moiety that interacts with an amino acid from a second pocket of said RTA chain.
17. The compound of claim 16, wherein said inhibitor interacts in a nonpolar fashion with side chains from amino acids selected from the group of Tyr 80, Val 82, Phe 57, Thr 77 and Arg 56.
18. The compound of claim 13, wherein said compound is a pteroic acid analog.
19. A compound of claim 18, wherein said analog is 2- amino-5- (4\'-carboxyphenylaminomethyl)-3, 4-dihydro-4-oxo-7H- pyrrolo [2, 3-d] pyrimidine.
20. The compound of claim 13, wherein said inhibitor is selected from the group consisting of 2, 4-diamino-6- hydroxypyrimidine, 2-amino-4, 6-dihydroxy-5-methylpyrimidine, 2- amino-4-hydroxy-6-methylpyrimidine, 4-amino-6-hydroxy-2- methylpyrimidine, 2-amino-6-hydroxy-8-mercaptopurine, 2, 6- diaminothiazolo [4, 5-d] pyrimidin-7 (6H)-one, 2-amino-5- (4\'- aminobenzylphenylamino) methyl-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3- d] pyrimidine and 2-amino-6-methyl-3, 4-dihydro-4-oxo-7H- pyrrolo [2, 3-d] pyrimidine.
21. A compound effective for inhibiting Shiga toxin, said compound able to act within an active site of STA1 and having a n aromatic heterocyclic molecular core, wherein said aromatic heterocyclic molecular core resembles an adenine moiety in size and shape and is derivatized with polar substituents such that said polar substituents interact in said active site of RTA.
22. The compound of claim 21, wherein said compound is a pteroic acid analog.
23. 2 3. The compound of claim 22, wherein said analog is 2- amino-5- [ (4\'-benzoyl) methylamino]-3, 4-dihydro-4-oxo-7H- pyrrolo [2, 3-d] pyrimidine.
Description:
RICIN INHIBITORS AND METHODS FOR USE THEREOF

BACKGROUND OF THE INVENTION Federal Funding Legend This invention was supported in part by federal funds, GM 30048 from the National Institutes of Health and by contract number DAMD17-94-C-4006 from the U. S. Army. The U. S. Government may have rights in this invention.

Field of the Invention The present invention relates generally to protein chemistry. More specifically, the present invention relates to the identification of inhibitors of the heterodimeric cytotoxin ricin, using computer modeling of the ricin active site as a template in structure- based drug design. Moreover, the present invention relates to the use of the identified inhibitors as antidotes to ricin or to facilitate immunotoxin treatment by controlling non-specific cytotoxicity. In addition, Shiga toxin, a compound essentially identical to the protein originally isolated from Shigella dysenteriae, has an active protein chain that is a homologue of the ricin active chain, and catalyzes the same depurination reaction. Thus, the present invention is drawn additionally to identifying inhibitors as antidotes to Shiga toxin.

Description of the Related Art Ricin is a potent heterodimeric cytotoxin easily isolated from the seeds of the castor plant, Ricinus communis. The protein consists of a lectin B chain, which can bind cell surfaces and is linked by disulfide bonds to an A chain (RTA), which enzymatically depurinates a key adenine residue in 28 S rRNA (Endo and Tsurugi, J.

Biol. Chem., 262 : 8128-30 (1987)). The molecular structure of ricin is represented in Figures 1A and 1B while the adenine depurination reaction mechanism is presented in Figure 2.

Ricin is an extraordinarily toxic molecule that attacks ribosomes, thereby inhibiting protein synthesis. Ricin has an LD5Q of approximately 1, ug/kg body weight for mice, rats, and dogs and is ten times more potent against rabbits (Olsnes, S. and Pihl, A.,"Toxic Lectins and Related Proteins,"The Molecular Action of Toxins and Viruses, pp 52-105 (1982)). The toxic dose for humans is likely to be in the Rg/kg range which ranks it among the most toxic substances known.

The protein has been used extensively in the design of therapeutic immunotoxins. In these constructs, ricin, RTA, or a related toxin is chemically or genetically linked to an antibody to form a so-called"magic bullet", which can target preferentially those cell lines carrying antigenic markers recognized by the antibody (Frankel, A. E., ed., Tmmunotoxins, Kluwer academic publishers, Boston (1988)).

Ricin also has been used as a poison agent. The protein gained notoriety when it was used in the famous"umbrella tip" assassination of Georgy Markov and was also used in an unsuccessful attempt to poison the Soviet dissident Alexander Solzhenitsyn. More recently, ricin was prepared by a militant anti-tax group which planned to poison IRS personnel.

Given the above, there is interest in identifying or designing potent inhibitors of ricin. These inhibitors could, in principle, be used to facilitate immunotoxin treatment by helping to control nonspecific cytotoxicity, or could be used as antidotes to poison attacks. Recently there has been interest in structure-based drug design--that is, using the knowledge of protein structure to identify enzyme inhibitors. The most common paradigm for this overall process has been called an"iterative protein crystallographic algorithm"by Appelt et al., J. Med. Chem. 34 : 1925-34 (1991), where the design of inhibitors for thymidylate synthetase is described. The main concept of the approach is that the protein active site is used as a template to design or to identify complementary ligands. The identified putative ligands are ranked and tested kinetically.

Promising inhibitor candidates are bound to the protein target and analyzed crystallographically for comparison with the proposed model. Additionally, alterations are made in the inhibitor to improve binding and a new round of tests is carried out on the altered compound.

A number of laboratories have used variations on this protocol to design efficacious inhibitors. For example, the search program DOCK (see Kuntz et al., J. Mol. Biol. 161 : 269-88 (1982)) was used to predict that the known anti-psychotic drug haloperidol would bind to the HIV protease (DesJarlais et al., Proc. Natl. Acad. Sci USA 87 : 6644-48 (1990)). Crystallographic studies together with computer aided search methods also were used in the design of inhibitors of purine nucleoside phosphorylase. The program GRID (Goodford, P. J., J. Med. Chem. 28 : 849-57 (1985)) has been used to design two very successful inhibitors of influenza virus. Certain chemical substitutions to the sialic acid substrate were predicted to b e

energetically favorable, based on interactions with the known X-ray structure of the enzyme. Subsequent binding assays revealed Ki values as low as 0. 2 nM. These compounds not only inhibited neuraminidase but retarded viral replication in cultured cell and animal models as well.

The X-ray structure of ricin has been solved (Montfort et al., J. Biol. Chem. 262 : 5398-03 (1987)), refined to 2. 5 A (Rutenber et al., Proteins 10 : 240-50 (1991)), and described in detail (Katzin et al., Proteins 10 : 251-59 (1991) ; and Rutenber and Robertus, Proteins 10 : 260-69 (1991)). The structure of RTA expressed in Escherichia coli has been resolved to 2. 3 A resolution for monoclinic crystals (Mlsna et al., Prot. Sci. 2 : 429-35 (1993)), and recently to 1. 8 A resolution for a tetragonal form (Weston et al., J. Mol. Biol. 244 : 410-422 (1994)).

The X-ray model of RTA allowed identification of a number of amino acids which were hypothesized to be important for substrate binding and for the depurination mechanism ; these residues include Glu 177, Arg 180, Trp 211, Tyr 80 and Tyr 123. Site-directed mutagenesis of the cloned RTA gene has been used to examine the relative significance of these residues (see, e. g., Schlossman et al., Mol. Cell. Biol. 9 : 5012-21 (1989) ; Frankel et al., Mol. Cell. Biol.

10 : 6257-63 (1990) ; Ready et al., Proteins 10 : 270-78 (1991) ; and Kim and Robertus, Protein Engineering 5 : 775-79 (1992)). In addition, Monzingo and Robertus, J. Mol. Biol. 227 : 1136-45 (1992), carried out an X-ray analysis of substrate analogs in the RTA active site, examining FMP, adenyl guanosine (ApG) and guanyl adenosine (GpA). The structure of important substrate and analog bases, together with the numbering scheme used in energy minimizations, is shown in Figure 8.

Several closely related mechanisms of action have been proposed which incorporate elements of both the structural and

kinetic analyses. It is likely that the susceptible adenine base binds between tyrosines 80 and 123 while forming specific hydrogen bonds with the backbone carbonyl and amido nitrogen of Val 81 and with the carbonyl of Gly 121. In the hydrolysis, the leaving adenine is at least partially protonated by Arg 180, and Glu 177 may stabilize a putative oxycarbonium transition state or, more likely, act as a base to polarize the attacking water.

Shiga toxin, a compound essentially identical to the protein originally isolated from Shigella dysenteriae, has an A chain which is activated by proteolysis, generating an active Al enzyme and an A2 fragment, which remains bound until a disulfide bond linking them is reduced. The active Al chain (STA1) is a homologue of RTA, and catalyzes the same depurination reaction. E. coli strains can carry the gene for Shiga toxin, and this renders them pathogenic. Human infection by Shiga toxin lead to hemorrhagic colitis and hemolytic- uremic syndrome--commonly referred to as"hamburger disease"--a severe and often fatal form of food poisoning. It is difficult to control outbreaks of hamburger disease because antibiotics tend to lyse the bacteria, releasing the destructive toxin into the system, aggravating tissue damage and internal bleeding. An effective inhibitor of STA1 would be a powerful adjunct to treatment of hemorrhagic colitis and hemolytic-uremic syndrome.

The Shiga toxin gene has been cloned and engineered to express the enzyme. A comparison of the amino acid sequences of RTA and STA show clearly that they are homologues. An energy- minimized model of STA1 was constructed (Deresiewicz et al, Biochemistry 31 : 3272-80 (1992)), and it was noted that key residues conserved among plant Ribosome Inactivating Protein (RIP) enzymes are conserved in STA1 (Katzin et al., Proteins 10 : 251-59 (1991)).

Until the present invention, no inhibitors for RTA or STA1 had been identified. Even FMP, which is known to bind RTA, is not an effective inhibitor of the RTA. Thus, the prior art is deficient in identifying compounds that are effective inhibitors of ricin. The present invention fulfills this long-standing need and desire in the art.

SUMMARY OF THE INVENTION The present invention provides compounds that are effective inhibitors of the cytotoxic proteins ricin or Shiga toxin.

In one object of the present invention, there is provided a compound effective for inhibiting ricin, the compound able to act within an active site of RTA and having an aromatic heterocyclic molecular core, wherein the aromatic heterocyclic molecular core resembles an adenine moiety in size and shape and is derivatized with at least one polar substituent such that the polar substituent interacts in the active site of RTA. In another embodiment of this object of the invention there is provided a compound wherein said polar substituent is an amine group able to donate hydrogen bonds to a carbonyl oxygen of Val 81 or Gly 121. As another embodiment of this object of the present invention, there is provided a compound wherein the inhibitor further comprises at least one pendant group which binds an amino acid adjacent to the active site. In addition, the inhibitor further may comprise at least one moiety which reacts with a shallow channel in the RTA chain.

In another object of the present invention, there is provided a compound effective for inhibiting ricin, the compound being able to act within an active site of RTA, having nonpolar

interactions with a side chain of an amino acid in the active site selected from the group Tyr 80, Ile 172 or Tyr 123, and having polar interactions with a side chain of an amino acid in the active site, selected from the group of carbonyl oxygens of Gly 121 or Val 81, backbone amides of Val 81 or Tyr 123, or atoms on side chains of Arg 180, Tyr 80, Tyr 123 or Asn 78. In addition, an embodiment of the present object provides an inhibitor further comprising at least one nonpolar moiety which interacts with side chains from Trp 211, Leu 45, Val 256, Tyr 257 or Thr 77 of said RTA chain. Alternatively, the inhibitor may possess at least one polar moiety which interacts with the carbonyl oxygens of Thr 77 or Tyr 257, or side chains from Asn 47 and Arg 258 of the RTA chain, or further comprise at least one polar moiety which interacts with an amino acid from a second pocket o f said RTA chain. In addition, the inhibitor may interact in a nonpolar fashion with side chains from amino acids selected from the group of Tyr 80, Val 82, Phe 57, Thr 77 and Arg 56.

In yet another object of the present invention, there is provided a compound effective for inhibiting ricin, wherein said compound is a pteroic acid analog.

In an additional object of the present invention, there is provided a compound effective for inhibiting Shiga toxin, the compound able to act within an active site of STA1 and having a n aromatic heterocyclic molecular core, wherein the aromatic heterocyclic molecular core resembles an adenine moiety in size and shape and is derivatized with polar substituents such that the polar substituents interact in the active site of RTA.

An additional object of the present invention, there is provided a compound effective for inhibiting Shiga toxin, wherein said compound is a pteroic acid analog.

Additionally, the present invention provides a method for identifying a compounds effective for inhibiting ricin, comprising the steps of performing at least one technique from the group of molecular modeling, crystallography, and energy minimization ; protein synthesis assay ; and phage display, and a method for identifying a compounds effective for inhibiting Shiga toxin, comprising the steps of performing at least one technique from the group of molecular modeling, crystallography, and energy minimization ; protein synthesis assay ; and phage display.

Another object of the present invention is to provide the pteroic acid analog as a pharmaceutical compound as an antidote to ricin or to facilitate immunotoxin treatment by controlling non- specific cytotoxicity.

Other and further aspects, features, and advantages of the present invention will be apparent from the following description of the presently preferred embodiments of the invention. These embodiments are given for the purpose of disclosure.

BRIEF DESCRIPTION OF THE DRAWINGS The appended drawings have been included herein so that the above-recited features, advantages and objects of the invention will become clear and can be understood in detail. These drawings form a part of the specification. It is to be noted, however, that the appended drawings illustrate preferred embodiments of the invention and should not be considered to limit the scope of the invention.

Figures 1A and 1B shows a ribbon representation of the ricin backbone. Figure 1A shows the entire ricin protein. The A

chain is in the upper right and the B chain is at the lower left. The two lactose moieties are represented as pairs of discs bound to the B chain. The chains have been separated to facilitate viewing. The disulfide bond linking the chains is indicated in the lower right portion of the molecule. The A chain has been divided into three domains as follows : residues 1-117 shown in small dot pattern, residues 118-210 shown as open ribbon, and residues 211-267 shown in large dot pattern. Figure 1B shows a close up representation of the ricin A chain.

Figure 2 shows the reaction mechanism by which ricin depurinates an adenine residue in 28 S rRNA.

Figure 3 shows the inhibition of ricin A chain (RTA) by pteroic acid (Figure 3A). Pteroic acid (squares) has a minor inhibitory effect on the protein synthesis using Artemia salinas ribosomes. Protein synthesis is reduced to about 20% in the presence of 1. 5 nM RTA (circles) but toxin action is inhibited by increasing doses of pteroic acid. Figure 3B replots the level of RTA activity as a function pteroic acid concentration.

Figures 4A and 4B show the electron density for pteroic acid complexed to RTA. This is OMIT density based on Fo-Fc coefficients and phases from the protein model in which Tyr 80 has been displaced to form the normal substrate-binding site. The difference density baskets are contoured at the 3s level. The refined position of PTA is shown superimposed on the difference density.

Figure 5 illustrates the interactions between pteroic acid and RTA. Hydrogen bonds are shown as dashed lines with the lengths indicated.

Figure 6 shows the superposition of the formycin (Figure 6A) and pterin rings (Figure 6B) in the RTA active site.

Carbon atoms are solid, nitrogen atoms have a light pattern and oxygens have a darker pattern. Hydrogen bonds between PTA and RTA are dashed lines and those between FMP and RTA (Monzingo and Robertus, J. Mol. Biol. 227 : 1136-45 (1992)) are dotted.

Figures 7A and 7B show the binding of neopterin in the RTA active site. Hydrogen bonds are shown as dashed lines and the lengths are indicated.

Figure 8 shows the structure and numbering of important substrate analogs for RTA. Nucleosides based on adenine, guanine, and formycin are linked to ribose at ring position 9. Pterin derivatives discussed in this specification are modified at position 6.

Figures 9A and 9B show the crystallographically observed complex of RTA and FMP. The figure shows the substrate recognition cleft occupied by FMP, an AMP analog. Hydrogen bonds are shown as dashed lines, bonds within RTA are lighter than those within FMP. Carbon and hydrogen atoms are black, oxygen is a light pattern, and nitrogen and phosphorous are a darker pattern.

Figures 10A to 10T show the chemical structures and partial charges for ring compounds examined in this study. Charges were derived from ab initio minimization and used in subsequent docking experiments with RTA.

Figures 11A to 11F show the bonding interactions predicted between RTA and various ring compounds. Partially charged compounds were docked into RTA and minimized in SYBYL. Hydrogen bonds are shown as dashed lines, bonds within RTA are lighter than those within the ligand ; carbon and hydrogen atoms are black, oxygen is a light pattern, and nitrogen is a darker pattern. The figure shows RTA interactions with 11A) formycin ring, 11B) adenine, 11C)

pterin (l) tautomer, 11D) pterin (3) tautomer, HE) guanine (4) tautomer, and 11F) AHA (3-amino-4-hydroxybenzoic acid).

Figure 12 shows the interactions between RTA and two forms of pterin. Figure 12A shows the hydrogen bonds and charge distribution of RTA interacting with pterin (l) and Figure 12B shows the interactions with pterin (3). Strong bonds (2. 6 to 2. 8 A) are dashed lines and weaker bonds (3. 0 to 3. 3 A) are dotted. Partial charges are indicated by the 8 symbol.

Figure 13 shows the correlation of ligand structure for the design of RTA inhibitors. Figure 13A shows the observed binding of nucleotides like FMP and AMP. Figure 13B shows the limitations of guanine nucleotides as ligands. Figure 13C and 13D show two representations for the binding of pterin based inhibitors. Figure 13 C and 13D differ in the hydrogen binding of pterin to Argl80.

Figure 14 illustrates the structures of two synthetic intermediates and of the pteroic acid analog 2-amino-5- [ (4\'- benzoyl) methylamino]-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 21.

Figure 15 depicts the RTA activity of 7-deazaguanine.

RTA activity is plotted versus inhibitor concentration.

Figure 16 shows a stereogram of the binding of 7- deazaguanine to the ricin A chain.

Figure 17 shows the synthesis scheme for 2, 4-diamino-6- hydroxypyrimidine 1.

Figure 18 shows the a-chloroketone based synthesis of 2-amino-6-methyl-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 6 from 2, 4-diamino-6-hydroxypyrimidine 1.

Figure 19 shows the a-chloroketone based synthesis of 2-amino-6-phenethyl-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 9.

Figure 20 shows the method of synthesis for 5-bromo- 2, 6-diamino-4-oxo-pyrimidine 13.

Figure 21 shows the method of synthesis for 2, 6- diaminothiazolo [4, 5-d] pyrimidin-7 (6H)-one 14 from 2, 4-diamino-6- hydroxypyrimidine 1.

Figure 22 shows the methods of synthesis for shows the method of synthesis of 2-amino-5- [ (4\'-benzoyl) methylamino]-3, 4- dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 21 (Figure 22A) and 2- amino-5- (4\'-aminobenzylphenylamino) methyl-3, 4-dihydro-4-oxo-7H- pyrrolo [2, 3-d] pyrimidine 24 (Figure 22B) from guanine. The synthesis of 24 used the same method as that of 21 except that 2-N- octanoylamino-5- (N, N-dibenzyl methyl amino)-3, 4-dihydro-4-oxo-7H- pyrrolo [2, 3-d] pyrimidine 18 is reacted with methyl 4-aminobenzoate 19 instead of 2-N-octanoylamino-5- (4\'-methylcarboxy phenylaminomethyl)-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 20.

Figure 23 shows the synthesis of 5-amino-2-methyl-6H- oxazolo [5, 4-d] pyrimidin-7-one (32).

Figure 24 shows the electron density for 8-methyl-9- oxoguanine complexed to RTA. Omit density was calculated with amplitudes Fobserved minus Fcalculated, and phases from the refined protein in which Tyr 80 has been displaced to form the normal substrate binding site.

Figure 25 compares the binding of pterin of PTA (Figure 25A) and 8-methyl-9-oxoguanine (Figure 25B) to the active site of RTA. Hydrogen bonds are indicated with dashed lines. The labeled N1

and N3 positions of the two heterocycles reflect the different nomenclature for pterin (Figure 25A) and guanine (Figure 25B) Figure 26 shows a stereogram of the 9-deazaguanine complex with RTA. Hydrogen bonds are indicated with dashed lines.

DETAILED DESCRIPTION OF THE INVENTION As used herein, the term"ricin"refers to a heterodimeric cytotoxic protein isolated from Ricinus communis, or a protein expressed in a heterologous system from the ricin gene. Ricin contains a lectin B chain (RTB) and an N-glycosidase A chain (RTA).

As used herein, the term"RTA substrate specificity pocket"refers to the portion of the RTA active site that binds adenine.

The term"RTA active site"refers to that portion of RTA which binds the rRNA substrate. The active site is large, generally polar, and interacts with the highly charged rRNA. The purine ring of the rRNA adenine rests between the side chains of tyrosines 80 and 123, and four polar protein groups in the active site (0121, 081, NH81 and NH1 of Arg 180) make hydrogen bonds with the substrate.

As used herein, the term"secondary pocket"refers to the portion of the RTA molecule that is a fairly large and generally nonpolar cavity of the opposite side of Tyr 80 from the adenine- binding site.

As used herein, the term"shallow channel"refers to the portion of the RTA molecule that forms a sort of rut, runs from the specificity pocket along its surface, and is flanked by Arg residues (213, 258, 48 and 56, moving away from the catalytic site).

As used herein, the term"inhibitor"refers to a small molecule which binds to RTA and retards its enzyme action.

As used herein the term"immunotoxin"refers to refers to a chemical or genetic linkage between an antibody with a binding preference for a target cell type, and a toxic protein such as ricin or RTA.

As used herein the term"non-specific cytotoxicity"refers to the intoxication of cells which are not targeted by an immunotoxin.

As used herein the term"pteroic acid"or"PTA"refers to 4- [[(2-amino-1, 4-dihydro-4-oxo-6-teridinyl) methyl] amino] benzoic acid.

As used herein, the term"formycin monophosphate"or "FMP"refers to an AMP analog in which the adenine ring is replaced by a pyrazolopyrimidine ring (that is, the N9 of adenine is replaced by C and C8 is replaced by N).

As used herein, the term"neopterin"refers to 2-amino-6- (1, 2, 3-trihydroxypropyl)-4 (3H)-pteridone.

As used herein, the term"guanine (4) tautomer"refers to the guanine tautomer in which hydrogen is bonded to N3 and not to N1.

As used herein, the term"inhibitor design"refers to a process of identifying inhibitors or improving the properties of a known inhibitor.

As used herein, the term"molecular modeling"o r "computer modeling"refers to fitting a model ligand compound, such as an inhibitor candidate, into the active site space of a protein model and adjusting the models according to chemical principles to simulate their interactions ; this is usually done with a computer display system.

As used herein, the term"X-ray structure"refers to a model of a protein or protein\'ligand complex which is based on experimental data obtained by X-ray diffraction from crystals of the protein or complex.

As used herein, the term"SYBYL"refers to a computer program distributed by TRIPOS Corporation which can be used to model proteins or complexes, and also can be used to calculate the energy of interaction between atoms or groups, based on various approximations of chemical laws.

As used herein, the term"6-31g** basis set"refers to a mathematical function used by the program SYBYLto approximate the distribution of electrons within atomic orbitals.

As used herein, the term"X-PLOR"refers to a suite of computer programs written by Axel Brunger and used routinely to adjust models to their lowest energy state, consistent with the observed experimental X-ray diffraction data.

As used herein, the term"energy minimization"refers to adjusting the atoms or chemical groups within a model so that the new model assumes the lowest energy--that is, the most likely conformation of the model, based on approximations to chemical laws.

As used herein the term"OMIT"refers to a kind of electron density map based on the observed X-ray amplitude data for a protein and phases calculated from a partial model.

As used herein the term"Fo-Fc"refers to the coefficients used to compute a kind of electron density map commonly called a "difference Fourier", and used to show those parts of a protein or complex which have changed as a result of some action, such as, but not limited to, inserting a ligand.

Ricin is a potent cytotoxin which has been used widely in the construction of therapeutic agents such as immunotoxins.

Recently it has been used by governments and underground groups as a poison. There is interest in identifying and designing effective inhibitors of the ricin A chain (RTA). The present invention utilizes computer-assisted searches to identify compounds that bind in the RTA active site, which normally recognizes a specific adenine base on rRNA. Kinetic assays indicated that pteroic acid (PTA) could inhibit RTA activity with an apparent Ki of 0. 6 mM, and a 2. 3 A crystal structure of the complex revealed the mode of binding. In the same way that PTA was identified as an inhibitor, other compounds have been assayed in the method of the present invention. The present invention is directed to compounds that are effective inhibitors of the cytotoxic protein ricin and methods for screening for additional inhibitors.

In addition, a molecular model of STA1 based on RTA has been proposed. The X-ray structure of the Shiga toxin B chain has been solved, as has the intact (AlA2) B5 toxin, which showed that the proposed model was sound. Respective matching of the RTA and STA1 active site finds Arg 180, the key catalytic residue in RTA, corresponds to Arg 170 in STA1, Glu 177 with Glu 167, Tyr 80 with Tyr 77, Tyr 123 with Tyr 114, and Trp 211 with Trp207. The X-ray structure of the Shiga toxin holotoxin (AlA2Bs) has been solved (Fraser, et al., Nature Struc. Biol. 1 : 59-69 (1994), confirming the general similarity of RTA and STA1.

Thus, the present invention provides a method of screening for a ricin inhibitor, comprising the steps of : identifying a test compound, wherein said test compound binds to the active site of the ricin A chain, has an aromatic heterocyclic molecular core resembling

an adenine moiety in size and shape and is derivatized with polar substituents ; and measuring the ability of said test compound to displace binding of a second compound known to bind to the active site of the ricin A chain, wherein said second compound is selected from the group of pteroic acid, pteroic acid analogs, neopterin, pterin (3) tautomer, guanine (4) tautomer, 2-amino-5- (4\'- carboxyphenylaminomethyl)-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3- d] pyrimidine 21, and 2-amino-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 15 wherein a competitive displacement of said second compound by the test compound indicates that the test compound is a ricin inhibitor. Preferably, the test compound acts within the active site of ricin A chain, has nonpolar interactions with a side chain of an amino acid in said active site selected from the group Tyr 80, Ile 172 or Tyr 123, and has polar interactions with a side chain of an amino acid in the active site selected from the group of carbonyl oxygens of Gly 121 or Val 81, backbone amides of Val 81 or Tyr 123 and atoms on side chains of Arg 180, Tyr 80, Tyr 123 or Asn 78. The test compound may further comprise at least one nonpolar moiety which interacts with side chains from Trp 211, Leu 45, Val 256, Tyr 257 or Thr 77 of said RTA chain. The test compound may further comprise at least one polar moiety which interacts with the carbonyl oxygens of Thr 77 or Tyr 257, or side chains from Asn 47 and Arg 258 of said RTA chain.

The test compound may further comprise at least one polar moiety which interacts with an amino acid from a second pocket of said RTA chain. Preferably, the identifying occurs by a technique selected from the group consisting of molecular modeling, crystallography, energy minimization protein synthesis assay and phage display.

The present invention also provides a method of screening for a Shiga toxin inhibitor, comprising the steps of : identifying a test

compound, wherein said test compound acts within an active site of Shiga toxin A chain, has an aromatic heterocyclic molecular core resembling an adenine moiety in size and shape and is derivatized with polar substituents such that said polar substituents interact in said active site of Shiga toxin A chain ; and measuring the ability of said test compound to displace binding of a second compound known to bind to the active site of the Shiga toxin, wherein said second compound is selected from the group of pteroic acid, pteroic acid analogs, neopterin, pterin (3) tautomer, guanine (4) tautomer, 2-amino-5- (4\'- carboxyphenylaminomethyl)-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3- d] pyrimidine 21, and 2-amino-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 15, wherein a competitive displacement of said second compound by the test compound indicates that the test compound is a Shiga toxin inhibitor. Preferably, the identifying occurs by a technique selected from the group consisting of molecular modeling, crystallography, energy minimization protein synthesis assay and phage display.

The present invention is also directed to a compound effective for inhibiting ricin, said compound able to act within an active site of RTA and having an aromatic heterocyclic molecular core, wherein said aromatic heterocyclic molecular core resembles an adenine moiety in size and shape and is derivatized with polar substituents such that said polar substituents interact in said active site of RTA. In one aspect, the polar substituent is an amine group able to donate hydrogen bonds to a carbonyl oxygen of Val 81 or Gly 121. In another aspect, inhibitor further comprises at least one pendant group which binds an amino acid adjacent to said active site.

The inhibitor may further comprise at least one moiety which reacts with a shallow channel in said RTA chain.

The present invention is also directed to a compound effective for inhibiting ricin, said compound being able to act within an active site of RTA, having nonpolar interactions with a side chain of an amino acid in said active site selected from the group Tyr 80, Ile 172 or Tyr 123, and having polar interactions with a side chain of an amino acid in said active site selected from the group of carbonyl oxygens of Gly 121 or Val 81, backbone amides of Val 81 or Tyr 123, or atoms on side chains of Arg 180, Tyr 80, Tyr 123 or Asn 78.

The inhibitor may further comprise at least one nonpolar moiety that interacts with side chains from Trp 211, Leu 45, Val 256, Tyr 257 or Thr 77 of said RTA chain. The inhibitor may further comprise at least one polar moiety which interacts with the carbonyl oxygens of Thr 77 or Tyr 257, or side chains from Asn 47 and Arg 258 of said RTA chain. The inhibitor may further comprise at least one polar moiety that interacts with an amino acid from a second pocket of said RTA chain. In one aspect, the inhibitor may interact in a nonpolar fashion with side chains from amino acids selected from the group of Tyr 80, Val 82, Phe 57, Thr 77 and Arg 56. The compound may be a pteroic acid analog such as compound #21. Representative examples of ricin inhibitor compounds include 2, 4-diamino-6- hydroxypyrimidine 1, 2-amino-6-methyl-3, 4-dihydro-4-oxo-7H- pyrrolo [2, 3-d] pyrimidine 6, 2, 6-diaminothiazolo [4, 5-d] pyrimidin- 7 (6H)-one 14, 2-amino-5- (4\'-carboxyphenylaminomethyl)-3, 4- dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 21, 2-amino-5- (4\'- aminobenzylphenylamino) methyl-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3- d] pyrimidine 24, 2-amino-4, 6-dihydroxy-5-methylpyrimidine 25, 2- amino-4-hydroxy-6-methylpyrimidine 26, 4-amino-6-hydroxy-2- methylpyrimidine 27, 2-amino-6-hydroxy-8-mercaptopurine 28 and 5-amino-2-methyl-6H-oxooxazolo [5, 4d] pyrimidin-7-one 32.

The present invention is also directed to a compound effective for inhibiting Shiga toxin, said compound able to act within an active site of STA1 and having an aromatic heterocyclic molecular core, wherein said aromatic heterocyclic molecular core resembles an adenine moiety in size and shape and is derivatized with polar substituents such that said polar substituents interact in said active site of RTA. In one aspect, the compound is a pteroic acid analog. A representative example of such an analog is 2-amino-5- [ (4\'- benzoyl) methylamino]-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3- d] pyrimidine.

The following examples are given for the purpose of illustrating various embodiments of the invention and are not meant to limit the present invention in any fashion.

EXAMPLE 1 Preparation Of Recomhinant RTA And Protein Synthesis Assay Recombinant RTA was prepared as described by Ready et al., Proteins 10 : 270-78 (1991). Recombinant RTA was prepared as described by Chaddock and Roberts, Protein Engineering 6 : 425-31 (1993). Briefly, 50 ml cultures of E. coli JM101 harboring plasmid pUTA (Ready et al, Proteins 10 : 270-78 (1991)) was grown overnight at 37°C, then used to inoculate 500 ml 2YT media and grown for 2 hours at 30°C. Expression of RTA was induced by addition of 50 gel 0. 5M isopropyl-l-thio-ß-galactoside (IPTG) and incubation for 4 hours at 30°C. Cells were harvested by centrifugation, resuspended in 30 ml 5 mM sodium phosphate (pH 6. 5), and broken with a French pressure cell. Following centrifugation to remove cell debris, the supernatant

was loaded onto a 100 ml carboxymethyl-sepharose (Pharmacia) column equilibrated in the same buffer. Unbound protein was washed from the column with 1500 ml buffer followed by 200 ml buffer containing 100 mM sodium chloride. RTA was eluted with a linear 500 ml 0. 1-0. 3 M sodium chloride gradient. Peak RTA-containing fractions were pooled and stored at 4°C at approximately 1 mg/ml.

Typical yields were 40-50 mg RTA/1 media for wild type RTA.

RTA activity was measured by its ability to inhibit protein synthesis in an in vitro protein synthesis assay using Artemia salinas ribosomes, as described previously by Ready et al, Biochem. Biophys.

Acta 740 : 19-28 (1983) and Ready et al, Proteins 10 : 270-78 (1991).

Thirty-four pmoles of Artemia salinas ribosomes were incubated in a volume of 120 u, l for 5 minutes at 25°C with inhibitor (pteroic acid) ranging in concentration from 0. 062 mM to 1. 67 mM, with or without 1. 5 mM spermine. The reaction was stopped by the addition of 400 nmoles of anti-RTA antibody. 100 ml of a protein synthesis mix containing 14C-Phe and wheat germ extract was added to all the tubes.

After a 15 minute incubation at 25 °C, the reaction was stopped by the addition of 5% TCA. The mixture was heated to 90°C for 10 minutes, then passed through a Millipore nitrocellulose filter. The filters were dried and counted by liquid scintillation counting. Inhibitor candidates were judged by their ability to disrupt RTA action against the Artemia ribosomes.

Pterin compounds including pterin-6 carboxylic acid, neopterin, pteroic acid, and folic acid were purchased from Sigma Chemical Company (St. Louis). Concentrations were evaluated spectroscopically, assuming the same ultraviolet spectrum as folic acid, that is e283 = 25, 000. Because of the limited solubility of the

pterins, the ribosome inhibition assays and controls were carried out at pH 9.

EXAMPLE 2 Crystallography And Molecular Modeling Crystals of RTA were grown in the monoclinic form (see Robertus et al., J. Biol. Chem. 262 : 19-20 (1987) ; and Mlsna et al., Prot. Sci. 2 : 429-35 (1993)) or in a tetragonal form (Weston et al., J.

Mol. Biol. 244 : 410-22 (1994)). Three-dimensional diffraction data were collected on a San Diego Multiwire Systems area detector (see Hamlin, Methods in Enzymology 114 : 416-52 (1985)) with a Rigaku RU-200 X-ray source operating at 50 kV, 110 mA with a graphite monochromator. Data was collected using the method of Xuong et al., J. Appl. Cryst. 18 : 342-50 (1985), and reduced and evaluated using the University of California, San Diego (UCSD) software system (Howard et al., Methods in Enzymology 114 : 453-72 (1985)). Rotation and translation searches, as well as crystallographic refinement of the RTA inhibitor complex (energy minimization and simulated annealing), were carried out using the X-PLOR package (Brunger, Crystallographic refinement by simulated annealing, Crystallographic Computing 4 : Techniques and New Technologies, Oxford : Clarendon Press, pp. 126- 140. (1988)). Molecular modeling was carried out using FRODO (see Jones, FRODO : A graphics fitting program for macromolecules, Computational Crystallography, Oxford : Oxford University Press, pp.

303-317 (1982)), running on an Evans and Sutherland PS390.

Molecular modeling programs CHEMX and SYBYL were purchased from Chemical Design Ltd. (Oxon, England) and Tripos Inc.

(St. Louis) respectively. Interaction energies were calculated with the Tripos force field. The energy minimization was terminated when a 0. 05 energy gradient value was reached using the Powell algorithm as implemented in SYBYL. A distance dependent and a constant function (dielectric constant = 1. 5) were used for dielectric effects, with the expectation that the distance dependent function can simulate the internal protein dielectric, while the constant function simulates the solvent effect on the protein. Non-bonded cutoff was 8. 0 C. The net atomic charges used in minimization were from the Kollman method (Weiner et al., J. Am. Chem. Soc., 106 : 765-84 (1984)) for RTA and from ab initio calculation (6-31g** basis set) for ligands ; these were calculated on CRAY YMP and COIL located at the National Institutes of Health.

EXAMPLE 3 Energy Minimization According to the ab initio minimization, the hybridization of an amino N atom connected to an aromatic ring is variable depending on its environment. Table 1 lists the deviation from the ring plane of two hydrogens on such N atoms. It shows that some N atoms are sp2 types, and some are combinations of sp2 and sp3 types.

If the deviation sum from the aromatic plane is larger than 3 0 degrees, the N atom was dealt with as sp3 type in the molecular mechanics calculation.

TABLE 1 The deviation sum (degree) of hydrogens of ammo N atom from aromatic ring plane according to ah mtio (6-31g**) rmmnnzation structuredeviation formycin 41. 1 adenine (l) 0. 7 adenine (2) 44. 0 Pterin (l) 37. 1 pterin (2) 0. 8 pterin (3) 36. 9 pterin (4 ion) 0. 1 guanine (l) 38. 1 guanine (2) 27. 8 guanine (3) 32. 1 guanine (4) 41. 3 2, 4-diamino-pteridine * 3-Amino-4-hydroxybenzoic acid (ion) 55. 7 S2 (ion) 0. 0 * Two amino N atoms, on which all hydrogens are in the aromatic ring plane.

In the energy minimization protocol, the energies for the crystal complexes of ricin with FMP (Monzingo, and Robertus, J. Mol.

Biol. 227 : 1136-45 (1992)) and pteroic acid were minimized first.

Next, these structures were reminimized after FMP and pteroic acid were replaced by formycin and pterin. Then the formycin or pterin was replaced by other structures one by one. The newly-formed complexes also were fully minimized by the criteria mentioned above for obtaining the interaction energies between ricin and ligands. The interaction energy was obtained by subtracting the unliganded ricin energy and ligand energy (in its binding conformation) from th e

complex energy. Some of the structures studied were small and not very different from formycin, adenine or pterin. It was assumed that these small structures causes only small changes in the RTA conformation and there is no requirement to perform molecular dynamics simulation for conformation searching.

EXAMPLE 4 Identification Of Pterrnr Acid As An Inhibitor A query was constructed using the geometric and bonding parameters of the observed FMP ligand, and a search of the NCI data base was made with CHEMX. Among the compounds predicted to have hydrogen bond donors and acceptors in an orientation favorable for RTA interaction was the pterin derivative called pteroic acid (PTA).

SYBYL calculations suggested that PTA might bind better than FMP.

Interaction enthalpies, calculated assuming one negative charge on each ligand, were-106 Kcal/mole for pteroic acid, as compared with- 89 Kcal/mole for FMP. It is important to note that the interaction enthalpies do not consider such terms as configurational entropy nor the effect of solvent interactions, and thus are not meant to represent the free energy of binding the ligand to the protein. Rather, the calculation is intended to serve as a rough guide to comparing the likely affinities of related compounds for the target enzyme.

To determine if pteroic acid was indeed an inhibitor of RTA, it was tested in the protein synthesis assay ; because of the limited solubility of PTA, the experiments and controls were carried out at pH 9. Figure 3A shows protein synthesis on 300 nM Artemia ribosomes in the presence of varying concentrations of pteroic acid, and in the presence or absence of 1. 5 nM RTA. The figure shows that

pteroic acid has a minor inhibitory effect on protein synthesis itself--a 1. 7 mM concentration of pteroic acid inhibits protein synthesis about 20%. Figure 3B replots the data as the fraction of possible RTA activity versus the concentration of pteroic acid. It shows that 0. 6 mM pteroic acid reduces RTA action by 50% and defines the effective Ki for the inhibitor.

Several efforts were made to observe crystallographically the binding of pteroic acid to RTA. The strongly diffracting tetragonal crystals are grown at pH 4. 5 and would not bind the inhibitor either when soaked in a saturated solution or when cocrystallized. However, monoclinic crystals at pH 9 were made 20 % saturated with pteroic acid and soaked for seven days. Although the crystals cracked initially, over the soaking time most of the cracks annealed.

Diffraction showed that the complex was not isomorphous with the native. Native cell parameters are a = 42. 6, b = 68. 1, c = 50. 2 A andp 112. 9° (Robertus et al., J. Biol. Chem. 262 : 19-20 (1987)). The pteroic acid complex crystals have the following parameters : a = 39. 0, b = 64. 4, c = 49. 4 A and ß= 108. 0°.

Three dimensional data were collected to 2. 3 A including 68, 262 observations of 10, 999 reflections ; reflection intensities scaled with R = 6. 5 %. Since the crystals were not isomorphous with the native, the RTA\'PTA complex was solved using molecular replacement methods. The monoclinic RTA was used as a model for a rotation search using the program X-PLOR. Data between 10 and 4 A showed no rotation. A translation search was then executed with X-PLOR. The translation search revealed a translation of 0. 65 A along x, and of 2. 8 2 A along z. The initial positioning was completed using the rigid body refinement option of X-PLOR ; the R factor to 3 A resolution was 0. 36.

A difference map with coefficients (Fo-Fc) was calculated and showed a mixture of positive and negative peaks within the active site. These peaks could be interpreted as the binding of PTA in the area previously seen to bind FMP. The pterin ring occupies the same site as the formycin ring of FMP or the adenine ring of the ApG dinucleotide. This requires that the Tyr 80 side chain rotate roughly 45° to stack with the aromatic surface of the pterin, in a way similar to that seen for FMP (see Monzingo and Robertus, J. Mol. Biol.

227 : 1136-45 (1992)). The Tyr 80 position in the model was adjusted and a second (Fo-Fc) map computed. This OMIT map is shown in Figures 4A and 4B. There is clear electron density for the pterin moiety of the inhibitor ; the benzoic acid density is weaker but still allows that group to be readily positioned. PTA was built into the density and the model subjected to four successive rounds of energy minimization and simulated annealing using X-PLOR. Between each round of refinement, the resolution was increased by 0. 2 A beginning at 3. 0 A and ending at 2. 3 A. Another difference map was calculated and 45 waters were added to the structure. After a final round o f energy minimization and isotropic temperature factor refinement, a final R-factor of 17. 9% was calculated for all data to 2. 3 A. Isotropic temperature factor refinement showed the pterin ring, with average B values of 11, was more rigidly held than the benzoic acid moiety where average B values were 26. The occupancy of the pteroic acid refined to 1. 0, as compared to 0. 5 seen previously for FMP (see Monzingo and Robertus, J. Mol. Biol. 227 : 1136-45 (1992)).

Figure 5 illustrates the interactions between PTA and the active site of RTA. The pterin ring binds in roughly the position of the substrate adenine, sandwiched between the rings of Tyr 80 and Tyr 123. For comparison, Figure 6A shows a superposition of the

formycin ring of FMP (Monzingo and Robertus, J. Mol. Biol. 227 : 1136- 45 (1992)) and Figure 6B shows the pterin ring of PTA. The 2. 6A distance between N1 of pterin and the carbonyl oxygen of Gly 121 suggests that the pterin is stabilized in a tautomeric form with a hydrogen on N1 (shown in Table 2) and this can be donated to the carbonyl group. Note that the 2-amino group of pterin also donates a hydrogen bond to the carbonyl oxygens of Gly 121 as well as to Val 81 ; this roughly mimics the role of the 6-amino group on an adenine substrate. In the same way, N3 of pteroic acid resembles N1 of adenine, receiving a hydrogen bond from the amido N of Val 81. The adenine substrate receives a hydrogen bond at N3 from Arg 180, and indeed this is the likely route of substrate protonation in the catalytic mechanism. Pteroic acid receives two hydrogen bonds from Arg 180, at the 4-oxo and N5 positions.

The benzoate moiety of PTA bends around the side chain of Tyr 80 and probably makes some nonpolar interactions with it. In fact, the benzoate ring appears to bind on the surface of a pocket--a site hypothesized by Monzingo and Robertus, J. Mol. Biol. 227 : 1136- 45 (1992), to be a second recognition site which may accommodate the guanine base of a natural rRNA substrate. Two hydrogen bonds are made to the benzoate, one between the carboxylate and Asn 78, and one to a water molecule which also bonds Arg 258. Although the pterin moiety of PTA makes strong and specific interaction within the adenine recognition site, it seems apparent that the benzoate moiety of PTA is not optimized for interactions within the second site.

EXAMPLE 5 Screening Additional Compounds Following the analysis of PTA, structural and kinetic studies were carried out on a number of other pterin-based compounds, as shown in Table 2. These compounds differ essentially in the size of the group attached to pterin at position 6.

TABLE 2 Pterin Based Compounds as Potential Inhibitors of RTA RTA Substituent Name Inhibitor X-ray Model OOC-6-carboxy pterin No No OH neopterin Ki=2mM 2. 5 A HO OH Pteric acid Ki = 0. 6 mM 2. 5 A NH IPTA) \'onc ~OOC) Folic acid No No NU \'ONC

Neopterin is pterin with propane triol at the six position.

It was observed in the protein synthesis assay to be a modest inhibitor of RTA, with a Ki > 2 mM (data not shown). Monoclinic crystals of an RTA complex were obtained by cocrystallization ; they were isomorphous with the native, and diffracted to 2. 5 A resolution.

58, 838 observations of 9, 030 reflections to 2. 5 A resolution were collected (Rmerge = 6. 6 %), the data reduced, and difference Fourier maps calculated. The pterin ring bound in the site previously seen to bind the formycin ring of FMP (Monzingo and Robertus, J. Mol. Biol.

227 : 1136-45 (1992) and the pterin ring of PTA as described above.

As with PTA and FMP, the side chain of Tyr 80 was displaced by the ligand. Neopterin was positioned by hand, including the propane triol moiety, and X-PLOR used to adjust position of the ligand and the protein in a simulated annealing refinement ; the R factor refined to 19. 5 % for all data. Coordinates for both the PTA and neopterin complexes have been submitted to the Brookhaven Protein Data Bank.

The binding of neopterin to RTA is shown in Figures 7A and 7B. The orientation of the pterin ring is similar to that seen for PTA. One difference is the bonding of Arg 180 to the inhibitor. It bonds to the 4-oxo group of neopterin, but does not bond to N5 as occurs in the PTA complex. Instead, a second bond is formed with the proximal hydroxyl of the propane triol moiety. Associated with this rearrangement is a 7° rotation of the pterin ring. The other atoms of the propane triol moiety of neopterin make no interaction with RTA.

Since the propane triol moiety of neopterin is much smaller than the corresponding substituents in PTA, it does not interact with Tyr 80 in the same way and appears to lack the van der Waals contribution to binding which might be expected in PTA. This is consistent with energy calculations using SYBYL, which show interruption enthalpies

between RTA and PTA to be-106 Kcal/mole and those between RTA and neopterin to be-73 Kcal/mole. Kinetic inhibition data, although not definitive, also suggest PTA is a better inhibitor and more likely to exhibit tighter binding. The stronger binding of PTA compared with neopterin argues that design of a pterin-based moiety derivatized at the 6 position which makes a strong and specific interaction with RTA is the preferred approach in creating an efficient inhibitor.

Crystals of RTA*neopterin were isomorphous with native RTA, but crystals of RTA*PTA were reproducibly not isomorphous ; RTA. PTA crystals exhibited a contraction of about 3 A along both the a and c axes. A least squares superpositioning of the two complexes showed very few significant structural differences between them. The largest changes were at the N and C termini, both far removed from the substrate (inhibitor) binding site. The chemical differences between the two inhibitors are centered on the size and charge of the moiety derivatizing pterin at position 6. However, there are no direct interactions of the bound inhibitor with any crystallographically- related molecules. Furthermore, the protein conformation itself is very similar in these two active sites, and so it does not appear that inhibitor binding induces protein changes which then are involved directly in crystal packing. It may be that inhibitor binding triggers a series of very subtle changes that propagate throughout the protein and cause a rearrangement of packing, although this seems unlikely.

It is also possible that saturating the solution with PTA has an unspecified solvent effect that tends to dehydrate the crystals and causes them to shrink. Except for this phenomenon, the binding of the two inhibitors seems to have roughly the same effect on RTA conformation, moving the Tyr 80 side chain. Further, inhibitor design can use this observed binding as a foundation and that novel inhibitor

models can be fit to this template without undo concern about predicting protein responses.

Neither pterin-6-carboxylic acid nor folic acid acted as inhibitors of RTA within the limits of their solubility. Efforts to soak the compounds into monoclinic crystals, even at saturating conditions, and cocrystallization efforts failed to produce stable complexes for X-ray analysis.

Given knowledge of the binding of PTA and neopterin it was predictable that neither pterin-6-carboxylic acid nor folic acid would bind to RTA. Refinement and energy minimization using SYBYL suggested that the interaction enthalpy for pterin-6-carboxylic acid, with one negative charge, was low (-64 Kcal/mole) compared with PTA or neopterin. Inspection of the hypothetical binding, assuming that the pterin moiety assumes the position seen in PTA, suggests that the 6 carboxylate is near the Glu 177 residue of RTA. Charge repulsion may therefore prevent binding of this simple compound.

Folic acid does not inhibit RTA, nor does it bind to RTA crystals. Folic acid resembles PTA accept that the benzoic acid moiety is derivatized with long and negatively charged glutamic acid.

Electrostatic mapping of the RTA surface shows the mouth of the second recognition site, beyond the area binding the PTA benzoate, is generally negative due to the presence of residues like aspartates 75, 96, and 100. This negative charge likely repels the folic acid group, driving it into solution and limiting binding compared with PTA.

Figures 10A-10T show the net atomic charges, derived from ab initio minimization with the 6-31g** basis set, on additional ligands examined by the methods of the present invention. Table 3 lists the molecular energies of these compounds in atomic units, au, where 1 au = 627. 5 Kcal. Table 3 also lists the interaction energies

between RTA and each ligand calculated from molecular mechanics methods (Tripos force field). The interaction energies are broken into Van der Waals and electrostatic components. The active site binding geometries for the energy minimized complexes with formycin, adenine, pterin (l), pterin (3) and guanine (4), computed by molecular mechanics, are presented in Figures 11A-11F.

TABLE 3 : molecular energy (au) from ab initio minimization and interaction energy (Kcal/mol) between ricin and these structures from molecular mechanics minimization molecular Interaction energy (Kcal/mol) structure energy (au) vdw electrostatic total formycin-464. 4888712-19. 2-24. 5-43. 7 adenine (l)-464. 5361520-18. 0-19. 4-37. 4 adenine (2)-464. 5211196-19. 0-20. 5-39. 5 Pterin (l)-577. 2721321-23. 1-22. 0-45. 1 pterin (2)-577. 2702506-22. 8-21. 8-44. 6 pterin (3)-577. 2624667-22. 6-31. 8-53. 4 pterin (4, ion)-577. 6574933-22. 9-24. 6-47. 5 guanine (l)-539. 4125625-18. 7-16. 1-34. 8 guanine (2)-539. 4129359 guanine (3)-539. 4055644-21. 8-20. 0-40. 8 guanine (4)-539. 4005964-20. 8-33. 0-53. 8 2-hydroxy formycin-539. 3674887-18. 7-29. 1-47. 8 2, 4-diamino-pteridine-557. 4495408-23. 2-18. 9-42. 1 aha (ion)-547. 6530542-19. 6-48. 5-68. 1 emodin-948. 2356142-30. 0-44. 0-74. 3 rhodizonic acid-677. 4887203-19. 6-25. 5-45. 1 lumazine-597. 0937652-17. 5-22. 5-40. 0 S1-597. 1110054-18. 0-25. 7-43. 7 S2 (ion)-466. 8918064-14. 5-41. 1-55. 6 S3-450. 2572027-15. 7-24. 1-39. 8 vdw---van der Waals interaction energy (kcal/mol) ; aha---3-Amino-4-hydroxybenzoic acid

Based on energy minimization, the formycin ring structure forms 7 hydrogen bonds in the RTA binding site (Fig 11 A).

The ligand has an interaction energy with RTA of-43. 7 Kcal/mol. The calculated orientation of the free formycin generally is similar to the ring seen crystallographically in the RTAFMP complex (Figures 9A/9B). However, note the free ring is shifted slightly so that an extra hydrogen bond is donated by the backbone NH of Tyr 123 to N8 of formycin. In addition, the guanidinium of Arg 180 donates two bonds instead of one to N3. The lack of the ribose ring allows the Arg 1 80 side chain to reorient in this way.

Adenine is the natural substrate for RTA and is recognized by the enzyme, at least when it is part of a GAGA loop context. For the two major adenine tautomers (see Figure 10), the molecular energy of adenine (l) is 9. 4 Kcal/mol (0. 015 au) lower than adenine (2), suggesting it is the normal solution form. Adenosine nucleotides are derivatized at N9 and so are frozen into the tautomer (l) form. Figure llb shows the interaction of adenine (l) with RTA. Table 3 shows that the free adenine (2) interacts slightly more strongly with RTA than does adenine (l). The main reason for this is that in adenine (2) the hydrogen on N7 can be donated to O of Gly 121. Weston et al., J. Mol. Biol. 244 : 410-422 (1994), found that AMP was hydrolyzed in a tetragonal form of RTA crystals leaving only adenine bound. The occupancy was low, however, and it was not possible to assess the tautomeric form of adenine. In the design of inhibitors, derivatives of either form would seem equivalent.

The interaction energies between RTA and the four underivatized pterin tautomers (one is in the ion form) were minimized and are listed in Table 3. Pterin (l), the most stable tautomer in aqueous solution, has an interaction energy around-45

Kcal/mol (Figure 11C), as strong as for formycin. The exocyclic 6- amino (N12 in Figure 1), donates hydrogen bonds to O of Gly 121 and Val 81, reminiscent of the amine group in formycin and adenine. The 4-oxo atom (Oll in Figure 8) receives bonds from NH of Val 81 and the Oyof Ser 176, while N8 receives a bond from the NH of Tyr 123 and N5 receives two bonds from Arg 180. Because the interaction energy is as strong as for formycin, it is reasonable to assume that pterin (l) may bind weakly in the RTA active site. However, pterin (l) forms fewer and generally weaker hydrogen bonds with RTA than does pterin (3), shown in Figure 11D. A comparison of the binding of the two pterin forms is also diagrammed in Figure 12, where strong hydrogen bonds are shown as dashed lines and weaker ones as dotted lines. It shows the pterin (l) complex forms three strong and four weak bonds, while pterin (3) complex forms six strong and two weak bonds. In particular, loss of the hydrogen from N3 of pterin (l) allows that atom to receive a bond from the NH of Val 81 (2. 8 C), and acquisition of a hydrogen at N1 allows it to donate a bond to the O of Gly 121 (2. 6 C). Consistent with this, the interaction energy for pterin (3) is about 8 Kcal/mole more favorable than for pterin (l). The X-ray data are consistent with the binding of pterin (3). The observed distance between pterin N3 and the backbone N of Val 81 is 3. 2 C, impossible if both have hydrogens on them. Also the distance between O of Gly 121 and N1 and the exocyclic N of pterin are 2. 6 and 2. 8 A respectively, suggesting two strong hydrogen bonds. It is apparent that the increased stability of binding tautomer (3) to RTA compensates for the cost of shifting from the more stable tautomer (l).

It is apparent from the interaction energies that the strong binding of PTA results not only from interactions with the pterin ring,

even its tautomer (3) form, but also from additional interactions with the benzoate-containing group linked at C6. Compared with PTA, the free pterin group also orients in a slightly different fashion based on molecular mechanics calculations. In particular, a hydrogen bond can be formed from Oy of Ser 176 to the 4-oxo of free pterin, whereas this is much weaker (3. 4 A) in the PTA complex. However, Ser 176 should be exploited in the design of inhibitors.

The binding of pterin (3) to RTA may result from simple equilibrium considerations. Pterin (3) is much less common in solution than form 1, but once bound the equilibrium shifts to produce more form 3. Alternatively, it may that the more stable pterin (l) initially recognizes the RTA active site. Then, a shift of the protons from N3 to N12 and a shift from N12 to N1 take place, creating the pterin (3) tautomer ; Figure 12 illustrates this process.

Tautomerization may be triggered by the repulsive forces between the positive charges on Arg 180, Ser 176, and Val 81 of the protein, and the positive charge on N3 of pterin (l). At the same time, there is an electrostatic attraction to move protons to N12 that is surrounded by negative charges from Gly 121 and Val 81.

The success of binding pterins to RTA raises questions about binding guanine derivatives, which have many of the same chemical features as pterins. Ab initio calculation shows that guanine (l), the form seen in nucleotides, is a low energy tautomer, while guanine (4) is the highest. From the interaction energies in Table 3, it can be seen that the guanine (4) interacts most strongly with RTA, as strongly as pterin (3). This putative binding to RTA is shown in Figure 11E, where the carbonyl (6-oxo) of guanine (4) accepts two hydrogen bonds from the guanidinium group of Arg 180, while other groups interact with RTA in a fashion similar to pterin (3). The

stability of guanine (4) compared to guanine (l) arises from the same considerations as were described for pterin (3) binding compared t o pterin (l). Interaction energy calculations suggest guanines are reasonable inhibitor candidates. Attempts to bind 2-amino-6, 8- dihydroxypurine (8-hydroxy guanine) to RTA were unsuccessful although the solubility of such compounds is limited and crystallization conditions may not allow strong binding to occur.

Since the interaction energies of guanine compounds is computed to be high compared to adenine, one might ask about the specificity of RTA for adenine. Figure 13 integrates information about the bonding of adenines, pterins, and guanines. It shows that guanine nucleotides cannot bind to the RTA active site for trivial steric reasons. Figure 13A shows the binding of adenine-like nucleotides such as FMP and ApG seen crystallographically. The R group would be a ribose linking the base into a larger rRNA molecule. Figures 13B and 13C show alternative binding schemes for the binding of pterin based ligands (Figures 13B and 13C differ in the hydrogen binding between the pterin ligand and Argl80) ; PTA and neopterin are derivatized at position 6, shown by the R. Figure 13D shows that the guanine (4) base can assume an electronic configuration like pterin, but compared with adenine, it is flipped over. Note that the R group at position 7 now clashes with the protein. Non-nucleoside guanine compounds might still be expected to bind well to RTA if they are modified at positions 8 or 9. Compounds like 2-amino-6, 8-dihydrooxypurine did not bind to RTA, but are poorly soluble and it simply may be that under crystal conditions these compounds cannot achieve high enough concentrations to bind well to the protein. On the other hand, it should be noted that, unlike pterin (l), the interaction energy of guanine (l) with RTA is quite low, about-34. 8 Kcal/mol, compared to-

45. 1 Kcal/mol. If a model as shown in Figure 12 is appropriate, the initial recognition of guanine may be very unlikely.

In addition to the purine and pterin bases, a number of other compounds were analyzed as potential RTA inhibitors. The interaction energy of 2, 4-diamino-pteridine with RTA is-42. 1 Kcal/mol (Table 3). This energy is substantially weaker than for pterin because the 4-amino hydrogens repulse the guanidinium group of Arg 180 and N1 repulses the carbonyl O of Gly 121. Rhodizonic acid has 4 carbonyl groups and 2 hydroxyl groups ; it might seem to b e flexible enough to form a variety of interactions with RTA. Its interaction energy is about-45 Kcal/mole, similar to formycin, and yet it failed to bind to RTA crystals. It should be noted however, that the interaction energy for FMP is-89 Kcal/mole and the additional interactions outside the specificity pocket may be important to judging the effectiveness of binding. In a similar fashion, lumazine, Sl, and S3 (Table 3 and Figure 10) interact with RTA at only modest strength, implying that they are unlikely to be good inhibitors of RTA themselves, but can serve as the parent compound for other inhibitor classes.

Table 3 shows the interaction energies of AHA and emodin with RTA to be relatively strong, based on molecular mechanics.

Although they are predicted to bind well to RTA, they have quite different interaction geometry compared with adenine, pterin, and formycin. The putative binding of AHA is shown in figure 11F.

Attempts to soak these highly colored compounds into RTA crystals failed, showing that a strong interaction energy is not sufficient to assure ligand binding.

In summary, the RTA specificity site is formed largely by the backbone of residues Gly 121, Asn 122, Tyrl23 and Val 81, and by

side chains of Arg 180, Tyr 80, Tyrl23, Glu177, and Trp 211. The interaction energy from hydrogen bonds is a major factor in determining RTA substrate specificity, and in particular at recognizing one adenine base in a particular rRNA context. The main donors and acceptors of RTA binding site to a base are the carbonyls of Gly 121 and Val 81, the amine of Val 81 and the side chain Arg 18 0.

Obviously, the carbonyls and backbone amides are not flexible. The Arg 180 chain is not flexible either, because it forms an aromatic stack with Trp 211 and is strongly fixed by hydrogen bonds with the carboxylate of Glu177 and the peptide carbonyl of Asn 78. The guanidinium side chain can move only within very limited space, and in the minimizations of the present invention, it always maintains its geometry. Taken together, this suggests that the RTA active site has evolved to act as a fairly rigid receptor to bind specific, complementary ring structures bearing a resemblance to adenine. It is possible to bind compounds, like pterins, which make additional contacts and have stronger interaction energies with the active site.

This is an important consideration in inhibitor design. RTA has evolved to recognize an extended rRNA substrate and its strong binding to ribosomes (Km ~ 1, uM) arises partly from interactions at sites remote from the catalytic adenine recognition site. It is these areas that may need to be targeted for optimum inhibitor design.

EXAMPLE 6 Phage Display And Inhibitor Design In addition to query-based search methods, phage display search methods are used to identify peptides which bind strongly to

RTA and Shiga toxin. The method is useful in inhibitor/drug design in two ways. First, peptides can be identified that inhibit the toxins directly. Second, the crystal structure of tight-binding peptides likely reveals novel enzyme\'ligand interactions since the peptides differ markedly from the heterocyclic systems observed crystallographically.

The novel interactions identified are then fed into the query-based searches.

The principles of the phage display method are generally well known (Clackson and Wells, TIBTech 12 : 173-84 (1994)). M 1 3 are filamentous phage with a major coat protein called P8 and a surface protein P3 involved in infectivity ; both proteins have been engineered for phage display. P3 is localized at one end of the phage and is present in only 3 to 5 copies ; the amino terminus of P3 can accommodate large peptide inserts without interfering with folding or infectivity. A commercial library was obtained from New England Biolabs containing all of the 1. 3 x 109 (=207) random heptapeptides displayed on P3. This is about the upper limit of library size, given current transformation efficiencies for E. coli, and the limits of convenient handling of cell volumes. The manufacturer\'s simplest protocol for screening suggests attaching RTA to the wells of a microtiter plate via nonspecific hydrophobic interactions.

Roughly 1011 phage (about 100 copies of each displayed heptamer) were applied and allowed to bind. Those with poor binding were simply rinsed away. The tighter binding phage were then eluted, as suggested, using free RTA. The eluted phage were amplified in E coli, titered, and used in a second screening round. A total of four rounds were carried out, presumably selecting for stronger binding peptides (the phage were eluted in rounds 3 and 4 using low p H instead of free RTA, as is discussed below). Ten phage from the last

round were sequenced over the insert region, using the commercially available sequencing primer. Seven phage showed a consensus tetrapeptide sequence : Ser- (Ser/Thr)- (Leu/Ile)-Pro.

Native RTA is presented to the phage library, as opposed to RTA which has been nonspecifically adsorbed to microtiter plates.

The standard RTA expression vector used in these experiments has been engineered to include coding for an amino terminal His6 peptide.

The His-tagged protein (his-RTA) is soluble and fully active and has been purified on standard Ni-NTA columns, indicating it binds Ni well.

Microtiter plates obtained from QIAGEN Corporation coated with a Ni- NTA matrix binds the His-RTA to the wells through the His leader, displaying a fully active enzyme with its active site facing the solvent.

Phage binding to the RTA plates could be eluted with free RTA in rounds 1 and 2, but free RTA could not release phage in round 3. This has been interpreted to indicate that free RTA is not a strong enough competitor to release the peptide ligands and denaturation of RTA with low pH is required to allow phage elution. The implication is that the very strongest ligands were probably not released in the early rounds and were not amplified. Ultimately, peptides were selected with a clear affinity for RTA. Stronger stripping methods, such as low pH are used to identify better ligands.

Once a strongly-binding peptide is identified, it must be synthesized, being acetylated at the N terminus and amidated at the C terminus. The neutralization of the terminal charges minimize the chances of repulsion between the peptide and charged groups in the RTA active site. Once synthesized, the peptide is purified and tested as an inhibitor of RTA (for example, in the protein synthesis assay described in Example 1), and is then diffused into RTA crystals or

cocrystallized with RTA for crystallographic examination as described in Example 2.

In addition to the linear-heptapeptide library, a cyclic peptide library may be utilized. A cyclic peptide is characterized by a fixed disulfide bond defining a loop of six random amino acids. The disulfide is flanked by at least one random amino acid on each end.

Alternatively, a similar library is used, having four randomized residues in the loop, or having a random hexapeptide loop and four randomized residues flanking the disulfide. Cyclic peptides are advantageous in that they are structurally more rigid and constrained, and have greatly reduced flexibility compared to linear peptides.

Binding of the corresponding free cyclic peptides to the enzyme target is entropically favorable compared with linear peptides and they tend to have stronger association constants. For example, it has been shown that the affinities of cyclic peptides for streptavidin exceed those of similar linear peptides by 100-1000 fold (Giebel, et al., Biochemistry 34 : 15430-35 (1995). Further, X-ray structures confirm the notion that the cyclic peptides are more compact and constrained in the active site.

EXAMPLE 7 Inhibitor Design For STA1 Inhibitor design for STA1 follows the paths described in the Examples above for RTA ; that is, computer searches, kinetic assays and phage display techniques are utilized. The plasmid pSC25 expresses STA. A stop codon is engineered into this plasmid at a position to produce the 32 kD active STA1 peptide, removing the

inhibiting STA2 peptide instead of cleaving it with proteases. The engineered 28 kD protein is expressed and isolated for crystallization experiments. Further, a polyHis sequence is engineered to both the N and C termini of STA1, which expedites purification and allows for isolation of useful crystals since each modified construct represents an independent protein, as far as crystallization is concerned. Such tags attached to RTA did not hamper full enzymatic activity, and it is likely that identical tags will not hamper STA1 enzymatic activity. The work described above has provided a model of STA that is used for structure-based drug design ; however, a workable crystal model of STA makes the best use of the iterative crystallographic algorithm of drug design.

Expression of the STA1 protein in the engineered pSC25 allows for kinetic tests to be carried out on various inhibitors.

Obviously, due to the similarities between STA1 and RTA, it makes sense to try inhibitors on STA1 that have been identified as being effective on RTA. For example, pteroic acid inhibits RTA, and is tested on STA1. There are, however, differences between RTA and STA1 that potentially are important. Gly 121 is a key residue in the RTA active site, since its carbonyl group receives a hydrogen bond from the 2 amino group of pteroic acid. The corresponding residue in STA1 is Ser 112, and the disposition of the Ser side chain suggests that it may also interact with pteroic acid. However, the Serll2 hydroxyl can form an additional bond, donating to N8 of the inhibitor, that is not formed by the RTA Gly, suggesting that pteroic acid may bind to STA1 more effectively that it does to RTA.

EXAMPLE 8 Synthetic Approach To Identifying Inhihit (lrs Once compounds are identified though the above- referenced computer modeling, crystallographic, energy minimization, phage display or protein synthesis assay means, the compounds are synthesized. A recent advance in synthetic organic chemistry involves the multiple simultaneous (multiplex) synthesis of an analogous series of compounds. Rather than performing one to three reactions in a given day, 30-100 reactions are run using an appropriate reaction manifold. The key is the purification of the individual products, which can be facilitated by supporting the fixed reagent on an insoluble polymer allowing simple filtration, followed by cleavage of the product from the support.

Alternatively, polymer-supported reagents are used t o affect coupling of the two reactants, or to scavenge side-products, thereby converting them to insoluble material, which is easily separated from the reaction mixture by filtration. In yet another alternative, solution-phase chemistry is employed, followed by standard aqueous work-up procedures or selective crystallization of the of the product from the reaction mixture.

A relatively inexpensive apparatus facilitating parallel syntheses and allowing temperature control, variable speed mixing, and the maintenance of an inert atmosphere over each reactor has been devised. Work-up and purification of the multiple reactions run in this apparatus is simplified by the use of a solid-supported scavenger which reacts with excess reagent to form an insoluble product that is easily separated from the reactions mixture by filtration. Between 10 and 100 mg of each product is synthesized. In

addition to providing ample sample for biological evaluation, this allows each product\'s identity and purity to be assayed using NMR techniques.

Scheme 1 below shows one class of parallel synthesis starting with pteroic acid and using carbonyldiimidazole or oxalyl chloride coupling agents or solid-phase synthesis with the pteroic acid group attached through the 2-amino substituent to a polymer support.

Scheme 2 shows a second example of multiplex synthesis using 2- pivaloyl-6-formylpterin as the starting material. Condensation of this material with a variety of amines leads to imine intermediates that are reduced to afford the 2-pivaloyl-protected analogs. Purification o f these protected analogs is simplified by use of a polymer-supported reducing agent, which is removed by filtration. The protected analogs are subjected to mild deprotection conditions and isolated by crystallization.

In parallel with the multiplex synthesis approach, traditional synthetic techniques are employed to prepare analogs of compounds identified as effective in inhibiting RTA or STA1. For example, analogs of pteroic acid in which the pterin moiety is replaced with an alternative heterocycle are prepared. Focus is o n heterocycles predicted to adopt a tautomeric structure with strong binding to the adenine specificity pocket. The hydrogen binding properties of pterin (3) are present, but do not have to undergo an energetically unfavorable tautomerization. The structure below shows the generalized form of a series of these compounds. In this series, the steric and hydrogen bonding potential about the sulfur is varied systematically by manipulating the oxidation state from the sulfide to sulfoxide to sulfone.

Additional changes in the hydrogen bonding potential are explored by examining both the 1, 2, 4-benzothiadiazine (I, X-CH) and 1, 2, 4-pyridothisadiazine 1, 1-dioxides ; the preferred tautomeric form in solution is the 4-H form shown in structure I. An additional advantage to these compounds when compared to pteroic acid is increased aqueous solubility. In addition, modeling within the RTA specificity site suggests that the polar interaction between the sulfone group (Y=SO2) and Arg 180 is substantially stronger than those made to the 4-oxo group of PTA. Finally, when developing inhibitors into drugs, it is important to note that compounds identified may bear a strong resemblance to known and approved compounds ; for example, the compounds described above resemble benzothiadiazine diuretics like Aldactazine and Diuril.

The synthesis of compounds of general structure I proceeds based upon literature precedence. A general route to 2- amino-substituted 1, 2, 4-benzothiadiazines involving displacement of a 3-methylthio group has been reported by Pirotte, et al, J. Med. Chem.

36 : 3211-13 (1993) ; and Tullio et al, Tetrahedron 51 : 3221-34 (1995) have demonstrated the synthesis of 3-amino-substituted [2, 3] pyrido- 1, 2, 4-thiadiazines from 3-amino-2-sulfamoylpyridine, as shown in Scheme 3. Although 1, 2, 4-pyridothisdiazine-1, 1-dioxides are known, there have not been reports of 1, 2, 4-pyridothisdiazines in which the sulfur is not oxidized, or is in the sulfoxide oxidation state. However, Finch et al., J. Org. Chem. 45 : 3416-21 (1980) reported routes to the corresponding 1, 2, 4-benzothiadiazines and 1, 2, 4,-benzotthisdiasine 1- oxides, Scheme 4. Adaptation of these rounds allows production of the corresponding pyrido-analogs.

Scheme 1 N N NH2 RNH N rant HIC H02C N N N-H 2 HO2Cw N NH EX : 1 NHR rHR Scheme N NHPiv I 1 ArNH= \'IHI N OHC N O N N NHPiv I) txane HCI ArNH N 0 N O \ANxNyNH2 ArNH zbN o N 0 Structure I X = N or CH Y = S, SO, or SO2 Scheme 3 gg NH, urea N SONT, H N U P, Ss N"S" NH N su O O H NH NaHC03 merl //nu o 0 H N SMe MHZ /N / N sI\'llN H N NHR / N . N U O

Scheme 4 H IV NH, N SMe PhSONSO R I AcOH SMe I//N reflux 11 . J SMe \. NTos N SMe \ N SMe sl---, N S S/N j j 11 0 |RNH |RNH HH NHR N NHR R SON n EXAMPLE 9 Binding interactions of small ring compounds with the RTA active site Computer-aided search methods, as described supra, indicated the compounds based on the pterin group act as inhibitors of ricin. For example, pteroic acid acts as a modest inhibitor with an

IC50 of about 0. 6 mM (Yan, X., et al., J. Mol. Biol, 266 : 1043-1049 (1997)). A series of computations on the binding interactions of small ring compounds with the RTA active site helped in the creation of a chemical template of the active site to aid in further drug design (Yan, X., et al., Proteins, 31 : 33-41 (1998)). The results are shown in Table 4 below. The data clearly demonstrate that the pterin (3) tautomer binds better than formycin, which is in agreement with kinetic and crystallographic results ; it is further suggested that guanine compounds, in the correct orientation, may be effective inhibitors.

TABLE 4 Interaction energy in kcal/mol, between RTA and small ring compounds RTA Interaction Energy Aqueous phase Solvation vdW Elec Total (IE) Net Compound energy energy (SE) (IE-SE) formycin 84. 6-20. 8-19. 2-24. 5-43. 7-22. 9 adenine (l) 63. 7-23. 1-18. 0-19. 4-37. 4-14. 3 adenine (2) 69. 9-23. 6-19. 0-20. 5-39. 5-15. 9 Pterin (l) 26. 6-23. 6-23. 1-22. 0-45. 1-21. 5 pterin (2) 33. 4-21. 8-22. 8-21. 8-44. 6-22. 8 pterin (3) 30. 5-25. 2-22. 6-31. 8-53. 4-28. 2 pterin (4, ion)-62. 4-75. 0-22. 9-24. 6-47. 5-27. 5 guanine (l) 22. 1-26. 8-18. 7-16. 1-34. 8-8. 0 guanine (2) 27. 4-22. 1-21. 8-20. 0-40. 8-18. 7 guanine (3) 30. 0-25. 3-20. 8-33. 0-53. 8-28. 5 vdw-van der Waals interaction energy (kcal/mol)

EXAMPLE 10 Synthesls of CTIlanine-hased compouncls Experimental All evaporations were carried out in vacuo with a rotary evaporator. Melting points were determined on a Thomas-Hoover Capillary melting point apparatus and are uncorrected. Unless otherwise noted, 1H and 13C NMR spectra were determined in DMSO- d6 on a Varian spectrometer operating at 300 and 75. 47 MHz, respectively. The chemical shift values are expressed in ppm relative to tetramethylsilane as internal standard. All mass spectra were obtained by chemical ionization using methane as the ionizing gas.

Thin-layer chromatography (TLC) was performed on Merck silica gel 60 F254plates with fluorescent indicator. Proportions of solvents used for TLC are by volume. All solvents and reagents were purchased from Aldrich Chemical Co. and Fisher Scientific and were used as received.

Table 5 shows a correlation of compound number, structure and compound name.

TABLE 5 Cmpd Structure Name 1. 0 2, 4-Diamino-6- hydroxypyrimidine H2N N NH2 2. H Guanidine H 2N NW 3.

a-Hydroxy acetone (7. 41 g, 0. 1 mol) was dissolved in dry Et3N (100 mL) and cooled to-20 °C. TsCl (21. 9 g, 0. 12 mol) in dry CH2C12 (50 mL) was added slowly in such a way that the temperature was maintained below-5 °C. After complete addition of TsCl, the reaction mixture was washed with water (200 mL), aq. 6N HC1 (100 mL), water (200 mL) and dried over MgS04. Filtration followed by evaporation of solvent afforded the crude product which was chromatographed on silica gel (20% EtOAc/Hexane) afforded the product as a colourless oil. Yield : 17. 7 g, 83% IH NMR (300 MHz, DMSO-d6) : 8 2. 20 (s, 3H), 2. 46 (s, 3H), 4. 49 (s, 2H), 7. 36-7. 82 (AA\'BB\' J=7. 8 Hz, 4H) 13C NMR (75 MHz, DMSO-d6) : # 21. 6, 26. 4, 71. 9, 127. 9, 130. 0, 132. 2, 145. 5, 201. 0 HRMS (CI) Calculated for CHXNO.

Zheng, L. ; Berridge, M. S. ; Ernsberger, P. J. Med. Chem., 1994, 37, 3219.

2-Amino-6-methyl-3,4-dihydro-4-oxo-7H-pyrrolo[2,3-d] pyrimidine (8-methyl-7-dea7. aguanine) (6) A mixture of 2, 6-diamino-4-hydroxypyrimidine (6. 30 g, 5 0 mmol) and NaOAc (4. 1 g, 50 mmol) in water (150 mL) was heated to 100 °C to form the homogeneous solution. a-tosyl acetone (4. 70 g, 50 mmol) was added all at once and stirred further at 100 °C for 5 hours.

The reaction mixture was cooled to 4 °C and filtered to afford the product as a pink coloured solid. The solid was heated to 100 °C under vacuum to afford buff colored solid. Yield : 5. 42 g, 66%, Rf : 0. 37 (20% MeOH/CHC13) IH NMR (300 MHz, DMSO-d6) : 8 2. 13 (d, J = 0. 9 Hz, 3H), 5. 82 (q, J = 0. 9 Hz, 1H), 5. 93 (br s, 1H), 10. 10 (br s, 1H), 10. 77 (br s, 1H). HRMS (CI) Calcd for C7H9N40 165. 0776 ; Found : 165. 0773. The spectral data are compatible with an authentic sample.

Mechanism :

The synthesis method is summarized in Figure 18. 4 and 5 were used as the ketone derivatives. Use of 5 resulted in a yield of 76% as opposed to a 66% when 4 was used as the ketone derivative.

2-amino-6-pher)ethy1-3,4-dihydro-4-oxo-7H-pyrroIo[2,3-d]p yrimidine Figure 19 presents the method used to synthesize 9.

Phenylpropanoyl chloride 7 was reacted with CH2N2 in the presence of gaseous hydrogen chloride and ether to produce 1-chloro-4-phenyl-2- butanone 8. 8 was reacted with 1 to produce 2-amino-6-phenethyl- 3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 9 at a 39% yield.

2, 4-diamino-4-noxy methylfuro [2, 3-d] pyrymidme (12) The synthesis of 12 is shown in Figure 20. Phenoxyacetyl chloride 10 was reacted with CH2N2 in the presence of gaseous

hydrogen chloride and ether to produce 1-chloro-3-phenoxy-2- propanone 11. 11 was reacted with 1 to produce 12 at a 39% yield.

2, 6-d1amniothia7.o1o[4,5-d]pynnndm-7(6H)-one(14) Figure 21 shows the method of synthesis of 14. 1 was reacted with Br2 in the presence of sodium acetate and acetic acid to produce 5-bromo-2, 6-diamino-4-oxo-pyrimidine 13. 13 was then reacted with CN2SH4 in ethanol to produce 14.

2-Amino-3,4-dihydro-4-oxo-7H-pyrrolo[2,3-d]-pyrimidine (7- deazaguanine) (1-5 (Davoll, 1960) was prepared by an improved, one-step procedure. Into a stirred solution of DMF (240 mL) and water (40 mL) at room temperature were added 2, 4-diamino-6- hydroxypyrimidine (12. 6 g, 0. 1 mol) and NaOAc (8. 2 g, 0. 1 mol).

After the reaction mixture became homogeneous, 50% aqueous a- chloroacetaldehyde (15. 7 g, 0. 1 mol) was added all at once. The reaction mixture was stirred at room temperature for 2 days and the solvent was evaporated under vacuum. The residue was triturated with 10 mL of water, cooled, and filtered. The solid was purified over silica gel column chromatography to afford 10. 8 g (72%) of a buff colored solid whose\'H and l3C NMR spectra matched those reported.

Synthesis of 2-amino-5-(4\'-carboxyphenylaminomethyl)-3,4-dihydro-



4-oxo-7H-pyrrolo [2, 3-d] Dyrimi dine (21)





2-amino-5- (4\'-carboxyphenylaminomethyl)-3, 4-dihydro- 4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 21 is an analog of PTA ; additionally it contains a 7-deazaquanine moiety. Figure 14A depicts

the structures of synthetic intermediates and of 21 that contain the 7- deazaguanine moiety.

Figure 22A details the synthesis of 21. 2-Amino-3, 4- dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 15 is reacted with octanoyl chloride in pyridine to produce 2-N, 7-dioctanoyl-3, 4-dihydro- 4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 16, which is then treated with ammonia in methanol to produce 2-N-octanoyl-3, 4-dihydro-4-oxo-7H- pyrrolo [2, 3-d] pyrimidine 17. 17 is reacted with (PhCH2) in formaldehyde to produce 2-N-octanoylamino-5- (N, N-dibenzyl methyl amino)-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 18. A suspension of 18 (220 mg, 0. 412 mmol) and methyl 4-aminobenzoate 19 (331. 2 mg, 2. 06 mmol) in 1. 5 ml of methanol : THF (1 : 1). The reaction was carried out in a sealed tube at 70° C for 24 hours with stirring. The reaction mixture was cooled to room temperature and evaporated to dryness. The residue was applied to a silica gel chromatography column and washed with 200 ml of chloroform : methanol (95 : 5). The product was eluted with chloroform : methanol (90 : 10) to yield 98 mg (54%) of 2-N- octanoylamino-5- (4\'-methylcarboxy phenylaminomethyl)-3, 4-dihydro- 4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 20 as an off-white solid.

20 can be further characterized chemically as follows : TLC Rf 0. 39 (CHCl3 : MeOH, 90 : 10) ; 1H NMR (DMSO-d6) __ 7. 66 (d, 10 Hz, 2H, C6H4), 6. 88 (s, 1H, 6-CH), 6. 77 (t, 6Hz, 1H, 2-NH), 6. 64 (d, 1 0 Hz, 2H, C6H4), 4. 39 (bd, 2H, 5-CCH2NH), 3. 72 (s, 3H, C02CH3), 2. 42 (t, 2H, C7HI5), 1. 56 (bt, 2H, C7Hl5), 1. 26 (bs, 8H, C7HI5), 0. 90 (bt, 3H, C7HIs) ; MS (CI+) 440. 14 (M+H+) ; HRMS (CI+) 439. 221475 (C23H29NsO4), Calculated Mass 439. 221955.

To a solution of 20 (103 mg, 0. 234 mmol) in 500 u ; l of EtOH and 2 ml of H2 O was added 5 M KOH (117. 25 1, 0. 586 mmol).

This mixture was stirred for 48 hours, at which time the reaction became homogeneous. The solution was evaporated and the residue taken up in H2O and the pH adjusted to 7. The H2O was removed in vacuo and the residue was taken up in MeOH. The mixture was centrifuged and the MeOH was removed. More MeOH was added ; this process was carried out three times. The MeOH fractions were combined, evaporated and the residue applied to a flash silica gel chromatography column and eluted with CHCl3 : MeOH : AcOH (90 : 9 : 1).

6. 2 mg (95) of 21 was isolated as a yellow solid. 21 can be further chemically characterized : TLC Rf 0. 09 (CHCl3 : MeOH : AcOH, 90 : 9 : 1) ; lH NMR (D2O) _ 7. 64 (d, 10 Hz, 2H, C6H4), 6. 79 (d, 10 Hz, 2H, C6H4), 6. 61 (s, 1H, 6-CH). and 2-amino-5- (4\'-armnobenzy1phenylamino) methy1-3, 4-d1hydro-4- oxo-7H-pyrrolo [2, 3-d] pyrimidine (24) The synthesis of 24 is shown in Figure 22B. The synthesis of 24 follows the same steps as the synthesis of 21 except that 18 is reacted with 4, 4\'-methylene dianiline 22 to form 2-N-octanoylamino- 5- (4\'-aminobenzyl phenylamino) methyl-3, 4-dihydro-4-oxo-7H- pyrrolo [2, 3-d] pyrimidine 23. 23 is then reacted with 5M KOH to form 2 4.

5-Amino-2-methyl-6H-oxazolo[5,4-d]pyrimidin-7-one (9-methyl- oxoguanine) (32) Figure 23 depicts the synthesis of 32. Sodium pieces (0. 575 g, 25 mmol) were added to ethanol (50 mL) under argon with stirring. After all the sodium had reacted, guanidine hydrocholride 2 (2. 38 g, 25 mmol) was added all at once, and the reaction mixture was stirred further for 10 m. The reaction mixture was heated to reflux

and during this time diethylacetamido malonate 30 (Zambito and Howe, 1973) (5. 43 g, 25 mmol) was added. After 2 h, ethanol (50 mL) was added, and the reaction mixture was refluxed for an additional 12 h. The reaction mixture was allowed to cool and was filtered. The residue was washed with ethanol and chloroform and was dried. The resulting solid was dissolved in water (25 mL) and precipitated by the addition of 50% hydrochloric acid. The solution was filtered, and the solid was washed with water and acetone and then dried under vacuum to afford 5-acetylamino-2-amino-4, 6- dihydroxy pyrimidine 31 (Tayor and Cain, 1949) as a white solid (4. 1 g, 89%) ; mp : >300 °C ;\'H NMR (300 MHz, D20) : 8 1. 88 (s) ; 13C NMR (75 MHz, DMSO-d6) : 8 21. 9 (q), 92. 6 (s), 154. 5 (s), 168. 2 (s), 174. 9 (s) ; HRMS (CI) Calcd for C6HgN403 185. 0675 ; Found : 185. 0668. The 5- acetylamino-2-amino-4, 6-dihydroxypyrimidine 31 (184 mg, 1 mmol) was taken in a sublimation apparatus and was heated slowly to 280 to 300 °C under vacuum for 14 h. The resulting sublimed solid was dissolved in methanol and purified by silica gel chromatography to afford the product 32 as a pale yellow solid (22. 7 mg, 14%). mp : > 300 °C ; IH NMR : 8 2. 38 (s, 3H), 6. 85 (br s, 2H), 11. 0 (br s, 1H) ; 13C NMR : 8 13. 8 (q), 111. 2 (s), 154. 6 (s), 155. 6 (s), 156. 1 (s), 166. 4 (s) ; HRMS Calcd for C6H7N402 167. 0569 ; Found : 167. 0573 EXAMPLE 11 Guanine-based compound x-ray data collection Diffraction data were collected on a Rigaku Raxis IV area detector. X-rays were generated by a Rigaku RU200 rotating anode generator operated at 50 mV, 100 mA. Data were collected at room

temperature for all crystals except RTA complexed with 8-methyl-9- oxoguanine. That crystal was soaked for 5 min in mother liquor containing 15 to 25% PEG 8000 and transferred to paratone. After removal of surface-bound water, the crystal was flash-frozen in liquid nitrogen, and mounted in a nitrogen cold stream. The program DENZO was used to process the diffraction images and SCALEPACKwas used to reduce the data (Otwinowski and Minor, 1997). The starting model used for molecular replacement was the monoclinic crystal form of recombinant RTA (Mlsna et al, 1993). The RTA2, 5-diamino- 4, 6-dihydroxypyrimidine co-crystal was triclinic with a= 42. 07A, b=65. 95 A, c=50. 81 A, a=93. 90, =112. 4, andy=84. 43. A molecular replacement solution was found with 2 monomers in the asymmetric unit using EPMR (Kissinger and Gehlhaar, Agouron Pharmaceuticals, Inc.). XPLOR was used for model refinement and generating Fourier maps with amplitudes Fobserved (Fo) minus Fcalculated (Fc) and 2Fo- Fc (Brunger, 1992). O was used for model building (Jones, 1991).

EXAMPLE 12 GnanmR-based compound structure deternnnatirm Soaking RTA crystals with 8-methyl-9-oxoguanine caused them to crack. Useful crystals were formed by crystallizing the enzyme in the presence of the inhibitor. In addition, only cryocooled co-crystals survived exposure to X-rays. The crystals were not perefectly isomorphous with native RTA crystals, and their diffraction was limited to 3.1 Å, suggesting poor crystal order. An electron density map is shown in Figure 24.

The binding orientation of 8-methyl-9-oxoguanine in the RTA active site was unexpected. Although many of the anticipated

pterin-like hydrogen bonds were formed, they were made in novel ways. This is illustrated in Figure 25 where the orientation of 8- methyl-9-oxoguanine is compared with the pterin ring of PTA. The orientation of 8-methyl-9-oxoguanine can be roughly generated from PTA by flipping the bicyclic ring system 180° around the horizontal axis and then rotating it 60° around an axis normal to the page. In both cases the exocyclic amine group sits in a pocket where it donates hydrogen bonds to the carbonyl groups of Val 81 and Gly 121. At the same time a ring nitrogen accepts a hydrogen bond from the amido group of Val 81. These are key interactions also preserved in the binding of adenine substrate and formycin (Monzingo and Robertus, 1992). As has been previously noted, the low energy tautomer of pterin has a hydrogen on N3, but this would clash with the amido hydrogen of Val 81. As a consequence, RTA binds the pterin (3) tautomer with the hydrogen on N1, which it donates in a hydrogen bond to the carbonyl oxygen of Gly 121, as shown in Figure 25 The cost of this rearrangement is more than compensated for by the formation of the hydrogen bond (Yan et al., 1998). In binding 8- methyl-9-oxoguanine no such promotion is required since the low energy form has the hydrogen already at N1 (using guanine nomenclature which differs from pterin nomenclature).

The novel orientation of 8-methyl-9-oxoguanine also brings the carbonyl oxygen into a position to accept a hydrogen bond from the amido group of Tyr 123, an interaction not seen in PTA. In PTA, Arg 180 donates two hydrogen bonds to the exocyclic oxygen at position 4 and to the ring nitrogen at position 5. The 8-methyl-9- oxoguanine orientation requires the side chain of Arg 180 to rotate slightly and allows only one bond to be formed to the ring oxygen a t position 9. It would seem that if 8-methyl-9-oxoguanine were

oriented like PTA, that the same number of hydrogen bonds could be made as are observed in the flipped and rotated mode. However, the strength of these interactions is difficult to estimate and it can only be assumed that the new orientation is energetically more favorable by a subtle amount.

The X-ray structures of the 7-deazaguanine and 9- deazaguanine complexes with RTA were determined to see if they would bind like the pterins or like 8-methyl-9-oxoguanine. Despite their differing pattern of hydrogen bonding groups, both compounds adopted a conformation like 8-methyl-9-oxoguanine in the crystal.

The binding is shown for 9-deazaguanine in Figure 26. In this structure two water molecules make additional interactions with the inhibitor which were not present in the 8-methyl-9-oxoguanine structure. In the 7-deazaguanine complex (not shown) N9 of 7- deazaguanine donates a hydrogen to the backbone carbonyl of Ala 79, an interaction not seen in other ligand complexes with RTA.

EXAMPLE 13 Crystallization of RTAinhihitor comnlexes RTA monoclinic crystals were grown for the purpose of studying inhibitor binding to the enzyme (Robertus et al., 1987 ; Mlsna et al., 1993). Briefly, 30 u, l of protein at 2. 2 to 2. 4 mg/ml in mother liquor (75 mM Tris-HCl pH 8. 9, 10 mM BME, 1 mM EDTA, 4. 1 % polyethylene glycol MW 8000 (PEG 8000)) was suspended over 200 RI of mother liquor at 4 °C. Microseeds were transferred to the sitting drop for nucleation. The compounds 5-amino-2-methyl-6H- oxooxazolo [5, 4d] pyrimidin-7-one (8-methyl-9-oxoguanine) and 2, 5-

diamino-4, 6-dihydroxypyrimidine were co-crystallized under native RTA monoclinic crystal conditions with a final ligand concentration of 4 mM and 10 mM respectively. The compound 7-deazaguanine was solubilized in 0. 033 N NaOH. Five gel of the ligand solution was added to a 25 ul RTA crystal drop for a 4 day soak. The final concentration of ligand was estimated at 4 mM. The inhibitor 9DG was co- crystallized under RTA monoclinic conditions. Five jil of 9DG in 0. 045 N NaOH was added to 25 gel of mother liquor to give an estimated 4 mM final concentration of ligand.

EXAMPLE 14 Inhibition of RTA by guanine-based compounds The potency of these guanine-based compounds to inhibit RTA action against ribosomes was tested. The Artemia salinas ribosome in vitro protein synthesis assay as described supra was used (see Example 1) to determine efficacy of these compounds as inhibitors of RTA.

Guanine-like compounds in the correct orientation can inhibit ricin. Energy calculations suggest that guanine, rotated 180° around the long center line of the purine base with respect to the adenine orientation, could bind and inhibit RTA (Yan et al, 1998). As shown in Table 6 guanine has an ICso of 0. 9 mM. Xanthine as an inhibitor has an IC50 of 3. 6 mM, confirming the importance of the exocyclic amine group at position 2 of the base. This amine group is expected to donate hydrogen bonds to the carbonyl oxygens of Val 8 1 and Gly 121. Subsequent inhibitor candidates would all preserve the

pyrimidine ring structure seen in guanine (that is 2-amino-6-hydroxy- pyridine).

Table 6 shows the ICso values for a number of single ring compounds and of bi-cyclic compounds that retained the 6-membered ring of guanine. These include 8-mercaptoguanine, thiazolpyrimidine, 9-deazagunaine, 7-deazaguanine, 8-methyl-7-deazaguanine, and 8- methyl-9-oxoguanine. A plot of RTA activity versus inhibitor concentration is shown in Figure 15. Figure 16 shows a stereogram representing the binding of 7-deazaguanine to RTA. The IC50 values for these compounds range from 2. 8 mM for 7-deazaguaninem 0. 4 mM for 8-methyl-9-oxoguanine and 0. 25 mM for. The latter compound is quite soluble in water and exhibits the best ICso of any compound we have measured. This compound has the potential to be the main molecular platform for other inhibitor designs.

21 was also tested and it was found to be an inhibitor of RTA. Preliminary tests suggest that the In50 for this inhibitor is about 0. 7 mM. This implies that a 7-deazaguanine moiety is a better framework than the pterins on which to construct RTA and STA inhibitor compounds. It is also possible to make longer pendent versions of 21 that will go into the"second pocket". That is, in addition to nonpolar interactions, one can make specific polar interactions.

TABLE 6 RTA Inhibitors. 2D Structure Name ICso (RTA) _ N Pteroic acid 0. 6 mM HOOO 9 J _ Fi HO 5-diamino-4, 6- 2. 2 mM _ dihydroxypyrimidine Y H 0 0 2, 4-diamino-6- 3. 4 mM _ hydroxypyrimidine H2N1N NH2 0 2-amino-4-hydroxy-6-2. 5 mM H-\'I methylpyrimidine H2N1NA OH 4-amino-6-hydroxy-2-3. 2 mM _ methylpyrimidine N NH2 3-amino-pyridothiadiazine No inhibition S., 2-amino-6-hydroxypurine 0. 9 mM (guanine) G9 I NH N NH N o 2, 6-dihydroxypurine 3. 6 mM _ (xanthine) r H T N N NH2 2-amino-6-hydroxy-8-0. 56 mM N HO mercaptopurine g (8-mercaptoguanine) . _ Thiazolopyrimidine 2. 0 mM _ tS\'\IH v 2-amino-3, 4-dihydro-4-oxo-1. 4 mM 9H-pyrrolo [2, 3-d]-pyrimidine q NH (9-deazaguanine) H ,,, 2 2-amino-3, 4-dihydro-4-oxo- 2. 8 mM 7H-pyrrolo [2, 3-d]-pyrimidine (7-deazaguanine) o 2-amino-6-methyl-3, 4-2. 1 mM dihydro-4-oxo-7H-pyrrolo [2, 3- WNH d]-pyrimidine (8-methyl-7-deazaguanine) 0 0 N NH2 5-amino-2-methyl-6H-0. 40 mM _ oxooxazolo [5, 4d] pyrimidin-7- u JtL NH one (8-methyl-9-oxoguanine) N N-4- [2- (2-amino-4 (3H)-oxo-0. 60 mM 7H-pyrrolo 7H-pyrrolo [2, 3-d] pyrimidin-6- yl) ethyl] benzoic acid 0 OOH 5- (4- 0. 70 mM _ carboxyphenylaminomethyl) 2- octanoylaminopyrrolo [2, 3- -i d] rimidin-4 (3H)-one 2-amino-5- (4\'- 0. 25 mM carboxyphenylaminomethyl)- H2N1SNJH 4-dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine

The results suggest that 2, 4-diamino-6-hydroxypyrimidine 19 2-amino-6-methyl-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 6, 2, 6-diaminothiazolo [4, 5-d] pyrimidin-7 (6H)-one 14, 2-amino-5- (4\'-carboxyphenylaminomethyl)-3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3- d] pyrimidine 21, 2-amino-5- (4\'-aminobenzylphenylamino) methyl- 3, 4-dihydro-4-oxo-7H-pyrrolo [2, 3-d] pyrimidine 24, 2-amino-4, 6-

dihydroxy-5-methylpyrimidine 2 5, 2-amino-4-hydroxy-6- methylpyrimidine 26, 4-amino-6-hydroxy-2-methylpyrimidine 27, 2- amino-6-hydroxy-8-mercaptopurine 28 and 5-amino-2-methyl-6H- oxooxazolo [5, 4d] pyrimidin-7-one 32 would be effective ricin A chain inhibitors.

EXAMPLE 15 Design of additional ricin A chain (RTA) inhibitors based on guanin or structurally similar to guanine/RTA complexes The structural information obtained from the pteroic acid/ricin A chain complex was used to design additional molecules that exploit the pterin binding pocket consisting in part of residues G121, V81, R180 to design additional RTA inhibitors. PTA bound RTA X-ray structure revealed that a more favorable interaction between RTA and the inhibitor could be made if the 6-membered pyrazine ring of the PTA were replaced with a 5-membered ring. Most of these molecules were derivative of guanine or structurally similar molecules.

It is contemplated that pendant groups to a 9-oxoguanine platform can be added. Based on the PTA structure, and earlier predictions of guanine binding to RTA, the site of the PABA moiety linkage to pterin would superimpose with position 8 of guanine. It would seem sterically favored to add a linker on position 8 rather than position 7, given the predicted close proximity of Arg 180 to position 7. However, 8-methyl-9-oxoguanine, 7-deazaguanine, and 9- deazaguanine all bind in an orientation that places position 8 in Van der Waals contact with Arg 180. Interestingly, the ICso values for the

near isomeric compounds N-4- [2- (2-amino-4 (3H)-oxo-7H- pyrrolo [2, 3-d] pyrimidin-6-yl) ethyl] benzoic acid and 5- (4- carboxyphenylaminomethyl) 2-octanoylaminopyrrolo [2, 3- d] pyrimidine-4 (3H)-one shown in TABLE ? are 0. 6 mM and 0. 7 mM. The fact that both molecules bind with essentially the same affinity suggests both position 7 and 8 are feasible for further design efforts. However, if one considers the ICso data, PTA binding, the predicted binding of guanine-based compounds, and the observed binding of 8- methyl-9-oxoguanine, 7-deazaguanine, and 9-deazaguanine to RTA as a whole, it suggests that guanine-based molecules appended at position 8 are likely to revert to the PTA binding mode.

Any patents or publications mentioned in this specification are indicative of the levels of those skilled in the art to which the invention pertains. Further, these patents and publications are incorporated by reference herein to the same extent as if each individual publication was specifically and individually indicated to be incorporated by reference.

One skilled in the art will appreciate readily that the present invention is well adapted to carry out the objects and obtain the ends and advantages mentioned, as well as those objects, ends and advantages inherent herein. The present examples, along with the methods, procedures, treatments, molecules, and specific compounds described herein are presently representative of preferred embodiments, are exemplary, and are not intended as limitations on the scope of the invention. Changes therein and other uses will occur to those skilled in the art which are encompassed within the spirit of the invention as defined by the scope of the claims.