# Patent

## Title
Methods of using the structure coordinates of molecules comprising an IMPDH-like binding pocket

## Abstract

The present invention relates to a data storage medium encoded with the corresponding structure coordinates of molecules and molecular complexes which comprise the active site binding pockets of IMPDH. Such data storage material is capable of displaying such molecules and molecular complexes, or their structural homologues, as a graphical three-dimensional representation on a computer screen. This invention also relates to methods of using the structure coordinates to solve the structure of homologous proteins or protein complexes. In addition, this invention relates to methods of using the structure coordinates to screen and design compounds, including inhibitory compounds, that bind to IMPDH or homologues thereof. This invention also relates to molecules and molecular complexes which comprise the active site binding pockets of IMPDH or close structural homologues of the active site binding pockets. This invention also relates to compounds and pharmaceutical compositions which are inhibitors of IMPDH.

## Background
This application is a divisional of U.S. application Ser. No. 08/640,164, filed Apr. 30, 1996, now U.S. Pat. No. 6,128,582.


TECHNICAL FIELD OF INVENTION
The present invention relates to a data storage medium encoded with the structural coordinates of crystallized molecules and molecular complexes which comprise the active site binding pockets of IMPDH. Such data storage material is capable of displaying such molecules and molecular complexes, or their structural homologues, as a graphical three-dimensional representation on a computer screen. This invention also relates to methods of using the structure coordinates to solve the structure of homologous proteins or protein complexes. In addition, this invention relates to methods of using the structure coordinates to screen and design compounds, including inhibitory compounds, that bind to IMPDH or homologues thereof. This invention also relates to molecules and molecular complexes which comprise the active site binding pockets of IMPDH or close structural homologues of the active site binding pockets. This invention also relates to compounds and pharmaceutical compositions which are inhibitors of IMPDH.
BACKGROUND OF THE INVENTION
The synthesis of nucleotides in organisms is required for the cells in those organisms to divide and replicate. Nucleotide synthesis in mammals may be achieved through one of two pathways: the de novo synthesis pathway or the salvage pathway. Different cell types use these pathways to a different extent.
Inosine-5′-monophosphate dehydrogenase (IMPDH; EC 1.1.1.205) is an enzyme involved in the de novo synthesis of guanosine nucleotides. IMPDH catalyzes the NAD-dependent oxidation of inosine-5′-monophosphate (IMP) to xanthosine-5′-monophosphate (XMP), which is the committed step in guanosine nucleotide synthesis. [R. C. Jackson et. al., Nature, 256, pp. 331–333 (1975)].
IMPDH is ubiquitous in eukaryotes, bacteria and protozoa [Y. Natsumeda and S. F. Carr, Ann. N.Y. Acad. Sci., 696, pp. 88–93 (1993). The prokaryotic forms share 30–40% sequence identity with the human enzyme. Regardless of species, the enzyme follows an ordered Bi—Bi reaction sequence, where IMP binding precedes that of NAD, and NADH is released prior to XMP [S. F. Carr et al., J. Biol. Chem. 268, pp. 27286–27290 (1993); E. W. Holmes; Biochim. Biophys. Acta 364, pp. 209–217 (1974)]. This mechanism differs from that of most other known NAD-dependent dehydrogenases, which have either a random order of substrate addition or require that NAD bind before substrate.
The de novo synthesis of guanosine nucleotides, and thus the activity of IMPDH, is particularly important in B- and T-lymphocytes. These cells depend on the de novo, rather than salvage pathway to generate sufficient levels of nucleotides necessary to initiate a proliferative response to mitogen or antigen [A. C. Allison et. al., Lancet II, 1179 (1975) and A. C. Allison et. al., Ciba Found. Symp., 48, 207 (1977)]. Thus, IMPDH is an attractive target for selectively inhibiting the immune system without also inhibiting the proliferation of other cells.
Immunosuppression has been achieved by inhibiting a variety of enzymes including for example, the phosphatase calcineurin (inhibited by cyclosporin A and FK-506); dihydroorotate dehydrogenase, an enzyme involved in the biosynthesis of pyrimidines (inhibited by leflunomide and brequinar); the kinase FRAP (inhibited by rapamycin); and the heat shock protein hsp70 (inhibited by deoxyspergualin). [See B. D. Kahan, Immunological Reviews, 136, pp. 29–49 (1993); R. E. Morris, The Journal of Heart and Lung Transplantation, 12(6), pp. S275–S286 (1993)].
Inhibitors of IMPDH are also known. U.S. Pat. Nos. 5,380,879 and 5,444,072 and PCT publications WO 94/01105 and WO 94/12184 describe mycophenolic acid (MPA) and some of its derivatives as potent, uncompetitive, reversible inhibitors of human IMPDH type I and type II. MPA has been demonstrated to block the response of B and T-cells to mitogen or antigen [A. C. Allison et. al., Ann. N.Y. Acad. Sci., 696, 63 (1993).




Immunosuppressants, such as MPA, are useful drugs in the treatment of transplant rejection and autoimmune diseases [R. E. Morris, Kidney Intel., 49, Suppl. 53, S-26 (1996)]. MPA, however, is characterized by undesirable pharmacological properties, such as gastrointestinal toxicity and poor bioavailability. [L. M. Shaw, et. al., Therapeutic Drug Monitoring, 17, pp. 690–699, (1995)].
Nucleoside analogs such as tiazofurin, ribavirin and mizoribine also inhibit IMPDH [L. Hedstrom et. al., Biochemistry, 29, pp. 849–854 (1990)]. These compounds, however, are not specific for IMPDH.
Mycophenolate mofetil, a prodrug which quickly liberates free MPA in vivo, was recently approved to prevent acute renal allograft rejection following kidney transplantation [L. M. Shaw et al., Therapeutic Drug Monitoring, 17, pp. 690–699 (1995); H. W. Sollinger, Transplantation, 60, pp. 225–232 (1995)]. However, because of gastrointestinal and other side-effects, the therapeutic potential of this drug appears limited [L. M. Shaw et al., Therapeutic Drug Monitoring, 17, pp. 690–699 (1995); A. C. Allison and E. M. Eugui Immunological Rev., 136, pp. 5–28 (1993)].
It is also known that IMPDH plays a role in other metabolic events. Increased IMPDH activity has been observed in rapidly proliferating human leukemic cell lines and other tumor cell lines, indicating IMPDH is a target for anti-cancer as well as immunosuppressive chemotherapy [M. Nagai et. al., Cancer Res., 51, pp. 3886–3890 (1991)].
IMPDH has also been shown to play a role in the proliferation of smooth muscle cells, indicating that inhibitors of IMPDH, such as MPA or rapamycin, may be useful in preventing restenosis or other hyperproliferative vascular diseases [C. R. Gregory et al., Transplantation, 59, pp. 655–61 (1995); PCT publication WO 94/12184; and PCT publication WO 94/01105].
Additionally, IMPDH has been shown to play a role in viral replication in some viral cell lines [S. F. Carr, J. Biol. Chem., 268, pp. 27286–27290 (1993)]. Analogous to lymphocyte and tumor cell lines, the implication is that the de novo, rather than the salvage, pathway is critical in the viral replication process.
Thus, there remains a need for potent IMPDH inhibitors with improved pharmacological properties. Such inhibitors would have therapeutic potential as immunosuppressants, anti-cancer agents, anti-vascular hyperproliferative agents and anti-viral agents. Specifically, such compounds may be used in the treatment of transplant rejection and autoimmune diseases, such as rheumatoid arthritis, multiple sclerosis, juvenile diabetes, asthma, inflammatory bowel disease, as well as in the treatment of cancer and tumors, such as lymphomas and leukemia, vascular diseases, such as restenosis, and viral replication diseases, such as retroviral diseases and herpes.
Two isoforms of human IMPDH, designated type I and type II, have been identified [F. R. Collart and E. Huberman, J. Biol. Chem., 263, pp. 15769–15772 (1988); Y. Natsumeda et. al., J. Biol. Chem., 265, pp. 5292–5295 (1990)]. Each is 514 amino acids, and they share 84% sequence identity. Both IMPDH type I and type II form active tetramers in solution, with subunit molecular weights of 56 kDa [Y. Yamada et. al., Biochemistry, 27, pp. 2737–2745 (1988)].
Both human IMPDH isoforms have been characterized by their cDNA and amino acid sequences [Y. Natsumeda et. al., J. Biol. Chem., 265, pp. 5222–5295 (1990)]. Chinese hamster IMPDH has been characterized by its cDNA and its amino acid sequence [F. R Collart and E. Hubermann, J. Biol. Chem., 263, pp. 15769–15772 (1988). Knowledge of the primary structure, i.e., amino acid sequence, of IMPDH, however, does not allow prediction of its tertiary structure. Nor does it afford an understanding of the structural, conformational and chemical interactions of IMPDH with MPA, IMP, or other compounds or inhibitors.
The crystal structure of IMPDH has not been reported. Nor has the crystal structure of a IMPDH homologue or a IMPDH co-complex been reported. The need, therefore, exists for determining the crystal structure of IMPDH to provide a more accurate description of the structure of IMPDH to aid in the design of improved IMPDH inhibitors. The crystal structure of a complex comprising IMPDH, IMP and MPA would provide such a description.
SUMMARY OF THE INVENTION
Applicants have solved this problem by achieving, for the first time, the crystallization of Chinese hamster type II IMPDH in complex with XMP* and MPA and have solved the three-dimensional structure of that complex. This isoform of IMPDH differs by only 6 amino acids from the human type II enzyme while retaining similar enzymatic characteristics and the ability to be inhibited by MPA. Solving this crystal structure has allowed applicants to determine the key structural features of IMPDH, particularly the shape of its XMP* binding pocket and its MPA binding pocket.
Thus, the present invention provides molecules or molecular complexes that comprise all or parts of either one or two of these binding pockets, or homologues of these binding pockets that have similar three-dimensional shapes.
The invention also provides a machine readable storage medium which comprises the structure coordinates of IMPDH, including all or any parts of the XMP* and MPA binding pockets. Such storage medium encoded with these data are capable of displaying on a computer screen or similar viewing device, a three-dimensional graphical representation of a molecule or molecular complex which comprises such binding pockets or similarly shaped homologous binding pockets.
The invention also provides methods for designing, evaluating and identifying compounds which bind to all or parts of the aforementioned binding pockets. Such compounds are potential inhibitors of IMPDH or its homologues.
The invention also provides novel classes of compounds, and pharmaceutical compositions thereof, that are useful as inhibitors of IMPDH or its homologues.
The invention also provides a method for determining at least a portion of the three-dimensional structure of molecules or molecular complexes which contain at least some structurally similar features to IMPDH. This is achieved by using at least some of the structural coordinates obtained for the Chinese hamster type II IMPDH/XMP*/MPA complex.
The invention also provides a method for crystallizing a Chinese hamster type II IMPDH/XMP*/MPA complex and related complexes.
The invention also provides a crystalline form of a Chinese hamster type II IMPDH/XMP*/MPA complex and related complexes.



BRIEF DESCRIPTION OF THE FIGURES
FIG. 1 lists the atomic structure coordinates for Chinese hamster IMPDH type II (hereafter referred to as IMPDH) in complex with XMP* and MPA as derived by X-ray diffraction from a crystal of that complex. The following abbreviations are used in FIG. 1:
“Atom type” refers to the element whose coordinates are measured. The first letter in the column defines the element.
“X, Y, Z” crystallographically define the atomic position of the element measured.
“B” is a thermal factor that measures movement of the atom around its atomic center.
“Occ” is an occupancy factor that refers to the fraction of the molecules in which each atom occupies the position specified by the coordinates. A value of “1” indicates that each atom has the same conformation, i.e., the same position, in all molecules of the crystal.
FIG. 2 is a stereo ribbon diagram depicting the fold and conformation of IMPDH in three-dimensions as determined by X-ray crystallography.
FIG. 3 is a topology diagram of the IMPDH fold.
FIG. 4 is a ribbon drawing of the IMPDH tetramer.
FIG. 5 depicts a stereo view of the conformation and electron densities of XMP* and MPA and some of the IMPDH active site amino acid residues surrounding these ligands.
FIG. 6 schematically depicts the XMP*-IMPDH interactions.
FIG. 7 schematically depicts the MPA-IMPDH interactions.
FIG. 8 shows the specific activity of IMPDH mutants for IMP substrates and the inhibition of the mutants by MPA.
FIG. 9 shows a diagram of a system used to carry out the instructions encoded by the storage medium of FIGS. 10 and 11.
FIG. 10 shows a cross section of a magnetic storage medium.
FIG. 11 shows a cross section of a optically-readable data storage medium.



DETAILED DESCRIPTION OF THE INVENTION
The following abbreviations are used throughout the application:















A = Ala = Alanine



V = Val = Valine



L = Leu = Leucine



I = Ile = Isoleucine



P = Pro = Proline



F = Phe = Phenylalanine



W = Trp = Tryptophan



M = Met = Methionine



G = Gly = Glycine



S = Ser = Serine



T = Thr = Threonine



C = Cys = Cysteine



Y = Tyr = Tyrosine



N = Asn = Asparagine



Q = Gln = Glutamine



D = Asp = Aspartic Acid



E = Glu = Glutamic Acid



K = Lys = Lysine



R = Arq = Arginine



H = His = Histidine










IMPDH=Chinese hamster type II inosine monophosphate dehydrogenase
IMP=inosine monophosphate
MPA=mycophenolic acid
NAD=nicotinamide adenine dinucleotide
KCl=potassium chloride
BME=2-mercaptoethanol
EDTA=ethylenediaminetetracetic acid
Tris=tris(hydroxymethyl)aminomethane
PEG=polyethylene glycol
LiCl=lithium chloride
MES=morpholinoethyl sulfonic acid
MeP=1-methyl-2-pyrrolidinone.
PMSF=phenylmethylsulfonyl fluoride
XMP*=a form of IMP in which the 2-position hydrogen has been replaced with a covalent bond to Cys 331 of IMPDH.
IMPDH/IMP/MPA=IMPDH in complex with IMP and MPA
IMPDH/XMP*/MPA=IMPDH in complex with XMP* and MPA
Additional definitions are set forth in the specification where necessary.
In order that the invention described herein may be more fully understood, the following detailed description is set forth.
Applicants have provided crystals comprising IMPDH/XMP*/MPA, which are suitable for X-ray crystallographic analysis. The IMPDH/XMP*/MPA complex is derived from an IMPDH/IMP complex in which a covalent bond has formed between the C2 carbon of IMP and the sulfur atom of Cys 331 (FIGS. 5–7), yielding an oxidized IMP thioimidate intermediate (XMP*) which is trapped by MPA.
Applicants have solved the three-dimensional structure of the IMPDH/XMP*/MPA complex using high resolution X-ray crystallography. Thus, in one embodiment of this invention is provided an IMPDH/XMP*/MPA crystal. Preferably, the crystal has tetragonal space group P4. More preferably, the crystal comprises rectangular shaped unit cells, each unit cell having the dimensions a=b=110.6±5 Å, and c=111.0±5 Å. Most preferably, the crystallized enzyme is a tetramer.
Importantly, applicants' invention has provided, for the first time, information about the shape and structure of the XMP* and MPA active site binding pockets of IMPDH.
Binding pockets are of significant utility in fields such as drug discovery. The association of natural ligands or substrates with the binding pockets of their corresponding receptors or enzymes is the basis of many biological mechanisms of action. Similarly, many drugs exert their biological effects through association with the binding pockets of receptors and enzymes. Such associations may occur with all or any parts of the binding pocket. An understanding of such associations will help lead to the design of drugs having more favorable associations with their target receptor or enzyme, and thus, improved biological effects. Therefore, this information is valuable in designing potential inhibitors of IMPDH-like binding pockets.
The term “binding pocket”, as used herein, refers to a region of a molecule or molecular complex, that, as a result of its shape, favorably associates with another chemical entity or compound.
The term “IMPDH-like binding pocket” refers to a portion of a molecule or molecular complex whose shape is sufficiently similar to all or any parts of the XMP* and MPA active site binding pockets of IMPDH as to bind common ligands. This commonality of shape is defined by a root mean square deviation from the structure coordinates of the backbone atoms of the amino acids that make up the binding pockets in IMPDH (as set forth in FIG. 1) of not more than 1.5 Å. How this calculation is obtained is described below.
The “active site binding pockets” or “active site” of IMPDH refers to the area on the IMPDH enzyme surface where conversion of IMP to XMP occurs, and where MPA exerts its inhibitory effect. In resolving the crystal structure of Chinese hamster IMPDH type II, applicants have determined that IMPDH amino acids 67, 68, 69, 70, 73, 93, 273, 274, 275, 276, 277, 303, 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 337, 339, 340, 364, 365, 366, 367, 368, 385, 386, 387, 388, 389, 391, 411, 412, 413, 414, 415, 416, 419, 420, 439, 440, 441, 442, 443, 469, 470, 500, 501, 502, 503, 504, 505, and 506, are situated close enough to either XMP* or MPA (within 7 Å) to interact with these ligands. It will be readily apparent to those of skill in the art that the numbering of amino acids in other isoforms of IMPDH may be different than that isolated from Chinese hamster.
Each of those amino acids is defined by a set of structure coordinates as set forth in FIG. 1. The term “structure coordinates” refers to Cartesian coordinates derived from mathematical equations related to the patterns obtained on diffraction of a monochromatic beam of X-rays by the atoms (scattering centers) of an IMPDH complex in crystal form. The diffraction data are used to calculate an electron density map of the repeating unit of the crystal. The electron density maps are then used to establish the positions of the individual atoms of the IMPDH enzyme or enzyme complex.
Those of skill in the art understand that a set of structure coordinates for an enzyme or an enzyme-complex or a portion thereof, is a relative set of points that define a shape in three dimensions. Thus, it is possible that an entirely different set of coordinates could define a similar or identical shape. Moreover, slight variations in the individual coordinates will have little effect on overall shape. In terms of binding pockets, these variations would not be expected to significantly alter the nature of ligands that could associate with those pockets.
The term “associating with” refers to a condition of proximity between a chemical entity or compound, or portions thereof, and an IMPDH molecule or portions thereof. The association may be non-covalent—wherein the juxtaposition is energetically favored by hydrogen bonding or van der Waals or electrostatic interactions—or it may be covalent.
The variations in coordinates discussed above may be generated because of mathematical manipulations of the IMPDH/XMP*/MPA structure coordinates. For example, the structure coordinates set forth in FIG. 1 could be manipulated by crystallographic permutations of the structure coordinates, fractionalization of the structure coordinates, integer additions or subtractions to sets of the structure coordinates, inversion of the structure coordinates or any combination of the above.
Alternatively, modifications in the crystal structure due to mutations, additions, substitutions, and/or deletions of amino acids, or other changes in any of the components that make up the crystal could also account for variations in structure coordinates. If such variations are within an acceptable standard error as compared to the original coordinates, the resulting three-dimensional shape is considered to be the same. Thus, for example, a ligand that bound to the active site binding pocket of IMPDH would also be expected to bind to another binding pocket whose structure coordinates defined a shape that fell within the acceptable error. Such modified complexes or the binding pocket(s) thereof are also within the scope of this invention.
Various computational analyses are therefore necessary to determine whether a molecule or the binding pocket portion thereof is sufficiently similar to all or parts of the IMPDH binding pockets described above. Such analyses may be carried out in current software applications, such as the Molecular Similarity application of QUANTA (Molecular Simulations Inc., San Diego, Calif.) version 4.1, and as described in the accompanying User's Guide.
The Molecular Similarity application permits comparisons between different structures, different conformations of the same structure, and different parts of the same structure. The procedure used in Molecular Similarity to compare structures is divided into four steps: 1) load the structures to be compared; 2) define the atom equivalences in these structures; 3) perform a fitting operation; and 4) analyze the results.
Each structure is identified by a name. One structure is identified as the target (i.e., the fixed structure); all remaining structures are working structures (i.e., moving structures). Since atom equivalency within QUANTA is defined by user input, for the purpose of this invention we will define equivalent atoms as protein backbone atoms (N, Ca, C and O) for all conserved residues between the two structures being compared. We will also consider only rigid fitting operations.
When a rigid fitting method is used, the working structure is translated and rotated to obtain an optimum fit with the target structure. The fitting operation uses an algorithm that computes the optimum translation and rotation to be applied to the moving structure, such that the root mean square difference of the fit over the specified pairs of equivalent atom is an absolute minimum. This number, given in angstroms, is reported by QUANTA.
For the purpose of this invention, any molecule or molecular complex or binding pocket thereof that has a root mean square deviation of conserved residue backbone atoms (N, Cα, C, O) of less than 1.5 Å when superimposed on the relevant backbone atoms described by structure coordinates listed in FIG. 1 are considered identical. More preferably, the root mean square deviation is less than 1.0 Å.
The term “root mean square deviation” means the square root of the arithmetic mean of the squares of the deviations from the mean. It is a way to express the deviation or variation from a trend or object. For purposes of this invention, the “root mean square deviation” defines the variation in the backbone of a protein from the backbone of IMPDH or a binding pocket portion thereof, as defined by the structure coordinates of IMPDH described herein.
Therefore, according to one embodiment, the present invention provides a molecule or molecular complex comprising all or any parts of the binding pocket defined by structure coordinates of IMPDH amino acids situated within about 7 Å of MPA in the crystal, i.e., amino acids 68, 69, 93, 273, 274, 275, 276, 277, 303, 322, 324, 325, 326, 327, 328, 330, 331, 332, 333, 334, 337, 339, 340, 364, 413, 414, 415, 416, 420, 439, 440, 441, 442, 469, and 470 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
Preferably, the binding pocket is defined by structure coordinates of IMPDH amino acids situated within about 5 Å of MPA in the crystal, i.e., amino acids 274, 275, 276, 277, 303, 322, 324, 325, 326, 331, 333, 414, 415, and 441 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
More preferably, the binding pocket is defined by structure coordinates of IMPDH amino acids situated within about 3.5 Å of MPA in the crystal, i.e., amino acids 275, 276, 303, 325, 326, 331, 333 and 441 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
According to another embodiment, the present invention provides a molecule or molecular complex comprising all or any parts of the binding pocket defined by structure coordinates of IMPDH amino acids situated within about 7 Å of XMP* in the crystal, i.e., amino acids 67, 68, 69, 70, 73, 274, 275, 276, 303, 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 364, 365, 366, 367, 368, 385, 386, 387, 388, 389, 391, 411, 412, 413, 414, 415, 416, 419, 440, 441, 442, 443, 500, 501, 502, 503, 504, 505, and 506 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
Preferably, the binding pocket is defined by structure coordinates of IMPDH amino acids situated within about 5 Å of XMP* in the crystal, i.e., amino acids 68, 69, 70, 303, 322, 326, 327, 328, 329, 330, 331, 332, 333, 335, 364, 365, 366, 367, 385, 386, 387, 388, 411, 413, 414, 415, 416, 419, 441, 442, 443, 501, 502, 503, and 504 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
More preferably, the binding pocket is defined by structure coordinates of IMPDH amino acids situated within about 3.5 Å of XMP* in the crystal, i.e., amino acids 68, 70, 322, 328, 329, 331, 332, 335, 364, 366, 387, 388, 411, 413, 414, 415, 441, 442, 501, and 502 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
According to another embodiment, the present invention provides a molecule or molecular complex comprising all or any parts of the binding pocket defined by structure coordinates of IMPDH amino acids situated within about 7 Å of either MPA or XMP* in the crystal, i.e., amino acids 67, 68, 69, 70, 73, 93, 273, 274, 275, 276, 277, 303, 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 337, 339, 340, 364, 365, 366, 367, 368, 385, 386, 387, 388, 389, 391, 411, 412, 413, 414, 415, 416, 419, 420, 439, 440, 441, 442, 443, 469, 470, 500, 501, 502, 503, 504, 505, and 506, according to FIG. 1, or a homologue of said molecule or molecular complex comprising a binding pocket or pockets that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
Preferably, the binding pocket that is defined by structure coordinates of those IMPDH amino acids situated within about 5 Å of either MPA or XMP* in the crystal, i.e., amino acids 68, 69, 70, 274, 275, 276, 277, 303, 322, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 335, 364, 365, 366, 367, 385, 386, 387, 388, 411, 413, 414, 415, 416, 441, 442, 443, 501, 502, 503, and 504 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of those amino acids of not more than 1.5 Å.
More preferably, the binding pocket is defined by structure coordinates of IMPDH amino acids situated within about 3.5 Å of MPA or XMP* in the crystal, i.e., amino acids 68, 70, 275, 276, 303, 322, 325, 326, 328, 329, 331, 332, 333, 335, 364, 366, 387, 388, 411, 413, 414, 415, 441, 442, 501, and 502 according to FIG. 1 +/− a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
Even more preferred are molecules or molecular complexes that are defined by the entire set of structure coordinates in FIG. 1 +/− a root mean square deviation from the conserved backbone atoms of those amino acids of not more than 1.5 Å. An alternative more preferred embodiment of this invention is a molecular complex that comprises amino acids 1–514 of IMPDH, XMP*, and MPA.
In order to use the structure coordinates generated for the IMPDH/XMP*/MPA complex or one of its binding pockets or homologues thereof, it is sometimes necessary to convert them into a three-dimensional shape. This is achieved through the use of commercially available software that is capable of generating three-dimensional graphical representations of molecules or portions thereof from a set of structure coordinates.
Therefore, according to another embodiment of this invention is provided a machine-readable storage medium comprising a data storage material encoded with machine readable data which, when using a machine programmed with instructions for using said data, is capable of displaying a graphical three-dimensional representation of any of the molecule or molecular complexes of this invention that have been described above.
According to another embodiment, the present invention provides a machine-readable data storage medium, comprising a data storage material encoded with machine readable data which, when using a machine programmed with instructions for using said data, is capable of displaying a graphical three-dimensional representation of a molecule or molecular complex comprising all or any parts of a binding pocket defined by structure coordinates of IMPDH amino acids situated within about 7 Å of MPA in the crystal, i.e., amino acids 68, 69, 93, 273, 274, 275, 276, 277, 303, 322, 324, 325, 326, 327, 328, 330, 331, 332, 333, 334, 337, 339, 340, 364, 413, 414, 415, 416, 420, 439, 440, 441, 442, 469, and 470 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
Preferably, the binding pocket defined by structure coordinates of IMPDH amino acids situated within about 5 Å of MPA in the crystal, i.e., amino acids 274, 275, 276, 277, 303, 322, 324, 325, 326, 331, 333, 414, 415, and 441 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
More preferably, the binding pocket is defined by structure coordinates of IMPDH amino acids situated within about 3.5 Å of MPA in the crystal, i.e., amino acids 275, 276, 303, 325, 326, 331, 333 and 441 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
According to another embodiment, the present invention provides a machine-readable data storage medium, comprising a data storage material encoded with machine readable data which, when using a machine programmed with instructions for using said data, is capable of displaying a graphical three-dimensional representation of a molecule or molecular complex comprising all or any parts of a binding pocket defined by structure coordinates of IMPDH amino acids situated within about 7 Å of XPM* in the crystal, i.e., amino acids 67, 68, 69, 70, 73, 274, 275, 276, 303, 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 364, 365, 366, 367, 368, 385, 386, 387, 388, 389, 391, 411, 412, 413, 414, 415, 416, 419, 440, 441, 442, 443, 500, 501, 502, 503, 504, 505, and 506 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
Preferably, the binding pocket is defined by structure coordinates of IMPDH amino acids situated within about 5 Å of XMP* in the crystal, i.e., amino acids 68, 69, 70, 303, 322, 326, 327, 328, 329, 330, 331, 332, 333, 335, 364, 365, 366, 367, 385, 386, 387, 388, 411, 413, 414, 415, 416, 419, 441, 442, 443, 501, 502, 503, and 504 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
More preferably, the binding pocket is defined by structure coordinates of IMPDH amino acids situated within about 3.5 Å of XMP* in the crystal, i.e., amino acids 68, 70, 322, 328, 329, 331, 332, 335, 364, 366, 387, 388, 411, 413, 414, 415, 441, 442, 501, and 502 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
According to another embodiment, the present invention provides a machine-readable storage medium is capable of displaying a graphical three-dimensional representation of a molecule or molecular complex comprising all or any parts of a binding pocket defined by structure coordinates of IMPDH amino acids situated within about 7 Å of MPA or XMP* in the crystal, i.e., amino acids 67, 68, 69, 70, 73, 93, 273, 274, 275, 276, 277, 303, 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 337, 339, 340, 364, 365, 366, 367, 368, 385, 386, 387, 388, 389, 391, 411, 412, 413, 414, 415, 416, 419, 420, 439, 440, 441, 442, 443, 469, 470, 500, 501, 502, 503, 504, 505, and 506 according to FIG. 1, or a homologue of said molecule or molecular complex, wherein said homologue comprises a binding pocket that has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
Preferably, the binding pocket is defined by structure coordinates of IMPDH amino acids situated within about 5 Å of MPA or XMP* in the crystal, i.e., amino acids 68, 69, 70, 274, 275, 276, 277, 303, 322, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 335, 364, 365, 366, 367, 385, 386, 387, 388, 411, 413, 414, 415, 416, 441, 442, 443, 501, 502, 503, and 504 according to FIG. 1 +/− a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
More preferably, the binding pocket is defined by structure coordinates of IMPDH amino acids situated within about 3.5 Å of MPA or XMP* in the crystal, i.e., amino acids 68, 70, 275, 276, 303, 322, 325, 326, 328, 329, 331, 332, 333, 335, 364, 366, 387, 388, 411, 413, 414, 415, 441, 442, 501, and 502 according to FIG. 1 +/− a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å.
Even more preferred is a machine-readable data storage medium that is capable of displaying a graphical three-dimensional representation of a molecule or molecular complex that is defined by the structure coordinates of all of the amino acids in FIG. 1 +/− a root mean square deviation from the backbone atoms of those amino acids of not more than 1.5 Å.
According to an alternate embodiment, the machine-readable data storage medium comprises a data storage material encoded with a first set of machine readable data which comprises the Fourier transform of the structure coordinates set forth in FIG. 1, and which, when using a machine programmed with instructions for using said data, can be combined with a second set of machine readable data comprising the X-ray diffraction pattern of a molecule or molecular complex to determine at least a portion of the structure coordinates corresponding to the second set of machine readable data.
FIG. 9 demonstrates one version of these embodiments. System 10 includes a computer 11 comprising a central processing unit (“CPU”) 20, a working memory 22 which may be, e.g, RAM (random-access memory) or “core” memory, mass storage memory 24 (such as one or more disk drives or CD-ROM drives), one or more cathode-ray tube (“CRT”) display terminals 26, one or more keyboards 28, one or more input lines 30, and one or more output lines 40, all of which are interconnected by a conventional bidirectional system bus 50.
Input hardware 36, coupled to computer 11 by input lines 30, may be implemented in a variety of ways. Machine-readable data of this invention may be inputted via the use of a modem or modems 32 connected by a telephone line or dedicated data line 34. Alternatively or additionally, the input hardware 36 may comprise CD-ROM drives or disk drives 24. In conjunction with display terminal 26, keyboard 28 may also be used as an input device.
Output hardware 46, coupled to computer 11 by output lines 40, may similarly be implemented by conventional devices. By way of example, output hardware 46 may include CRT display terminal 26 for displaying a graphical representation of a binding pocket of this invention using a program such as QUANTA as described herein. Output hardware might also include a printer 42, so that hard copy output may be produced, or a disk drive 24, to store system output for later use.
In operation, CPU 20 coordinates the use of the various input and output devices 36, 46, coordinates data accesses from mass storage 24 and accesses to and from working memory 22, and determines the sequence of data processing steps. A number of programs may be used to process the machine-readable data of this invention. Such programs are discussed in reference to the computational methods of drug discovery as described herein. Specific references to components of the hardware system 10 are included as appropriate throughout the following description of the data storage medium.
FIG. 10 shows a cross section of a magnetic data storage medium 100 which can be encoded with a machine-readable data that can be carried out by a system such as system 10 of FIG. 9. Medium 100 can be a conventional floppy diskette or hard disk, having a suitable substrate 101, which may be conventional, and a suitable coating 102, which may be conventional, on one or both sides, containing magnetic domains (not visible) whose polarity or orientation can be altered magnetically. Medium 100 may also have an opening (not shown) for receiving the spindle of a disk drive or other data storage device 24.
The magnetic domains of coating 102 of medium 100 are polarized or oriented so as to encode in manner which may be conventional, machine readable data such as that described herein, for execution by a system such as system 10 of FIG. 9.
FIG. 11 shows a cross section of an optically-readable data storage medium 110 which also can be encoded with such a machine-readable data, or set of instructions, which can be carried out by a system such as system 10 of FIG. 9. Medium 110 can be a conventional compact disk read only memory (CD-ROM) or a rewritable medium such as a magneto-optical disk which is optically readable and magneto-optically writable. Medium 100 preferably has a suitable substrate 111, which may be conventional, and a suitable coating 112, which may be conventional, usually of one side of substrate 111.
In the case of CD-ROM, as is well known, coating 112 is reflective and is impressed with a plurality of pits 113 to encode the machine-readable data. The arrangement of pits is read by reflecting laser light off the surface of coating 112. A protective coating 114, which preferably is substantially transparent, is provided on top of coating 112.
In the case of a magneto-optical disk, as is well known, coating 112 has no pits 113, but has a plurality of magnetic domains whose polarity or orientation can be changed magnetically when heated above a certain temperature, as by a laser (not shown). The orientation of the domains can be read by measuring the polarization of laser light reflected from coating 112. The arrangement of the domains encodes the data as described above.
Thus, in accordance with the present invention, data capable of displaying the three dimensional structure of IMPDH and portions thereof and their structurally similar homologues is stored in a machine-readable storage medium, which is capable of displaying a graphical three-dimensional representation of the structure. Such data may be used for a variety of purposes, such as drug discovery.
For example, the structure encoded by the data may be computationally evaluated for its ability to associate with chemical entities. Chemical entities that associate with IMPDH may inhibit IMPDH, and are potential drug candidates. Alternatively, the structure encoded by the data may be displayed in a graphical three-dimensional representation on a computer screen. This allows visual inspection of the structure, as well as visual inspection of the structure's association with chemical entities.
Thus, according to another embodiment, the invention relates to a method for evaluating the potential of a chemical entity to associate with any of the molecules or molecular complexes set forth above. This method comprises the steps of: a) employing computational means to perform a fitting operation between the chemical entity and a binding pocket of the molecule or molecular complex; and b) analyzing the results of said fitting operation to quantify the association between the chemical entity and the binding pocket. The term “chemical entity”, as used herein, refers to chemical compounds, complexes of at least two chemical compounds, and fragments of such compounds or complexes.
For the first time, the present invention permits the use of molecular design techniques to identify, select and design chemical entities, including inhibitory compounds, capable of binding to IMPDH-like binding pockets.
Applicants' elucidation of the XMP* and MPA binding sites on IMPDH provides the necessary information for designing new chemical entities and compounds that may interact with either or both IMPDH-like binding pockets, in whole or in part. This elucidation also enables the evaluation of structure-activity data for analogs of MPA or other compounds which bind to IMPDH-like binding pockets.
Throughout this section, discussions about the ability of an entity to bind to, associate with or inhibit a IMPDH-like binding pocket refers to features of the entity alone. Assays to determine if a compound binds to IMPDH are well known in the art [B. Magasanik et al., J. Biol. Chem., 226, p. 339 (1957)].
The design of compounds that bind to or inhibit IMPDH-like binding pockets according to this invention generally involves consideration of two factors. First, the entity must be capable of physically and structurally associating with parts or all of the IMPDH-like binding pockets. Non-covalent molecular interactions important in this association include hydrogen bonding, van der Waals interactions, hydrophobic interactions and electrostatic interactions.
Second, the entity must be able to assume a conformation that allows it to associate with the IMPDH-like binding pocket directly. Although certain portions of the entity will not directly participate in these associations, those portions of the entity may still influence the overall conformation of the molecule. This, in turn, may have a significant impact on potency. Such conformational requirements include the overall three-dimensional structure and orientation of the chemical entity in relation to all or a portion of the binding pocket, or the spacing between functional groups of an entity comprising several chemical entities that directly interact with the IMPDH-like binding pocket or homologues thereof.
The potential inhibitory or binding effect of a chemical entity on a IMPDH-like binding pocket may be analyzed prior to its actual synthesis and testing by the use of computer modeling techniques. If the theoretical structure of the given entity suggests insufficient interaction and association between it and the IMPDH-like binding pocket, testing of the entity is obviated. However, if computer modeling indicates a strong interaction, the molecule may then be synthesized and tested for its ability to bind to a IMPDH-like binding pocket. This may be achieved by testing the ability of the molecule to inhibit IMPDH using the assays described in Examples 7 and 8. In this manner, synthesis of inoperative compounds may be avoided.
A potential inhibitor of a IMPDH-like binding pocket may be computationally evaluated by means of a series of steps in which chemical entities or fragments are screened and selected for their ability to associate with the IMPDH-like binding pockets.
One skilled in the art may use one of several methods to screen chemical entities or fragments for their ability to associate with a IMPDH-like binding pocket. This process may begin by visual inspection of, for example, a IMPDH-like binding pocket on the computer screen based on the IMPDH structure coordinates in FIG. 1 or other coordinates which define a similar shape generated from the machine-readable storage medium. Selected fragments or chemical entities may then be positioned in a variety of orientations, or docked, within that binding pocket as defined supra. Docking may be accomplished using software such as Quanta and Sybyl, followed by energy minimization and molecular dynamics with standard molecular mechanics force fields, such as CHARMM and AMBER.
Specialized computer programs may also assist in the process of selecting fragments or chemical entities. These include:

    1. GRID (P. J. Goodford, “A Computational Procedure for Determining Energetically Favorable Binding Sites on Biologically Important Macromolecules”, J. Med. Chem., 28, pp. 849–857 (1985)). GRID is available from Oxford University, Oxford, UK.
    2. MCSS (A. Miranker et al., “Functionality Maps of Binding Sites: A Multiple Copy Simultaneous Search Method.” Proteins: Structure, Function and Genetics, 11, pp. 29–34 (1991)). MCSS is available from Molecular Simulations, San Diego, Calif.
    3. AUTODOCK (D. S. Goodsell et al., “Automated Docking of Substrates to Proteins by Simulated Annealing”, Proteins: Structure, Function, and Genetics, 8, pp. 195–202 (1990)). AUTODOCK is available from Scripps Research Institute, La Jolla, Calif.
    4. DOCK (I. D. Kuntz et al., “A Geometric Approach to Macromolecule-Ligand Interactions”, J. Mol. Biol., 161, pp. 269–288 (1982)). DOCK is available from University of California, San Francisco, Calif.


Once suitable chemical entities or fragments have been selected, they can be assembled into a single compound or complex. Assembly may be preceded by visual inspection of the relationship of the fragments to each other on the three-dimensional image displayed on a computer screen in relation to the structure coordinates of IMPDH. This would be followed by manual model building using software such as Quanta or Sybyl [Tripos Associates, St. Louis, Mo.].
Useful programs to aid one of skill in the art in connecting the individual chemical entities or fragments include:

    1. CAVEAT (P. A. Bartlett et al, “CAVEAT: A Program to Facilitate the Structure-Derived Design of Biologically Active Molecules”, in Molecular Recognition in Chemical and Biological Problems”, Special Pub., Royal Chem. Soc., 78, pp. 182–196 (1989); G. Lauri and P. A. Bartlett, “CAVEAT: a Program to Facilitate the Design of Organic Molecules”, J. Comput. Aided Mol. Des., 8, pp. 51–66 (1994)). CAVEAT is available from the University of California, Berkeley, Calif.
    2. 3D Database systems such as ISIS (MDL Information Systems, San Leandro, Calif.). This area is reviewed in Y. C. Martin, “3D Database Searching in Drug Design”, J. Med. Chem., 35, pp. 2145–2154 (1992).
    3. HOOK (M. B. Eisen et al, “HOOK: A Program for Finding Novel Molecular Architectures that Satisfy the Chemical and Steric Requirements of a Macromolecule Binding Site”, Proteins: Struct., Funct., Genet., 19, pp. 199–221 (1994). HOOK is available from Molecular Simulations, San Diego, Calif.


Instead of proceeding to build an inhibitor of a IMPDH-like binding pocket in a step-wise fashion one fragment or chemical entity at a time as described above, inhibitory or other IMPDH binding compounds may be designed as a whole or “de novo” using either an empty binding site or optionally including some portion(s) of a known inhibitor(s). There are many de novo ligand design methods including:

    1. LUDI (H.-J. Bohm, “The Computer Program LUDI: A New Method for the De Novo Design of Enzyme Inhibitors”, J. Comp. Aid. Molec. Design, 6, pp. 61–78 (1992)). LUDI is available from Molecular Simulations Incorporated, San Diego, Calif.
    2. LEGEND (Y. Nishibata et al., Tetrahedron, 47, p. 8985 (1991)). LEGEND is available from Molecular Simulations Incorporated, San Diego, Calif.
    3. LeapFrog (available from Tripos Associates, St. Louis, Mo.).
    4. SPROUT (V. Gillet et al, “SPROUT: A Program for Structure Generation)”, J. Comput. Aided Mol. Design, 7, pp. 127–153 (1993)). SPROUT is available from the University of Leeds, UK.


Other molecular modeling techniques may also be employed in accordance with this invention [see, e.g., N. C. Cohen et al., “Molecular Modeling Software and Methods for Medicinal Chemistry, J. Med. Chem., 33, pp. 883–894 (1990); see also, M. A. Navia and M. A. Murcko, “The Use of Structural Information in Drug Design”, Current Opinions in Structural Biology, 2, pp. 202–210 (1992); L. M. Balbes et al., “A Perspective of Modern Methods in Computer-Aided Drug Design”, in Reviews in Computational Chemistry, Vol. 5, K. B. Lipkowitz and D. B. Boyd, Eds., VCH, New York, pp. 337–380 (1994); see also, W. C. Guida, “Software For Structure-Based Drug Design”, Curr. Opin. Struct. Biology, 4, pp. 777–781 (1994)].
Once a compound has been designed or selected by the above methods, the efficiency with which that entity may bind to an IMPDH binding pocket may be tested and optimized by computational evaluation. For example, an effective IMPDH binding pocket inhibitor must preferably demonstrate a relatively small difference in energy between its bound and free states (i.e., a small deformation energy of binding). Thus, the most efficient IMPDH binding pocket inhibitors should preferably be designed with a deformation energy of binding of not greater than about 10 kcal/mole, more preferably, not greater than 7 kcal/mole. IMPDH binding pocket inhibitors may interact with the binding pocket in more than one conformation that is similar in overall binding energy. In those cases, the deformation energy of binding is taken to be the difference between the energy of the free entity and the average energy of the conformations observed when the inhibitor binds to the protein.
An entity designed or selected as binding to an IMPDH binding pocket may be further computationally optimized so that in its bound state it would preferably lack repulsive electrostatic interaction with the target enzyme and with the surrounding water molecules. Such non-complementary electrostatic interactions include repulsive charge-charge, dipole-dipole and charge-dipole interactions.
Specific computer software is available in the art to evaluate compound deformation energy and electrostatic interactions. Examples of programs designed for such uses include: Gaussian 94, revision C (M. J. Frisch, Gaussian, Inc., Pittsburgh, Pa. ©1995); AMBER, version 4.1 (P. A. Kollman, University of California at San Francisco, ©1995); QUANTA/CHARMM (Molecular Simulations, Inc., San Diego, Calif. ©1995); Insight II/Discover (Molecular Simulations, Inc., San Diego, Calif. ©1995); DelPhi (Molecular Simulations, Inc., San Diego, Calif. ©1995); and AMSOL (Quantum Chemistry Program Exchange, Indiana University). These programs may be implemented, for instance, using a Silicon Graphics workstation such as an Indigo2 with “IMPACT” graphics. Other hardware systems and software packages will be known to those skilled in the art.
Another approach enabled by this invention, is the computational screening of small molecule databases for chemical entities or compounds that can bind in whole, or in part, to a IMPDH binding pocket. In this screening, the quality of fit of such entities to the binding site may be judged either by shape complementarity or by estimated interaction energy [E. C. Meng et al., J. Comp. Chem., 13, pp. 505–524 (1992)].
Thus, enabled by this invention are compounds that inhibit IMPDH by associating directly with the MPA binding site. In one embodiment, such compounds comprise a surrogate for the bicyclic ring portion of MPA and additional functionality which imparts affinity for IMPDH. Preferably, the compounds comprise a substantially hydrophobic core capable of making van der Waals contact with the bound XMP* intermediate and one or more of the following IMPDH residues: Asp 274, Ser 275, Ser 276, Gln 277, Asn 303, Arg 322, Gly 324, Met 325, Gly 326, Thr 333, Met 414, Gly 415 and Gln 441, and further comprise one or more atoms that are substituted via linkers onto the substantially hydrophobic core, wherein the atoms are capable of forming one or more hydrogen bonds with residues such as Asp 274, Thr 333, Gln 441, Gly 326, or nearby residues on IMPDH. Preferably, these linkers form additional van der Waals interactions with residues His 93, Gly 251, Thr 252, His 253, Asp 256, Arg 259, Leu 273, Phe 282, and Gln 283, or nearby residues. More preferably, such compounds have a strain energy of 10 kcal/mol or less. Even more preferably, these compounds contain fewer than three secondary amide bonds. Most preferably, these compounds have a molecular weight of less than 1000.
For example, according to another embodiment, these compounds are those of formula (I):




wherein:
A is —C(O)— or —S(O)2—;
B is —O— or —N(R1)—;
D is —CH2— or —CH2CH2—;
E is selected from the group consisting of C(OH), C(SH), C(NH2), C(NHR1), C(F), N and —N(O−)—;
G is a C1-5 straight or branched alkyl group or a C2-5 straight or branched alkenyl or alkynyl group, wherein any CH2 group is optionally replaced by —O— or —N(H)—;
J is selected from the group consisting of CH2, O, N, S, and a 5–6 membered monocyclic ring system, containing from 0–4 heteroatoms selected from the group consisting of —O—, —N—, and —S—, and optionally containing one or more double bonds, wherein when a heteroatom is attached to at least one —CH2—, the —CH2— is optionally substituted with ═O;
K is a C1-5 straight or branched alkyl group or a C2-5 straight or branched alkenyl or alkynyl group, wherein any CH2 of the alkyl, alkenyl, or alkynyl groups is optionally replaced by —O—, —N(H)—, or —S—;
L is selected from the group consisting of —R1, —OH, —OR1, —NH2, and —N(H) (R1);
M is a C0-3 straight or branched alkyl group or a C2-5 straight or branched alkenyl group;
N is selected from the group consisting of —OH, —OR1, —NH2, —N(H)(R1), —CO2H, —F, —Cl, —S(O2)NH2, —S(O2)N(H)(R1), —NO2, and —CN;
provided that when M is —CH2—, then N is not —OH;
U is a substituent bound to an alkyl carbon atom of K and selected from the group consisting of —OH, —OR1, —NH2, —N(H) (R1), —CO2H, —F, —S(O2) NH2, —SO2N(H) (R1), —NO2, and —CN; and
R1 is a C1-4 straight or branched alkyl.
In another embodiment, this invention provides compounds of formula (II):




wherein,
A, B, E, G, J, K and U are as defined above in formula (I);
each D is independently selected from —CH2— or —CH2CH2—;
Q is selected from the group consisting of —O—, —N(H)—, —N(R1)—, and —CH2—; and
T is selected from the group consisting of —CH2—, —C(O)— and —S(O2)—.
In another embodiment, this invention provides compounds of formula (III):




wherein:
A, B, D, E, G, J, K and U are as defined above in formula (I);
each D is as defined above in formula (II);
Q is selected from the group consisting of —O—, —N(H)—, —N(R1)—, and —CH2—; and
T is selected from the group consisting of —CH2—, —C(O)— and —S(O2)—.
Unless expressly stated to the contrary, the terms “—SO2-” and “—S(O)2—” as used herein refer to a sulfone or sulfone derivative (i.e., both appended groups linked to the S), and not a sulfinate ester.
The term “immunosuppressant” refers to a compound or drug which possesses immune response inhibitory activity. Examples of such agents include cyclosporin A, FK506, rapamycin, leflunomide, deoxyspergualin, prednisone, azathioprine, mycophenolate mofetil, OKT3, ATAG and mizoribine.
IMPDH-mediated disease refers to any disease state in which the IMPDH enzyme plays a regulatory role in the metabolic pathway of that disease. Examples of IMPDH-mediated disease include transplant rejection and autoimmune diseases, such as rheumatoid arthritis, multiple sclerosis, juvenile diabetes, asthma, and inflammatory bowel disease, as well as cancer, viral replication diseases and vascular diseases.
Pharmaceutical compositions of this invention comprise a compound of formulae I–III or a pharmaceutically acceptable salt thereof; an additional agent selected from an immunosuppressant, an anti-cancer agent, an anti-viral agent, or an anti-vascular hyperproliferation compound; and any pharmaceutically acceptable carrier, adjuvant or vehicle. Alternate compositions of this invention comprise a compound of formulae I–III or a pharmaceutically acceptable salt thereof; and a pharmaceutically acceptable carrier, adjuvant or vehicle. Such composition may optionally comprise an additional agent selected from an immunosuppressant, an anti-cancer agent, an anti-viral agent, or an anti-vascular hyperproliferation compound.
The pharmaceutical compositions of this invention may be administered orally, parenterally, by inhalation spray, topically, rectally, nasally, buccally, vaginally or via an implanted reservoir. We prefer oral administration or administration by injection. The pharmaceutical compositions of this invention may contain any conventional non-toxic pharmaceutically-acceptable carriers, adjuvants or vehicles. In some cases, the pH of the formulation may be adjusted with pharmaceutically acceptable acids, bases or buffers to enhance the stability of the formulated compound or its delivery form. The term parenteral as used herein includes subcutaneous, intracutaneous, intravenous, intramuscular, intra-articular, intrasynovial, intrasternal, intrathecal, intralesional and intracranial injection or infusion techniques.
Dosage levels of between about 0.01 and about 100 mg/kg body weight per day, preferably between about 0.5 and about 75 mg/kg body weight per day of the IMPDH inhibitory compounds described herein are useful in a monotherapy for the prevention and treatment of IMPDH mediated disease. Typically, the pharmaceutical compositions of this invention will be administered from about 1 to about 5 times per day or alternatively, as a continuous infusion. Such administration can be used as a chronic or acute therapy. The amount of active ingredient that may be combined with the carrier materials to produce a single dosage form will vary depending upon the host treated and the particular mode of administration. A typical preparation will contain from about 5% to about 95% active compound (w/w). Preferably, such preparations contain from about 20% to about 80% active compound.
When the compositions of this invention comprise a combination of an IMPDH inhibitor of formulae I–III and one or more additional therapeutic or prophylactic agents, both the IMPDH inhibitor and the additional agent should be present at dosage levels of between about 10 to 100%, and more preferably between about 10 to 80% of the dosage normally administered in a monotherapy regimen.
The structure coordinates set forth in FIG. 1 can also be used to aid in obtaining structural information about another crystallized molecule or molecular complex. This may be achieved by any of a number of well-known techniques, including molecular replacement.
Therefore, in another embodiment this invention provides a method of utilizing molecular replacement to obtain structural information about a molecule or molecular complex whose structure is unknown comprising the steps of:
a) crystallizing said molecule or molecular complex of unknown structure;
b) generating an X-ray diffraction pattern from said crystallized molecule or molecular complex; and
c) applying at least a portion of the structure coordinates set forth in FIG. 1 to the X-ray diffraction pattern to generate a three-dimensional electron density map of the molecule or molecular complex whose structure is unknown.
By using molecular replacement, all or part of the structure coordinates of the IMPDH/XMP*/MPA complex as provided by this invention (and set forth in FIG. 1) can be used to determine the structure of a crystallized molecule or molecular complex whose structure is unknown more quickly and efficiently than attempting to determine such information ab initio.
Molecular replacement provides an accurate estimation of the phases for an unknown structure. Phases are a factor in equations used to solve crystal structures that can not be determined directly. Obtaining accurate values for the phases, by methods other than molecular replacement, is a time-consuming process that involves iterative cycles of approximations and refinements and greatly hinders the solution of crystal structures. However, when the crystal structure of a protein containing at least a homologous portion has been solved, the phases from the known structure provide a satisfactory estimate of the phases for the unknown structure.
Thus, this method involves generating a preliminary model of a molecule or molecular complex whose structure coordinates are unknown, by orienting and positioning the relevant portion of the IMPDH/XMP*/MPA complex according to FIG. 1 within the unit cell of the crystal of the unknown molecule or molecular complex so as best to account for the observed X-ray diffraction pattern of the crystal of the molecule or molecular complex whose structure is unknown. Phases can then be calculated from this model and combined with the observed X-ray diffraction pattern amplitudes to generate an electron density map of the structure whose coordinates are unknown. This, in turn, can be subjected to any well-known model building and structure refinement techniques to provide a final, accurate structure of the unknown crystallized molecule or molecular complex [E. Lattman, “Use of the Rotation and Translation Functions”, in Meth. Enzymol., 115, pp. 55–77 (1985); M. G. Rossmann, ed., “The Molecular Replacement Method”, Int. Sci. Rev. Ser., No. 13, Gordon & Breach, New York (1972)].
The structure of any portion of any crystallized molecule or molecular complex that is sufficiently homologous to any portion of the IMPDH/XMP*/MPA complex can be resolved by this method.
In a preferred embodiment, the method of molecular replacement is utilized to obtain structural information about a molecule or molecular complex, wherein the complex comprises at least one IMPDH subunit or homologue.
The structure coordinates of IMPDH as provided by this invention are particularly useful in solving the structure of other crystal forms of IMPDH or IMPDH complexes.
Furthermore, the structure coordinates of IMPDH as provided by this invention are useful in solving the structure of IMPDH mutants, which may optionally be crystallized in co-complex with a chemical entity. The crystal structures of a series of such complexes may then be solved by molecular replacement and compared with that of wild-type IMPDH. Potential sites for modification within the various binding sites of the enzyme may thus be identified. This information provides an additional tool for determining the most efficient binding interactions, for example, increased hydrophobic interactions, between IMPDH and a chemical entity or compound.
The structure coordinates are also particularly useful to solve the structure of crystals of IMPDH or IMPDH homologues co-complexed with a variety of chemical entities. This approach enables the determination of the optimal sites for interaction between chemical entities, including candidate IMPDH inhibitors and IMPDH. For example, high resolution X-ray diffraction data collected from crystals exposed to different types of solvent allows the determination of where each type of solvent molecule resides. Small molecules that bind tightly to those sites can then be designed and synthesized and tested for their IMPDH inhibition activity.
All of the complexes referred to above may be studied using well-known X-ray diffraction techniques and may be refined versus 1.5–3 Å resolution X-ray data to an R value of about 0.20 or less using computer software, such as X-PLOR [Yale University, ©1992, distributed by Molecular Simulations, Inc.; see, e.g., Blundell & Johnson, supra; Meth. Enzymol., vol. 114 & 115, H. W. Wyckoff et al., eds., Academic Press (1985)]. This information may thus be used to optimize known IMPDH inhibitors, and more importantly, to design new IMPDH inhibitors.
In another embodiment of this invention is provided a method for preparing a IMPDH/XMP*/MPA crystal comprising the steps of:
a. forming a complex between IMPDH and IMP;
b. adding NAD and MPA to the complex formed in step a;
c. monitoring the accumulation of the IMPDH/XMP*/MPA complex; and
d. crystallizing the complex formed in step c in the presence of MeP.
In order that this invention be more fully understood, the following examples are set forth. These examples are for the purpose of illustration only and are not to be construed as limiting the scope of the invention in any way.
EXAMPLE 1
Purification of IMPDH and Preparation of an IMPDH/XMP*/MPA Complex for Crystallization
Chinese hamster IMPDH type II was purified by published methods [H. J. Gilbert et. al., Biochem. J., 183, pp. 481–494 (1979); and T. Ikegami et al., Life Sci., 40, pp. 2277–2282 (1987)] with the following modifications: Typically 100 g of E. coli cell paste containing the over-expressed protein was suspended in 5 volumes of Buffer a (50 mM Tris-HCl, 300 mM KCl, 2 mM EDTA, 10 mM BME, 1.5 M urea, pH 8.0 at 4° C.). Protease inhibitors were added (0.2 mM PMSF, and 1 mg/liter each of Pepstatin, Leupeptin, and E-64) and the cells were lysed in a microfluidizer (Microfluidics Corporation) at about 4° C. Cellular debris was removed by centrifugation at 45,000×g for 40 minutes at 4° C. Crystalline ammonium sulfate was added slowly to the supernatant, while stirring at 4° C., to a final concentration of 25 g per 100 ml supernatant. The ammonium sulfate solution was allowed to stabilize for 1 hour at 4° C., whereupon the precipitated IMPDH was harvested by centrifugation, as above. The ammonium sulfate pellets were resuspended in Buffer B (50 mM Tris-HCl, 300 mM KCl, 2 mM EDTA, 10 mM BME, 10% glycerol, pH 8.0 at 4° C.) and loaded (20% at a time while the remainder was stored at −70° C.) onto a column (5×12.5 cm) of IMP-Sepharose, prepared from epoxy-activated Sepharose (Pharmacia) essentially as described by Gilbert et al. (1979). The affinity column was washed with more than 4 column volumes of Buffer B, then the enzyme was eluted from the column with the same buffer containing 10 mM IMP (pH readjusted to 8.0 at 4° C., with KOH).
At this point the protein appeared to be >95% pure, as judged by SDS PAGE, but combined size-exclusion chromatography and light-scattering analyses indicated about 30% of the protein existed as high molecular weight aggregates. The IMPDH-containing fractions were pooled and concentrated by ultra-filtration to about 5 mg/ml (50 ml), then fractionated on a Sephacryl S-300 (Pharmacia) size-exclusion column (5×90 cm, eluting at 3 ml/min in Buffer B, at 4° C.). This step removed the aggregated protein and virtually all of the IMP used to elute the enzyme from the affinity matrix. Aggregated IMPDH did not re-appear upon re-concentration, but some purified samples still contained small quantities of IMP. The residual IMP was removed by exhaustive dialysis, when necessary. SDS PAGE and N-terminal sequence analyses indicated that the IMPDH samples were greater than 99% pure at this point.
Samples for crystallization studies were made as follows. One hundred milligrams of purified apo-IMPDH (2–4 mg/ml) was placed in a stirred ultrafiltration cell with a 30 kDa MW cut-off membrane. A stock of 100 mM MPA in MeP was added to a final concentration of 2 mM MPA. IMP and NAD were then added in a 2-fold molar excess over IMPDH protomers and the solutions were allowed to equilibrate at room temperature for 30 minutes. Samples were concentrated and exchanged into fresh Buffer B containing 2 mM MPA, until no IMP, NAD or NADH was detectable by high-performance capillary electrophoresis (HPCE). The fully inhibited samples were concentrated to 20–60 mg/ml, harvested, and centrifuged at 45,000×g for 20 minutes at 4° C. 100 ml aliquots were stored at −70° C. until use in crystallization experiments.
EXAMPLE 2
Crystallization of the IMPDH/XMP*/MPA Complex
Crystals of inhibited IMPDH were grown by vapor diffusion when protein at 20 mg/ml in Buffer B was mixed with reservoir (10% polyethylene glycol 6000, 1 M LiCl, 100 mM morpholino ethyl sulfonic (MES) acid, 5.4% MeP (v/v), 36 mM BME, pH 5.88) at a 4:2 ratio, and allowed to stand over the reservoir solution at 22° C. Crystals grew within 72 hours to form blocks of approximate dimension 0.15 mm×0.15 mm×0.5 mm. The crystals belong to space group P4, with unit cell dimensions a=b=110.6 Å, c=111.0 Å, and angles α=β=γ=90°. There are two IMDPH/XMP*/MPA complexes per asymmetric unit of the crystal. Analysis of the dissolved crystals showed identical protein composition as the original complex solution.
The addition of MeP was critical in obtaining well-diffracting crystals, and was subsequently identified in electron density maps positioned between two phenylalanine residues in the crystal. It was also critical to use a Chinese hamster form of IMPDH; human type II IMPDH did not yield well-diffracting crystals under conditions which did yield well-diffracting crystals of Chinese hamster type II IMPDH.
Those of skill in the art will appreciate that the aforesaid crystallization conditions can be varied and still produce crystals of IMPDH or IMPDH complexes or homologues suitable for structural analysis. Such variations may be used alone or in combination, and include final protein complex concentrations between 1 mg/ml and 100 mg/ml; any combination of IMPDH/XMP*/MPA complex to precipitant ratios; pH buffer concentrations between 1 mM and 500 mM; any concentration of BME or other sulphur reducing agent between 0 mM and 100 mM; pH ranges between 4.0 and 9.0; polyethylene glycol (PEG) concentrations between 1% and 25% (g/100 ml); PEG weights between 2000 and 20000; LiCl or other salt concentration between 50 and 2000 mM; any concentration or type of detergent; any temperature between −5° C. and 30° C.; and crystallization of IMPDH/XMP*/MPA complexes by batch, liquid bridge, or dialysis method using these conditions or variations thereof.
EXAMPLE 3
Crystal Structure Determination of the IMPDH/XMP*/MPA Complex
Crystals were equilibrated with 10% PEG 6000, 1M LiCl, 100 mM MES, 9% glycerol, 3% MeP, 2 mM MPA, pH 5.88 prior to heavy atom derivitization or native data collection, and then transferred to a 15% glycerol version of this stock just before X-ray data collection at −165° C. Native and derivative data sets were collected on frozen crystals by oscillation photography on a Rigaku R-AXIS IIC phosphor imaging area detector mounted on a Rigaku RU200 rotating anode generator (Molecular Structure Corp., Houston, Tex.), operating at 50 kV and 100 mA. Measured intensities were integrated, scaled, and merged using software supplied by the manufacturer (Molecular Structure Corp., Houston, Tex.). Thirty-five heavy atom reagents were tested, and K2WO4, PCMBS, PbCl, EuCl3, tet-HgCl furan, bis-Hg bithiophene, and bis-Hg benzofuran, were identified as compounds that would bind to IMPDH and be useful to determine initial phase angles for the calculation of electron density maps of the IMPDH/XMP*/MPA complex.
Heavy atom positions were located and confirmed with Patterson or difference Fourier syntheses using PHASES [W. Furey et al., S. Am. Cryst. Assoc. Mtg. Summ. 18, p. 73 (1990)]. Heavy atom parameters were refined with PHASES, and used to compute multiple isomorphous replacement (MIR) phases. The mean figure of merit was 0.65 to 4.0 Å. MIR phases were improved and extended by cycles of solvent flattening (B. C. Wang, Meth. Enzym. 115, pp. 90–112 (1985)), and phase combination using SIGMAA (R. J. Reed, Acta Crystallogr., A42, pp. 140–149 (1986). The initial electron density map revealed the alpha-helices and beta-strands of the barrel cores, but the insertion domain (residues 110–224) in each subunit was poorly seen. The molecular model for IMPDH/XMP*/MPA was built into electron density maps using QUANTA (Quanta version 4.1, Molecular Simulations Inc., Burlington Mass., 1995). Cycles of model building, positional refinement, and simulated annealing using XPLOR-3.1 (A. T. Brunger, X-PLOR (Version 3.1), Yale Univ., New Haven (1993)), and phase combination were done until the switch to phases calculated from the model could be made. The current model is consistent with the derived amino-acid sequence of Chinese hamster IMPDH, and the chemical nature of the heavy atom substitutions. The R-factor is 21.7% against X-ray data between 8 and 2.6 Å resolution, with root-mean-square deviation from ideal bond lengths and angles of 0.007 Å and 1.63° respectively. Definitions for Patterson synthesis, Fourier synthesis, R-factor, figure of merit, multiple isomorphous replacement, and other phrases in this section not previously defined can be found in T. L. Blundell and L. N. Johnson, Protein Crystallography, Academic Press (1976). Additional data collection and structure determination statistics is listed in Table 1.
EXAMPLE 4
IMP and MPA Binding to IMPDH
The crystal structure reveals the interactions between XMP*, MPA and IMPDH that lead to tight-binding of these ligands. This structure reveals that a covalent bond has been formed between the C2 carbon of IMP and the sulfur atom of Cys 331 (FIGS. 5–7), to yield an oxidized IMP thioimidate intermediate (XMP*). Cys 331 has also been shown to form a covalent bond with [8-14C] IMP when NAD is present, [J. A. Huete-Perez et al., Biochemistry, 34, pp. 13889–13894 (1995)] and with the purine ring of 6-Cl-IMP, which inactivates IMPDH [L. C. Antonino et al., Biochemistry, 33, pp. 1760–1765 (1994)]. Together these observations confirm an important role for this residue in catalysis.
Many additional interactions between XMP* and the enzyme are observed (FIGS. 5–7). The XMP* phosphate forms hydrogen bonds to the amide nitrogen atoms of Ser 329, Gly 366, Gly 387 Ser 388, and the side chain hydroxyl groups of Ser 329 and Tyr 411. Two water molecules near Gly 387 are also within hydrogen bonding distance of the phosphate moiety. The hypoxanthine ring makes three hydrogen bonds with IMPDH. The first is between the amide nitrogen of Met 414 and the N7 nitrogen, the second between the C6 carbonyl oxygen and the amide nitrogen of Gly 415, and the third is between the N1 nitrogen and the carbonyl oxygen of Gln 441. The hypoxanthine ring may also be stabilized by an interaction between the N3 nitrogen and a water molecule.
The ribose ring, which adopts a C3′-endo conformation, also contributes significantly to binding. The structure shows that the O2′ and O3′ hydroxyl groups form a hydrogen-bonding network to Ser 68 and Asp 364. The ribose O3′ hydroxyl group accepts a proton from the Ser 68 side chain, and donates a proton to the carboxylate group of Asp 364. Van der Waals contacts between XMP*, Met 70, and Ile 330 are also observed.
Mycophenolic acid is a potent, uncompetitive, inhibitor of Chinese hamster IMPDH, and the structure reveals many interactions between MPA and IMPDH active site residues (FIGS. 5–7). One face of the bicyclic ring system is stacked on the XMP* hypoxanthine ring, while the other makes contact with the main-chain atoms of Ser 276. Together, the hexenoic acid tail, methyl substituent, and methoxy group of MPA make van der Waals contacts with the side chain atoms of Asp 274, Ser 275, Ser 276, Asn 303, Arg 322, and Gln 441. Six hydrogen bonds between the drug and IMPDH are also observed. These include hydrogen bonds between the O2 lactone oxygen and the amide nitrogen of Gly 326, and the Cl carbonyl oxygen and hydroxyl group of Thr 333. The hexenoic acid tail of MPA adopts a bent conformation, unlike the extended conformation seen in NMR studies [G. M. Makara et al., J. Med. Chem., 39, pp. 1236–1242 (1996)] and the crystal structure of free MPA [W. Harrison et al., J. Chem. Soc. Perkin Trans. II, pp. 1542–1544 (1972)] allowing the carboxylate group to form hydrogen bonds with the amide nitrogen and side-chain hydroxyl groups of Ser 276 [F. H. Allen et al., J. Chem. Info. Comp. Sci. 31, p. 187 (1991)]. Additionally, the C7 phenolic oxygen forms hydrogen bonds to the side chain hydroxyl group of Thr 333 and the side chain amide of Gln 441.
Enzymatic Mechanism of IMPDH and Inhibition by Mycophenolic Acid
Aspects of the IMPDH catalyzed reaction can be addressed with the crystal structure. The IMPDH catalyzed oxidation of IMP results in transfer of hydrogen to the nicotinamide ring of NAD, forming NADH and XMP. Since direct transfer of the hydride is energetically unfavorable, two mechanisms involving activation of IMP at the inosine C2 position have been proposed [L. Hedstrom and C. C. Wang, Biochemistry, 29, 849–854 (1990)]. In the first mechanism, water, aided by an active site base to provide OH− attack, is added in an initial step at C2. Hydride transfer to NAD then occurs from the tetrahedral intermediate thus formed, producing the enol tautomer of XMP. In the second mechanism, nucleophilic attack on IMP occurs from an active site cysteine thiol. This is followed by hydride transfer to NAD, to yield a covalently bound thioimidate intermediate which is hydrolyzed to XMP in a subsequent step. The crystal structure, in combination with recent results that also demonstrate the formation of an IMPDH-substrate covalent adduct [J. A. Huete-Perez et al., Biochemistry, 34, pp. 13889–13894 (1995)]; J. O. Link and K. Straub, J. Am. Chem. Soc., 118, pp. 2091–2092 (1996)], strongly supports the second mechanism. The direct observation of the covalently bound thioimidate as the oxidized IMP intermediate confirms that enzyme-catalyzed oxidation of IMP occurs via attack of Cys 331 at the C2 position, and excludes a general base mechanism where water is added to the inosine ring in an early step.
Other aspects of the IMPDH catalyzed reaction can be addressed with the crystal structure. Although the MPA-inhibited complex does not contain NAD or NADH, a combination of structural and chemical evidence allows the nicotinamide ring to be modeled into the active site. The nicotinamide ring must be oriented to allow hydride transfer from the C2 position of IMP to the C4 position of NAD. Further, hydride transfer occurs more readily if the nicotinamide and hypoxanthine rings are nearly parallel [Y.-D. Wu et al., J. Am. Chem. Soc., 117, pp. 4100–4108 (1995)], consistent with the favorable interactions provided by stacking between the nicotinamide ring and the bound substrate, and as observed in the structures of glutathione reductase (1GET) and NADH peroxidase (2NPX). It is also known that hydride transfer occurs on the beta face of NAD [D. Cooney et al., Biochim. Biophys. Acta, pp. 89–93 (1987)]. If no large conformational changes occur between the time NADH leaves and MPA binds, these structural considerations support an earlier prediction [L. Hedstrom and C. C. Wang, Biochemistry, 29, pp. 849–854 (1990)] that during hydride transfer the nicotinamide of NAD occupies a position similar to the 6,5 ring system of MPA. In this orientation, the nicotinamide amide moiety would form hydrogen bonds with Gly-324, Thr 333, Gly-326 and Asn-303.
The structure of the inhibited complex also indicates that the phenolic hydroxyl group of MPA, which forms hydrogen bonds to Thr 333 and Gin 441, may be a replacement for the catalytic water that hydrolyzes the thioimidate intermediate to produce XMP. In the absence of MPA, a water molecule in the vicinity of the MPA hydroxyl would be stabilized by hydrogen bonds with Thr 333 and Gln 441, and would be properly positioned for nucleophilic attack at the C2 carbon of the thioimidate intermediate. Therefore, structural features of the bound orientation of MPA indicate it is both a nicotinamide ring and a catalytic water mimic. This hypothesis is consistent with a report that des-hydroxy-MPA is at least 1,000-fold less potent in a cellular assay than MPA [Y. S. Or et al., ACS Meeting, Chicago, Poster No. 112 (1995)]. It has also been shown that the presence of a hydroxyl group which is able to mimic a catalytic water molecule can lead to as much as a 10 kcal/mol improvement in binding [R. Wolfenden and W. M. Kati, Acc. Chem. Res., 24, pp. 209–215 (1991)].
Similarity of IMPDH to Other Enzymes
Insight into the structure and mechanism of related enzymes is revealed in part by the IMPDH structure. Sequence database searching identified GMP reductase as the closest IMPDH homologue, with 63% similarity and 37% identity over a region of 150 amino acids that includes the active site cysteine as well as the phosphate binding site. The high level of sequence conservation around the active site suggests that GMP reductase has a similar fold and active site geometry to IMPDH. This is supported by the observation that GMP reductase, like IMPDH, binds substrate before co-factor [T. Spector et al., J. Biol. Chem., 254, pp. 2308–2315 (1979)].
The structure of the IMPDH α/β barrel is similar to that of other flavin and nicotinamide dependent oxidoreductases, including glycolate oxidase, the NADPH-dependent aldo-keto reductases, and triethylamine dehydrogenase. The C-terminally sequence encoded phosphate binding site is also conserved, but in IMPDH this site is occupied by the IMP ribose phosphate rather than by the phosphates of the flavin or NAD(P) cofactors as seen in the other enzymes. This suggests that the NAD binding site in IMPDH is novel, and may help explain the specificity seen with Inhibitors that bind in this site, such as MPA and thiazole adenine dinucleotide [H. J. Lee et al., Cancer Res., 45, pp. 5512–5520 (1985); L. Hedstrom and C. C. Wang, Biochemistry, 29, pp. 849–854 (1990)]. In contrast, nucleoside analog inhibitors that are competitive with IMP, such as mizoribine and ribavirin phosphate, are more likely to recognize the consensus nucleotide phosphate binding site found in other enzymes, and are thus likely to be less specific than MPA.
FIGS. 2–7 further depict the structure of the IMPDH/XMP*/MPA complex. Thus, FIG. 2 depicts the fold and conformation of IMPDH in three dimensions as determined by x-ray crystallography. The structure is viewed from the C-terminal end of the β-strands that form the α/β barrel. The α-helices on the outside of the barrel are labeled α1 through α8. The portion of the sub-domain that is ordered is shown. Cys 331 is labeled and sits over one end of the barrel. An arrow marks the location of the flap (residues 400–450) that, together with the active site loop, helps form the active site pocket.
FIG. 3 depicts a topological diagram of the IMPDH fold. Secondary structure was assigned using the Kabsch and Sander algorithm, along with visual inspection [W. Kabsch and C. Sander, Biopolymers, 22, pp. 2577–2637 (1983)]. The β-strands and α-helices that form the α/β barrel core are labeled β1 through β8 and α1 through α8. The strands and helices that are not part of the barrel are labeled βA through βP and αA through αG. Cys 331 is located on the loop between strand β6 and helix αD. The sub-domain starts after helix α2 and ends at strand α3, and includes strands βF through βI and helices βB and αC. Parts of the structure that were not visible in electron density maps are marked “???”. The flap (residues 400–450) between strand β8 and helix α8 includes strands βJ through βN and helix αE. Residue assignments for the strands and helices that form the α/β barrel are as follows: β1, residues 65–68; α1, 76–85; β2, 88–91; α2, 98–109; β3, 245–250; α3, 256–266; β4, 270–274; α4, 281–293; β5, 298–302; α5, 307–316; β6, 320–324; α6, 343–355; β7, 360–364; α7, 370–378; β8, 382–386; and α8, 456–469. Residue assignments for the remaining secondary structural elements are as follows: βA, residues 18–19; βB, 35–38; βC, 40–42; βD, 53–55; βE, 59–61; βF, 114–116; βG, 186–189; βH, 206–212; βI, 220–223; βJ, 400–402; βK, 406–409; βL, 411–413; βM, 438–440; βN, 443–448; βO, 489–492; βP, 509–513; αA, 21–23; αB, 194–200; αC, 225–232; αD, 333–337; αE, 416–420; αF, 476–484; and αG, 495–501.
FIG. 4 depicts a ribbon drawing of the IMPDH tetramer, viewed down the crystallographic four-fold axis [M. Carson, J. Appl. Cryst., 24, pp. 958–961 (1991)]. Most of the tetramer-related contacts are made between adjacent barrels, and the surface area buried at each subunit interface is approximately 4000 Å2. Bound potassium ions are seen adjacent to Cys 331 at each subunit interface. Several other contacts are noteworthy. Residues 41–43 form a beta-strand parallel with residues 279–281 in an adjacent subunit. Residues 502–503 make van der Waals contact with Cys 331 in an adjacent subunit. Residues 507–510 form an anti-parallel beta-strand with residues 444–447 of the active site flap of an adjacent subunit, and these two strands can be seen above the potassium ions.
FIG. 5 depicts a stereo view of the IMPDH active site electron density [M. Carson, J. Appl. Cryst., 24, pp. 958–961 (1991)]. Refined coordinates of the XMP* thioimidate intermediate and MPA are shown in thick bonds superimposed on the SigmaA-weighted [R. J. Read, Acta Cryst., A42, pp. 140–149 (1986)] 2Fobs−Fcalc electron density map contoured at 2.0σ. The hypoxanthine ring makes a covalent bond to the sulfur of Cys 331. Some of the side chains that interact with substrate or inhibitor to form the active site pocket are shown using thin bonds. There are only six differences in amino acid sequence between Chinese hamster and human type II IMPDH (R173C, N215D, L265Q, M290I, E292D, and C327S). These are all at least 15 Å away from the active site, except C327S, which is a conservative mutation and points away from the active site pocket. Thus, these interactions should also apply to the human type II form of IMPDH.
FIG. 6 depicts a schematic representation of the XMP*-IMPDH interactions. There is a covalent bond between the sulfur atom of Cys 331 and the C2 carbon of the hypoxanthine ring.
FIG. 7 depicts a schematic representation of the MPA-IMPDH interactions. All proximal water molecules observed in electron density difference maps are labeled “H2O”. All distances pertain to non-hydrogen atoms.
EXAMPLE 5
Mutational and Kinetic Analysis of the Active Site
A series of mutants was generated to determine the roles of human type II IMPDH active site residues in catalysis and inhibitor binding (FIG. 8).
Knowledge of the three-dimensional structure which we have determined permits rationalization of the observed phenotype. Mutation of either Cys 331 or Asp 364 to Ala effectively abolished IMPDH activity relative to wild type. The crystal structure shows that these two residues are in direct contact with XMP*. Cys 331 forms a covalent bond to XMP*, confirming its role in catalysis, while the side-chain of Asp 364 forms a hydrogen bond to the ribose moiety of XMP*. Changing Ser 329 to Ala reduces enzyme activity to 13% of the wild type. This side chain forms a hydrogen bond to the phosphate of XMP*. Residues which make direct contact with MPA also were modified. Substituting Thr 333 with Ile, and Gln 441 with Ala, increased the Ki app. of MPA 300-fold and 25-fold respectively. The Thr 333 to Ile mutation is of particular interest, since it has been observed in murine blastoma cells selected for 10,000-fold increased resistance to MPA [S. D. Hodges et al., J. Biol. Chem. 264, p. 18137 (1989); T. Lightfoot and F. F. Synder, Biochem. Biophys. Acta, 1217, p. 156 (1994)]. The crystal structure shows a hydrogen bond network between Thr 333, Gin 441 and MPA, with the phenolic oxygen in particularly good hydrogen bonding distance and geometry with the side chain hydroxyl of Thr 333 (FIGS. 5–7). The carboxylate group of the hexenoic acid tail of MPA forms two hydrogen bonds with Ser 276. Mutating this residue to Ala disrupts this interaction and leads to a 7-fold increase in Ki app. In contrast, mutations of other active site residues, such as Ser 275, Ser 327 and Gln 368, have little effect on catalytic activity or drug inhibition. The structure indicates these residues do not contact substrate or inhibitor directly.
Kinetic experiments with human IMPDH have revealed that, although IMP can bind to the enzyme in the absence of potassium, this ion is required for the reaction to proceed [B. Xiang et al., J. Biol. Chem., 271, pp. 1435–1440 (1996)]. The crystal structure of this invention explains this observation. We observed a large peak in difference electron density maps that was surrounded by a water molecule and five main-chain carbonyl oxygens, including that of Cys 331 (data not shown). The average distance of the peak center to each of the surrounding ligands was 3.1 Å, and the peak was hexagonally coordinated. The nature and positioning of the ligands suggested a potassium binding site. Thus, potassium may organize protein conformation around the active site, and could help position Cys 331 for catalysis. In addition, three of the carbonyl oxygen ligands reside on residues near the C-terminus of an adjacent IMPDH subunit, suggesting that potassium may also stabilize the tetramer form of IMPDH (FIG. 4).
FIG. 8 depicts the specific activity of IMPDH mutants for IMP substrate. Mutations were made in the human type II IMPDH cDNA cloned into a pT7 blue vector (Novagen) by a four primer PCR method using Pfu DNA polymerase (Stratagen) [A. Rashchian et al., PCR Methods and Applications, 2, pp. 124–130 (1992)]. PCR products were digested with appropriate restriction enzymes and cloned into unique sites within IMPDH cDNA. Mutants were sequenced in the area containing PCR products and the surrounding restriction sites. The full length IMPDH cDNA carrying the confirmed mutation was then subcloned into a pSPC27 vector in the IMPDH-deficient E. coli strain H712 [H. J. J. Nijkamp and P. G. Haan, Biochim. Biophys. Acta 145, pp. 31–40 (1967)]. Cultures (500 ml) were grown at 37° C. for 14 to 16 hours after IPTG induction, and typically yielded 2 grams of cell paste. Cells were resuspended in 50 mM Tris, pH 8.0, 150 mM KCl, 3 mM EDTA, 2 mM DTT buffer containing 10% urea, and then lysed by addition of lysozyme (1 mg/g cell paste) and sonication. IMPDH wild-type and mutant proteins were precipitated from crude lysates by 25% w/v ammonium sulfate. Up to 150 mgs of at least 65%–70% pure IMPDH was obtained by this single purification step. The partially purified IMPDH was resuspended in 50 mM Tris, pH 8.0, 100 mM KCl, 3 mM EDTA, 2 mM DTT, 10% glycerol buffer and used for kinetic analysis. Second order rate constants (k cat/Km, M−1 sec−1) were calculated from IMP titration data at saturating (400 mM). NAD concentration obtained by monitoring the rate of NADH production at 340 nm at 37° C. FIG. 8 also presents the inhibition of various mutant IMPDH proteins by MPA. Ki app. values for MPA were obtained from the rate vs. inhibitor data of IMPDH mutants at saturating IMP and NAD concentrations. The data were fit to the equation for tight-binding uncompetitive inhibition using the program KineTic 3.0 [D. W. Marquardt, J. Soc. Ind. Appl. Math., 11, pp. 431–441 (1963)].
EXAMPLE 6
Use of IMPDH/XMP*/MPA Coordinates for Inhibitor Design
The coordinates in FIG. 1 are used to design compounds, including inhibitory compounds, that associate with IMPDH or homologues of IMPDH. This process may be aided by using a machine-readable data storage medium encoded with a set of machine-executable instructions, wherein the recorded instructions are capable of displaying a three-dimensional representation of the IMPDH/XMP*/MPA complex or a portion thereof. The graphical representation is used according to the methods described herein to design compounds, including inhibitory compounds, that bind to IMPDH. Such compounds may associate with IMPDH at the active site, in the XMP* binding site, in the MPA binding site, or in parts or all of both sites.
Compounds Which Bind at the MPA Binding Site
The process outlined above is used to design a compound that inhibits IMPDH by associating directly with the MPA binding site, and is expected to bind in the MPA binding site. Such compounds will have a surrogate for the bicyclic ring portion of MPA and additional functionality that imparts affinity for IMPDH. The compound comprises a substantially hydrophobic core capable of making van der Waals contact with the bound XMP* intermediate and one or more of the following residues on IMPDH: Asp 274, Ser 275, Ser 276, Gln 277, Asn 303, Arg 322, Gly 324, Met 325, Gly 326, Thr 333, Met 414, Gly 415 and Gln 441. Substituted onto this hydrophobic core is one or more atoms capable of forming one or more hydrogen bonds with residues such as Asp 274, Thr 333, Gln 441, Gly 326, or nearby residues on IMPDH. These hydrogen bonding atoms are connected by linkers to the substantially hydrophobic core. These linkers may form additional van der Waals interactions with residues His 93, Gly 251, Thr 252, His 253, Asp 256, Arg 259, Leu 273, Phe 282, and Gln 283, or nearby residues. In the lowest-energy bound conformation, this molecule will have a strain energy of 10 kcal/mol or less. Further, this molecule will contain fewer than three secondary amide bonds and will have a molecular weight of less than 1000.
EXAMPLE 7
IMPDH Activity Inhibition Assay
IMP dehydrogenase activity was assayed following an adaptation of the method first reported [B. Magasanik et al., J. Biol. Chem. 226, p. 339 (1957)]. Enzyme activity was measured spectrophotometrically, by monitoring the increase in absorbance at 340 nm due to the formation of NADH (e340 is 6220 M−1 cm−1). The reaction mixture contained 0.1 M Tris pH 8.0, 0.1 M KCl, 3 mM EDTA, 2 mM DTT, 0.1 M IMP and enzyme (IMPDH human type II) at a concentration of 15 to 50 nM. This solution is incubated at 37° C. for 10 minutes. The reaction is started by adding NAD to a final concentration of 0.1M and the initial rate is measured by following the linear increase in absorbance at 340 nm for 10 minutes. For reading in a standard spectrophotometer (path length 1 cm) the final volume in the cuvette is 1.0 ml. The assay has also been adapted to a 96 well microtiter plate format; in this case the concentrations of all the reagents remain the same and the final volume is decreased to 200 μl.
For the analysis of inhibitors, the compound in question is dissolved in DMSO to a final concentration of 20 mM and added to the initial assay mixture for preincubation with the enzyme at a final volume of 2–5% (v/v). The reaction is started by the addition of NAD, and the initial rates measured as above. Ki determinations are made by measuring the initial velocities in the presence of varying amounts of inhibitor and fitting the data using the tight-binding equations of Henderson [P. J. F. Henderson, Biochem. J., 127, p. 321 (1972)].
EXAMPLE 8
Immunosuppression (Mitogenesis) Assays Cell Source and Culture
Fresh peripheral blood lymphocytes (PBLs) from LeukoPak cells or whole blood from random normal blood donors (tested HIV-negative and hepatitis negative) are isolated and separated by density centrifugation over Histopaque 1077 (Sigma Chemical Co., St. Louis, Mo.). The murine CTLL cytotoxic T-cell line and the human Jurkat T-cell line are available from ATCC (CTLL-2 ATCC TIB214, JURKAT CLONE E6-1 ATCC TIB152). The human allogeneic B-cell lines used for activation of the fresh PBLs are EBV-transformed lymphocytes from normal healthy adult donors with two completely different HLA haplotypes. All cell lines are routinely tested for the presence of Mycoplasma contamination using the Gibco Mycotect test kit and found to be Mycoplasma-free. Culture medium consisted of RPMI 1640 (Gibco, Grand Island, N.Y.) containing penicillin (50 U/ml) and streptomycin (50 μg/ml), L-glutamine 2 mM, BME (5×10−5), 10% heat-inactivated FCS and 10 mM HEPES.
Compound Solutions and Titrations
All chemical stocks are dissolved in DMSO. Titrations of compounds are made into the medium the individual assay are carried out in, i.e., complete RPMI or HB 104 for final diluted concentrations, using multiple three-fold dilutions from 1 μM or 10 μM stock solutions.
Mitogenesis Assays (“PMA” and “OKT3”)
The inhibitory effect of test compounds on the proliferation of human PBLs in response to mitogens [W. K. Waithe et al., Handbook of Experimental Immunology, 3d Ed., Blackwell Scientific Publications, Oxford (1978); B. B. Mishell et al., Selected Methods in Cellular Immunology, W. H. Freeman and Co., San Francisco, Calif. (1980)] are assessed by stimulation of 5×104 cells with OKT3 (10−4 dilution final) or PMA (10 ng/ml) plus ionomycin (250 ng/ml) in the presence or absence of different concentrations of test compounds and control drugs (CsA, FK506, rapamycin) in final volume of 200 μl per well in 96 well round bottomed plates. After 48 h incubation (37° C., 5% CO2), cells are pulsed with 1 μCi of 3H-Leucine, harvested 24 h later with a Tom Tek cell harvester, and counted in LKB β-scintillation counter. Results (cpm) are compared with controls with medium alone, and concentrations causing 50% reduction in counts (IC50) are calculated.
While we have described a number of embodiments of this invention, it is apparent that our basic examples may be altered to provide other embodiments which utilize the products and processes of this invention. Therefore, it will be appreciated that the scope of this invention is to be defined by the appended claims rather than by the specific embodiments which have been represented by way of example.

## Claims

1. A method for selecting a chemical entity that associates with all or part of a binding pocket of a molecule or molecular complex, or a homologue of said binding pocket, with a deformation energy not greater than about 10 kcal/mole, wherein the structure coordinates of inosine monophosphate dehydrogenase (IMPDH) amino acids 68, 69, 93, 273, 274, 275, 276, 277, 303, 322, 324, 325, 326, 327, 328, 330, 331, 332, 333, 334, 337, 339, 340, 364, 413, 414, 415, 416, 420, 439, 440, 441, 442, 469, and 470 according to FIG. 1 characterize the binding pocket, and wherein the homologue of said binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) determining structure coordinates of IMPDH amino acids 68, 69, 93, 273, 274, 275, 276, 277, 303, 322, 324, 325, 326, 327, 328, 330, 331, 332, 333, 334, 337, 339, 340, 364, 413, 414, 415, 416, 420, 439, 440, 441, 442, 469, and 470 according to FIG. 1 to characterize the binding pocket;
b) employing computational means which utilize all or part of said structure coordinates and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
c) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
d) outputting said quantified deformation energy to a suitable output hardware; and
e) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

2. A method for selecting a chemical entity that associates with all or part of a binding pocket of a molecule or molecular complex, or a homologue of said binding pocket, with a deformation energy not greater than about 10 kcal/mole, wherein structure coordinates of inosine monophosphate dehydrogenase (IMPDH) amino acids 275, 276, 303, 325, 326, 331, 333 and 441 according to FIG. 1 characterize the binding pocket, and wherein the homologue of said binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) determining structure coordinates of IMPDH amino acids 275, 276, 303, 325, 326, 331, 333 and 441 according to FIG. 1 to characterize the binding pocket;
b) employing computational means which utilize all or part of said structure coordinates and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
c) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
d) outputting said quantified deformation energy to a suitable output hardware; and
e) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

3. A method for selecting a chemical entity that associates with all or part of a binding pocket of a molecule or molecular complex, or a homologue of said binding pocket, with a deformation energy not greater than about 10 kcal/mole, wherein structure coordinates of inosine monophosphate dehydrogenase (IMPDH) amino acids 274, 275, 276, 277, 303, 322, 324, 325, 326, 331, 333, 414, 415, and 441 according to FIG. 1 characterize the binding pocket, and wherein the homologue of said binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) determining structure coordinates of IMPDH amino acids 274, 275, 276, 277, 303, 322, 324, 325, 326, 331, 333, 414, 415, and 441 according to FIG. 1 to characterize the binding pocket;
b) employing computational means which utilize all or part of said structure coordinates and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
c) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
d) outputting said quantified deformation energy to a suitable output hardware; and
e) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

4. The method according to claim 1, 2 or 3, further comprising the steps of:
f) repeating steps a) to e) with a second chemical entity that associates with all or another part of said binding pocket, or homologue thereof;
g) optionally, visually inspecting the relationship of the selected first and second chemical entity to each other in relation to the binding pocket or homologue thereof on a computer screen using the three-dimensional graphical representation of the binding pocket or homologue thereof and said selected first and second chemical entity; and
h) assembling the selected first and second chemical entity into a compound or complex that associates with all or part of said binding pocket or homologue thereof by model building.

5. A method for selecting a chemical entity that associates with all or part of a binding pocket of a molecule or molecular complex, or a homologue of said binding pocket, with a deformation energy not greater than about 10 kcal/mole, wherein the structure coordinates of inosine monophosphate dehydrogenase (IMPDH) amino acids 67, 68, 69, 70, 73, 274, 275, 276, 303, 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 364, 365, 366, 367, 368, 385, 386, 387, 388, 389, 391, 411, 412, 413, 414, 415, 416, 419, 440, 441, 442, 443, 500, 501, 502, 503, 504, 505, and 506 according to FIG. 1 characterize the binding pocket, and wherein the homologue of said binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) determining structure coordinates of IMPDH amino acids 67, 68, 69, 70, 73, 274, 275, 276, 303, 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 364, 365, 366, 367, 368, 385, 386, 387, 388, 389, 391, 411, 412, 413, 414, 415, 416, 419, 440, 441, 442, 443, 500, 501, 502, 503, 504, 505, and 506 according to FIG. 1 to characterize the binding pocket;
b) employing computational means which utilize all or part of said structure coordinates and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
c) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
d) outputting said quantified deformation energy to a suitable output hardware; and
e) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

6. A method for selecting a chemical entity that associates with all or part of a binding pocket of a molecule or molecular complex, or a homologue of said binding pocket, with a deformation energy not greater than about 10 kcal/mole, wherein structure coordinates of inosine monophosphate dehydrogenase (IMPDH) amino acids 68, 70, 322, 328, 329, 331, 332, 335, 364, 366, 387, 388, 411, 413, 414, 415, 441, 442, 501, and 502 according to FIG. 1 characterize the binding pocket, and wherein the homologue of said binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) determining structure coordinates of IMPDH amino acids 68, 70, 322, 328, 329, 331, 332, 335, 364, 366, 387, 388, 411, 413, 414, 415, 441, 442, 501, and 502 according to FIG. 1 to characterize the binding pocket;
b) employing computational means which utilize all or part of said structure coordinates and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
c) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
d) outputting said quantified deformation energy to a suitable output hardware; and
e) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

7. A method for selecting a chemical entity that associates with all or part of a binding pocket of a molecule or molecular complex, or a homologue of said binding pocket, with a deformation energy not greater than about 10 kcal/mole, wherein structure coordinates of inosine monophosphate dehydrogenase (IMPDH) amino acids 68, 69, 70, 303, 322, 326, 327, 328, 329, 330, 331, 332, 333, 335, 364, 365, 366, 367, 385, 386, 387, 388, 411, 413, 414, 415, 416, 419, 441, 442, 443, 501, 502, 503, and 504 according to FIG. 1 characterize the binding pocket, and wherein the homologue of said binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å comprising the steps of:
a) determining structure coordinates of IMPDH amino acids 68, 69, 70, 303, 322, 326, 327, 328, 329, 330, 331, 332, 333, 335, 364, 365, 366, 367, 385, 386, 387, 388, 411, 413, 414, 415, 416, 419, 441, 442, 443, 501, 502, 503, and 504 according to FIG. 1 to characterize the binding pocket;
b) employing computational means which utilize all or part of said structure coordinates and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
c) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
d) outputting said quantified deformation energy to a suitable output hardware; and
e) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

8. The method according to claim 5, 6 or 7, further comprising the steps of:
f) repeating steps a) to e) with a second chemical entity that associates with all or another part of said binding pocket, or homologue thereof;
g) optionally, visually inspecting the relationship of the selected first and second chemical entity to each other in relation to the binding pocket or homologue thereof on a computer screen using the three-dimensional graphical representation of the binding pocket or homologue thereof and said selected first and second chemical entity; and
h) assembling the selected first and second chemical entity into a compound or complex that associates with all or part of said binding pocket or homologue thereof by model building.

9. A method for selecting a chemical entity that associates with all or part of a binding pocket of a molecule or molecular complex, or a homologue of said binding pocket, with a deformation energy not greater than about 10 kcal/mole, wherein the structure coordinates of inosine monophosphate dehydrogenase (IMPDH) amino acids 67, 68, 69, 70, 73, 93, 273, 274, 275, 276, 277, 303, 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 337, 339, 340, 364, 365, 366, 367, 368, 385, 386, 387, 388, 389, 391, 411, 412, 413, 414, 415, 416, 419, 420, 439, 440, 441, 442, 443, 469, 470, 500, 501, 502, 503, 504, 505, and 506 according to FIG. 1 characterize the binding pocket, and wherein the homologue of said binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) determining structure coordinates of IMPDH amino acids 67, 68, 69, 70, 73, 93, 273, 274, 275, 276, 277, 303, 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 337, 339, 340, 364, 365, 366, 367, 368, 385, 386, 387, 388, 389, 391, 411, 412, 413, 414, 415, 416, 419, 420, 439, 440, 441, 442, 443, 469, 470, 500, 501, 502, 503, 504, 505, and 506 according to FIG. 1 to characterize the binding pocket;
b) employing computational means which utilize all or part of said structure coordinates and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
c) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
d) outputting said quantified deformation energy to a suitable output hardware; and
e) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

10. The method of claim 9, wherein the docking of the chemical entity with all or part of the binding pocket utilizes shape complementarity or is followed by molecular dynamics or energy minimization.

11. The method of claim 9, wherein the docking of the chemical entity with all or part of the binding pocket includes visual inspection on a computer screen using a computer program capable of generating a three-dimensional graphical representation of said structure coordinates and structure coordinates of said chemical entity.

12. The method according to any one of claims 1, 5 or 9, further comprising the steps of:
f) contacting the selected chemical entity with the molecule or molecular complex; and
g) selecting the chemical entity that inhibits the catalytic activity of the molecule or molecular complex.

13. The method of any one of claims 1, 5 and 9 wherein the method is for selecting a chemical entity that associates with all or part of a binding pocket of a molecule or molecular complex with a deformation energy not greater than about 7 kcal/mole and step e) comprises selecting the chemical entity if said deformation energy is not greater than 7 kcal/mole.

14. A method for selecting a chemical entity that associates with all or part of a binding pocket of a molecule or molecular complex, or a homologue of said binding pocket, with a deformation energy not greater than about 10 kcal/mole, wherein structure coordinates of inosine monophosphate dehydrogenase (IMPDH) amino acids 68, 70, 275, 276, 303, 322, 325, 326, 328, 329, 331, 332, 333, 335, 364, 366, 387, 388, 411, 413, 414, 415, 441, 442, 501, and 502 according to FIG. 1 characterize the binding pocket, and wherein the homologue of said binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å comprising the steps of:
a) determining structure coordinates of IMPDH amino acids 68, 70, 275, 276, 303, 322, 325, 326, 328, 329, 331, 332, 333, 335, 364, 366, 387, 388, 411, 413, 414, 415, 441, 442, 501, and 502 according to FIG. 1 to characterize the binding pocket;
b) employing computational means which utilize all or part of said structure coordinates and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
c) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
d) outputting said quantified deformation energy to a suitable output hardware; and
e) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

15. A method for selecting a chemical entity that associates with all or part of a binding pocket of a molecule or molecular complex, or a homologue of said binding pocket, with a deformation energy not greater than about 10 kcal/mole, wherein structure coordinates of inosine monophosphate dehydrogenase (IMPDH) amino acids 68, 69, 70, 274, 275, 276, 277, 303, 322, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 335, 364, 365, 366, 367, 385, 386, 387, 388, 411, 413, 414, 415, 416, 441, 442, 443, 501, 502, 503, and 504 according to FIG. 1 characterize the binding pocket, and wherein the homologue of said binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å comprising the steps of:
a) determining structure coordinates of IMPDH amino acids 68, 69, 70, 274, 275, 276, 277, 303, 322, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 335, 364, 365, 366, 367, 385, 386, 387, 388, 411, 413, 414, 415, 416, 441, 442, 443, 501, 502, 503, and 504 according to FIG. 1 to characterize the binding pocket;
b) employing computational means which utilize all or part of said structure coordinates and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
c) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
d) outputting said quantified deformation energy to a suitable output hardware; and
e) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

16. The method according to claim 9, 14 or 15, further comprising the steps of:
f) repeating steps a) to e) with a second chemical entity that associates with all or another part of said binding pocket, or homologue thereof;
g) optionally, visually inspecting the relationship of the selected first and second chemical entity to each other in relation to the binding pocket or homologue thereof on a computer screen using the three-dimensional graphical representation of the binding pocket or homologue thereof and said selected first and second chemical entity; and
h) assembling the selected first and second chemical entity into a compound or complex that associates with all or part of said binding pocket or homologue thereof by model building.

17. A method for selecting a chemical entity that associates with all or part a binding pocket of a molecule or molecular complex of inosine monophosphate dehydrogenase (IMPDH), or a homologue thereof, with a deformation energy not greater than about 10 kcal/mole, wherein said homologue of said IMPDH binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) employing computational means which utilize all or part of structure coordinates of said binding pocket defined by IMPDH amino acids 68, 69, 93, 273, 274, 275, 276, 277, 303, 322, 324, 325, 326, 327, 328, 330, 331, 332, 333, 334, 337, 339, 340, 364, 413, 414, 415, 416, 420, 439, 440, 441, 442, 469, and 470 according to FIG. 1, and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
b) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
c) outputting said quantified deformation energy to a suitable output hardware; and
d) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

18. A method for selecting a chemical entity that associates with all or part a binding pocket of a molecule or molecular complex of inosine monophosphate dehydrogenase (IMPDH), or a homologue thereof, with a deformation energy not greater than about 10 kcal/mole, wherein said homologue of said IMPDH binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) employing computational means which utilize all or part of structure coordinates of said binding pocket defined by IMPDH amino acids 275, 276, 303, 325, 326, 331, 333 and 441 according to FIG. 1, and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
b) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
c) outputting said quantified deformation energy to a suitable output hardware; and
d) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

19. A method for selecting a chemical entity that associates with all or part a binding pocket of a molecule or molecular complex of inosine monophosphate dehydrogenase (IMPDH), or a homologue thereof, with a deformation energy not greater than about 10 kcal/mole, wherein said homologue of said IMPDH binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) employing computational means which utilize all or part of said binding pocket structure coordinates defined by IMPDH amino acids 274, 275, 276, 277, 303, 322, 324, 325, 326, 331, 333, 414, 415, and 441 according to FIG. 1, and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
b) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
c) outputting said quantified deformation energy to a suitable output hardware; and
d) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

20. A method for selecting a chemical entity that associates with all or part a binding pocket of a molecule or molecular complex of inosine monophosphate dehydrogenase (IMPDH), or a homologue thereof, with a deformation energy not greater than about 10 kcal/mole, wherein said homologue of said IMPDH binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) employing computational means which utilize all or part of structure coordinates of said binding pocket defined by IMPDH amino acids 67, 68, 69, 70, 73, 274, 275, 276, 303, 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 364, 365, 366, 367, 368, 385, 386, 387, 388, 389, 391, 411, 412, 413, 414, 415, 416, 419, 440, 441, 442, 443, 500, 501, 502, 503, 504, 505, and 506 according to FIG. 1, and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
b) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
c) outputting said quantified deformation energy to a suitable output hardware; and
d) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

21. A method for selecting a chemical entity that associates with all or part a binding pocket of a molecule or molecular complex of inosine monophosphate dehydrogenase (IMPDH), or a homologue thereof, with a deformation energy not greater than about 10 kcal/mole, wherein said homologue of said IMPDH binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) employing computational means which utilize all or part of structure coordinates of said binding pocket defined by IMPDH amino acids 68, 70, 322, 328, 329, 331, 332, 335, 364, 366, 387, 388, 411, 413, 414, 415, 441, 442, 501, and 502 according to FIG. 1, and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
b) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
c) outputting said quantified deformation energy to a suitable output hardware; and
d) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

22. A method for selecting a chemical entity that associates with all or part a binding pocket of a molecule or molecular complex of inosine monophosphate dehydrogenase (IMPDH), or a homologue thereof, with a deformation energy not greater than about 10 kcal/mole, wherein said homologue of said IMPDH binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) employing computational means which utilize all or part of structure coordinates of said binding pocket defined by IMPDH amino acids 68, 69, 70, 303, 322, 326, 327, 328, 329, 330, 331, 332, 333, 335, 364, 365, 366, 367, 385, 386, 387, 388, 411, 413, 414, 415, 416, 419, 441, 442, 443, 501, 502, 503, and 504 according to FIG. 1, and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
b) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
c) outputting said quantified deformation energy to a suitable output hardware; and
d) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

23. A method for selecting a chemical entity that associates with all or part of a binding pocket of a molecule or molecular complex of inosine monophosphate dehydrogenase (IMPDH), or a homologue thereof, with a deformation energy not greater than about 10 kcal/mole, wherein said homologue of said IMPDH binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) employing computational means which utilize all or part of structure coordinates of said binding pocket defined by IMPDH amino acids 67, 68, 69, 70, 73, 93, 273, 274, 275, 276, 277, 303, 322, 323, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 334, 335, 337, 339, 340, 364, 365, 366, 367, 368, 385, 386, 387, 388, 389, 391, 411, 412, 413, 414, 415, 416, 419, 420, 439, 440, 441, 442, 443, 469, 470, 500, 501, 502, 503, 504, 505, and 506 according to FIG. 1, and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
b) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
c) outputting said quantified deformation energy to a suitable output hardware; and
d) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

24. The method of any one of claims 17, 20 and 23, wherein the method is for selecting a chemical entity that associates with all or part of a binding pocket of a molecule or molecular complex of IMPDH, or homologue thereof, with a deformation energy not greater than about 7 kcal/mole and step d) comprises selecting the chemical entity if said deformation energy is not greater than 7 kcal/mole.

25. The method of claim 23, wherein the docking of the chemical entity with all or part of the binding pocket utilizes shape complementarity or is followed by molecular dynamics or energy minimization.

26. The method of claim 25, wherein the method is for selecting a chemical entity wherein the chemical entity is an inosine monophosphate dehydrogenase (IMPDH) inhibitor.

27. The method according to any one of claims 17, 20 or 23, further comprising the steps of:
e) contacting the selected chemical entity with the molecule or molecular complex; and
f) selecting the chemical entity that inhibits the catalytic activity of the molecule or molecular complex.

28. The method of claim 23, wherein the docking of the chemical entity with all or part of the binding pocket includes visual inspection on a computer screen using a computer program capable of generating a three-dimensional graphical representation of said structure coordinates and structure coordinates of said chemical entity.

29. A method for selecting a chemical entity that associates with all or part a binding pocket of a molecule or molecular complex of inosine monophosphate dehydrogenase (IMPDH), or a homologue thereof, with a deformation energy not greater than about 10 kcal/mole, wherein said homologue of said IMPDH binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) employing computational means which utilize all or part of structure coordinates of said binding pocket defined by IMPDH amino acids 68, 70, 275, 276, 303, 322, 325, 326, 328, 329, 331, 332, 333, 335, 364, 366, 387, 388, 411, 413, 414, 415, 441, 442, 501, and 502 according to FIG. 1, and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
b) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
c) outputting said quantified deformation energy to a suitable output hardware; and
d) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.

30. A method for selecting a chemical entity that associates with all or part a binding pocket of a molecule or molecular complex of inosine monophosphate dehydrogenase (IMPDH), or a homologue thereof, with a deformation energy not greater than about 10 kcal/mole, wherein said homologue of said IMPDH binding pocket has a root mean square deviation from the backbone atoms of said amino acids of not more than 1.5 Å, comprising the steps of:
a) employing computational means which utilize all or part of structure coordinates of said binding pocket defined by IMPDH amino acids 68, 69, 70, 274, 275, 276, 277, 303, 322, 324, 325, 326, 327, 328, 329, 330, 331, 332, 333, 335, 364, 365, 366, 367, 385, 386, 387, 388, 411, 413, 414, 415, 416, 441, 442, 443, 501, 502, 503, and 504 according to FIG. 1, and structure coordinates of the chemical entity, to dock the chemical entity with all or part of said binding pocket or homologue thereof;
b) quantifying the deformation energy between the chemical entity and all or part of the binding pocket or homologue thereof;
c) outputting said quantified deformation energy to a suitable output hardware; and
d) selecting the chemical entity if said deformation energy is not greater than about 10 kcal/mole.