# Patent

## Title
Crystal structure of FMS-like tyrosine kinase

## Abstract

The invention relates to FMS-like tyrosine kinase (FLT3), FLT3 binding pockets or FLT3-like binding pockets. The invention relates to a computer comprising a data storage medium encoded with the structure coordinates of such binding pockets. The invention also relates to methods of using the structure coordinates to solve the structure of homologous proteins or protein complexes. The invention relates to methods of using the structure coordinates to screen for and design compounds that bind to FMS-like tyrosine kinase protein, complexes of FMS-like tyrosine kinase protein, homologues thereof, or FLT-3-like protein or protein complexes. The invention also relates to crystallizable compositions and crystals comprising an FMS-like tyrosine kinase cytoplasmic domain or homologues thereof. The invention also relates to methods of identifying inhibitors of the cytoplasmic domain of FMS-like tyrosine kinase protein.

## Background
CROSS-REFERENCE TO RELATED APPLICATIONS
The present application claims the benefit under 35 U.S.C. § 119 of U.S. Provisional patent application No. 60/503,270, filed Sep. 15, 2003, and U.S. Provisional patent application No. 60/540,391, filed Jan. 29, 2004, the entire contents of both applications being incorporated herein by reference.


TECHNICAL FIELD OF INVENTION
The present invention relates to human FMS-like tyrosine kinase (FLT3), FLT3 binding pockets or FLT3-like binding pockets. The present invention provides a computer comprising a data storage medium encoded with the structure coordinates of such binding pockets. 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 for and design compounds, including inhibitory compounds, that bind to FLT3 protein, FLT3 protein complexes, homologues thereof, or FLT3-like protein or FLT3-like protein complexes. The invention also relates to crystallizable compositions and crystals comprising FLT3 cytoplasmic domain.
BACKGROUND OF THE INVENTION
The FMS-like tyrosine kinase (FLT3) is a type III receptor tyrosine kinase that is thought to play a key role in hematopoiesis. Certain classes of FLT3 mutations cause constitutively activated forms of the receptor that are found in significant numbers of patients with acute myeloid leukemia (AML). The mutations occur in either the activation loop, for example, as point mutations of Asp835, or as internal tandem duplication (ITD) sequences in the juxtamembrane (JM) domain.
FLT3 (FMS-like tyrosine kinase 3) (Gilliland, D. G., and Griffin, J. D., Blood 100: 1532-1542 (2002); D. Kottaridis, P., et. al., Br. J. Haematol. 122, 523-538 (2003); Stirewalt, D. L. and Radich, J. P. Nat. Rev. Cancer 3: 650-665 (2003)), also known as FLK-2 (fetal liver kinase 2) and STK-1 (human stem cell kinase 1), belongs to a family of type III receptor tyrosine kinases (RTKs) (Rosnet, O., et. al., Oncogene 6: 1641-1650 (1991); Rosnet, O., et. al., Genomics 9: 380-385 (1991); Small, D., et. al., Proc. Natl. Acad. Sci. USA 91, 459-463 (1994); Matthews, W., et. al., Cell 65: 1143-1152 (1991)). Members of a subset of this family include FLT3, platelet-derived growth factor receptors α and β (PDGFRα and PDGFRβ) (Yarden, Y., et. al., Nature 323: 226-232 (1986); Claesson-Welsh, L., et. al., Methods Enzymol. 198: 72-77 (1991); Claesson-Welsh, L et. al., Proc. Natl. Acad. Sci. USA 86: 4917-4921 (1989); Matsui, T., et. al., Science 243: 800-804 (1989)), FMS (Stanley, E. R., et. al., J. Cell. Biochem. 21:151-159 (1983)) and cKIT (Yarden, Y et. al., EMBO J. 6: 3341-3351 (1987); Mol, C. D., et. al., J. Biol. Chem. 278: 31461-31464 (2003)) and are characterized by an extracellular domain consisting of five immunoglobulin-like (Ig-like) domains, a single transmembrane region, a cytoplasmic juxtamembrane domain (JM) and a cytoplasmic tyrosine kinase domain interrupted by a kinase insert domain (KID) (Agnes, F., et. al., Gene 145: 283-288 (1994); Rosnet, O., and Birnbaum, D. (1993) Crit. Rev. Oncog. 4, 595-61; Scheijen, B., and Griffin, J. D. Oncogene 21: 3314-3333 (2002)). Two groups independently reported the cloning of the flt3 gene (Rosnet, O., et. al., Oncogene, supra; Rosnet, O., et. al., Genomics, supra; Matthews, W et. al., supra).
Subsequently, FL, the ligand for FLT3, and a type I transmembrane protein was cloned from mouse (Lyman, S. D., et. al., Stem Cells 12 Suppl 1: 99-107; discussion 108-110 (1994); Lyman, S. D et. al., Oncogene 11: 1165-1172 (1995); Hannum, C., et. al., Nature 368: 643-648 (1994); Savvides, S. N., et. al., Nat. Struct. Biol. 7: 486-491 (2000)). The binding of FL leads to dimerization, activation and autophosphorylation of the receptor and subsequent activation of several signaling pathways including STAT5 (Zhang, S., et. al., J. Exp. Med. 192: 719-728 (2000)), Ras/mitogen-activated protein kinase (MAPK) and phosphatidylinositol 3 kinase (PI3K)/AKT pathways. The human flt3 gene encodes a 993 amino acid protein of which residues 572-603 and 604-958 represent the JM and tyrosine kinase domains, respectively (Rosnet, O., et. al., Blood 82: 1110-1119 (1993)).
FLT3 is primarily expressed in immature hematopoietic cells (Rosnet, O., et. al., Genomics, supra; deLapeyriere, O., et. al., Differentiation 58: 351-359 (1995)) and is essential for the normal function of stem cells and the immune system (deLapeyriere, O., et. al., supra; Brasel, K., et. al., Leukemia 9: 1212-1218 (1995); Turner, A. M., et. al., Blood 88: 3383-3390 (1996)). FLT3 is also found in placenta, gonads and brain (Maroc, N., et. al., Oncogene 8: 909-918 (1993)) and is expressed in high levels in a wide range of hematopoietic malignancies including 70-100% of acute myelogenous leukemia (AML), acute lymphoblastic leukemia (ALL) and chronic myelogenous leukemia (Rosnet, O., et. al., Acta Haematol. 95: 218-223 (1996); Drexler, H. G. Leukemia 10: 588-599 (1996)).
Two distinct types of FLT3 mutations have been identified in up to 41% of AML patients. Internal tandem duplication (ITD) mutations within the JM domain contribute to about 17-34% of FLT3 activating mutations in AML (Nakao, M., et. al., Leukemia 10: 1911-1918 (1996); Thiede, C., et. al., Blood 99: 4326-4335 (2002)). FLT3-ITD has also been detected at low frequency in myelodysplastic syndrome (MDS) (Yokota, S., et. al., Leukemia 11: 1605-1609 (1997)); Horiike, S., et. al., Leukemia 11: 1442-1446 (1997)). The ITDs are always in-frame, and are limited to the JM domain. However, they vary in length and position from patient to patient. These repeat sequences may serve to disrupt the autoinhibitory activity of the JM domain resulting in the constitutive activation of FLT3. Point mutations at aspartate 835 within the activation loop of the FLT3 kinase domain represent a second class of activating mutations (Yamamoto, Y., et. al., Blood 97: 2434-2439 (2001); Abu-Duhier, F. M., et. al., Br. J. Haematol. 113: 983-988 (2001)). FLT3-Asp835 mutations also lead to constitutive activation of the receptor and have been reported in 7% of AML, 3% of MDS and 3% of all cases. The most common substitution is Asp835Tyr, but other substitutions including Asp835Val, Asp835His, Asp835Glu and Asp835Asn have also been reported (Yamamoto, Y., et. al., supra). Both FLT3-ITD and FLT3-Asp835 mutations are associated with FLT3 autophosphorylation and phosphorylation of downstream targets (Yamamoto, Y., et. al., supra; Mizuki, M., et. al., Blood 96: 3907-3914 (2000); Mizuki, M., et. al., Blood 101: 3164-3173 (2003); Hayakawa, F., et. al., Oncogene 19: 624-631 (2000)).
A novel class of constitutively activated FLT3 mutants has been recently identified in AML patients in which isoleucine 836 is either deleted (FLT3-Ile836del) or substituted with methionine and arginine (FLT3-Ile836Met+Arg) (Thiede, C., et. al., supra). In mice, injection of FLT3-ITD transformed cells results in leukemia-like syndrome (Mizuki, M., et. al., (2000), supra). Several FLT3 inhibitors, such as PKC412 (N-benzoyl staurosporine) (Fabbro, D., et. al., Anticancer Drug Des. 15: 17-28 (2000); Weisberg, E., et. al., Cancer Cell 1: 433-443 (2002)), CT53518 (also known as MLN518) (Kelly, L. M., et. al., Cancer Cell 1: 421-432 (2002)), SU11248 (O'Farrell, A. M., et. al., Blood 101: 3597-3605 (2003)), SU5614 (Spiekermann, K., et. al., Blood 101: 1494-1504 (2003)), and SU5416 (Giles, F. J., et. al., Blood 102: 795-801 (2003)), have been shown to have antitumor activity. Collectively, these data suggest that FLT3 is an attractive therapeutic target for the development of kinase inhibitors for AML and other associated diseases.
SUMMARY OF THE INVENTION
The present invention provides the first time the crystal structure of the autoinhibited, unphosphorylated form of FLT3 comprised of the kinase domain (minus the KID) and the JM domain. This structure shows, for the first time, the autoinhibitory conformation of a complete JM domain in type III class receptor tyrosine kinases. The structure of FLT3 provides direct insight into the mechanism by which the JM domain exerts its autoinhibitory effect on the catalytic activity of the kinase domain, which is likely utilized by other members of the type III receptor tyrosine kinases. The role of the internal tandem duplication (ITD) in the constitutive activation of FLT3 in acute myeloid leukemia (AML) patients is also revealed. This structure provides a framework to explain the aberrant behavior of FLT3 in disease and to understand the possible mechanisms by which the ITD can switch FLT3 from an inactive to a catalytically active form. An analysis of the sequence alignment of other members of the PDGFR family suggests that the FLT3 mechanism of autoinhibition may be a general one applicable to all members of this family. The structure also presents a rationale for the structure-based design of small molecule FLT3 inhibitors as therapeutic agents, thus addressing the need for novel drugs for the treatment of AML and related diseases.
The present invention also provides molecules comprising FLT3 binding pockets, or FLT3-like binding pockets that have similar three-dimensional shapes. In one embodiment, the molecules are FLT3 or FLT3-like proteins, protein complexes, or homologues thereof. In another embodiment, the molecules are FLT3 cytoplasmic domains or homologues thereof. In another embodiment, the molecules are in crystalline form.
The invention provides crystallizable compositions and crystal compositions comprising the cytoplasmic domain of human FLT3 or a homologue thereof with or without a chemical entity.
The invention provides a computer comprising a machine-readable storage medium, comprising a data storage material encoded with machine-readable data, wherein the data defines the binding pockets or domains according to the structure coordinates of molecules or molecular complexes of FLT3 or FLT3-like proteins, protein complexes or homologues thereof. The invention also provides a computer comprising the data storage medium. Such storage medium when read and utilized by a computer programmed with appropriate software can display, on a computer screen or similar viewing device, a three-dimensional graphical representation of such binding pockets or domains. In one embodiment, the structure coordinates of said molecules or molecular complexes are produced by homology modeling of the coordinates of FIG. 1A.
The invention also provides methods for designing, selecting, evaluating and identifying and/or optimizing compounds which bind to the molecules or molecular complexes or their binding pockets. Such compounds are potential inhibitors of FLT3, FLT3-like proteins or their 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 FLT3, particularly FLT3 homologues. This is achieved by using at least some of the structure coordinates obtained from a FLT3 cytoplasmic domain.
The invention provides a crystal comprising a cytoplasmic domain of an FMS-like tyrosine kinase protein or a homologue thereof, wherein the cytoplasmic domain of the FMS-like tyrosine kinase protein is selected from the group consisting of amino acid residues 564-993 of SEQ ID NO:1, amino acid residues 564-958 of SEQ ID NO:1, amino acid residues 564-710 and 762-958 of SEQ ID NO:1, amino acid residues 570-710 and 783-947 of SEQ ID NO:1, amino acid residues 570-958 of SEQ ID NO:1, and amino acid residues 570-710 and 762-958 of SEQ ID NO:1.
The invention also provides a crystal according to paragraph 14, wherein the cytoplasmic domain of the FMS-like tyrosine kinase protein comprises amino acid residues 564-710 and 762-958 of SEQ ID NO:1.
The invention provides a crystallizable composition comprising a cytoplasmic domain of an FMS-like tyrosine kinase protein or a homologue thereof, wherein the cytoplasmic domain of the FMS-like tyrosine kinase protein is selected from the group consisting of amino acid residues 564-993 of SEQ ID NO:1, amino acid residues 564-958 of SEQ ID NO:1, amino acid residues 564-710 and 762-958 of SEQ ID NO:1, amino acid residues 570-710 and 783-947 of SEQ ID NO:1, amino acid residues 570-958 of SEQ ID NO:1, and amino acid residues 570-710 and 762-958 of SEQ ID NO:1.
The invention also provides a crystallizable composition according to paragraph 16, wherein the cytoplasmic domain of the FMS-like tyrosine kinase protein comprises amino acid residues 564-710 and 762-958 of SEQ ID NO:1.
The invention provides a computer comprising:

    
    
        (a) a machine-readable data storage medium, comprising a data storage material encoded with machine-readable data, wherein the data defines a binding pocket or domain selected from the group consisting of:
        
            (i) a set of amino acid residues which are identical to human FMS-like tyrosine kinase (FLTs) amino acid residues H809, R810, and D811 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the set of amino acid residues and the FLT3 amino acid residues is not greater than about 2.0 Å;
            (ii) a set of amino acid residues comprising at least five amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the at least five amino acid residues and the FLT3 amino acid residues which are identical is not greater than about 2.0 Å;
            (iii) a set of amino acid residues comprising at least six amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, K644, A657, L658, E661, M664, L802, K805, S806, C807, V808, H809, R810, D811, C828, D829, F830, G831, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the at least six amino acid residues and the FLT3 amino acid residues which are identical is not greater than about 2.0 Å; and
            (iv) a set of amino acid residues that are identical to FLT3 amino acid residues according to FIG. 1A, wherein the root mean square deviation between the set of amino acid residues and the FLT3 amino acid residues is not more than about 3.0 Å;
        
        
        (b) a working memory for storing instructions for processing the machine-readable data;
        (c) a central processing unit coupled to the working memory and to the machine-readable data storage medium for processing the machine-readable data and a means for generating three-dimensional structural information of the binding pocket or domain; and
        (d) output hardware coupled to the central processing unit for outputting said three-dimensional structural information of the binding pocket or domain, or information produced using the three-dimensional structural information of the binding pocket or domain.
    
    


The invention also provides the computer according to paragraph 18, wherein the binding pocket is produced by homology modeling of the structure coordinates of the FMS-like tyrosine kinase amino acid residues according to FIG. 1A. The invention further provides the computer according to paragraph 18, wherein means for generating three-dimensional structural information is provided by means for generating a three-dimensional graphical representation of the binding pocket or domain.
The invention provides the computer according to paragraph 18, wherein the output hardware is a display terminal, a printer, CD or DVD recorder, ZIPTM or JAZTM drive, a disk drive, or other machine-readable data storage device.
The invention provides a method of using a computer for selecting an orientation of a chemical entity that interacts favorably with a binding pocket or domain selected from the group consisting of:

    
    
        (i) a set of amino acid residues which are identical to human FMS-like tyrosine kinase (FLT3) amino acid residues H809, R810, and D811 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the set of amino acid residues and the FLT3 amino acid residues is not greater than about 2.0 Å;
        (ii) a set of amino acid residues comprising at least five amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the at least five amino acid residues and the FLT3 amino acid residues which are identical is not greater than about 2.0 Å;
        (iii) a set of amino acid residues comprising at least six amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, K644, A657, L658, E661, M664, L802, K805, S806, C807, V808, H809, R810, D811, C828, D829, F830, G831, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the at least six amino acid residues and the FLT3 amino acid residues which are identical is not greater than about 2.0 Å; and
        (iv) a set of amino acid residues that are identical to FLT3 amino acid residues according to FIG. 1A, wherein the root mean square deviation between the set of amino acid residues and the FLT3 amino acid residues is not more than about 3.0 Å;
        the method comprising the steps of:
        (a) providing the structure coordinates of the binding pocket or domain on a computer comprising means for generating three-dimensional structural information from the structure coordinates;
        (b) employing computational means to dock a first chemical entity in the binding pocket or domain;
        (c) quantifying the association between the chemical entity and all or part of the binding pocket or domain for different orientations of the chemical entity; and
        (d) selecting the orientation of the chemical entity with the most favorable interaction based on the quantified association.
    
    


The invention also provides the method according to paragraph 21, further comprising the step of (e) generating a three-dimensional graphical representation of the binding pocket or domain prior to step (b). The invention further provides the method according to paragraph 21, wherein energy minimization, molecular dynamics simulations, or rigid-body minimizations are performed simultaneously with or following step (b). The invention provides the method according to paragraph 21, further comprising the steps of:

    
    
        (e) repeating steps (b) through (d) with a second chemical entity; and
        (f) selecting at least one of the first or second chemical entity that interacts more favorably with the binding pocket or domain based on the quantified association of the first or second chemical entity.
    
    


The invention provides a method of using a computer for selecting an orientation of a chemical entity with a favorable shape complementarity in a binding pocket selected from the group consisting of:

    
    
        (i) a set of amino acid residues which are identical to human FMS-like tyrosine kinase (FLT3) amino acid residues H809, R810, and D811 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the set of amino acid residues and the FLT3 amino acid residues is not greater than about 2.0 Å;
        (ii) a set of amino acid residues comprising at least five amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the at least five amino acid residues and the FLT3 amino acid residues which are identical is not greater than about 2.0 Å;
        (iii) a set of amino acid residues comprising at least six amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, K644, A657, L658, E661, M664, L802, K805, S806, C807, V808, H809, R810, D811, C828, D829, F830, G831, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the at least six amino acid residues and the FLT3 amino acid residues which are identical is not greater than about 2.0 Å; and
        (iv) a set of amino acid residues that are identical to FLT3 amino acid residues according to FIG. 1A, wherein the root mean square deviation between the set of amino acid residues and the FLT3 amino acid residues is not more than about 3.0 Å;
        the method comprising the steps of:
        (a) providing the structure coordinates of the binding pocket and all or part of the JM-B binding motif bound therein on a computer comprising means for generating three-dimensional structural information from the structure coordinates;
        (b) employing computational means to dock a first chemical entity in the binding pocket;
        (c) quantitating the contact score of the chemical entity in different orientations; and
        (d) selecting the orientation with the highest contact score.
    
    


The invention also provides method according to paragraph 23, further comprising the step of:

    
    
        (e) generating a three-dimensional graphical representation of the binding pocket and all or part of the JM-B binding motif bound therein prior to step (b). The invention also provides the method according to paragraph 23, further comprising the steps of:
        (e) repeating steps (b) through (d) with a second chemical entity; and
        (f) selecting at least one of the first or second chemical entity that has a higher contact score based on the quantitated contact score of the first or second chemical entity.
    
    


The invention provides a method for identifying a candidate inhibitor of a molecule or molecular complex comprising a binding pocket or domain selected from the group consisting of:

    
    
        (i) a set of amino acid residues which are identical to human FMS-like tyrosine kinase (FLT3) amino acid residues H809, R810, and D811 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the set of amino acid residues and the FLT3 amino acid residues is not greater than about 2.0 Å;
        (ii) a set of amino acid residues comprising at least five amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the at least five amino acid residues and the FLT3 amino acid residues which are identical is not greater than about 2.0 Å;
        (iii) a set of amino acid residues comprising at least six amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, K644, A657, L658, E661, M664, L802, K805, S806, C807, V808, H809, R810, D811, C828, D829, F830, G831, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the at least six amino acid residues and the FLT3 amino acid residues which are identical is not greater than about 2.0 Å; and
        (iv) a set of amino acid residues that are identical to FLT3 amino acid residues according to FIG. 1A, wherein the root mean square deviation between the set of amino acid residues and the FLT3 amino acid residues is not more than about 3.0 Å;
    
    


comprising the steps of:

    
    
        (a) using a three-dimensional structure of the binding pocket or domain to design, select or optimize a plurality of chemical entities;
        (b) contacting each chemical entity with the molecule or the molecular complex;
        (c) monitoring an inhibitory effect of the catalytic activity of the molecule or molecular complex by each chemical entity; and
        (d) selecting a chemical entity based on the inhibitory effect of the chemical entity on the catalytic activity of the molecule or molecular complex.
    
    


The invention provides a method of designing a compound or complex that interacts with a binding pocket or domain selected from the group consisting of:

    
    
        (i) a set of amino acid residues which are identical to human FMS-like tyrosine kinase (FLT3) amino acid residues H809, R810, and D811 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the set of amino acid residues and the FLT3 amino acid residues is not greater than about 2.0 Å;
        (ii) a set of amino acid residues comprising at least five amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the at least five amino acid residues and the FLT3 amino acid residues which are identical is not greater than about 2.0 Å;
        (iii) a set of amino acid residues comprising at least six amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, K644, A657, L658, E661, M664, L802, K805, S806, C807, V808, H809, R810, D811, C828, D829, F830, G831, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between the at least six amino acid residues and the FLT3 amino acid residues which are identical is not greater than about 2.0 Å; and
        (iv) a set of amino acid residues that are identical to FLT3 amino acid residues according to FIG. 1A, wherein the root mean square deviation between the set of amino acid residues and the FLT3 amino acid residues is not more than about 3.0 Å;
        comprising the steps of:
        (a) providing the structure coordinates of the binding pocket or domain on a computer comprising means for generating three-dimensional structural information from the structure coordinates;
        (b) using the computer to dock a first chemical entity in part of the binding pocket or domain;
        (c) docking at least a second chemical entity in another part of the binding pocket or domain;
        (d) quantifying the association between the first or second chemical entity and part of the binding pocket or domain;
        (e) repeating steps (b) to (d) with another first and second chemical entity,
        (f) selecting a first and a second chemical entity based on the quantified association of both the first and second chemical entity;
        (g) optionally, visually inspecting the relationship of the selected first and second chemical entity to each other in relation to the binding pocket or domain on a computer screen using the three-dimensional graphical representation of the binding pocket or domain and the first and second chemical entity; and
        (h) assembling the selected first and second chemical entity into a compound or complex that interacts with said binding pocket or domain by model building.
    
    


The method provides a method of utilizing molecular replacement to obtain structural information about a molecule or a molecular complex of unknown structure, wherein the molecule is sufficiently homologous to a cytoplasmic domain of an FLT3 protein, comprising the steps of:

    
    
        (a) crystallizing the molecule or molecular complex;
        (b) generating an X-ray diffraction pattern from the crystallized molecule or molecular complex; and
        (c) applying at least a portion of the structure coordinates set forth in FIG. 1A or a homology model thereof to the X-ray diffraction pattern to generate a three-dimensional electron density map of at least a portion of the molecule or molecular complex of unknown structure; and
        (d) generating a structural model of the molecule or molecular complex from the three-dimensional electron density map.
    
    


The invention also provides a method according to paragraph 27, wherein the molecule is selected from the group consisting of the FMS-like tyrosine kinase protein, and a homologue of a cytoplasmic domain of the FMS-like tyrosine kinase protein.
The invention provides the method according to paragraph 27, wherein the molecular complex is selected from the group consisting of the FMS-like tyrosine kinase protein complex and a homologue of the cytoplasmic domain of FMS-like tyrosine cytoplasmic complex.
The invention provides a method for identifying a candidate inhibitor that interacts with a binding site of an FMS-like tyrosine kinase protein or a homologue thereof, comprising the steps of:

    
    
        (a) obtaining a crystal comprising a cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof, wherein the crystal is characterized with space group P43212 and has unit cell parameters of a=b=80.67 Å, c=150.16 Å;
        (b) obtaining the structure coordinates of amino acids of the crystal of step (a), wherein the structure coordinates are set forth in FIG. 1A-1 to 1A-50;
        (c) generating a three-dimensional model of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof using the structure coordinates of the amino acids obtained in step (b), a root mean square deviation from backbone atoms of said amino acids of not more than ±2.0 Å;
        (d) determining a binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof from said three-dimensional model; and
        (e) performing computer fitting analysis to identify the candidate inhibitor which interacts with said binding site.
    
    


The invention provides the method according to paragraph 30, further comprising the step of: (f) contacting the identified candidate inhibitor with the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof in order to determine the effect of the inhibitor on FMS-like tyrosine kinase protein activity.
The invention provides the method according to paragraph 30, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined in step (d) comprises the structure coordinates according to FIG. 1A-1 to 1A-50 of amino acid residues H809, R810, and D811, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
The invention provides the method according to paragraph 30, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined in step (d) comprises the structure coordinates according to FIG. 1A-1 to 1A-50 of amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
The invention provides the method according to paragraph 30, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined in step (d) comprises the structure coordinates according to FIG. 1A-1 to 1A-50 of amino acid residues F621, K644, A657, E661, M664, L802, S806, C807, V808, H809, R810, D811, D829, and L832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
The invention provides a method for identifying a candidate inhibitor that interacts with a binding site of a cytoplasmic domain of an FMS-like tyrosine kinase protein or a homologue thereof, comprising the steps of:

    
    
        (a) obtaining a crystal comprising the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof, wherein the crystal is characterized with space group P43212 and has unit cell parameters of a=b=80.67 Å, c=150.16 Å;
        (b) obtaining the structure coordinates of amino acids of the crystal of step (a);
        (c) generating a three-dimensional model of said FMS-like tyrosine kinase protein or said homologue thereof using the structure coordinates of the amino acids generated in step (b), a root mean square deviation from backbone atoms of said amino acids of not more than ±2.0 Å;
        (d) determining a binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof from said three-dimensional model; and
        (e) performing computer fitting analysis to identify the candidate inhibitor which interacts with said binding site.
    
    


The invention provides the method according to paragraph 35, further comprising the step of:

    
    
        (f) contacting the identified candidate inhibitor with the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof in order to determine the effect of the inhibitor on FMS-like tyrosine kinase protein activity.
    
    


The invention provides the method according to paragraph 35, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined in step (d) comprises the structure coordinates according to FIG. 1A-1 to 1A-50 of amino acid residues H809, R810, and D811, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
The invention provides the method according to paragraph 35, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined in step (d) comprises the structure coordinates according to FIG. 1A-1 to 1A-50 of amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
The method according to paragraph 35, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined in step (d) comprises the structure coordinates according to FIG. 1A-1 to 1A-50 of amino acid residues F621, K644, A657, E661, M664, L802, S806, C807, V808, H809, R810, D811, D829, and L832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
The invention provides a method for identifying a candidate inhibitor that interacts with a binding site of a cytoplasmic domain of an FMS-like tyrosine kinase protein or a homologue thereof, comprising the step of determining a binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or the homologue thereof from a three-dimensional model to design or identify the candidate inhibitor which interacts with said binding site.
The invention provides the method according to paragraph 40, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined comprises the structure coordinates according to FIG. 1A-1 to 1A-50 of amino acid residues H809, R810, and D811, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
The invention provides the method according to paragraph 40, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined comprises the structure coordinates according to FIG. 1A-1 to 1A-50 of amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
The invention provides the method according to paragraph 40, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined comprises the structure coordinates according to FIG. 1A-1 to 1A-50 of amino acid residues F621, K644, A657, L658, E661, M664, L802, K805, S806, C807, V808, H809, R810, D811, C828, D829, F830, G831, and L832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
The invention provides a method for identifying a candidate inhibitor of a molecule or molecular complex comprising a binding pocket or domain selected from the group consisting of:

    
    
        (i) a set of amino acid residues which are identical to human FMS-like tyrosine kinase (FLT3) amino acid residues H809, R810, and D811 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between said set of amino acid residues and said FMS-like tyrosine kinase amino acid residues is not greater than about 2.0 Å;
        (ii) a set of amino acid residues comprising at least five amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between said at least five amino acid residues and said FMS-like tyrosine kinase amino acid residues which are identical is not greater than about 2.0 Å;
        (iii) a set of amino acid residues comprising at least six amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, K644, A657, L658, E661, M664, L802, K805, S806, C807, V808, H809, R810, D811, C828, D829, F830, G831, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between said at least six amino acid residues and said FMS-like tyrosine kinase amino acid residues which are identical is not greater than about 2.0 Å; and
        (iv) a set of amino acid residues that are identical to FMS-like tyrosine kinase amino acid residues according to FIG. 1A, wherein the root mean square deviation between said set of amino acid residues and said FMS-like tyrosine kinase amino acid residues is not more than about 3.0 Å;
    
    


comprising the steps of:

    
    
        (a) using a three-dimensional structure of the binding pocket or domain to design, select or optimize a plurality of chemical entities; and
        (b) selecting said candidate inhibitor based on the inhibitory effect of said chemical entities on a cytoplasmic domain of a FMS-like tyrosine kinase protein or a cytoplasmic domain of a FMS-like tyrosine kinase protein homologue on the catalytic activity of the molecule or molecular complex.
    
    





BRIEF DESCRIPTION OF THE FIGURES
The following abbreviations are used in FIG. 1A:
“Atom type” refers to the element whose coordinates are measured. The first letter in the column defines the element.
“Resid” refers to the amino acid residue in the molecular model.
“X, Y, Z” 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 the molecules.
FIG. 1A (1A-1 to 1A-50) lists the atomic coordinates for human FLT3 (amino acid residues 570-710 and 783-947 of the cytoplasmic domain of human FLT3 protein (GenBank accession no. NP—004110; SEQ ID NO:1)) as derived from X-ray diffraction. Residues 649-654 were not included in the final model. The coordinates are shown in Protein Data Bank (PDB) format. Residues “PO4 B”, “CAP C” and “TIP D” represent phosphate, CAPS (3-cyclohexylamino-1-propanesulfonic acid) and water molecules, respectively.
FIG. 2A depicts the structure of autoinhibited FLT3 using ribbon diagram (shown on left) and surface diagram (shown on right) and highlights the spatial arrangement of the various structural elements of the FLT3 molecule. The N-terminal kinase domain (labeled and shown in dark gray) and the C-terminal kinase domain (labeled and shown in medium gray) comprise the standard kinase fold. The activation loop (labeled in the right diagram) is folded up between the two kinase domains. The JM domain (labeled and shown in light gray) nearly spans the length of the molecule. All tyrosines in the JM domain and the activation loop are displayed as ‘stick’ representations.
FIG. 2B is the same as FIG. 2A except the molecule is rotated 90 degrees clockwise when viewed down the vertical axis.
FIG. 3A depicts a ribbon representation of FLT3 highlighting the JM domain. A molecular surface representation of the JM-B shows its central position relative to the activation loop, αC, glycine-rich loop, JM-S and JM-Z.
FIG. 3B depicts a stereoview of a ribbon representation of FLT3. The JM domain is displayed as a ‘stick’ model.
FIG. 4 shows a closeup view of the active site region of FLT3. Tyr572 and Tyr842 have hydrogen bonds to Glu661 and Asp811, respectively, which in turn are involved in salt bridges.
FIG. 5A shows a superposition of the closed activation loop from FLT3 on the closed IRK-I activation loop (dark gray).
FIG. 5B shows a superposition of the closed activation loop from FLT3 on the partially open FGFR activation loop (dark gray).
FIG. 5C shows a superposition of the closed activation loop from FLT3 on the fully open IRK-A activation loop (dark gray).
FIG. 6A shows the secondary structure assignments for FLT3 (amino acid residues 572-958 of the cytoplasmic domain (full length protein: GenBank accession no. NP—004110; SEQ ID NO:1)) in which α-helices and β-strands are represented by striped and solid bars, respectively. Dotted bars delineate secondary structures that are present in IRK3P (Hubbard et. al., 1997, supra), but not in autoinhibited FLT3.
FIG. 6B show a sequence alignment of the JM domain for FLT3 (amino acid residues 564-609 of the cytoplasmic domain (full length protein: GenBank accession no. NP—004110; SEQ ID NO:1)), PDGFRα (SEQ ID NO:2), PDGFRβ (SEQ ID NO:3), cKIT (SEQ ID NO:4), and FMS (SEQ ID NO:5). Identical and homologous residues are shaded in gray. Gray arrows above the sequence alignment represent β-strands. Asterisks denote location of residues in PDGFRβ which when mutated to alanine give rise to constitutive activation. These PDGFRβ mutations include: Y530A, W534A, V536A, 1537A, L555A, Y557A and the double mutant Y547A/Y549A (PDGFRβ numbers are given).
FIG. 7 shows a diagram of a system used to carry out the instructions encoded by the storage medium of FIGS. 8 and 9.
FIG. 8 shows a cross section of a magnetic storage medium.
FIG. 9 shows a cross section of a optically-readable data storage medium.



DETAILED DESCRIPTION OF THE INVENTION
In order that the invention described herein may be more fully understood, the following detailed description is set forth.
Throughout the specification, the word “comprise” or variations such as “comprises” or “comprising” will be understood to imply the inclusion of a stated integer or groups of integers but not the exclusion of any other integer or groups of integers.
The following abbreviations are used throughout the application:

















A =
Ala =
Alanine
T =
Thr =
Threonine


V =
Val =
Valine
C =
Cys =
Cysteine


L =
Leu =
Leucine
Y =
Tyr =
Tyrosine


I =
Ile =
Isoleucine
N =
Asn =
Asparagine


P =
Pro =
Proline
Q =
Gln =
Glutamine


F =
Phe =
Phenylalanine
D =
Asp =
Aspartic Acid


W =
Trp =
Tryptophan
E =
Glu =
Glutamic Acid


M =
Met =
Methionine
K =
Lys =
Lysine


G =
Gly =
Glycine
R =
Arg =
Arginine


S =
Ser =
Serine
H =
His =
Histidine









As used herein, the following definitions shall apply unless otherwise indicated.
The term “about” when used in the context of root mean square deviation (RMSD) values takes into consideration the standard error of the RMSD value, which is ±0.1 Å.
The term “associating with” refers to a condition of proximity between a chemical entity or compound, or portions thereof, and a binding pocket or binding site on a protein. The association may be non-covalent—wherein the juxtaposition is energetically favored by hydrogen bonding, hydrophobic, van der Waals or electrostatic interactions—or it may be covalent.
The term “binding pocket” refers to a region of a molecule or molecular complex, that, as a result of its shape, favorably associates with a chemical entity. The term “pocket” includes, but is not limited to, cleft, channel or site. FLT3, FLT3-like molecules or homologues thereof may have binding pockets which include, but are not limited to, peptide or substrate binding, JM-B binding and ATP-binding sites. The shape of a first binding pocket may be largely pre-formed before binding of a chemical entity, may be formed simultaneously with binding of a chemical entity, or may be formed by the binding of another chemical entity to a different binding pocket of the molecule, which in turn induces a change in shape of the first binding pocket.
The term “catalytic active site” or “active site” refers to the portion of the protein kinase to which nucleotide substrates bind. For example, the catalytic active site of FLT3 is at the interface between the N-lobe and the C-lobe.
The term “chemical entity” refers to chemical compounds, complexes of at least two chemical compounds, and fragments of such compounds or complexes. The chemical entity can be, for example, a ligand, substrate, nucleotide triphosphate, nucleotide diphosphate, phosphate, nucleotide, agonist, antagonist, inhibitor, antibody, peptide, protein or drug. In one embodiment, the chemical entity is an inhibitor or substrate for the active site of FLT3 proteins or protein complexes, or homologues thereof. The first and second chemical entities referred to in the present invention may be identical or distinct from each other. When iterative steps of using first and second chemical entities are carried out, taken as a pair, the first and second chemical entities used in repeated steps should be different from the first and second chemical entities of the steps.
The term “complex” or “molecular complex” refers to a protein associated with a chemical entity.
The term “conservative substitutions” refers to residues that are physically or functionally similar to the corresponding reference residues. That is, a conservative substitution and its reference residue have similar size, shape, electric charge, chemical properties including the ability to form covalent or hydrogen bonds, or the like. Preferred conservative substitutions are those fulfilling the criteria defined for an accepted point mutation in Dayhoff et. al., Atlas of Protein Sequence and Structure, 5: 345-352 (1978 & Supp.), which is incorporated herein by reference. Examples of conservative substitutions are substitutions including but not limited to the following groups: (a) valine, glycine; (b) glycine, alanine; (c) valine, isoleucine, leucine; (d) aspartic acid, glutamic acid; (e) asparagine, glutamine; (f) serine, threonine; (g) lysine, arginine, methionine; and (h) phenylalanine, tyrosine.
The term “contact score” refers to a measure of shape complementarity between the chemical entity and binding pocket, which is correlated with an RMSD value obtained from a least square superimposition between all or part of the atoms of the chemical entity and all or part of the atoms of the ligand bound (for example, AMP-PNP or an inhibitor) in the binding pocket according to FIG. 1 or 2. The docking process may be facilitated by the contact score or RMSD values. For example, if the chemical entity moves to an orientation with high RMSD, the system will resist the motion. A set of orientations of a chemical entity can be ranked by contact score. A lower RMSD value will give a higher contact score. See Meng et. al. J. Comp. Chem., 4, 505-524 (1992).
The term “correspond to” or “corresponding amino acids”, when used in the context of amino acid residues that correspond to FLT3 amino acids, refers to particular amino acids or analogues thereof in a tyrosine kinase that correspond to amino acids in the human FLT3 protein. The corresponding amino acid may be an identical, mutated, chemically modified, conserved, conservatively substituted, functionally equivalent or homologous amino acid, when compared to the FLT3 amino acid to which it corresponds. For example, the following are examples of FLT3 amino acid residues that correspond to cKIT amino acid residues: F570:P551 and L576:W557 (the identity of the FLT3 residue is listed first; its position is indicated using FLT3 sequence numbering; and the identity of cKIT residue is given at the end).
Methods for identifying a corresponding amino acid are known in the art and are based upon sequence, structural alignment, its functional position or a combination thereof, as compared to the FLT3 protein. For example, corresponding amino acids may be identified by superimposing the backbone atoms of the amino acids in FLT3 and another protein using well known software applications, such as QUANTA (Molecular Simulations, Inc., San Diego, Calif. ©1998, 2000; Accelrys ©2001, 2002). The corresponding amino acids may also be identified using sequence alignment programs such as the “bestfit” program or CLUSTAL W Alignment Tool (Higgins D. G., et. al., Methods Enzymol., 266: 383-402 (1996)).
The term “crystallization solution” refers to a solution which promotes crystallization comprising at least one agent, including a buffer, one or more salts, a precipitating agent, one or more detergents, sugars or organic compounds, lanthanide ions, a poly-ionic compound, and/or stabilizer.
The term “docking” refers to orienting, rotating, or translating a chemical entity in the binding pocket, domain, molecule or molecular complex or portion thereof based on distance geometry or energy. Docking may be performed by distance geometry methods that find sets of atoms of a chemical entity that match sets of sphere centers of the binding pocket, domain, molecule or molecular complex or portion thereof. See Meng et. al. J. Comp. Chem., 4, 505-524 (1992). Sphere centers are generated by providing an extra radius of given length from the atoms (excluding hydrogen atoms) in the binding pocket, domain, molecule or molecular complex or portion thereof. Real-time interaction energy calculations, energy minimizations or rigid-body minimizations (Gschwend, et. al., J. Mol. Recognition, 9:175-186 (1996)) can be performed during or after orientation of the chemical entity to facilitate docking. For example, interactive docking experiments can be designed to follow the path of least resistance. If the user in an interactive docking experiment makes a move to increase the energy, the system will resist that move. However, if that user makes a move to decrease energy, the system will favor that move by increased responsiveness. (Cohen, et. al., J. Med. Chem. 33:889-894 (1990)). Docking can also be performed by combining a Monte Carlo search technique with rapid energy evaluation using molecular affinity potentials. See Goodsell and Olson, Proteins: Structure, Function and Genetics 8:195-202 (1990). Software programs that carry out docking functions include but are not limited to MATCHMOL (Cory et. al., J Mol. Graphics, 2, 39 (1984); MOLFIT (Redington, Comput. Chem., 16, 217 (1992)) and DOCK (Meng et. al., supra).
The term “domain” refers to a structural unit of the FLT3 protein or homologue. The domain can comprise a binding pocket, a sequence or structural motif.
The term “full-length FLT3” refers to the complete human FLT3 protein, which includes an extracellular domain consisting of five immunoglobulin-like (Ig-like) domains, a single transmembrane region, and a cytoplasmic domain (amino acid residues 1 to 993; GenBank accession no. NP—004110; SEQ ID NO:1). The cytoplasmic domain includes a cytoplasmic juxtamembrane domain (JM) and a cytoplasmic kinase domain interrupted by a kinase insert domain (KID).
The term “FLT3-like” refers to all or a portion of a molecule or molecular complex that has a commonality of shape with all or a portion of the FLT3 protein. For example, in the FLT3-like JM-B binding pocket, the commonality of shape is defined by a root mean square deviation of the structure coordinates of the backbone atoms between the amino acids in the FLT3-like JM-B binding pocket and the FLT3 amino acids in the FLT3 JM-B binding pocket (as set forth in FIG. 1A). Compared to the amino acids of the FLT3 binding pocket, the corresponding amino acid residues in the FLT3-like binding pocket may or may not be identical. Depending on the set of FLT3 amino acid residues that define the FLT3 JM-B binding pocket, one skilled in the art would be able to locate the corresponding amino acids that define an FLT3-like binding pocket in a protein based on sequence or structural homology.
The term “FLT3 protein complex” or “FLT3 homologue complex” refers to a molecular complex formed by associating the FLT3 protein or FLT3 homologue with a chemical entity, for example, a ligand, a substrate, nucleotide triphosphate, nucleotide diphosphate, phosphate, an agonist or antagonist, inhibitor, antibody, drug or compound.
The term “generating a three-dimensional structure” or “generating a three-dimensional representation” refers to converting the lists of structure coordinates into structural models or graphical representations in three-dimensional space. This can be achieved through commercially or publicly available software. A model of a three-dimensional structure of a molecule or molecular complex can thus be constructed on a computer screen by a computer that is given the structure coordinates and that comprises the correct software. The three-dimensional structure may be displayed or used to perform computer modeling or fitting operations. In addition, the structure coordinates themselves, without the displayed model, may be used to perform computer-based modeling and fitting operations.
The term “homologue of FLT3 cytoplasmic domain” or “FLT3 cytoplasmic domain homologue” refers to the cytoplasmic domain of a protein in the type III class of receptor tyrosine kinases that is at least 70%, 80%, 90%, 95%, 96%, 97%, 98%, 99% or greater than 99% identical in sequence to the cytoplasmic domain (approximately amino acid residues 564-993 of SEQ ID NO:1) of human FLT3 protein and retains FLT3 kinase activity. In one embodiment, the homologue is at least 95%, 96%, 97%, 98% or 99% identical in sequence to human FLT3 cytoplasmic domain, and has conservative mutations as compared to human FLT3 cytoplasmic domain. The homologue can be an FLT3 cytoplasmic domain from another species, or the foregoing human FLT3 cytoplasmic domain with mutations, conservative substitutions, additions, deletions or a combination thereof. Such animal species include, but are not limited to, mouse, rat, a primate such as monkey or other primates.
The term “homology model” refers to a structural model derived from known three-dimensional structure(s). Generation of the homology model, termed “homology modeling”, can include sequence alignment, residue replacement, residue conformation adjustment through energy minimization, or a combination thereof.
The term “interaction energy” refers to the energy determined for the interaction of a chemical entity and a binding pocket, domain, molecule or molecular complex or portion thereof. Interactions include but are not limited to one or more of covalent interactions, non-covalent interactions such as hydrogen bond, electrostatic, hydrophobic, aromatic, van der Waals interactions, and non-complementary electrostatic interactions such as repulsive charge-charge, dipole-dipole and charge-dipole interactions. As interaction energies are measured in negative values, the lower the value the more favorable the interaction.
The term “JM-B binding motif” or “JM-B” refers to a short segment in the FLT3 cytoplasmic domain or homologue thereof including amino acid residues Tyr572-Met578.
The term “JM-B binding pocket” refers to the binding pocket for the JM-B binding motif. The JM-B binding pocket comprises the kinase domain amino acid residues found within the near vicinity of the heptapeptide section of the autoinhibitory JM-B binding motif of FLT3 cytoplasmic domain or homologue thereof (amino acid residues YESQLQM (residues 572-578 of SEQ ID NO:1)).
The term “motif” refers to a group of amino acid residues in the FLT3 protein or homologue that defines a structural compartment or carries out a function in the protein or homologue, for example, catalysis or structural stabilization, or phosphorylation. The motif may be conserved in sequence, structure and function. The motif can be contiguous in primary sequence or three-dimensional space. An example of a motif includes but is not limited to the activation loop.
The term “part of a binding pocket” refers to less than all of the amino acid residues that define the binding pocket. The structure coordinates of amino acid residues that constitute part of a binding pocket may be specific for defining the chemical environment of the binding pocket, or useful in designing fragments of an inhibitor that may interact with those residues. For example, the portion of amino acid residues may be key residues that play a role in ligand binding, or may be residues that are spatially related and define a three-dimensional compartment of the binding pocket. The amino acid residues may be contiguous or non-contiguous in primary sequence. In one embodiment, part of the binding pocket has at least two amino acid residues, preferably at least three, eight, fourteen or fifteen amino acid residues.
The term “part of an FLT3 protein” or “part of an FLT3 homologue” refers to less than all of the amino acid residues of an FLT3 protein or homologue. In one embodiment, part of the FLT3 protein or homologue defines the binding pockets, domains, sub-domains, and motifs of the protein or homologue. The structure coordinates of amino acid residues that constitute part of an FLT3 protein or homologue may be specific for defining the chemical environment of the protein, or useful in designing fragments of an inhibitor that interact with those residues. The portion of amino acid residues may also be residues that are spatially related and define a three-dimensional compartment of the binding pocket, motif or domain. The amino acid residues may be contiguous or non-contiguous in primary sequence. For example, the portion of amino acid residues may be key residues that play a role in ligand or substrate binding, peptide binding, antibody binding, catalysis, structural stabilization or degradation.
The term “quantified association” refers to calculations of distance geometry and energy. Energy can include but is not limited to interaction energy, free energy and deformation energy. See Cohen, supra.
The term “root mean square deviation” or “RMSD” refers to 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 FLT3, a binding pocket, a motif, a domain, or portion thereof, as defined by the structure coordinates of FLT3 described herein. It would be readily apparent to those skilled in the art that the calculation of RMSD involves standard error of ±0.1 Å.
The term “soaked” refers to a process in which a crystal is transferred to a solution containing a compound of interest.
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 a protein or protein 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 molecule or molecular complex.
The term “sub-domain” refers to a portion of a domain.
The term “substantially all of an FLT3 binding pocket” or “substantially all of an FLT3 protein” refers to all or almost all of the amino acids in the FLT3 binding pocket or protein. For example, substantially all of an FLT3 binding pocket can be 100%, 95%, 90%, 80%, or 70% of the residues defining the FLT3 binding pocket or protein.
The term “substrate binding pocket” refers to the binding pocket for a substrate of FLT3 or homologue thereof. A substrate is generally defined as the molecule upon which an enzyme performs catalysis. Natural substrates, synthetic substrates or peptides, or mimics of a natural substrates of FLT3 or homologue thereof may associate with the substrate binding pocket.
The term “sufficiently homologous to FLT3” refers to a protein that has a sequence identity of at least 25% compared to FLT3 protein. In other embodiments, the sequence identity is at least 40%. In other embodiments, the sequence identity is at least 50%, 60%, 70%, 80%, 90%, 95%, 96%, 97%, 98% or 99%.
The term “three-dimensional structural information” refers to information obtained from the structure coordinates. Structural information generated can include the three-dimensional structure or graphical representation of the structure. Structural information can also be generated when subtracting distances between atoms in the structure coordinates, calculating chemical energies for an FLT3 molecule or molecular complex or homologues thereof, calculating or minimizing energies for an association of an FLT3 molecule or molecular complex, or homologues thereof to a chemical entity.
Crystallizable Compositions and Crystals of an FLT3 Cytoplasmic Domain and Complexes Thereof
In one embodiment, the invention provides a crystallizable composition comprising an FLT3 cytoplasmic domain or its homologue. In another embodiment, the crystallizable composition further comprises a buffer that maintains pH between about 8.0 and 12.0, 0.1-5 M sodium phosphate, 0.1-5 M potassium phosphate, and 50-500 mM lithium sulfate. In certain embodiments, the crystallizable composition comprises equal volumes of a solution of an FLT3 cytoplasmic domain or a homologue thereof (10 mg/ml) in the presence of 2 mM AMP-PNP and 4 mM MgCl2, and 1.2 M sodium phosphate, 1.2 M potassium phosphate, 0.1 M CAPS (pH 10.5) and 0.2 M lithium sulfate.
According to another embodiment, the invention provides a crystal comprising an FLT3 cytoplasmic domain or its homologue. Preferably, the native crystal has a unit cell dimension of a=b=80.65 Å c=150.13 Å and belongs to space group P43212. It will be readily apparent to those skilled in the art that the unit cells of such a crystal composition may deviate ±1-2 Å from the above cell dimensions depending on the deviation in the unit cell calculations.
As used herein, the FLT3 cytoplasmic domain in the crystallizable compositions or crystals can be amino acids 564-993 of SEQ ID NO:1; amino acids 564-958 of SEQ ID NO:1; amino acids 564-947 of SEQ ID NO:1; (removing residues spanning H711-H761) amino acids 564-710 and 762-958 of SEQ ID NO:1; (removing residues spanning 711-782) amino acids 564-710 and 783-958 of SEQ ID NO:1; (removing residues spanning H711-H761) amino acids 570-710 and 762-958 of SEQ ID NO:1; (removing residues spanning 711-782) amino acids 570-710 and 783-958 of SEQ ID NO:1; amino acid residues 570-958 of SEQ ID NO:1; amino acid residues 570-947 of SEQ ID NO:1; (removing residues spanning H711-H761) amino acid residues 570-710 and 762-947 of SEQ ID NO:1; (removing residues spanning 711-782) amino acid residues 570-710 and 783-947 of SEQ ID NO:1; (removing residues spanning H711-H761) amino acid residues 564-710 and 762-947 of SEQ ID NO:1; and (removing residues spanning 711-782) amino acid residues 564-710 and 783-947 of SEQ ID NO:1. The homologue thereof can be any of the aforementioned amino acids with conservative substitutions, deletions or additions, to the extent that any substitutions, deletions or additions maintains an FLT3 kinase activity in the homologue; preferably the homologue with substitutions, deletions or additions is at least 70%, 80%, 90%, 95%, 96%, 97%, 98%, or 99% identical to one of the aforementioned. Preferably, the homologue with substitutions, deletions or additions is at least 80%, 90%, 95%, 96%, 97%, 98%, or 99% identical to one of the aforementioned.









1




MPALARDAGT VPLLVVFSAM IFGTITNQDL PVIKCVLINH KNNDSSVGKS
SEQ ID NO:1


 


51


SSYPMVSESP EDLGCALRPQ SSGTVYEAAA VEVDVSASIT LQVLVDAPGN


 


101


ISCLWVFKHS SLNCQPHFDL QNRGVVSMVI LKMTETQAGE YLLFIQSEAT


 


151


NYTILFTVSI RNTLLYTLRR PYFRKMENQD ALVCISESVP EPIVEWVLCD


 


201


SQGESCKEES PAVVKKEEKV LHELFGTDIR CCARNELGRE CTRLFTIDLN


 


251


QTPQTTLPQL FLKVGEPLWI RCKAVHVNHG FGLTWELENK ALEEGNYFEM


 


301


STYSTNRTMI RILFAFVSSV ARNDTGYYTC SSSKHPSQSA LVTIVGKGFI


 


351


NATNSSEDYE IDQYEEFCFS VRFKAYPQIR CTWTFSRKSF PCEQKGLDNG


 


401


YSISKFCNHK HQPGEYIFHA ENDDAQFTKM FTLNIRRKPQ VLAEASASQA


 


451


SCFSDGYPLP SWTWKKCSDK SPNCTEEITE GVWNRKANRK VFGQWVSSST


 


501


LNMSEAIKGF LVKCCAYNSL GTSCETILLN SPGPFPFIQD NISFYATIGV


 


551


CLLFIVVLTL LICHKYKKQF RYESQLQMVQ VTGSSDNEYF YVDFREYEYD


 


601


LKWEFPRENL EFGKVLGSGA FGKVMNATAY GISKTGVSIQ VAVKMLKEKA


 


651


DSSEREALMS ELKMMTQLGS HENIVNLLGA CTLSGPIYLI FEYCCYGDLL


 


701


NYLRSKREKF HRTWTEIFKE HNFSFYPTFQ SHPNSSMPGS REVQIHPDSD


 


751


QISGLHGNSF HSEDEIEYEN QKRLEEEEDL NVLTFEDLLC FAYQVAKGME


 


801


FLEFKSCVHR DLAARNVLVT HGKVVKICDF GLARDIMSDS NYVVRGNARL


 


851


PVKWMAPESL FEGIYTIKSD VWSYGILLWE IFSLGVNPYP GIPVDANFYK


 


901


LIQNGFKMDQ PFYATEEIYI IMQSCWAFDS RKRPSFPNLT SFLGCQLADA


 


951


EEAMYQNVDG RVSECPHTYQ NRRPFSREMD LGLLSPQAQV EDS






The FLT3 protein or its homologue may be produced by any well-known method, including synthetic methods, such as solid phase, liquid phase and combination solid phase/liquid phase syntheses; recombinant DNA methods, including cDNA cloning, optionally combined with site directed mutagenesis; and/or purification of the natural products.
Methods of Obtaining Crystals of an FLT3 Cytoplasmic Domain or Its Homologues
The invention also relates to a method of obtaining a crystal of an FLT3 cytoplasmic domain or homologue thereof, comprising the steps of:

    
    
        a) optionally producing and purifying an FLT3 cytoplasmic domain or homologue thereof;
        b) combining a crystallization solution with said FLT3 cytoplasmic domain or homologue thereof to produce a crystallizable composition; and
        c) subjecting the composition to conditions which promote crystallization and obtaining said crystal.
    
    


In another embodiment, the invention provides methods of obtaining crystals of an FLT3 cytoplasmic domain protein, a homologue thereof, or complexes thereof using the steps set forth above. In one embodiment, step (b) is performed with an FLT3 cytoplasmic domain or homologue thereof bound to a chemical entity. In another embodiment, the above method further comprises the step of soaking said crystal in a solution comprising a chemical entity that binds to the FLT3 cytoplasmic domain or homologue thereof.
In certain embodiments, the method of making crystals of an FLT3 cytoplasmic domain, a homologue, or an FLT3 cytoplasmic domain protein or homologue complex includes the use of a device for promoting crystallizations. Devices for promoting crystallization can include but are not limited to the hanging-drop, sitting-drop, sandwich-drop, dialysis, microbatch or microtube batch devices (U.S. Pat. Nos. 4,886,646, 5,096,676, 5,130,105, 5,221,410 and 5,400,741; Pav, S., et. al., Proteins Struct. Funct. Genet., 20: 98-102 (1994); Chayen, Acta. Cryst., D54: 8-15 (1998), Chayen, Structure, 5: 1269-1274 (1997), D'Arcy et. al., J. Cryst. Growth, 168: 175-180 (1996) and Chayen, J. Appl. Cryst., 30: 198-202 (1997), incorporated herein by reference). The hanging-drop, sitting-drop and some adaptations of the microbatch methods (D'Arcy et. al., J. Cryst. Growth, 168: 175-180 (1996) and Chayen, J. Appl. Cryst., 30: 198-202 (1997)) produce crystals by vapor diffusion. The hanging drop and sitting drop containing the crystallizable composition is equilibrated against a reservoir containing a higher or lower concentration of precipitant. As the drop approaches equilibrium with the reservoir, the saturation of protein in the solution leads to the formation of crystals.
Microseeding may be used to increase the size and quality of crystals. In this instance, microcrystals are crushed to yield a stock seed solution. The stock seed solution is diluted in series. Using a needle, glass rod, micro-pipet, micro-loop or strand of hair, a small sample from each diluted solution is added to a set of equilibrated drops containing a protein concentration equal to or less than a concentration needed to create crystals without the presence of seeds. The aim is to end up with a single seed crystal that will act to nucleate crystal growth in the drop.
It would be readily apparent to one of skill in the art to vary the crystallization conditions disclosed above to identify other crystallization conditions that would produce crystals of FLT3 protein, FLT3 protein complex, FLT3 cytoplasmic domain protein complex or homologue thereof, or FLT3 cytoplasmic domain homologue. Such variations include, but are not limited to, adjusting pH, protein concentration and/or crystallization temperature, changing the identity or concentration of salt and/or precipitant used, using a different method for crystallization, or introducing additives such as detergents (e.g., TWEEN 20 (monolaurate), LDOA, Brji 30 (4 lauryl ether)), sugars (e.g., glucose, maltose), organic compounds (e.g., dioxane, dimethylformamide), lanthanide ions, or poly-ionic compounds that aid in crystallizations. High throughput crystallization assays may also be used to assist in finding or optimizing the crystallization condition.
In certain embodiments, the crystal comprising a cytoplasmic domain of an FMS-like tyrosine kinase protein or a homologue thereof diffract X-rays to a resolution of at least 2.1 Å. In other embodiments, the crystal comprising a cytoplasmic domain of an FLT3 cytoplasmic domain, a homologue, or an FLT3 cytoplasmic domain protein or homologue complex diffract X-rays to a resolution of at least 5.0 Å, at least 3.5 Å, at least 3.0 Å, at least 2.5 Å, or at least 2.3 Å.
In certain embodiments, the crystal comprising a cytoplasmic domain of an FMS-like tyrosine kinase protein, a homologue thereof, or complexes thereof can produce an electron density map having resolution of at least 2.1 Å. In other embodiments, the crystal comprising a cytoplasmic domain of an FLT3 cytoplasmic domain, a homologue, or an FLT3 cytoplasmic domain protein or homologue complex can produce an electron density map having resolution of at least 5.0 Å, at least 3.5 Å, at least 3.0 Å, at least 2.5 Å, or at least 2.3 Å.
In certain embodiments, the electron density map produced above are sufficient to determine the atomic coordinates a cytoplasmic domain of an FMS-like tyrosine kinase protein or a homologue thereof.
Binding Pockets of FLT3 Protein or Its Homologues
As disclosed herein, applicants have provided the three-dimensional X-ray structure of FLT3 containing the complete JM domain in autoinhibited form. The atomic coordinate data is presented in FIG. 1A.
To use the structure coordinates generated for the FLT3 cytoplasmic domain or one of its binding pockets or an FLT3-like binding pocket, it may be necessary to convert the structure coordinates, or portions thereof, into a three-dimensional shape (i.e., a three-dimensional representation of these proteins and binding pockets). This is achieved through the use of a computer comprising commercially available software that is capable of generating three-dimensional representations or structures of molecules or molecular complexes, or portions thereof, from a set of structure coordinates. These three-dimensional representations may be displayed on a computer screen.
Binding pockets, also referred to as binding sites in the present invention, 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 part 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 the binding pockets of biologically important targets. The binding pockets of this invention are useful and important for drug design.
The conformations of FLT3 and other proteins at a particular amino acid site, along the polypeptide backbone, can be compared using well-known procedures for performing sequence alignments of the amino acids. Such sequence alignments allow for the equivalent sites on these proteins to be compared. Such methods for performing sequence alignment include, but are not limited to, the “bestfit” program and CLUSTAL W Alignment Tool, Higgins et. al., supra.
The JM-B binding pocket comprises the kinase domain amino acid residues found within the near vicinity of the heptapeptide section of the autoinhibitory JM-B binding motif of FLT3 (amino acid residues YESQLQM (residues 572-578 of SEQ ID NO:1)).
In one embodiment, the JM-B binding pocket comprises amino acid residues F621, K644, A657, E661, M664, L802, S806, C807, V808, H809, R810, D811, D829, and L832 according to the structure of the autoinhibited conformation of FLT3 in FIG. 1A. The above-identified amino acid residues were within 5 Å (“5 Å sphere amino acids”) of the heptapeptide section of the autoinhibitory JM-B binding motif of FLT3 (amino acid residues YESQLQM (residues 572-578 of SEQ ID NO:1)) in the binding pocket. These residues were identified using the program InsightII, which allow the display of the structure, and a software program to calculate the residues within 5 Å of the heptapeptide section of the autoinhibitory JM-B binding motif of FLT3 (amino acid residues YESQLQM (residues 572-578 of SEQ ID NO:1)). QUANTA (Molecular Simulations, Inc., San Diego, Calif. ©1998, 2000; Accelrys ©2001, 2002), O (T. A. Jones et. al., Acta Cryst., A47: 110-119 (1991)) and RIBBONS (Carson, J. Appl. Cryst., 24: 958-961 (1991)) may also be used to obtain the above residues.
In another embodiment, the JM-B binding pocket comprises amino acids F621, K644, A657, L658, E661, M664, L802, K805, S806, C807, V808, H809, R810, D811, C828, D829, F830, G831, and L832 according to the structure of FLT3 protein in FIG. 1A. These amino acid residues were within 8 Å (“8 Å sphere amino acids”) of heptapeptide section of the autoinhibitory JM-B binding motif of FLT3 (amino acid residues YESQLQM (residues 572-578 of SEQ ID NO:1)). These residues were identified using the above-mentioned programs. QUANTA, O and RIBBONS, supra may also be used to obtain the above residues.
In another embodiment, the JM-B binding pocket comprises amino acids F621, K644, A657, E661, M664, L802, S806, C807, V808, H809, R810, D811, D829, and L832 according to the structure of FLT3 protein in FIG. 1A. These amino acid residues within 3.8 Å of the JM-B binding motif (amino acid residues 572 to 578). These residues were identified using the program QUANTA (Molecular Simulations, Inc., San Diego, Calif. ©1998, 2000; Accelrys ©2001, 2002).
In another embodiment, the JM-B binding pocket comprises amino acids F621, E661, M664, L802, V808, R810, D829, and L832 according to the structure of FLT3 protein in FIG. 1A. These amino acid residues make contacts less than 3.8 Å in length with residues of the JM-B binding motif (F621, M664, L802 and L832 have hydrophobic interactions or van der Waals contacts; E661, V808, R810 and D829 form hydrogen bonds; and R810 is also involved in an ionic pairing with Glu573). These residues were identified using the program QUANTA (Molecular Simulations, Inc., San Diego, Calif. ©1998, 2000; Accelrys ©2001, 2002).
In another embodiment, the JM-B binding pocket comprises amino acids H809, R810, D811, D829, F830 and G831 according to the structure of FLT3 protein in FIG. 1A.
In another embodiment, the JM-B binding pocket comprises amino acids H809, R810, and D811 according to the structure of FLT3 protein in FIG. 1A.
It will be readily apparent to those of skill in the art that the numbering of amino acid residues in homologues of human FLT3 may be different than that set forth for human FLT3. Corresponding amino acid residues in homologues of FLT3 are easily identified by visual inspection of the amino acid sequences or by using commercially available homology software programs. Homologues of FLT3 include, for example, FLT3 from other species, such as non-humans primates, mouse, rat, etc.
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 variations in coordinates discussed above may be generated because of mathematical manipulations of the FLT3 structure coordinates. For example, the structure coordinates set forth in FIG. 1A could undergo crystallographic permutations, fractionalization, integer additions or subtractions, inversion, 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 a certain root mean square deviation as compared to the original coordinates, the resulting three-dimensional shape is considered encompassed by this invention. Thus, for example, a ligand that bound to the binding pocket of FLT3 would also be expected to bind to another binding pocket whose structure coordinates defined a shape that fell within the acceptable root mean square deviation.
Various computational analyses may be necessary to determine whether a molecule or the binding pocket or portion thereof is sufficiently similar to the FLT3 binding pockets described above. Such analyses may be carried out using well known software applications, such as ProFit (A. C. R. Martin, SciTech Software, ProFit version 1.8, University College London, www.bioinf.org.uk/software), Swiss-Pdb Viewer (Guex et. al., Electrophoresis, 18: 2714-2723 (1997)), the Molecular Similarity application of QUANTA (Molecular Simulations, Inc., San Diego, Calif. ©1998, 2000; Accelrys ©2001, 2002) and as described in the accompanying User's Guide, which are incorporated herein by reference.
The above programs permit comparisons between different structures, different conformations of the same structure, and different parts of the same structure. The procedure used in QUANTA (Molecular Simulations, Inc., San Diego, Calif. ©1998, 2000; Accelrys ©2001, 2002) and Swiss-Pdb Viewer 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 on the structures; and 4) analyze the results.
The procedure used in ProFit to compare structures includes the following steps: 1) load the structures to be compared; 2) specify selected residues of interest; 3) define the atom equivalences in the selected residues; 4) perform a fitting operation on the selected residues; and 5) analyze the results.
Each structure in the comparison 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 (Molecular Simulations, Inc., San Diego, Calif. ©1998, 2000; Accelrys ©2001, 2002) is defined by user input, for the purpose of this invention we will define equivalent atoms as protein backbone atoms N, C, O and Cα for all corresponding amino acids between the two structures being compared.
The corresponding amino acids may be identified by sequence alignment programs such as the “bestfit” program available from the Genetics Computer Group which uses the local homology algorithm described by Smith and Waterman in Advances in Applied Mathematics 2, 482-489 (1981), which is incorporated herein by reference. A suitable amino acid sequence alignment will require that the proteins being aligned share a minimum percentage of identical amino acids. Generally, a first protein being aligned with a second protein should share in excess of about 35% identical amino acids (Hanks, S. K., et. al., Science, 241, 42-52 (1988); Hanks, S. K. and Quinn, A. M. Methods in Enzymology, 200: 38-62 (1991)). The identification of equivalent residues can also be assisted by secondary structure alignment, for example, aligning the α-helices, β-sheets in the structure. The program Swiss-Pdb Viewer has its own best fit algorithm that is based on secondary sequence alignment.
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 the above programs. The Swiss-Pdb Viewer program sets an RMSD cutoff for eliminating pairs of equivalent atoms that have high RMSD values. An RMSD cutoff value can be used to exclude pairs of equivalent atoms with extreme individual RMSD values. In the program ProFit, the RMSD cutoff value can be specified by the user.
For the purpose of this invention, any molecule, molecular complex, binding pocket, motif, domain thereof or portion thereof that is within a root mean square deviation for backbone atoms (N, Cα, C, O) when superimposed on the relevant backbone atoms described by structure coordinates listed in FIG. 1A are encompassed by this invention.
One embodiment of this invention provides a crystalline molecule comprising a protein defined by structure coordinates of a set of amino acid residues that are identical to FLT3 amino acid residues according to FIG. 1A, wherein the RMSD between said set of amino acid residues and said FLT3 amino acid residues is not more than about 5.0 Å. In other embodiments, the RMSD between said set of amino acid residues and said FLT3 amino acid residues is not greater than about 4.0 Å, not greater than about 3.0 Å, not greater than about 2.0 Å, not greater than about 1.5 Å, not greater than about 1.0 Å, or not greater than about 0.5 Å.
In one embodiment, the present invention provides a crystalline molecule comprising all or part of a binding pocket defined by a set of amino acid residues comprising at least six amino acid residues which are identical to human FLT3 amino acid residues F621, K644, A657, L658, E661, M664, L802, K805, S806, C807, V808, H809, R810, D811, C828, D829, F830, G831, and L832 according to FIG. 1A, wherein the RMSD of the backbone atoms between said FLT3 amino acid residues and said at least six amino acid residues which are identical is not greater than about 3.0 Å. In other embodiments, the RMSD is not greater than about 2.0 Å, 1.0 Å, 0.8, 0.5 Å, 0.3 Å, or 0.2 Å. In other embodiments, the binding pocket is defined by a set of amino acid residues comprising at least four, six, eight, twelve or fifteen amino acid residues which are identical to said FLT3 amino acid residues.
In one embodiment, the present invention provides a crystalline molecule comprising all or part of a binding pocket defined by a set of amino acid residues which are identical to human FLT3 amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832 according to FIG. 1A, wherein the RMSD of the backbone atoms between said FLT3 amino acid residues and said set of amino acid residues which are identical is not greater than about 3.0 Å. In other embodiments, the RMSD is not greater than about 2.0 Å, 1.0 Å, 0.8, 0.5 Å, 0.3 Å, or 0.2 Å. In other embodiments, the binding pocket is defined by a set of amino acid residues comprising at least four, five, six or seven amino acid residues which are identical to said FLT3 amino acid residues.
In one embodiment, the present invention provides a crystalline molecule comprising all or part of a binding pocket defined by a set of amino acid residues comprising a set of amino acid residues which are identical to human FLT3 amino acid residues H809, R810, and D811 according to FIG. 1A, wherein the RMSD of the backbone atoms between said FLT3 amino acid residues and said set of amino acid residues which are identical is not greater than about 3.0 Å. In other embodiments, the RMSD is not greater than about 2.0 Å, 1.0 Å, 0.8, 0.5 Å, 0.3 Å, or 0.2 Å.
In one embodiment, the above molecule is FLT3 protein, FLT3 cytoplasmic domain or homologues thereof. In another embodiment, the above molecules are in crystalline form. An FLT3 protein may be human FLT3. Homologues of human FLT3 can be FLT3 from another species, such as a mouse, a rat or a non-human primate.
Computer Systems
According to another embodiment, this invention provides a machine-readable data storage medium, comprising a data storage material encoded with machine-readable data, wherein said data defines the above-mentioned molecules or molecular complexes or binding pockets thereof. In one embodiment, the data defines the above-mentioned binding pockets by comprising the structure coordinates of said amino acid residues according to FIG. 1A. To use the structure coordinates generated for FLT3, homologues thereof, or one of its binding pockets, it is at times necessary to convert them into a three-dimensional shape or to extract three-dimensional structural information from them. This is achieved through the use of commercially or publicly available software that is capable of generating a three-dimensional structure or a three-dimensional representation of molecules or portions thereof from a set of structure coordinates. In one embodiment, three-dimensional structure or representation may be displayed graphically.
Therefore, according to another embodiment, this invention provides a machine-readable data storage medium comprising a data storage material encoded with machine readable data. In one embodiment, a machine programmed with instructions for using said data is capable of generating a three-dimensional structure or three-dimensional representation of any of the molecules, or molecular complexes or binding pockets thereof, that are described herein.
This invention also provides a computer comprising:

    
    
        (a) a machine-readable data storage medium, comprising a data storage material encoded with machine-readable data, wherein said data defines any one of the above molecules or molecular complexes;
        (b) a working memory for storing instructions for processing said machine-readable data;
        (c) a central processing unit (CPU) coupled to said working memory and to said machine-readable data storage medium for processing said machine readable data and means for generating three-dimensional structural information of said molecule or molecular complex; and
        (d) output hardware coupled to said central processing unit for outputting three-dimensional structural information of said molecule or molecular complex, or information produced by using said three-dimensional structural information of said molecule or molecular complex.
    
    


In one embodiment, the data defines the binding pocket of the molecule or molecular complex.
Three-dimensional data generation may be provided by an instruction or set of instructions, such as a computer program or commands for generating a three-dimensional structure or graphical representation from structure coordinates, or by subtracting distances between atoms, calculating chemical energies for an FLT3 molecule or molecular complex or homologues thereof, or calculating or minimizing energies for an association of an FLT3 molecule or molecular complex or homologues thereof to a chemical entity. The graphical representation can be generated or displayed by commercially available software programs. Examples of software programs include but are not limited to QUANTA (Accelrys ©2001, 2002), O (Jones et. al., Acta Crystallogr. A47: 110-119 (1991)) and RIBBONS (Carson, J. Appl. Crystallogr., 24: 9589-961 (1991)), which are incorporated herein by reference. Certain software programs may imbue this representation with physico-chemical attributes which are known from the chemical composition of the molecule, such as residue charge, hydrophobicity, torsional and rotational degrees of freedom for the residue or segment, etc. Examples of software programs for calculating chemical energies are described in the Rational Drug Design section.
Information about said binding pocket or information produced by using said binding pocket can be outputted through display terminals, touchscreens, facsimile machines, modems, CD-ROMs, printers, a CD or DVD recorder, ZIPTM or JAZTM drives or disk drives. The information can be in graphical or alphanumeric form.
In one embodiment, the computer is executing an instruction such as a computer program for generating three-dimensional structure or docking. In another embodiment, the computer further comprises a commercially available software program to display the information as a graphical representation. Examples of software programs include but as not limited to, QUANTA (Accelrys ©2001, 2002), O (Jones et. al., Acta Crystallogr. A47: 110-119 (1991)) and RIBBONS (Carson, J. Appl. Crystallogr., 24: 9589-961 (1991)), all of which are incorporated herein by reference.
FIG. 7 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, CD-ROM drives or DVD-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 bi-directional system bus (50).
Input hardware (35), 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 (35) may comprise CD-ROM or DVD-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 (Molecular Simulations, Inc., San Diego, Calif. ©1998, 2000; Accelrys ©2001, 2002) as described herein. Output hardware may 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. Output hardware may also include a display terminal, touchscreens, facsimile machines, modems, a CD or DVD recorder, ZIPTM or JAZTM drives, disk drives, or other machine-readable data storage device.
In operation, CPU (20) coordinates the use of the various input and output devices (35), (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. 8 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. 7. 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. 7.
FIG. 9 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. 7. 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.
In one embodiment, the structure coordinates of said molecules or molecular complexes or binding pockets are produced by homology modeling of at least a portion of the structure coordinates of FIG. 1A. Homology modeling can be used to generate structural models of FLT3 homologues or other homologous proteins based on the known structure of FLT3 cytoplasmic domain. This can be achieved by performing one or more of the following steps: performing sequence alignment between the amino acid sequence of a molecule (possibly an unknown molecule) against the amino acid sequence of FLT3; identifying conserved and variable regions by sequence or structure; generating structure coordinates for structurally conserved residues of the unknown structure from those of FLT3; generating conformations for the structurally variable residues in the unknown structure; replacing the non-conserved residues of FLT3 with residues in the unknown structure; building side chain conformations; and refining and/or evaluating the unknown structure.
Software programs that are useful in homology modeling include XALIGN (Wishart, D. S., et. al., Comput. Appl. Biosci., 10: 687-88 (1994)) and CLUSTAL W Alignment Tool, Higgins et. al., supra. See also, U.S. Pat. No. 5,884,230. These references are incorporated herein by reference.
To perform the sequence alignment, programs such as the “bestfit” program available from the Genetics Computer Group (Waterman in Advances in Applied Mathematics 2, 482 (1981), which is incorporated herein by reference) and CLUSTAL W Alignment Tool (Higgins et. al., supra, which is incorporated by reference) can be used. To model the amino acid side chains of homologous molecules, the amino acid residues in FLT3 can be replaced, using a computer graphics program such as “O” (Jones et al, (1991) Acta Cryst. Sect. A, 47: 110-119), by those of the homologous protein, where they differ. The same orientation or a different orientation of the amino acid can be used. Insertions and deletions of amino acid residues may be necessary where gaps occur in the sequence alignment. However, certain portions of the active site of FLT3 and its homologues are highly conserved with essentially no insertions and deletions.
Homology modeling can be performed using, for example, the computer programs SWISS-MODEL available through Glaxo Wellcome Experimental Research in Geneva, Switzerland; WHATIF available on EMBL servers; Schnare et. al., J. Mol. Biol, 256: 701-719 (1996); Blundell et. al., Nature 326: 347-352 (1987); Fetrow and Bryant, Bio/Technology 11:479-484 (1993); Greer, Methods in Enzymology 202: 239-252 (1991); and Johnson et al, Crit. Rev. Biochem. Mol. Biol. 29:1-68 (1994). An example of homology modeling can be found, for example, in Szklarz G. D., Life Sci. 61: 2507-2520 (1997). These references are incorporated herein by reference.
Thus, in accordance with the present invention, data capable of generating the three-dimensional structure or three-dimensional representation of the above molecules or molecular complexes, or binding pockets thereof, can be stored in a machine-readable storage medium, which is capable of displaying structural information or a graphical three-dimensional representation of the structure. In one embodiment, means of generating three-dimensional information is provided by means for generating a three-dimensional structural representation of the binding pocket or protein or protein complex.
Rational Drug Design
The FLT3 structure coordinates or the three-dimensional graphical representation generated from these coordinates may be used in conjunction with a computer for a variety of purposes, including 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 FLT3 may inhibit or activate FLT3 or its homologues, 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.
In one embodiment, the invention provides a method of using a computer for selecting an orientation of a chemical entity that interacts favorably with a binding pocket or domain comprising the steps of:

    
    
        (a) providing the structure coordinates of said binding pocket or domain on a computer comprising means for generating three-dimensional structural information from said structure coordinates;
        (b) employing computational means to dock a first chemical entity in the binding pocket or domain;
        (c) quantifying the association between said chemical entity and all or part of the binding pocket or domain for different orientations of the chemical entity; and
        (d) selecting the orientation of the chemical entity with the most favorable interaction based on said quantified association.
    
    


In one embodiment, the docking is facilitated by said quantified association.
In one embodiment, the above method further comprises the following steps before step (a):

    
    
        (e) producing a crystal of a molecule or molecular complex comprising an FLT3 cytoplasmic domain or homologue thereof;
        (f) determining the three-dimensional structure coordinates of the molecule or molecular complex by X-ray diffraction of the crystal; and
        (g) identifying all or part of a binding pocket that corresponds to said binding pocket.
    
    


Three-dimensional structural information in step (a) may be generated by instructions such as a computer program or commands that can generate a three-dimensional representation; subtract distances between atoms; calculate chemical energies for an FLT3 molecule, molecular complex or homologues thereof; or calculate or minimize the chemical energies of an association of FLT3 molecule, molecular complex or homologues thereof to a chemical entity. These types of computer programs are known in the art. The graphical representation can be generated or displayed by commercially available software programs. Examples of software programs include but are not limited to QUANTA (Molecular Simulations, Inc., San Diego, Calif. ©1998, 2000; Accelrys ©2001, 2002), O (Jones et. al., Acta Crystallogr. A47: 110-119 (1991)) and RIBBONS (Carson, J. Appl. Crystallogr., 24: 9589-961 (1991)), which are incorporated herein by reference. Certain software programs may imbue this representation with physico-chemical attributes which are known from the chemical composition of the molecule, such as residue charge, hydrophobicity, torsional and rotational degrees of freedom for the residue or segment, etc. Examples of software programs for calculating chemical energies are described below.
The above method of paragraphs 160 and 161 may further comprise the following step after step (d): outputting said quantified association to a suitable output hardware, such as a CRT display terminal, a CD or DVD recorder, ZIPTM or JAZTM drive, a disk drive, or other machine-readable data storage device, as described previously. The method may further comprise generating a three-dimensional structure, graphical representation thereof, or both, of the protein, binding pocket, molecule or molecular complex prior to step (b).
One embodiment of this invention provides the above method, wherein energy minimization, molecular dynamics simulations, or rigid body minimizations are performed simultaneously with or following step (b).
The above method may further comprise the steps of:

    
    
        (e) repeating steps (b) through (d) with a second chemical entity; and
        (f) selecting at least one of said first or second chemical entity that interacts more favorably with said binding pocket or domain based on said quantified association of said first or second chemical entity.
    
    


In another embodiment, the invention provides the method of using a computer for selecting an orientation of a chemical entity with a favorable shape complementarity in a binding pocket comprising the steps of:

    
    
        (a) providing the structure coordinates of said binding pocket and all or part of the JM-B binding motif bound therein on a computer comprising means for generating three-dimensional structural information from said structure coordinates;
        (b) employing computational means to dock a first chemical entity in the binding pocket;
        (c) quantitating the contact score of said chemical entity in different orientations in the binding pocket; and
        (d) selecting an orientation with the highest contact score.
    
    


In one embodiment, the docking is monitored and directed or facilitated by the contact score.
The method above may further comprise the step of generating a three-dimensional graphical representation of the binding pocket and all or part of the JM-B binding motif bound therein prior to step (b).
The method above may further comprise the steps of:

    
    
        (e) repeating steps (b) through (d) with a second chemical entity; and
        (f) selecting at least one of said first or second chemical entity that has a higher contact score based on said quantitated contact score of said first or second chemical entity.
    
    


In another embodiment, the invention provides a method for screening a plurality of chemical entities to associate at a deformation energy of binding of no greater than 7 kcal/mol with said binding pocket:

    
    
        (a) employing computational means, which utilize said structure coordinates to dock one of said chemical entities from the plurality of chemical entities and said binding pocket;
        (b) quantifying the deformation energy of binding between the chemical entity and the binding pocket;
        (c) repeating steps (a) and (b) for each remaining chemical entity; and
        (d) outputting a set of chemical entities that associate with the binding pocket at a deformation energy of binding of not greater than 7 kcal/mol to a suitable output hardware.
    
    


In another embodiment, the method comprises the steps of:

    
    
        (a) constructing a computer model of a binding pocket of a molecule or molecular complex;
        (b) selecting a chemical entity to be evaluated by a method selected from the group consisting of assembling said chemical entity; selecting a chemical entity from a small molecule database; de novo ligand design of said chemical entity; and modifying a known agonist or inhibitor, or a portion thereof, of an FLT3 protein, or homologue thereof to produce said chemical entity;
        (c) employing computational means to dock said chemical entity to be evaluated in said binding pocket in order to provide an energy-minimized configuration of said chemical entity in the binding pocket; and
        (d) evaluating the results of said docking to quantify the association between said chemical entity and the binding pocket.
    
    


Alternatively, the structure coordinates of the FLT3 binding pockets may be utilized in a method for identifying a candidate inhibitor of a molecule or molecular complex comprising a binding pocket of FLT3. This method comprises the steps of:

    
    
        (a) using a three-dimensional structure of the binding pocket or domain of FLT3 to design, select or optimize a plurality of chemical entities;
        (b) contacting each chemical entity with the molecule and molecular complex;
        (c) monitoring the inhibition to the catalytic activity of the molecule or molecular complex by the chemical entity; and
        (d) selecting a chemical entity based on the effect of the chemical entity on the activity of the molecule or molecular complex.
    
    


In one embodiment, step (a) is carried out using a three-dimensional structure of the binding pocket or domain or portion thereof of the molecule or molecular complex. In another embodiment, the three-dimensional structure is displayed as a graphical representation.
In another embodiment, the method comprises the steps of:

    
    
        (a) constructing a computer model of a binding pocket of the molecule or molecular complex;
        (b) selecting a chemical entity to be evaluated by a method selected from the group consisting of assembling said chemical entity; selecting a chemical entity from a small molecule database; de novo ligand design of said chemical entity; and modifying a known agonist or inhibitor, or a portion thereof, of an FLT3 protein or homologue thereof to produce said chemical entity;
        (c) employing computational means to dock said chemical entity to be evaluated and said binding pocket in order to provide an energy-minimized configuration of said chemical entity in the binding pocket; and
        (d) evaluating the results of said docking to quantify the association between said chemical entity and the binding pocket;
        (e) synthesizing said chemical entity; and
        (f) contacting said chemical entity with said molecule or molecular complex to determine the ability of said chemical entity to activate or inhibit said molecule.
    
    


In one embodiment, the invention provides a method of designing a compound or complex that associates with all or part of the binding pocket of a cytoplasmic domain of an FLT3 protein comprising the steps of:

    
    
        (a) providing the structure coordinates of said binding pocket or domain on a computer comprising means for generating three-dimensional structural information from said structure coordinates;
        (b) using the computer to dock a first chemical entity in part of the binding pocket or domain;
        (c) docking a second chemical entity in another part of the binding pocket or domain;
        (d) quantifying the association between the first and second chemical entity and part of the binding pocket or domain;
        (e) repeating steps (b) to (d) with another first and second chemical entity and selecting a first and a second chemical entity based on said quantified association of all of said first and second chemical entity;
        (f) optionally, visually inspecting the relationship of the first and second chemical entity to each other in relation to the binding pocket or domain on a computer screen using the three-dimensional graphical representation of the binding pocket or domain and said first and second chemical entity; and
        (g) assembling the first and second chemical entity into a compound or complex that interacts with said binding pocket by model building.
    
    


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 FLT3 or FLT3-like binding pockets and domains.
Applicants' elucidation of binding pockets of FLT3 provides the necessary information for designing new chemical entities and compounds that may interact with FLT3 substrate, active site, JM-B binding pockets or FLT3-like substrate, active site or JM-B binding pockets, in whole or in part.
Throughout this section, discussions about the ability of a chemical entity to bind to, interact with or inhibit FLT3 binding pockets refer to features of the entity alone.
The design of compounds that bind to or inhibit FLT3 binding pockets according to this invention generally involves consideration of two factors. First, the chemical entity must be capable of physically and structurally associating with parts or all of the FLT3 binding pockets. Non-covalent molecular interactions important in this association include hydrogen bonding, van der Waals interactions, hydrophobic interactions and electrostatic interactions.
Second, the chemical entity must be able to assume a conformation that allows it to associate with the FLT3 binding pockets directly. Although certain portions of the chemical entity will not directly participate in these associations, those portions of the chemical 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 a chemical entity comprising several chemical entities that directly interact with the FLT3 or FLT3-like binding pockets.
The potential inhibitory or binding effect of a chemical entity on FLT3 binding pockets 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 FLT3 binding pockets, 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 an FLT3 binding pocket. This may be achieved by testing the ability of the molecule to inhibit FLT3 using the assays described in Chan, P. M., et. al., Mol. Cell. Biol. 23: 3067-3078 (2003), which is incorporated herewith by reference.
A potential inhibitor of an FLT3 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 FLT3 binding pockets.
One skilled in the art may use one of several methods to screen chemical entities or fragments or moieties thereof for their ability to associate with the binding pockets described herein. This process may begin by visual inspection of, for example, any of the binding pockets on the computer screen based on the FLT3 structure coordinates FIG. 1A, or other coordinates which define a similar shape generated from the machine-readable storage medium. Selected chemical entities, or fragments or moieties thereof 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 (Accelrys ©2001, 2002) and Sybyl (Tripos Associates, St. Louis, Mo.), followed by, or performed simultaneously with, energy minimization, rigid-body minimization (Gshwend, supra) 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 (Goodford, P. J., “A Computational Procedure for Determining Energetically Favorable Binding Sites on Biologically Important Macromolecules”, J. Med. Chem., 28: 849-857 (1985)). GRID is available from Oxford University, Oxford, UK.
2. MCSS (Miranker, A., et. al., “Functionality Maps of Binding Sites: A Multiple Copy Simultaneous Search Method.” Proteins Struct. Funct. Genet., 11: 29-34 (1991)). MCSS is available from Molecular Simulations, San Diego, Calif.
3. AUTODOCK (Goodsell, D. S., et. al., “Automated Docking of Substrates to Proteins by Simulated Annealing”, Proteins Struct., Funct., and Genet., 8: 195-202 (1990)). AUTODOCK is available from Scripps Research Institute, La Jolla, Calif.
4. DOCK (Kuntz, I. D., et. al., “A Geometric Approach to Macromolecule-Ligand Interactions”, J. Mol. Biol., 161: 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 FLT3. This would be followed by manual model building using software such as QUANTA (Accelrys ©2001, 2002) 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 (Bartlett, P. A., et. al., “CAVEAT: A Program to Facilitate the Structure-Derived Design of Biologically Active Molecules”, in Molecular Recognition in Chemical and Biological Problems, S. M. Roberts, Ed., Royal Society of Chemistry, Special Publication No. 78: pp. 182-196 (1989); Lauri, G. and Bartlett, P. A., “CAVEAT: A Program to Facilitate the Design of Organic Molecules”, J. Comp. Aid. Molec. Design, 8: 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 Martin, Y. C., “3D Database Searching in Drug Design”, J. Med. Chem., 35: 2145-2154 (1992).
3. HOOK (Eisen, M. B., 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: 199-221 (1994)). HOOK is available from Molecular Simulations, San Diego, Calif.
Instead of proceeding to build an inhibitor of an FLT3 binding pocket in a step-wise fashion one fragment or chemical entity at a time as described above, inhibitory or other FLT3 binding compounds may be designed as a whole or “de novo” using either an empty binding pocket or optionally including some portion(s) of a known inhibitor(s). There are many de novo ligand design methods including:
1. LUDI (Böhm, H.-J., “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 (Nishibata, Y., et. al., Tetrahedron, 47: 8985-8990 (1991)). LEGEND is available from Molecular Simulations Incorporated, San Diego, Calif.
3. LeapFrog (available from Tripos Associates, St. Louis, Mo.).
4. SPROUT (Gillet, V., et. al., “SPROUT: A Program for Structure Generation)”, J. Comp. Aid. Molec. Design, 7: 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., Cohen, N. C., et. al., “Molecular Modeling Software and Methods for Medicinal Chemistry, J. Med. Chem., 33: 883-894 (1990); see also, Navia, M. A. and Murcko, M. A., “The Use of Structural Information in Drug Design”, Current Opinions in Structural Biology, 2: 202-210 (1992); Balbes, L. M., et. al., “A Perspective of Modern Methods in Computer-Aided Drug Design”, in Reviews in Computational Chemistry, K. B. Lipkowitz and D. B. Boyd, Eds., VCH Publishers, New York, 5: pp. 337-379 (1994); see also, Guida, W. C., “Software For Structure-Based Drug Design”, Curr. Opin. Struct. Biology, 4: 777-781 (1994)).
Once a chemical entity has been designed or selected by the above methods, the efficiency with which that entity may bind to any of the above binding pockets may be tested and optimized by computational evaluation. For example, an effective 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 binding pocket inhibitors should preferably be designed with a magnitude of deformation energy of binding of not greater than about 10 kcal/mole, more preferably, not greater than 7 kcal/mole. 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.
A chemical entity designed or selected as binding to any one of the above binding pockets 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 (Accelrys ©2001, 2002); Insight II/Discover (Molecular Simulations, Inc., San Diego, Calif. ©1998); DelPhi (Molecular Simulations, Inc., San Diego, Calif. ©1998); 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 any of the above binding pocket. In this screening, the quality of fit of such entities to the binding pocket may be judged either by shape complementarity or by estimated interaction energy (Meng, E. C., et. al., J. Comp. Chem., 13: 505-524 (1992)).
According to another embodiment, the invention provides chemical entities which associate with an FLT3 binding pocket produced or identified by the method set forth above.
Another particularly useful drug design technique enabled by this invention is iterative drug design. Iterative drug design is a method for optimizing associations between a protein and a chemical entity by determining and evaluating the three-dimensional structures of successive sets of protein/chemical entity complexes.
In iterative drug design, crystals of a series of protein or protein complexes are obtained and then the three-dimensional structures of each crystal is solved. Such an approach provides insight into the association between the proteins and compounds of each complex. This is accomplished by selecting compounds with inhibitory activity, obtaining crystals of this new protein/compound complex, solving the three-dimensional structure of the complex, and comparing the associations between the new protein/compound complex and previously solved protein/compound complexes. By observing how changes in the compound affected the protein/compound associations, these associations may be optimized.
In some cases, iterative drug design is carried out by forming successive protein-compound complexes and then crystallizing each new complex. High throughput crystallization assays may be used to find a new crystallization condition or to optimize the original protein crystallization condition for the new complex. Alternatively, a pre-formed protein crystal may be soaked in the presence of an inhibitor, thereby forming a protein/compound complex and obviating the need to crystallize each individual protein/compound complex.
Any of the above methods may be used to design peptide or small molecule mimics of the JM-B binding motif which may have inhibitory effects on full-length FLT3 protein or fragments thereof, or on full-length FLT3 protein which is mutated in or fragments of the mutated protein thereof.
In one embodiment, the present invention provides a method for identifying a candidate inhibitor that interacts with a binding site of an FMS-like tyrosine kinase protein or a homologue thereof, comprising the steps of:

    
    
        (a) obtaining a crystal comprising a cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof, wherein the crystal is characterized with space group P43212 and has unit cell parameters of a=b=80.67 Å, c=150.16 Å;
        (b) obtaining the structure coordinates of amino acids of the crystal of step (a), wherein the structure coordinates are set forth in FIG. 1A-1 to 1A-50;
        (c) generating a three-dimensional model of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof using the structure coordinates of the amino acids generated in step (b), a root mean square deviation from backbone atoms of said amino acids of not more than ±2.0 Å;
        (d) determining a binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof from said three-dimensional model; and
        (e) performing computer fitting analysis to identify the candidate inhibitor which interacts with said binding site.
    
    


In one embodiment, the present invention provides the method for identifying a candidate inhibitor that interacts with a binding site of an FMS-like tyrosine kinase protein or a homologue thereof, further comprising the step of: (f) contacting the identified candidate inhibitor with the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof in order to determine the effect of the inhibitor on FMS-like tyrosine kinase protein activity.
In one embodiment, the present invention provides the method for identifying a candidate inhibitor that interacts with a binding site of an FMS-like tyrosine kinase protein or a homologue thereof, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined in step (d) comprises the structure coordinates according to FIGS. 1A-1 to 1A-50 of amino acid residues H809, R810, and D811, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
In one embodiment, the present invention provides the method for identifying a candidate inhibitor that interacts with a binding site of an FMS-like tyrosine kinase protein or a homologue thereof, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined in step (d) comprises the structure coordinates according to FIGS. 1A-1 to 1A-50 of amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
In one embodiment, the present invention provides the method for identifying a candidate inhibitor that interacts with a binding site of an FMS-like tyrosine kinase protein or a homologue thereof, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined in step (d) comprises the structure coordinates according to FIGS. 1A-1 to 1A-50 of amino acid residues F621, K644, A657, E661, M664, L802, S806, C807, V808, H809, R810, D811, D829, and L832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
In one embodiment, the present invention provides a method for identifying a candidate inhibitor that interacts with a binding site of a cytoplasmic domain of an FMS-like tyrosine kinase protein or a homologue thereof, comprising the steps of:

    
    
        (a) obtaining a crystal comprising the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof, wherein the crystal is characterized with space group P43212 and has unit cell parameters of a=b=80.67 Å, c=150.16 Å;
        (b) obtaining the structure coordinates of amino acids of the crystal of step (a);
        (c) generating a three-dimensional model of said FMS-like tyrosine kinase protein or said homologue thereof using the structure coordinates of the amino acids generated in step (b), a root mean square deviation from backbone atoms of said amino acids of not more than ±2.0 Å;
        (d) determining a binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof from said three-dimensional model; and
        (e) performing computer fitting analysis to identify the candidate inhibitor which interacts with said binding site. In one embodiment, the step of obtaining a crystal is optional.
    
    


In one embodiment, the present invention provides the method for identifying a candidate inhibitor that interacts with a binding site, further comprising the step of:

    
    
        (f) contacting the identified candidate inhibitor with the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof in order to determine the effect of the inhibitor on FMS-like tyrosine kinase protein activity.
    
    


One embodiment of this invention provides the method for identifying a candidate inhibitor that interacts with a binding site, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined in step (d) comprises the structure coordinates according to FIGS. 1A-1 to 1A-50 of amino acid residues H809, R810, and D811, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
One embodiment of this invention provides the method for identifying a candidate inhibitor that interacts with a binding site, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined in step (d) comprises the structure coordinates according to FIGS. 1A-1 to 1A-50 of amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
In one embodiment, the present invention provides the method for identifying a candidate inhibitor that interacts with a binding site, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined in step (d) comprises the structure coordinates according to FIGS. 1A-1 to 1A-50 of amino acid residues F621, K644, A657, E661, M664, L802, S806, C807, V808, H809, R810, D811, D829, and L832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
In one embodiment, the present invention provides a method for identifying a candidate inhibitor that interacts with a binding site of a cytoplasmic domain of an FMS-like tyrosine kinase protein or a homologue thereof, comprising the step of determining a binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or the homologue thereof from a three-dimensional model to design or identify the candidate inhibitor which interacts with said binding site.
In one embodiment, the present invention provides the method for identifying a candidate inhibitor that interacts with a binding site of a cytoplasmic domain of an FMS-like tyrosine kinase protein or a homologue thereof, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined comprises the structure coordinates according to FIGS. 1A-1 to 1A-50 of amino acid residues H809, R810, and D811, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
In one embodiment, the present invention provides the method for identifying a candidate inhibitor that interacts with a binding site of a cytoplasmic domain of an FMS-like tyrosine kinase protein or a homologue thereof, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined comprises the structure coordinates according to FIGS. 1A-1 to 1A-50 of amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
In one embodiment, the present invention provides the method for identifying a candidate inhibitor that interacts with a binding site of a cytoplasmic domain of an FMS-like tyrosine kinase protein or a homologue thereof, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase protein or said homologue thereof determined comprises the structure coordinates according to FIGS. 1A-1 to 1A-50 of amino acid residues F621, K644, A657, L658, E661, M664, L802, K805, S806, C807, V808, H809, R810, D811, C828, D829, F830, G831, and L832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.
One embodiment of this invention provides a method for identifying a candidate inhibitor of a molecule or molecular complex comprising a binding pocket or domain selected from the group consisting of:

    
    
        (i) a set of amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues H809, R810, and D811 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between said set of amino acid residues and said FMS-like tyrosine kinase amino acid residues is not greater than about 2.0 Å;
        (ii) a set of amino acid residues comprising at least five amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, E661, M664, L802, V808, R810, D829, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between said at least five amino acid residues and said FMS-like tyrosine kinase amino acid residues which are identical is not greater than about 2.0 Å;
        (iii) a set of amino acid residues comprising at least six amino acid residues which are identical to human FMS-like tyrosine kinase amino acid residues F621, K644, A657, L658, E661, M664, L802, K805, S806, C807, V808, H809, R810, D811, C828, D829, F830, G831, and L832 according to FIG. 1A, wherein the root mean square deviation of the backbone atoms between said at least six amino acid residues and said FMS-like tyrosine kinase amino acid residues which are identical is not greater than about 2.0 Å; and
        (iv) a set of amino acid residues that are identical to FMS-like tyrosine kinase amino acid residues according to FIG. 1A, wherein the root mean square deviation between said set of amino acid residues and said FMS-like tyrosine kinase amino acid residues is not more than about 3.0 Å;
    
    


comprising the steps of:

    
    
        (a) using a three-dimensional structure of the binding pocket or domain to design, select or optimize a plurality of chemical entities; and
        (b) selecting said candidate inhibitor based on the inhibitory effect of said chemical entities on said cytoplasmic domain of said FMS-like tyrosine kinase protein or said cytoplasmic domain of said FMS-like tyrosine kinase protein homologue on the catalytic activity of the molecule.
    
    


In one embodiment, the present invention provides a method of using a crystal of a cytoplasmic domain of said FMS-like tyrosine kinase protein or a homologue in an inhibitor screening assay comprising:

    
    
        (a) selecting a potential inhibitor by performing rational drug design with a three-dimensional structure determined for the crystal, wherein said selecting is performed in conjunction with computer modeling;
        (b) contacting the potential inhibitor with a kinase; and
        (c) detecting the ability of the potential inhibitor for inhibiting the kinase.
    
    


In certain embodiments, the ability of the potential inhibitor for inhibiting the kinase is assesed using an enzyme inhibition assay. In other embodiments, the ability of the potential inhibitor for inhibiting the kinase is performed using a cellular-based assay.
In one embodiment, the present invention provides a method for identifying a candidate inhibitor that interacts with a binding site of an FMS-like tyrosine kinase protein or a homologue thereof comprising:

    
    
        (a) obtaining a crystal of an FMS-like tyrosine kinase protein or a homologue thereof;
        (b) obtaining the atomic coordinates of the crystal; and
        (c) using the atomic coordinates and one or more molecular modeling techniques to identify the candidate inhibitor that interacts with a binding site of an FMS-like tyrosine kinase protein or a homologue thereof. In certain embodiments, the crystal comprising a cytoplasmic domain of an FMS-like tyrosine kinase protein or a homologue thereof. In one embodiment, the step of obtaining a crystal is optional.
    
    


In one embodiment, the present invention provides the method for identifying a candidate inhibitor that interacts with a binding site of an FMS-like tyrosine kinase protein or a homologue thereof, wherein the one or more molecular modeling techniques are selected from the group consisting of graphic molecular modeling and computational chemistry.
In one embodiment, the present invention provides the method for identifying a candidate inhibitor that interacts with a binding site of an FMS-like tyrosine kinase protein or a homologue thereof, further comprising the candidate inhibitor with the FMS-like tyrosine kinase protein or the homologue and detecting binding of the candidate inhibitor to the FMS-like tyrosine kinase protein or the homologue.
In one embodiment, the present invention provides a method of struture-based identification of candidate compounds for binding to an FMS-like tyrosine kinase protein or a homologue thereof, comprising:

    
    
        (a) constructing a three-dimensional structure of the FMS-like tyrosine kinase protein or a homologue thereof;
        (b) performing computer-assisted structure-based drug design with said structure of the FMS-like tyrosine kinase protein or a homologue; and
        (c) identifying at least one candidate inhibitor that is predicted to have a compatible conformation with a binding site of the structure of the FMS-like tyrosine kinase protein or a homologue.
    
    


In certain embodiments, the present invention provides for methods wherein the three-dimensional structure is visualized as a computer image generated when said atomic coordinates determined by X-ray diffraction are analyzed on a computer using a graphical display software program to create an electronic file of the image and visualizing the electronic file on a computer capable of representing the electronic file as a three-dimensional image.
Structure Determination of Other Molecules
The structure coordinates set forth in FIG. 1A can also be used in obtaining structural information about other crystallized molecules or molecular complexes. This may be achieved by any of a number of well-known techniques, including molecular replacement.
According to one 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 at least a portion of the structure coordinates set forth in FIG. 1A or homology model thereof, 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.
In another embodiment, the invention provides a computer for determining at least a portion of the structure coordinates corresponding to X-ray diffraction data obtained from a molecule or molecular complex having an unknown structure, wherein said computer comprises:

    
    
        (a) a machine-readable data storage medium comprising a data storage material encoded with machine-readable data, wherein said data comprises at least a portion of the structure coordinates of FLT3 according to FIG. 1A or a homology model thereof;
        (b) a machine-readable data storage medium comprising a data storage material encoded with machine-readable data, wherein said data comprises X-ray diffraction data obtained from said molecule or molecular complex having an unknown structure; and
        (c) instructions for performing a Fourier transform of the machine-readable data of (a) and for processing said machine-readable data of (b) into structure coordinates.
    
    


For example, the Fourier transform of at least a portion of the structure coordinates set forth in FIG. 1A or homology model thereof may be used to determine at least a portion of the structure coordinates of the molecule or molecular complex.
Therefore, another embodiment this invention provides a method of utilizing molecular replacement to obtain structural information about a molecule or a molecular complex of unknown structure wherein the molecule or molecular complex is sufficiently homologous to FLT3, 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;
        (c) applying at least a portion of the FLT3 structure coordinates set forth in one of FIG. 1A or a homology model thereof to the X-ray diffraction pattern to generate a three-dimensional electron density map of at least a portion of the molecule or molecular complex whose structure is unknown; and
        (d) generating a structural model of the molecule or molecular complex from the three-dimensional electron density map.
    
    


In one embodiment, the method is performed using a computer. In another embodiment, the molecule is selected from the group consisting of FLT3 protein and FLT3 cytoplasmic domain homologues. In another embodiment, the molecular complex is FLT3 cytoplasmic domain complex or homologue thereof.
By using molecular replacement, all or part of the structure coordinates of FLT3 as provided by this invention (and set forth in FIG. 1A) 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 may 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 FLT3 protein according to FIG. 1A 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: 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 structure of human FLT3 protein can be resolved by this method.
In one embodiment, the method of molecular replacement is utilized to obtain structural information about an FLT3 homologue. The structure coordinates of FLT3 as provided by this invention are particularly useful in solving the structure of FLT3 complexes that are bound by ligands, substrates and inhibitors.
Furthermore, the structure coordinates of FLT3 as provided by this invention are useful in solving the structure of FLT3 proteins that have amino acid substitutions, additions and/or deletions (referred to collectively as “FLT3 mutants”, as compared to naturally occurring FLT3). These FLT3 mutants 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 FLT3. Potential sites for modification within the various binding pockets 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 FLT3 and a chemical entity or compound.
The structure coordinates are also particularly useful in solving the structure of crystals of the cytoplasmic domain of FLT3 or 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 FLT3 inhibitors. 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 FLT3 inhibition activity.
All of the molecules and complexes referred to above may be studied using well-known X-ray diffraction techniques and may be refined using 1.5-3.4 Å resolution X-ray data to an R value of about 0.30 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)) or CNS (Brunger et. al., Acta Cryst., D54: 905-921, (1998)).
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
FLT3 Expression and Purification
The cytoplasmic domain (H564-S993) of FLT3 (full length sequence, amino acid residues 1 to 993; GenBank accession no. NP—004110; SEQ ID NO:1) was cloned from cDNA bone marrow library (Clonetech, CA, USA). The expressed protein (H564-S958) was engineered with an internal deletion, removing the residues spanning H711-H761, and containing an N-terminal hexa-histidine tagged protein with a thrombin cleavage site. Cloned into the baculoviral transfer vector pBEV10, for insect cell expression, the recombinant virus generated was plaque purified and amplified to obtain a high-titer clonal viral stock. For production, High-5 insect cells grown to 2×106 cell/ml in Excell-405 medium (JRH Bioscience, KS, US) were infected with virus at a multiplicity of infection (MOI) of 2.5 and incubated 72-96 h at 27° C.
Frozen cell paste was thawed in 8 volumes of Buffer A (50 mM Hepes pH 8.2, 500 mM NaCl, 20% (v/v) glycerol, 0.1% (v/v) Triton X-100, 3 mM β-mercaptoethanol, 5 mM imidazole, 50 μM diisopropyl fluorophosphate, 2 μg/ml Pepstatin, 1 μg/ml leupeptin, 1 μg/ml E64 (Roche Diagnostics Corp., Indianapolis, Ind.), and 25 μl/L Benzonase) and mechanically lysed in a microfluidizer (Microfluidics, Newton, Mass.). The lysate was centrifuged at 54,000×g for 1 h, and the supernatant incubated with metal affinity resin (Ni-NTA agarose, Qiagen Inc. Valencia, Calif.) overnight at 4° C. After extensive washing with 20 column volumes of Buffer A, the protein was eluted with Buffer A containing 200 mM imidazole with the pH readjusted to 8.2. The elution pool was concentrated by ultrafiltration (30 KDa MWCO) in an Amicon stirred cell concentrator (Millipore, Billerica, Mass.) and loaded onto a HR 16/60 Superdex-200 size-exclusion column (Amersham Biosciences, Uppsala, Sweden) equilibrated in Buffer B (50 mM Hepes 8.2, 200 mM NaCl, 15% (v/v) glycerol, 5 mM DTT, and 0.1% (w/v) β-octylglucopyranoside). The FLT3 protein was pooled based on SDS-PAGE analysis. The protein was diluted two-fold with buffer containing 50 mM Hepes pH 8.2, 20% (v/v) glycerol, 5 mM DTT to reduce the salt concentration to 100 mM NaCl and loaded onto a MonoQ HR (5/5) column (Amersham Biosciences Corp, Uppsala, Sweden) pre-equilibrated in Buffer C (50 mM Hepes 8.2, 100 mM NaCl, 20% (v/v) glycerol, 5 mM DTT). A gradient was developed from Buffer C to Buffer C containing 1 M NaCl over 60 column volumes and FLT3 was eluted at  ̃250 mM NaCl. The protein was pooled based on SDS-PAGE analysis and the buffer exchanged to 25 mM Tris pH 8.1, 500 mM NaCl, 5% (v/v) glycerol and 5 mM DTT using a G25 column (Amersham Biosciences, Uppsala, Sweden). The protein containing fractions were pooled and concentrated to 10 mg/ml using a 30 KDa MWCO Vivaspin concentrator (Vivascience, Hanover, Germany) in the presence of 2 mM AMP-PNP and 4 mM MgCl2.
EXAMPLE 2
Protein Crystallization for Native FLT3
Samples of the protein from Example 1 were subjected to ultracentrifugation at 90,000×g for 10 minutes prior to crystallization.
Crystals were grown in hanging drops by the vapor diffusion method where 0.5 μl each of FLT3 sample and reservoir solutions were combined. The reservoir condition was 1.2 M sodium phosphate, 1.2 M potassium phosphate, 0.1 M CAPS (3-cyclohexylamino-1-propanesulfonic acid) (pH 10.5) and 0.2 M lithium sulfate. Typically, the crystals took 1-3 weeks to appear and grew as hexagonal bipyramids to full size 2-4 days later. The morphology is representative of the space group P43212 with a unit cell of dimensions a=b=80.67 Å, c=150.16 Å containing one molecule per asymmetric unit.
EXAMPLE 3
X-ray Diffraction and Structure Determination of FLT3
Prior to data collection the crystals were transferred stepwise to reservoir solutions containing 30% (v/v) glycerol then flash frozen in liquid nitrogen and stored in liquid nitrogen. Diffraction data were collected at −180° C. on a Rigaku R-AXIS IV++ imaging system mounted on a Rigaku RU-H3R rotating anode x-ray generator (CuKα) operated at 50 kV and 100 mA. Diffraction intensities were integrated and scaled with CrystalClear (Rigaku Intl. Corp). Table 1 summarizes information about the data collection.
The structure of FLT3 was solved by the molecular replacement method using AMoRe (Navaza, J., Acta Cryst. A50: 157-163 (1994)) in the CCP4 suite of programs (Collaborative Computational Project, Number 4, “The CCP4 Suite: Programs for Protein Crystallography,” Acta Cryst. D50: 760-763 (1994)). The starting model for FLT3 was derived from the human vascular endothelial growth factor receptor 2 kinase domain (KDR) atomic coordinates (PDB entry 1VR2) where KDR residues identical to those in FLT3 were included in the search model. For those residues that were not identical, the side chains were pared back to the beta carbon atom.
Both the rotation and translation searches yielded a single consistent solution. The model corresponding to this solution was subjected to rigid body refinement, followed by torsional dynamics using CNX (Brünger, A. T., et al, Acta Cryst. D54: 905-921 (1998); (Accelrys, ©2000, 2001). All model building was performed using QUANTA (Molecular Simulations, Inc., San Diego, Calif. ©1998, 2000; Accelrys ©2001, 2002).
The electron density corresponding to side chains absent from the search model was generally clear and unambiguous in the kinase domain. Although, the electron density for certain regions of the activation loop and JM domain was visible, the entire course of each segment could not be traced in the preliminary electron density maps. The model was subjected to cycles of building/refinement. Composite difference electron maps were calculated which allowed the main chains of both the active site loop and the JM domain to be traced.
The final FLT3 structure contains the JM domain (residues 570 to 603), the kinase fold and activation loop (residues 604 to 710, 783 to 947), and 329 water molecules, one phosphate ion and one CAPS buffer molecule. During the course of the refinement, the electron density corresponding to residues 649-654 was poor and did not improve. Consequently, these residues that reside on the loop connecting β3 and αC were removed from the final model. Crystallographic refinement statistics are provided in Table 1.
EXAMPLE 4
Overview of FLT3 Structure
The principal features of the autoinhibited FLT3 structure include the bilobal kinase fold, the activation loop and the JM domain (FIGS. 2A and 2B). The kinase fold in FLT3 is that typically found in most protein kinases and consists of N- and C-terminal domains (N- and C-lobes) (Knighton, D. R., et. al., Science 253: 407-414 (1991)). The N-lobe contains a twisted five-stranded anti-parallel beta sheet adjacent to an α-helix.
According to the standard kinase nomenclature (FIG. 6A), the five β-strands in the N-lobe are termed β1-β5 and the sole N-lobe α-helix is termed αC. The secondary structure of the FLT3 C-lobe is predominantly α-helical, containing seven α-helices (αD, αE, αEF, αF-αI) and three β-strands (β6, β7 and β8). The activation loop is attached to the C-lobe and contains two additional β-strands (β10 and β11). These two β-strands comprise a small β-sheet that is characteristic of inactive kinases. This β-sheet is not present in activated kinase structures such as IRK (Hubbard et. al., 1997, supra) and cKIT (Mol et. al., supra) where the activation loop is fully open and the two β-strands are separated. However, in both activated structures, the unfolded activation loop generates β9 that pairs with β6. Even though β9 does not exist in the inactive FLT3 structure, β6 survives by forming a novel β-pair with βJ1 on the JM domain. The JM domain contains an additional β-strand termed βJ2.
Kinase N- and C-lobes are typically connected by a single flexible polypeptide strand that allows considerable rotational movement of the two domains relative to each other. This conformational diversity is observed in the multitude of available kinase structures (Huse, M., and Kuriyan, J. Cell 109: 275-282 (2002)) where the interlobe angular orientation can vary by approximately 20°. When the N-lobe is rotated away from the C-terminal domain, the kinase is in the catalytically ‘inactive’ form. Conversely, if the N-lobe is rotated towards the C-lobe, allowing key catalytic residues from both lobes to align, the kinase adopts the catalytically ‘active’ conformation. The crystal structure of autoinhibited FLT3 conforms to the prototypical conformation common to other inactive kinases that have a ‘closed’ activation loop folded between the two lobes of the ‘inactive’ kinase fold. A novel feature of the FLT3 structure is the presence of the complete JM domain that assumes its autoinhibited conformation and interacts with all key features of FLT3.
Activation Loop
Analogous to other kinases, the activation loop on FLT3 can be characterized as a long flexible peptide segment that connects β8 to αEF. The ends of the loop are delimited by the conserved three residue sequences Asp829-Phe830-Gly831 and Trp854-Met855-Ala856 located at the N- and C-terminal hinges of the loop, respectively. Asp829 is invariant in kinases and serves as the catalytic base in the phosphotransfer reaction (Bossemeyer, D., et. al., EMBO J. 12: 849-859 (1993)).
The crystal structures of FGFR (Mohammadi, M., et. al., supra), activated IRK (IRK-A) (Hubbard, S. R., supra) and inactive IRK (IRK-I) (Hubbard, S. R., et. al. (1994), supra) demonstrate the wide range of conformations available to the activation loop. In IRK-I, the activation loop is completely closed, allowing Tyr1143 to position itself in the active site, whereas the FGFR activation loop is partially unfolded and the IRK-A activation loop is fully open (FIG. 5A). These large conformational changes have been documented for other kinases as they transition from inactive to active states (Huse, M., and Kuriyan, supra). Kinase activation loops usually contain one to three tyrosine residues that can serve as phosphorylation sites. When these tyrosines are unphosphorylated, the activation loop typically assumes the closed conformation by folding into the cleft between the N- and C-lobes, thereby blocking access to the peptide substrate and ATP binding sites. In the activated state, which is associated with tyrosine phosphorylation, the activation loop adopts the open conformation and no longer restricts ATP and protein substrates from binding (Huse, M., and Kuriyan, supra).
The FLT3 activation loop in the crystal structure is similar to that observed in other inactive kinases such as IRK-I. The FLT3 loop contains a single unphosphorylated tyrosine residue (Tyr842) that adopts an orientation identical to its equivalent tyrosine in IRK-I (Tyr1143) and it is located in the same type of environment as its IRK-I analog. In particular, the phenolic ring points directly into the active site and is held in place by a hydrogen bond to an Asp side chain (Asp811) that, in turn, forms an ion pair with an arginine side chain (Arg834). The electron density corresponding to the Tyr842 side chain was one of the dominant features of the activation loop observed in the initial difference electron density maps. Except for the β10/β11 β-sheet, the electron density corresponding to the rest of the activation loop is significantly weaker, particularly for the residue segment running from Ile836 to Arg841. This is consistent with the IRK-I structure where the activation loop contains flexible peptide segments that are connected firmly to the C-lobe at each end and are held in place by the tightly bound tyrosine.
Point mutations of residues located on the activation loop generate constitutively activated forms of FLT3 that are implicated in AML (Yamamoto, Y., et. al., supra; Abu-Duhier, F. M., et. al., supra). The most common mutation is Asp835Tyr, however other mutations such as Asp835Thr and Asp835His have been reported (Yamamoto, Y., et. al., supra). Constitutive activation has also been reported for Ile836 mutations as well (Thiede, C., et. al., supra; Yamamoto, Y., et. al., supra). It has been suggested that Asp835 might provide essential stability to the closed form of the loop that is lost in the FLT3 mutants. Looking at Asp835 in the FLT3 structure, the side chain of Asp835 does indeed make a hydrogen bond with the main chain of Ser838. However, the significance of this, with respect to stability of the closed loop conformation, is not clear because the electron density in this region is relatively diffuse. Likewise, the possible role of Ile836 in stabilizing the closed conformation remains ambiguous due to the disorder and implied mobility of the activation loop in this region. Therefore, the electron density of the structure does not support the suggestion that Asp835 and/or Ile836 provides intrinsic stability to the closed form of the activation loop which, when lost in oncogenic mutants, prevents the closure of the activation loop.
Juxtamembrane Domain
The structure of the JM domain and its relation to the rest of the FLT3 structure is illustrated in FIGS. 2A, 2B, 3A and 3B. The JM domain can be divided into three distinct topological components: the JM binding motif (JM-B), the JM switch motif (JM-S) and the zipper or linker peptide segment (JM-Z) (FIGS. 3A and 3B). The JM-B (Tyr572-Met578) is termed the binding motif because it is nearly buried in the FLT3 structure. It is a short finger-like segment beginning with an invariant tyrosine (Tyr572) that resembles a ‘hook’. Even though the JM-B consists of only seven residues, it makes contacts with virtually every structural component implicated in the activation/inactivation cycle of the FLT3 cytoplasmic domain (FIGS. 3A and 3B). The JM-S (Val579-Val592) is a two-stranded anti-parallel β-twist that is situated in an external location on the C-lobe. It is attached to the JM-B by a continuous extension of a β-strand termed βJ1. The JM-S is defined as a switching motif because it contains two key tyrosine residues whose state of phosphorylation is implicated in the activation and regulation of receptor enzymatic activity (Mol. C. D., et. al., supra). Lastly, the JM-Z (Asp593-Trp603) is located at the C-terminus of the JM domain and is associated primarily with the N-lobe as it loops around the outside of αC forming fairly extensive contacts with this key α-helix.
Looking at the FLT3 JM domain in more detail, particularly at the N-terminal region of the peptide sequence, the segment from His564 to Arg571, now defined as part of the linker connecting the JM and the transmembrane (TM) domains, is disordered and not present in the FLT3 structure. The first residue identified in the electron density map is Tyr572 located at the N-terminus of the JM-B. The side chain of Tyr572 sits in a tight pocket where the aromatic ring makes extensive hydrophobic contacts with surrounding residues. The Tyr572 phenolic hydroxyloxygen forms a hydrogen bond with a side chain oxygen on the highly conserved Glu661 which, together with sequence-invariant Lys644, forms an ion pair which is critical for nucleotide binding. It is interesting to note that Tyr842, which serves as an anchoring point on the activation loop, interacts with another ion pair (Asp811:Arg834) in a similar manner (FIG. 4). Proceeding in the C-terminal direction from Tyr572, the polypeptide chain runs inward three residues before making a perpendicular turn at Gln575. The Gln575 and the next three residues in the polypeptide chain interact with β6 to form an antiparallel β-sheet. The electron density in this segment (Tyr572-Met578) is excellent and indicative of a well-ordered conformation. The JM-B binds in a central strategic area of FLT3 and makes contacts with a number of key structural components (Knighton, et. al., supra) including the glycine-rich loop, the activation loop, and the catalytically important αC.
The JM-B is nearly buried in the crystal structure of autoinhibited FLT3. It loses 350 Å2 of solvent accessible surface area as it attaches to its binding site. This accounts for 76% of the total accessible surface area available to the JM-B when it is not in contact with its neighboring structural elements in the crystal structure. This proportionally large area of contact suggests a stable interaction between the JM-B and the rest of the molecule. The strength of the interaction is likely enhanced by a moderate amount of interdigitation between side chains on the JM-B and those in the surrounding pocket as it is positioned between β6 on one side and αC on the other. Therefore, the JM-B can be thought of as a ‘wedge’ that stabilizes the inactive kinase conformation of FLT3 by preventing the N-lobe from rotating towards the C-lobe to generate the activated kinase fold.
An additional role of the JM-B in stabilizing the inactive kinase conformation is revealed by its proximity to the N-terminal hinge of the activation loop. This is shown in FIGS. 5A and 5C which depicts the relationship of the JM-B to the open and closed forms of the IRK activation loops. The autoinhibited FLT3 activation loop is in the closed conformation and, as long as the JM-B is in place, the activation loop cannot unfold in a manner similar to the fully open IRK-A form (FIG. 5C) or even to the partially open FGFR form (FIG. 5B). In fact, a segment of the IRK activation loop forms a β-strand (β9) that overlaps with the superposed βJ1 and similarly generates a β-pair with β6. This relationship between β9 and β6 exists in the activated structure of cKIT. Therefore, β6 can stabilize both the active and the inactive (autoinhibited) kinase conformation, depending on whether it interacts with the activation loop or JM domain, respectively.
The JM-S lies immediately adjacent to the JM-B and consists of fourteen residues that form a two-stranded β-twist motif. The JM-S starts as a four-residue extension to βJ1 where the transition from JM-B to JM-S occurs at Met578-Val579. The extended conformation of the βJ1 likely provides a rigid connection between the JM-B and JM-S that could directly transmit perturbations and movements from one motif to the other. βJ1 terminates at residue Val581 where a five-residue loop connects to a short four-residue β-strand (βJ2) that contains two conserved tyrosines (Tyr589 and Tyr591). The side chains on these tyrosines are adjacent to each other and lay between the JM-S and the C-lobe.
Unlike the JM-B, which is immersed in extensive intra-molecular contacts, the JM-S interacts with the rest of the molecule in a less extensive, less complementary way. As shown in FIGS. 3A and 4, the JM-S protrudes from FLT3 as an extension of the βJ1. The JM-S appears to be held in place by the rigidity of the β-strand network consisting of β6, βJ1, and βJ1. Any significant change in the orientation of the JM-S relative to the C-lobe would likely disrupt the β-network and de-stabilize binding of the JM-B. A dominant feature of the JM-S is the tyrosine pair (Tyr589 and Tyr591) pointed directly at the C-lobe. In fact, this tyrosine pair is the principal contact between the JM-S β-sheet and the C-lobe. Neither of the tyrosine side chains binds in a tight complementary pocket similar to Tyr572 and Tyr599. As such, the JM-S does not provide a significant amount of binding interface between it and the C-lobe that can directly enhance the attachment of the JM-B. Rather, the spatial arrangement of the JM-S relative to the C-lobe suggests the role of the JM-S is to provide a rigid and properly oriented framework requisite for the interposition of Tyr589 and Tyr591 between it and the C-lobe.
The remaining eleven residues of the JM domain comprise the JM-Z that folds up alongside the N-lobe. It is a simple peptide segment that loops over αC and terminates in the hinge region of the JM domain. Recent crystallographic studies on IRK, cKIT and other RTKs show that the JM-Z region can undergo large amplitude rotations away from the N-lobe by pivoting about its attachment point. This is demonstrated in the crystal structure of cKIT (PDB entry, 1PKG) that contains two independent copies of the activated molecule. The JM-Z on one molecule adopts the ‘attached’ conformation as it lays alongside the N-lobe, while the JM-Z on the other molecule in the asymmetric unit is ‘unattached’ or ‘unzipped’ from the N-lobe. The conformation of the FLT3 JM-Z is very similar to the ‘attached’ conformation observed in activated cKIT and inactivated IRK crystal structures. However, the JM-B and the N-terminal half of the JM-S in these structures are either missing or disordered.
FLT3 Tyr599 is conserved across the PDGFR family and is located near the JM hinge region. Like its cKIT equivalent on the ‘attached’ JM-Z, the tyrosine side chain is buried in a pocket where the hydroxyl group is hydrogen bonded to the equally conserved Glu604. In the ‘unattached’ cKIT JM-Z, the tyrosine is pulled out of its insertion pocket and is exposed to solvent. There is no biochemical or crystallographic evidence to suggest that Tyr599, or its equivalents in cKIT or IRK, are autophosphorylation sites. Thus, the principal function of Tyr599 could be one of providing the proper orientation and guidance to the JM hinge region as the JM domain transitions from the ‘unattached’ to the ‘attached’ conformation.
The JM Domain Autoinhibitory Mechanism
The activation of FLT3 is a stepwise process that is described for all the type III RTKs (Rosnet, O., and Birnbaum, D. Crit. Rev. Oncog. 4: 595-613 (1993)). It is initiated when FL, its specific endogenous ligand, binds to the extracellular domain of FLT3. This promotes dimerization of the receptor and concomitant juxtapositioning of the cytoplasmic domains. Once the dimer is formed, transphosphorylaton of specific tyrosine residues on the JM domain can take place. This activates full kinase activity which induces multiple signaling pathways that are involved in cell proliferation and activation (Turner, A. M., et. al., supra; Heldin, C. H. Cell 80: 213-223 (1995)). The kinase activity is negatively modulated by tyrosine phosphatases that dephosphorylate the tyrosines on the unbound JM domain. This allows the JM domain to adopt its autoinhibitory conformation.
There has been considerable speculation on how the JM domain regulates the activity of type III RTKs. Certainly, the phosphorylation state of key tyrosine residues plays a central role, but the precise nature of regulation has not been understood. However, the structure of autoinhibited FLT3 suggests a simple mechanism for kinase modulation by the JM domain. As shown earlier, the tyrosine pair (Tyr589 and Tyr591) is positioned between the JM-S β-sheet and the C-lobe. As such, it cannot accommodate the charged, bulky phosphate moieties of phosphorylated tyrosines and still maintain the bound conformation of the JM-B/JM-S complex. Without being bound by theory, when one or both of these tyrosines is phosphorylated, the JM-S cannot fold up properly and/or position itself next to the C-lobe in a manner conducive to the autoinhibited state. Conversely, when the phosphates are removed from the tyrosines by regulatory phosphatases, the JM-S can position itself next to the C-lobe and allow the JM-B to insert into its autoinhibitory binding site. This role of the JM domain is supported by Chan et. al. (Chan, P. M., et. al., Mol. Cell. Biol. 23: 3067-3078 (2003)) who synthesized the JM domain of cKIT, approximately 39 amino acid long, and showed that this peptide alone inhibits a cKIT construct that had the JM domain entirely deleted. Furthermore, other JM constructs that were either mutated or phosphorylated did not ‘autoinhibit’ the JM-minus cKIT construct.
According to this suggested mechanism, the role of the JM-Z is to correctly align and maintain the JM-S in the proper register during and after the transition between activated and inactive states of FLT3. Therefore the length of the JM-S is critical and should be conserved across all members of the PDGFR family. This is confirmed by analysis of the aligned sequences of the PDGFR family (FIG. 6B) where the lengths of all structural elements in the JM domain are conserved. The only exception occurs in the loop connecting the two JM-S β-strands. The variability in this loop region is inconsequential because it is located between the two aligned strands of the β-sheet that does not affect the position of the JM-S relative to the JM-B and the JM-Z.
Autoinhibitory Mechanism of the PDGRF Family
The amino acid sequence of the JM domain shows a high degree of homology across the PDGFR family (FIG. 6B). This homology likely correlates with conserved structural and mechanistic features of all the members. Irusta et. al. (Irusta, P. M., et al, J. Biol. Chem. 277, pp 38627-38634 (2002)) carried out a systematic alanine-scanning mutational analysis of the JM domain of murine PDGFRβ and identified a number of residues, namely Y530A, W534A, V536A, 1537A, L555A, Y557A, and the double mutant Y547A/Y549A that resulted in constitutive receptor activation (PDGFRβ numbers are given). These residues in PDGFRβ when mutated to alanine, result in constitutive activation, and therefore, correspond to key residues in the JM domain observed in the FLT3 structure. For example, substitution of alanine for residues corresponding to Tyr572 and Leu576, two key anchoring points for the JM-B, resulted in receptor activation. Conversely, mutation of amino acid residues corresponding to FLT3 amino acid residues that are not involved in interactions with the kinase domain, such as Gln580, Gly583 and Glu588, did not result in activation of PDGFRβ. In another experiment, an alanine was inserted between positions Ser574 and Gln575 and between positions Asp586 and Asn587 in the PDGFRβ sequence. These insertions fall within the JM-B region and within the JM-S loop region, respectively. Not surprisingly, the JM-B insertion resulted in constitutive activation while the loop insertion did not. The results of these mutational experiments, together with own structural analysis, strongly suggests that the JM sequence of the PDGFR family constitutes a common inhibitory domain that utilizes the autoinhibitory mechanism described for FLT3.
The common autoinhibitory mechanism for the PDGFR kinase family is the result of a highly conserved amino acid sequence in the JM region. Other kinases, with different JM amino acid sequences, utilize distinctly different autoregulatory mechanisms. For example, the JM region of Type I TGFβ receptor, in association with FKBP12, forms an autoinhibitory structure that interacts with the kinase domain to block catalytic activity, most significantly by displacing the critical αC helix in the N-lobe of the kinase domain and disrupting the conformation of the ATP-binding site (Huse M., et. al., Cell 96: 425-436 (1999)). In a similar manner, the autoinhibitory mechanism of EphB2 is facilitated by the JM domain as it interacts with the αC helix. Specifically, a helix on the JM domain introduces a kink in αC that stabilizes the inactive conformation (Wybenga-Groot, L. E., et. al., supra). All three of these distinct autoinhibitory mechanisms are examples of the conformational plasticity of protein kinases described by Huse and Kuriyan, supra.
Structural Implications for ITDs Found in AML
The internal tandem duplications (ITDs) found in AML patients (Nakao, M., et. al., supra; Thiede, C., et. al., supra) range from four to about sixty-eight amino acid residues (Kottaridis, P. D., et. al., supra). This insert domain made FLT3 autoinhibition ‘leaky’, as described by Schlessinger et. al. (Irusta, P. M., et. al., supra), allowing it to switch FLT3 from the inactive to the catalytically active conformation without the presence of its endogenous ligand FL. The ITD insertions generally occur in the JM-Z near the JM hinge region which, among other things, offsets the position of the JM-S in the FLT3 structure. This likely disturbs or prevents the optimal orientation of JM-S as it tries to position the JM-B in its binding site. In addition to this, the presence of an ITD could simply obviate the complementary interaction that exists between the normal JM-Z and the N-lobe. As such, the aberrant JM-S could be stuck in the ‘unattached’ conformation permanently.
While we have described a number of embodiments of this invention, it is apparent that our basic constructions may be altered to provide other embodiments which utilize the products, processes and methods of this invention.
All cited documents are incorporated herein by reference.







TABLE 1







Data Collection and Refinement Statistics










Data set
Native







Data collection




X-ray source
Rigaku RU-H3R



Space group
P43212



Unit cell parameters (Å)
a = b = 80.65;




c = 150.13



Resolution (Å)
30 − 2.1



Unique reflections
29690



Redundancy
6.82



Completeness (%)*
99.3 (94.1)



Rmerge*
0.056 (0.265)



<I/σ>*
18.6 (57)



Refinement



Reflections used
29428



Test reflections
1706



R-factor
0.209



free R-factor (% data)
0.246 (5.0)



RMS deviation



Bond lengths (Å)
0.009



Bond angles (°)
1.59



Protein atoms
2439



Solvent atoms
329



Phosphate atoms
5



CAPS atoms
13







*Values for the highest resolution shell are shown in parentheses.







  
    R
    merge
  
  =
  
    
      ∑
      hkl
    
    ⁢
    
      
        ∑
        i
      
      ⁢
      
        
          
          
            
              
                I
                ⁡
                
                  (
                  hkl
                  )
                
              
              i
            
            -
            
              〈
              
                I
                ⁡
                
                  (
                  hkl
                  )
                
              
              〉
            
          
          
        
        /
        
          
            ∑
            hkl
          
          ⁢
          
            
              ∑
              i
            
            ⁢
            
              〈
              
                
                  I
                  ⁡
                  
                    (
                    hkl
                    )
                  
                
                i
              
              〉
            
          
        
      
    
  



over i observations of reflection hkl.







  
    R
    ⁢
    
      -
    
    ⁢
    factor
  
  =
  
    ∑
    
      
        
        
          
            
            
              F
              obs
            
            
          
          -
          
            
            
              F
              calc
            
            
          
        
        
      
      /
      
        ∑
        
          
          
            F
            obs
          
          
        
      
    
  



where Fobs and Fcalc are the observed and calculated structure factors, respectively.



Free R-factor is calculated from a randomly chosen subset of reflections not used for refinement.

## Claims

1. A crystal of a cytoplasmic domain of an FMS-like tyrosine kinase 3 (FLT3) protein wherein said cytoplasmic domain is unphosphorylated and consists of, in the following order: an N-terminal hexa-histidine tag containing a thrombin cleavage site fused to the N-terminus of FLT3 amino acids 564-958 of SEQ ID NO:1, except for an internal deletion of amino acids 711-761, and wherein said crystal is characterized as having space group P43212 with unit cell parameters of a=b=80.67 Å and c=150.16 Å.

2. A crystallizable composition of a cytoplasmic domain of an FMS-like tyrosine kinase 3 (FLT3) protein wherein said cytoplasmic domain consists of, in the following order, an N-terminal hexa-histidine tag containing a thrombin cleavage site fused to the N-terminus of the polypeptide of amino acids 564-958 of SEQ ID NO:1, except for an internal deletion of amino acids 711-761 and wherein said crystallizable composition forms isomorphous crystals.

3. A method for identifying a candidate inhibitor that interacts with a binding site of an FMS-like tyrosine kinase 3 (FLT3) protein, comprising the steps of:
(a) diffracting the crystal according to claim 1 with X-rays to obtain a set of diffraction patterns;
(b) using the diffraction patterns in step a) to solve the three-dimensional protein structure and obtain structure coordinates for the amino acids of the crystal of step (a), wherein the structure coordinates are set forth in FIG. 1A-1 to 1A-50;
(c) generating a three-dimensional model of the cytoplasmic domain of said FLT3 protein using the structure coordinates of the amino acids obtained in step (b), wherein said structure coordinates have a root mean square deviation from the backbone atoms of said amino acids of no more than ±2.0 Å;
(d) determining a binding site of the cytoplasmic domain of said FLT3 protein from said three-dimensional model; and
(e) performing computer fitting analysis to identify a candidate inhibitor or inhibitors which interact with said binding site, wherein the candidate inhibitor or inhibitors that interact with said binding site are identified as a candidate inhibitor or inhibitors of FLT3.

4. The method according to claim 3, further comprising the step of:
(f) contacting the identified candidate inhibitor or inhibitors with the cytoplasmic domain of said FLT3 protein in vitro in order to determine the effect of the inhibitor on FLT3 protein activity.

5. The method according to claim 3, wherein the binding site of the cytoplasmic domain of said FLT3 protein determined in step (d) comprises the structure coordinates according to FIGS. 1A-1 to 1A-50 of amino acid residues His 809, Arg 810, and Asp 811, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.

6. The method according to claim 3, wherein the binding site of the cytoplasmic domain of said FLT3 protein determined in step (d) comprises the structure coordinates according to FIGS. 1A-1 to 1A-50 of amino acid residues Phe 621, Glu 661, Met 664, Leu 802, Val 808, Arg 810, Asp 829, and Leu 832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.

7. The method according to claim 3, wherein the binding site of the cytoplasmic domain of said FLT3 protein determined in step (d) comprises the structure coordinates according to FIGS. 1A-1 to 1A-50 of amino acid residues Phe 621, Lys 644, Ala 657, Glu 661, Met 664, Leu 802, Ser 806, Cys 807, Val 808, His 809, Arg 810, Asp 811, Asp 829, and Leu 832, wherein the root mean square deviation from the backbone atoms of said amino acids is not more than ±2.0 Å.

8. A method for identifying a candidate inhibitor that interacts with a binding site of the cytoplasmic domain of an FLT3 protein, comprising the steps of:
(a) diffracting the crystal according to claim 1 with X-rays to obtain a set of diffraction patterns;
(b) using the diffraction patterns in step a) to solve the three-dimensional protein structure and obtain structure coordinates for the amino acids of the crystal of step (a);
(c) generating a three-dimensional model of said FLT3 protein using the structure coordinates of the amino acids generated in step (b), wherein the structure coordinates have a root mean square deviation from backbone atoms of said amino acids of not more than ±2.0 Å;
(d) determining a binding site of the cytoplasmic domain of said FLT3 protein from said three-dimensional model; and
(e) performing computer fitting analysis to identify a candidate inhibitor which interacts with said binding site, wherein a candidate inhibitor that interacts with said binding site is identified as a candidate inhibitor of FLT3.

9. The method according to claim 8, further comprising the step of:
(f) contacting the identified candidate inhibitor with the cytoplasmic domain of said FLT3 protein in vitro in order to determine the effect of the inhibitor on FLT3 protein activity.

10. The method according to claim 8, wherein the binding site of the cytoplasmic domain of said FMS-like tyrosine kinase 3 protein determined in step (d) comprises amino acid residues His 809, Arg 810, and Asp 811.

11. The method according to claim 8, wherein the binding site of the cytoplasmic domain of said FLT3 protein determined in step (d) comprises amino acid residues Phe 621, Glu 661, Met 664, Leu 802, Val 808, Arg 810, Asp 829, and Leu 832.

12. The method according to claim 8, wherein the binding site of the cytoplasmic domain of said FLT3 protein determined in step (d) comprises amino acid residues Phe 621, Lys 644, Ala 657, Glu 661, Met 664, Leu 802, Ser 806, Cys 807, Val 808, His 809, Arg 810, Asp 811, Asp 829, and Leu 832.