The present invention relates to the design of crystallisable Polo-like Kinase 3 (PLK3) complexes and the X-ray analysis of crystalline molecules or molecular complexes of this protein. The present invention provides for the first time the crystal structure of PLK3 protein bound to adenosine. The present invention also provides crystalline molecules or molecular complexes that comprise binding pockets of PLK3 kinase (PLK3) and/or its structural homologues, the structure of these molecules or molecular complexes. The present invention further provides crystals of PLK3 complexed with adenosine and methods for producing these crystals. This invention also relates to crystallizable compositions from which the protein-ligand complexes may be obtained. The present invention also relates to a data storage medium encoded with the structural coordinates of molecules and molecular complexes that comprise the ATP-binding pockets of PLK3 or their structural homologues. The present invention also relates to a computer comprising such data storage material. The computer may generate a three-dimensional structure or graphical three-dimensional representation of such molecules or molecular complexes. This invention also relates to methods of using the structure coordinates to solve the structure of homologous proteins or protein complexes. This invention also relates to computational methods of using structure coordinates of the PLK3 complex(es) to screen for, identify and design compounds, including inhibitory compounds and antibodies, that interact with PLK3 or homologues thereof.
The search for new therapeutic agents has been greatly aided in recent years by a better understanding of the structure of enzymes and other biomolecules associated with diseases. One important class of enzymes that has been the subject of extensive study is protein kinases.
Protein kinases constitute a large family of structurally related enzymes that are responsible for the control of a variety of signal transduction processes within the cell. (See, Hardie, G. and Hanks, S. The Protein Kinase Facts Book, I and II, Academic Press, San Diego, Calif.: 1995). Protein kinases are thought to have evolved from a common ancestral gene due to the conservation of their structure and catalytic function. Almost all kinases contain a similar 250-300 amino acid catalytic domain. The kinases may be categorized into families by the substrates they phosphorylate (e.g., protein-tyrosine, protein-serine/threonine, lipids, etc.). Sequence motifs have been identified that generally correspond to each of these kinase families (See, for example, Hanks, S. K., Hunter, T., FASEB J., 9:576-596 (1995); Knighton et al., Science, 253:407-414 (1991); Hiles et al., Cell, 70:419-429 (1992); Kunz et al., Cell, 73:585-596 (1993); Garcia-Bustos et al., EMBO J., 13:2352-2361 (1994)).
In general, protein kinases mediate intracellular signaling by effecting a phosphoryl transfer from a nucleoside triphosphate to a protein acceptor that is involved in a signaling pathway. These phosphorylation events act as molecular on/off switches that can modulate or regulate the target protein biological function. These phosphorylation events are ultimately triggered in response to a variety of extracellular and other stimuli. Examples of such stimuli include environmental and chemical stress signals (e.g., osmotic shock, heat shock, ultraviolet radiation, bacterial endotoxin, and H2O2), cytokines (e.g., interleukin-1 (IL-1) and tumor necrosis factor α (TNF-α)), and growth factors (e.g., granulocyte macrophage-colony-stimulating factor (GM-CSF), and fibroblast growth factor (FGF)). An extracellular stimulus may affect one or more cellular responses related to cell growth, migration, differentiation, secretion of hormones, activation of transcription factors, muscle contraction, glucose metabolism, control of protein synthesis, and regulation of the cell cycle.
Many diseases are associated with abnormal cellular responses triggered by protein kinase-mediated events as described above. These diseases include, but are not limited to, autoimmune diseases, inflammatory diseases, bone diseases, metabolic diseases, neurological and neurodegenerative diseases, cancer, cardiovascular diseases, allergies and asthma, Alzheimer's disease, and hormone-related diseases. Accordingly, there has been a substantial effort in medicinal chemistry to find protein kinase inhibitors that are effective as therapeutic agents.
Among medically important serine/threonine kinases is the family of polo-like kinases that have an important role in several stages of mitosis. This family includes the mammalian PLK1, PLK2 (serum-inducible kinase, Snk), PLK3 (fibroblast growth factor-inducible kinase, Fnk/proliferation-related kinase, Prk) and PLK4 (Sak). Polo-like kinase 1 (PLK1) has been suggested to promote mitotic entry by activating cyclin B/cdk1 by phosphorylating cyclin B itself, by phosphorylating the Cdk1-activating phosphatase Cdc25c and by phosphorylating the Cdk1-inhibiting kinases Myt1/Wee1 (Abrieu, A., et al. Cell Sci 111, 1751-1757 (1998); Jackman, M. et al Nat Cell Biol. 5, 143-148 (2003); Nakajima, H. et al. J. Biol. Chem. 278, 25277-25280 (2003); Qian, Y. W. et al. Mol. Cell. Biol. 18, 4262-4271 (1998); Toyoshima-Morimoto, F. et al. Nature 410, 215-220 (2001); Watanabe, N. et al. Proc. Nat. Acad. Sci. USA 101, 4419-4424 (2004)). Furthermore, PLK1 has been implicated in centrosome maturation and separation, with defects giving rise to monopolar spindles (Lane, H. A., et al. J. Cell Biol. 135, 1701-1713 (1996); Sunkel, C. E., et al. J. Cell Sci. 89, 25-38 (1988)). Studies on human cells and in vitro have indicated a role for PLK1 in activating the anaphase-promoting complex/cyclosome (APC/C) (Golan, A. et al J. Biol. Chem. 277, 15552-15557 (2002); Kotani, S. et al. Mol. Cell 1, 371-380 (1998). PLK1 is located mostly in the cytoplasm during interphase and translocates to the nucleus in early mitosis. PLK1 also has been shown to associate with centrosomes from G2 up to metaphase, to translocate to kinetochores at metaphase, and to locate at the midbody from anaphase to telophase (Golsteyn, R. M. et al. J. Cell Biol. 129, 1617-1628 (1995); Lee et al. Mol. Cell Biol. 15, 7143-7151 (1995)).
In contrast to PLK1, PLK2 and PLK3 appear to perform different functions, and PLK3 plays some functions that might directly antagonize PLK1 function (Smits et al. Nat. Cell Biol. 2 672-676 (2000); Xie et al. J. Biol. Chem. 276, 43305-43312 (2001). Small interfering RNA (siRNA)-mediated inhibition of PLK3 expression suggests that PLK3 is essential for superoxide-induced cell death, and deletion analysis of PLK3 showed that the N-terminal domain (amino-acids 1-26) is essential for induction of delayed onset of apoptosis (Li Z. et al. J. Biol. Chem. 280, 16843-16850 (2005)). Over-expression of PLK1 has been observed to cause oncogenic transformation in NIH 3T3 cells (Smith et al. Biochem. Biophys. Res. Commun. 234, 397-405 (1997), while overexpression of PLK3 induces apoptosis (Conn et al. Cancer Res. 60, 6826-6831 (2000).
Polo-like kinases are characterized by an N-terminal catalytic Ser/Thr kinase domain and a C-terminal non-catalytic region containing two tandem Polo-boxes (Sonnhammer et al. Nucleic Acids Res. 26, 320-322 (1998). The Polo-box domain is required for correct subcellular localization (Lee et al., Proc. Natl. Acad. Sci. USA 95, 9301-9306 (1998); Ma et al., Mol. Cancer. Res. 1, 376-384 (2003); May et al., J. Cell Biol. 156, 23-28 (2002); Reynolds & Ohkura J. Cell Sci., 116, 1377-1387 (2003)), and the crystal structure of the polo-box domain in PLK1 has been determined (Elia et al. Cell 115, 83-95 (2003); Cheng et al., EMBO J. 22, 5757-5768 (2003)). The Polo-box domain recognizes a conserved phosphothreonine/serine motif in the substrate which binds along a positively-charged cleft located at the edge of the Polo-box interface. Mutations of key residues to disrupt the phosphodependent interactions abolishes cell-cycle dependent localization but does not affect its kinase activity (Lee et al., Proc. Natl. Acad. Sci. USA 95, 9301-9306 (1998), and suggests that substrate binding to the Polo-box domain is required for proper mitotic progression (Elia et al. Cell 115, 83-95 (2003)).
Because of their link to cellular proliferation, PLK family proteins have been associated with cancer development and progression. In fact, the over-expression of murine PLK1 in NIH3T3 cells has been observed to be transforming (Smith, M. R. et al. Biochem. Biophys. Res. Commun. 234, 397-405 (1997)). In addition, PLK1 has been seen to be over-expressed in a wide variety of primary tumor cell lines, including breast cancer (Wolf G. et al. Pathol. Res. Pract. 196, 753-759 (2000)), colorectal cancer (Takahashi T. et al. Cancer Sci. 94, 148-152 (2003)), non-small-cell lung cells (Wolf G. et al. Oncogene 14, 543-549 (1997)), head and neck squamous cell carcinomas (Knecht R. et al. Cancer Res. 59, 479-480 (1999)), ovarian carcinomas (Takai N. et al. Cancer Lett. 164, 41-49 (2001)), prostate (Weichert W. et al Prostate 60, 240-245 (2004)) and pancreatic cancer (Gray P. J. et al. Mol. Cancer. Ther. 3, 641-646 (2004)) and is positively correlated with aggressiveness and prognosis (Takai N., et al Oncogene 24, 287-291 (2005)).
Whilst the precise function of PLK2 is uncertain, the expression of PLK3 is negatively correlated with the development of certain cancers. PLK3 mRNA is either undetectable or markedly downregulated in lung carcinomas (Li B. et al. J. Biol. Chem. 271, 19402-19408 (1996)) head and neck squamous cell carcinomas (Dai W. et al. Genes Chromosomes Cancer 27, 332-336 (2000)) and in colon tumors (Dai W. et al. Int. J. Oncol. 20, 121-126 (2002)). In support of these observations, the expression of PLK3 in vitro reduces the rate of fibroblast cell proliferation (Dai W. et al. Genes Chromosomes Cancer 27, 332-336 (2000)) and the enforced expression of constructs expressing kinase-active PLK3 induces rapid cell cycle arrest and apoptosis (Conn C. W. et al. Cancer Res. 60, 6826-6831 (2000); Wang Q. et al. Mol. Cell. Biol. 22, 3450-3459 (2002)). PLK3 has been shown to be a stress response protein involved in DNA damage checkpoint (Bahassi E. M. et al. Oncogene 21, 6633-6640 (2002); Xie S. et al. J. Biol. Chem. 276, 43305-43312 (2001)). In this role, PLK3 only becomes phosphorylated following DNA damage or mitotic spindle disruption, whereupon its kinase activity is enhanced, otherwise, PLK3 is not phosphorylated during normal cell cycle progression. Consistent with this, possible substrates of PLK3 are thought to include both the tumor suppressors Chk2 and p53, which have both been shown to bind to PLK3 in vivo, and be phosphorylated by PLK3 in vitro.
The role of PLK3 in DNA damage checkpoint leads to the interest in designing PLK3 inhibitors that are effective as therapeutic agents. A challenge has been to find protein kinase inhibitors that act in a selective manner, targeting only PLK3. Since there are numerous protein kinases that are involved in a variety of cellular responses, non-selective inhibitors may lead to unwanted side effects. In this regard, the three-dimensional structure of the kinase would assist in the rational design of inhibitors. The determination of the amino acid residues in PLK3 binding pockets and the determination of the shape of those binding pockets would allow one to design selective inhibitors that bind favorably to this class of enzymes. The determination of the amino acid residues in PLK3 binding pockets and the determination of the shape of those binding pockets would also allow one to design inhibitors that can bind to PLK3.
For example, a general approach to designing inhibitors that are selective for an enzyme target is to determine how a putative inhibitor interacts with the three dimensional structure of the enzyme. For this reason it is useful to obtain the enzyme protein in crystal form and perform X-ray diffraction techniques to determine its three dimensional structure coordinates. If the enzyme is crystallized as a complex with a ligand, one can interactively elucidate the binding pocket, and determine both the shape of the enzyme binding pocket when bound to the ligand, as well as the amino acid residues that are capable of close contact with the ligand. By knowing the shape and amino acid residues in the binding pocket, one may design new ligands that will interact favorably with the enzyme. With such structural information, available computational methods may be used to predict how strong the ligand binding interaction will be. Such methods thus enable the design of inhibitors that bind strongly, as well as selectively to the target enzyme.
Despite the fact that the genes for PLK3 has been isolated and the amino acid sequence of PLK3 is known, no one has described X-ray crystal structural coordinate information of PLK3 protein. As discussed above, such information would be extremely useful in identifying and designing potential inhibitors of the PLK3 kinase or homologues thereof, which, in turn, could have therapeutic utility.
The structures of several serine/threonine kinases have been solved by X-ray diffraction and analyzed. Specifically, the crystal structures of P38 kinase (Wilson et al., J. Biol. Chem., 271, pp. 27696-27700 (1996)) and MAPKAP Kinase 2 (U.S. Provisional application 60/337,513) have been studied in detail.
To date, no crystal structures of PLK3 kinase have been reported. Thus the crystal structure of unphosphorylated PLK3 kinase domain complexes with inhibitors are of great importance for defining the active conformation of PLK3 kinase. This information is essential for the rational design of selective and potent inhibitors of PLK3.
The present invention provides for the first time, crystallizable compositions, crystals, and the crystal structures of a PLK3—inhibitor complex. The PLK3 protein used in these studies corresponds to a single polypeptide chain, which encompasses the complete catalytic kinase domain, amino acids 48 to 332. Solving this crystal structure has allowed the applicants to determine the key structural features of PLK3, particularly the shape of its substrate and ATP-binding pockets.
Thus, in one aspect, the present invention provides molecules or molecular complexes comprising all or parts of these binding pockets, or homologues of these binding pockets that have similar three-dimensional shapes.
In another aspect, the present invention further provides crystals of PLK3 complexed with adenosine and methods for producing these crystals. In this embodiment, PLK3 is unphosphorylated.
In a further aspect, the present invention provides crystallizable compositions from which PLK3-ligand complexes may be obtained.
In another aspect, the present invention provides for a general strategy for the design of protein constructs for producing crystallisable kinase domains.
In another aspect, the invention provides a data storage medium that comprises the structure coordinates of molecules and molecular complexes that comprise all or part of the PLK3 binding pockets. Such storage medium encoded with these data when read and utilized by a computer programmed with appropriate software displays, on a computer screen or similar viewing device, a three-dimensional graphical representation of a molecule or molecular complex comprising such binding pockets or similarly shaped homologous binding pockets.
In yet another aspect, the invention provides computational methods of using structure coordinates of the PLK3 complex to screen for, identify and design compounds, including inhibitory compounds and antibodies that interact with PLK3, PLK1, PLK2 or other homologues thereof. In certain embodiments, the invention provides methods for designing, evaluating and identifying compounds, which bind to the aforementioned binding pockets. In certain embodiments, such compounds are potential inhibitors of PLK3, PLK1, PLK2 or their homologues.
The following abbreviations are used in
“Atom type” refers to the identity of the element whose coordinates are measured. The first letter in the column defines the actual type of atom (i.e., C for carbon, N for nitrogen, O for oxygen, etc.).
“Resid” refers to the amino acid residue or molecule identity in the molecular model. The amino acid residues or molecule identity are designated A, B and W for PLK3, adenosine (ADN) or water (HOH), respectively.
“#” refers to the amino acid residue number.
“X, Y, Z” crystallographically define the atomic position of the element measured.
“B” is a thermal factor that measures movement of the atom around its atomic center.
“Occ” is an occupancy factor that refers to the fraction of the molecules in which each atom occupies the position specified by the coordinates. A value of “1” indicates that each atom has the same conformation, i.e., the same position, in all molecules of the crystal.
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 exclusion of any other integer or groups of integers.
The following abbreviations are used throughout the application:
Additional definitions are set forth below.
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 hydrogen bonding or Van der Waals or electrostatic interactions energetically favor the juxtaposition—or it may be covalent.
The term “binding pocket”, as used herein, refers to a region of a molecule or molecular complex, that, as a result of its shape and charge, favorably associates with another chemical entity or compound. The term “pocket” includes, but is not limited to, cleft, channel or site. PLK or PLK-like molecules may have binding pockets which include, but are not limited to, peptide or substrate binding, ATP-binding and antibody binding sites.
The term “chemical entity”, as used herein, refers to chemical compounds, complexes of at least two chemical compounds, and fragments of such compounds or complexes. The chemical entity may be, for example, a ligand, a substrate, a nucleotide triphosphate, a nucleotide diphosphate, phosphate, a nucleotide, an agonist, antagonist, inhibitor, antibody, drug, peptide, protein or compound.
“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, pp. 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 “corresponding amino acid” or “residue that corresponds to” refers to a particular amino acid or analogue thereof in a PLK3 homologue that corresponds to an amino acid in the PLK3 structure. The corresponding amino acid may be an identical, mutated, chemically modified, conserved, conservatively substituted, functionally equivalent or homologous amino acid when compared to the PLK3 amino acid to which it corresponds.
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 PLK3 structure. For example, corresponding amino acids may be identified by superimposing the backbone atoms of the amino acids in PLK3 and the PLK3 homologue using well known software applications, such as QUANTA® [Molecular Simulations, Inc., San Diego, Calif.; now part of Accelrys, San Diego, Calif.]. The corresponding amino acids may also be identified using 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 Adv. Appl. Math., 2, 482 (1981), which is incorporated herein by reference.
The term “domain” refers to a portion of the PLK3 protein or homologue that can be separated according to its biological function, for example, catalysis. The domain is usually conserved in sequence or structure when compared to other kinases or related proteins. The domain can comprise a binding pocket, or a sequence or structural motif.
The term “sub-domain” refers to a portion of the domain as defined above in the PLK3 protein or homologue. The catalytic kinase domain (amino acid residues 48 to 332) of PLK3 is a bi-lobal structure consisting of an N-terminal, β-strand sub-domain (residues 48 to 139) and a C-terminal, α-helical sub-domain (residues 144 to 332).
The term “catalytic active site” refers to the area of the protein kinase to which nucleotide substrates bind. The catalytic active site of plk3 is at the interface between the N-terminal, β-strand sub-domain and the C-terminal, α-helical sub-domain.
The “PLK3 ATP-binding pocket” of a molecule or molecular complex is defined by the structure coordinates of a certain set of amino acid residues present in the PLK3 structure, as described below. In general, the ligand for the ATP-binding pocket is a nucleotide such as ATP. This binding pocket is in the catalytic active site of the kinase domain. In the protein kinase family, the ATP-binding pocket is generally located at the interface of the α-helical and β-strand sub-domains, and is bordered by the glycine rich loop and the hinge [See, Xie et al., Structure, 6, pp. 983-991 (1998), incorporated herein by reference].
The term “PLK3-like” refers to all or a portion of a molecule or molecular complex that has a commonality of shape to all or a portion of the PLK3 protein. In the PLK3-like ATP-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 PLK3-like ATP-binding pocket and the amino acids in the PLK3 ATP-binding pocket (as set forth in
The term “part of an PLK3 ATP-binding pocket” or “part of an PLK3-like ATP-binding pocket” refers to less than all of the amino acid residues that define the PLK3 or PLK3-like ATP-binding pocket. The structure coordinates of residues that constitute part of an PLK3 or PLK3-like ATP-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 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 residues may be contiguous or non-contiguous in primary sequence. In one embodiment, part of the PLK3 or PLK3-like ATP-binding pocket is at least two amino acid residues, preferably, E140 and C142. In another embodiment, the amino acids are selected from the group consisting of L68, H114, V123, L139, E140, C142, K145 and F192.
The term “PLK3 kinase domain” refers to the catalytic domain of PLK3. The kinase domain includes, for example, the catalytic active site, which comprises the catalytic residues, the activation loop or phosphorylation lip, the DFG motif, and the glycine-rich phosphate anchor or glycine-rich loop [See, Xie et al., Structure, 6, pp. 983-991 (1998); R. Giet and C. Prigent, J. Cell Sci., 112, pp. 3591-3601 (1999), incorporated herein by reference]. The kinase domain in the PLK3 protein comprises residues from about 48 to 332.
The term “part of a PLK3 kinase domain” or “part of a PLK3-like kinase domain” refers to a portion of the PLK3 or PLK3-like catalytic domain. The structure coordinates of residues that constitute part of a PLK3 or PLK3-like kinase domain may be specific for defining the chemical environment of the domain, or useful in designing fragments of an inhibitor that may interact with those residues. For example, the portion of 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 domain. The residues may be contiguous or non-contiguous in primary sequence. For example, part of a PLK3 kinase domain can be the active site, the glycine-rich loop, the activation loop, or the catalytic loop [see Xie et al., supra].
The term “homologue of PLK3” or “PLK3 homologue” refers to a molecule or molecular complex that is homologous to PLK3 by three-dimensional structure or sequence and retains the kinase activity of PLK3. Examples of homologues include but are not limited to the following: human PLK3 with mutations, conservative substitutions, additions, deletions or a combination thereof; PLK3 from a species other than human; a protein comprising an PLK3-like ATP-binding pocket, a kinase domain; another member of the protein kinase family, preferably the Polo-like kinase family or the CDK kinase family.
The term “part of a PLK3 protein” or “part of a PLK3 homologue” refers to a portion of the amino acid residues of a PLK3 protein or homologue. In one embodiment, part of a PLK3 protein or homologue defines the binding pockets, domains, sub-domains, and motifs of the protein or homologue. The structure coordinates of residues that constitute part of a PLK3 protein or homologue may be specific for defining the chemical environment of the protein, or useful in designing fragments of an inhibitor that may interact with those residues. The portion of residues may also be residues that are spatially related and define a three-dimensional compartment of a binding pocket, motif or domain. The residues may be contiguous or non-contiguous in primary sequence. For example, the portion of 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 “PLK3 protein complex” or “PLK3 homologue complex” refers to a molecular complex formed by associating the PLK3 protein or PLK3 homologue with a chemical entity, for example, a ligand, a substrate, nucleotide triphosphate, an agonist or antagonist, inhibitor, drug or compound. In one embodiment, the chemical entity is selected from the group consisting of an ATP, a nucleotide triphosphate and an inhibitor for the ATP-binding pocket. In another embodiment, the inhibitor is an ATP analog such as MgAMP-PNP (adenylyl imidodiphosphate) or adenosine.
The term “motif” refers to a portion of the PLK3 protein or homologue that defines a structural compartment or carries out a function in the protein, for example, catalysis, structural stabilization, or phosphorylation. The motif may be conserved in sequence, structure and function when compared to other kinases or related proteins. The motif can be contiguous in primary sequence or three-dimensional space. The motif can comprise α-helices and/or β-sheets. Examples of a motif include but are not limited to a binding pocket, active site, phosphorylation lip or activation loop, the glycine-rich phosphate anchor loop and the catalytic loop [See, Xie et al., Structure, 6, pp. 983-991 (1998); R. Giet and C. Prigent, J. Cell Sci., 112, pp. 3591-3601 (1999)].
The term “root mean square deviation” or “RMSD” means the square root of the arithmetic mean of the squares of the deviations from the mean. It is a way to express the deviation or variation from a trend or object. For purposes of this invention, the “root mean square deviation” defines the variation in the backbone of a protein from the backbone of PLK3, a binding pocket, a motif, a domain, or portion thereof, as defined by the structure coordinates of PLK3 described herein.
The term “sufficiently homologous to PLK3” refers to a protein that has a sequence homology of at least 35% compared to PLK3 protein. In one embodiment, the sequence homology is at least 40%, at least 60%, at least 80%, at least 90% or at least 95%.
The term “soaked” refers to a process in which the crystal is transferred to a solution containing the compound of interest. In certain embodiments, the compound is diffused into the crystal.
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. It would be readily apparent to those skilled in the art that all or part of the structure coordinates of
The term “about” when used in the context of RMSD values takes into consideration the standard error of the RMSD value, which is ±0.1 Å.
The term “crystallization solution” refers to a solution that promotes crystallization. The solution comprises at least one agent, and may include 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 a stabilizer.
The term “generating a three-dimensional structure” or “generating a three-dimensional graphical representation” refers to converting the lists of structure coordinates into structural models in three-dimensional space. This can be achieved through commercially or publicly available software. The three-dimensional structure may be displayed as a graphical representation or used to perform computer modeling or fitting operations. In addition, the structure coordinates themselves may be used to perform computer modeling and fitting operations.
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 “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 a PLK3 molecule or molecular complex or homologues thereof, calculating or minimizing energies for an association of a PLK3 molecule or molecular complex or homologues thereof to a chemical entity.
Crystallizable Compositions and Crystals of PLK3 Complexes
According to one embodiment, the invention provides a crystal or crystallizable composition comprising unphosphorylated PLK3 protein. In one embodiment, the PLK3 protein is human PLK3 protein or a catalytic domain thereof. In another embodiment, the invention provides a crystal or crystallizable composition comprising unphosphorylated PLK3 catalytic domain. In another embodiment, the invention provides a crystal or a crystallizable composition comprising unphosphorylated PLK3 protein, or catalytic domain thereof, and a substrate analogue, such as, but not limited, to adenosine. In one embodiment, the invention provides a crystal or a crystallizable composition comprising unphosphorylated PLK3 protein, or catalytic domain thereof, and adenosine. In one embodiment, the aforementioned crystallizable composition further comprises a precipitant, 900-1100 mM sodium dihydrogen phosphate/1100-900 mM potassium hydrogen phosphate (i.e., the concentration of sodium di-hydrogen phosphate plus the concentration of di-potassium hydrogen phosphate is 2M) and a buffer that maintains pH at between about 6.6 and 8.6. The composition may further comprise a reducing agent, such as dithiothreitol (DTT) at between about 1 to 10 mM. The unphosphorylated PLK3 protein or complex is preferably 85-100% pure prior to forming the composition.
According to another embodiment, the invention provides a crystal composition comprising PLK3 protein complex. In one embodiment, the crystal has a unit cell dimension of a=146 Å, b=59 Å, c=47 Å, α=γ=90°, β=95° and belongs to space group C2. It will be readily apparent to those skilled in the art that the unit cells of the crystal compositions may deviate ±1-2 Å from the above cell dimensions and the beta angle may vary ±1-2 degrees depending on the deviation in the unit cell calculations.
As used herein, the PLK3 protein in the crystal or crystallizable compositions can be a truncated protein with amino acids 48 to 332 of SEQ ID NO:1; and the truncated protein with conservative substitutions.
The PLK3 protein 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. Preferably, the protein is overexpressed from a baculovirus system.
The invention also relates to a method of making crystals of PLK3 complexes or PLK3 homologue complexes. Such methods comprise the steps of:
In one embodiment, the invention provides a crystal or a crystallizable composition comprising unphosphorylated PLK3 protein, or catalytic domain thereof, and a chemical entity. In one embodiment, the chemical entity is selected from the group consisting of an ATP analogue, nucleotide triphosphate, nucleotide diphosphate, phosphate, adenosine, or an active site inhibitor. In another embodiment, the chemical entity is selected from the group consisting of adenosine, ATP, an ATP analogue, AMP-PNP, a nucleotide triphosphate, a nucleotide diphosphate, phosphate and an active site inhibitor. In another embodiment, the chemical entity is an ATP analogue. In yet another embodiment, the crystallization solution is as described previously. In another embodiment, the composition is treated with micro-crystals of PLK3 or PLK3 complexes or homologues thereof. In another embodiment, the composition is treated with micro-crystals of PLK3 complexes or homologues thereof.
In certain embodiments, the invention provides a method of making PLK3 crystals, the method comprising steps of:
In one embodiment, the crystallizable composition of step b) is made according to the conditions discussed above.
Devices for promoting crystallization can include but are not limited to the hanging-drop, sitting-drop, dialysis 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 et al., Proteins: Structure, Function, and Genetics, 20, pp. 98-102 (1994), incorporated herein by reference]. The hanging-drop or sitting-drop methods produce crystals by vapor diffusion. The hanging-drop, sitting-drop, and some adaptations of the microbatch methods [D'Arcy et al., J. Cryst. Growth, 168, pp. 175-180 (1996) and Chayen, J. Appl. Cryst., 30, pp. 198-202 (1997)] produce crystals by vapor diffusion. The hanging drop and sitting drop containing the crystallizable composition is equilibrated in a reservoir containing a higher or lower concentration of the precipitant. As the drop approaches equilibrium with the reservoir, the saturation of protein in the solution leads to the formation of crystals.
Microseeding or seeding may be used to obtain larger, or better quality (i.e., crystals with higher resolution diffraction or single crystals) crystals from initial micro-crystals. Microseeding involves the use of crystalline particles to provide nucleation under controlled crystallization conditions. Microseeding is used to increase the size and quality of crystals. In this instance, micro-crystals are crushed to yield a stock seed solution. The stock seed solution is diluted in series. Using a needle, glass rod 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 following the teachings of the specification to vary the crystallization conditions disclosed herein to identify other crystallization conditions that would produce crystals of PLK3 homologue, PLK3 homologue complex, PLK3 protein or other PLK3 protein complexes. 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 of crystallization, or introducing additives such as detergents (e.g., TWEEN 20 (monolaurate), LDAO, Brij 30 (4 lauryl ether)), sugars (e.g., glucose, maltose), organic compounds (e.g., dioxane, dimethylformamide), lanthanide ions or polyionic compounds that aid in crystallization. High throughput crystallization assays may also be used to assist in finding or optimizing the crystallization conditions.
Binding Pockets of PLK3 Protein or Homologues thereof
As disclosed above, applicants have provided for the first time the three-dimensional X-ray crystal structures of a PLK3-inhibitor complex. The crystal structure of PLK3 presented here is the first reported for PLK3. The invention is useful for inhibitor design to study the role of PLK3 in cell signaling. Atomic or structure coordinates are presented in
In order to use the structure coordinates generated for PLK3, their complexes, one of their binding pockets, or a PLK3-like binding pocket thereof, it is often at times necessary to convert the coordinates into a three-dimensional shape. This is achieved through the use of commercially available software that is capable of generating three-dimensional graphical representations (e.g., three-dimensional structures) of molecules or portions thereof from a set of structure coordinates.
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 ATP and substrate binding pockets of this invention will be important for drug design.
In one embodiment, the ATP-binding pocket comprises amino acids L68, G69, K70, G71, G72, F73, A74, R75, C76, A89, K91, E110, H114, V123, L139, E140, L141, C142, S143, R144, K145, S146, H149, G189, N190, F192, D203 and F204 using the structure of the PLK3-adenosine complex according to
In another embodiment, the ATP-binding pocket comprises amino acids R66, L67, L68, G69, K70, G71, G72, F73, A74, R75, C76, Y77, E78, Y88, A89, V90, K91, V92, I93, E110, H114, V123, R124, I137, L139, E140, L141, C142, S143, R144, K145, S146, L147, A148, H149, L188, G189, N190, F191, F192, I193, T194, E195, K200, V201, G202 and D203 using the structure of the PLK3—adenosine complex according to
In another embodiment, the ATP-binding pocket comprises amino acids R66, L67, L68, G69, K70, G71, G72, F73, A74, R75, C76, Y77, E78, Y88, A89, V90, K91, V92, I93, Q103, K106, I107, E110, I111, H114, V123, R124, F125, I137, F138, L139, E140, L141, C142, S143, R144, K145, S146, L147, A148, H149, W151, K152, H183, D185, K187, L188, G189, N190, F191, F192, I193, T194, E195, K200, V201, G202, D203, F204, G205, L206, T223 and Y226 using the structure of the PLK3—adenosine complex according to
In one embodiment, the ATP-binding pocket comprises the amino acids of L68, H114, V123, L139, E140, C142, K145 and F192 according to
Those of skill in the art understand that a set of structure coordinates for a molecule or a molecular-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 PLK3 structure coordinates. For example, the structure coordinates set forth in
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 PLK3 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 maybe necessary to determine whether a binding pocket, motif, domain or portion thereof of a molecule or molecular complex is sufficiently similar to the binding pocket, motif, domain or portion thereof of PLK3. Such analyses may be carried out in well known software applications, such as ProFit [A. C. R. Martin, SciTech Software, ProFit version 1.8, University College London, http://www.bioinf.org.uk/software], Swiss-Pdb Viewer [Guex et al., Electrophoresis, 18, pp. 2714-2723 (1997)], the Molecular Similarity application of QUANTA® [Molecular Simulations, Inc., San Diego, Calif.; now part of Accelrys, San Diego, Calif.] 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.; now part of Accelrys, San Diego, Calif.] 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 the above programs is defined by user input, for the purpose of this invention we will define equivalent atoms as protein backbone atoms (N, Cα, C and O) for PLK3 amino acids and corresponding amino acids in the 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 (1981), which is incorporated herein by reference. A suitable amino acid sequence alignment will require that the proteins being aligned share 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 with the second protein [Hanks et al., Science, 241, 42 (1988); Hanks and Quinn, Meth. Enzymol., 200, 38 (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. For programs that calculate an average of the individual RMSD values of the backbone atoms, 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.
The RMSD values between other protein kinases the PLK3 protein complex (
aWee1: Squire et al. Structure 13, 541-550 (2005); PDB accession number 1X8B
bp38: Wilson et al., J. Biol. Chem., 271, pp. 27696-27700 (1996); Z. Wang et al., Proc. Natl. Acad. Sci. U.S.A., 94, pp. 2327-2332 (1997); PDB Accession number 1WFC
cphosphoinositide dependent protein kinase 1 (PDK1): Komander et al. Biochem J., 375, 255-261 (2003); PDB Accession code 10KY
dInterleukin-2 Tyrosine Kinase: Brown et al., J. Biol. Chem. 279, pp. 18727 (2004); PDB Accession number 1SM2.
eMitogen activated protein kinase activated protein (MAPKAP) kinase 2: U.S. application Ser. No. 10/497,767 which is a national stage entry of PCT/US02/39070 which claims priority to U.S. Provisional Application No. 60/337,513, filed Dec. 5, 2001.
fHuman cyclin-dependent kinase 2 (cdk2): Brown et al. J. Biol. Chem. 274, 8746-8752 (1999)
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
Therefore, one embodiment of this invention provides a molecule or molecular complex comprising all or part of a PLK3 ATP-binding pocket defined by structure coordinates of PLK3 amino acids L68, G69, K70, G71, G72, F73, A74, R75, C76, A89, K91, E110, H114, V123, L139, E140, L141, C142, S143, R144, K145, S146, H149, G189, N190, F192, D203 and F204 according to
Another embodiment of this invention provides a molecule or molecular complex comprising all or part of an PLK3 ATP-binding pocket defined by structure coordinates of PLK3 amino acids R66, L67, L68, G69, K70, G71, G72, F73, A74, R75, C76, Y77, E78, Y88, A89, V90, K91, V92, I93, E110, H114, V123, R124, I137, L139, E140, L141, C142, S143, R144, K145, S146, L147, A148, H149, L188, G189, N190, F191, F192, I193, T194, E195, K200, V201, G202 and D203 according to
Yet other embodiments of this invention provide a molecule or molecular complex comprising all or part of other PLK3 ATP-binding pockets defined in this invention and by the structure coordinates of PLK3 amino acids according to
In one embodiment, the above molecules or molecular complexes are in crystalline form.
Computer Systems
According to another embodiment of this invention is provided a machine-readable data storage medium, comprising a data storage material encoded with machine-readable data, wherein said data comprises all or part of an PLK3 ATP-binding pocket defined by structure coordinates of PLK3 amino acids L68, G69, K70, G71, G72, F73, A74, R75, C76, A89, K91, E110, H114, V123, L139, E140, L141, C142, S143, R144, K145, S146, H149, G189, N190, F192, D203 and F204, according to
In other embodiments of this invention is provided a machine-readable data storage medium, comprising a data storage material encoded with machine-readable data, wherein said data comprises all or part of any molecule or molecular complex discussed in the above paragraphs.
In one embodiment of this invention is provided a computer comprising:
In one embodiment, a computer according to this invention comprises a working memory for storing instructions for processing the machine-readable data, a central-processing unit coupled to the working memory and to said machine-readable data storage medium for processing said machine-readable data into the three-dimensional structure. In one embodiment, the computer further comprises a display for displaying the three-dimensional structure as a graphical representation. In another embodiment, the computer further comprises commercially available software program to display the graphical representation. Examples of software programs include but are not limited to QUANTA® [Molecular Simulations, Inc., San Diego, Calif.; now part of Accelrys, San Diego, Calif.], O [Jones et al., Acta Cryst. A, 47, pp. 110-119 (1991)] and RIBBONS [M. Carson, J. Appl. Cryst., 24, pp. 958-961 (1991)], which are incorporated herein by reference.
In another embodiment, a computer according to this invention comprises executable code for:
This invention also provides a computer comprising:
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 a PLK3 molecule or molecular complex or homologues thereof, or calculating or minimizing energies for an association of a PLK3 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® [Molecular Simulations, Inc., San Diego, Calif.; now part of Accelrys, San Diego, Calif.], O [Jones et al., Acta Crystallogr. A47, pp. 110-119 (1991)] and RIBBONS [Carson, J. Appl. Crystallogr., 24, pp. 9589-961 (1991)], Coot [Emsley, P. and Cowtan, K. Acta Crystallogr. D60, pp. 2126-2132 Part 12 Sp. Iss. 1 (December 2004)], and PyMol™ [DeLano, W. L. The PyMOL Molecular Graphics System (2002) DeLano Scientific™, Palo Alto, Calif., USA.], 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, printers, modems, facsimile machines, CD-ROMs or disk drives. The information can be in graphical or alphanumeric form.
Input hardware 36, coupled to computer 11 by input lines 30, may be implemented in a variety of ways. Machine-readable data of this invention may be inputted via the use of a modem or modems 32 connected by a telephone line or dedicated data line 34. Alternatively or additionally, the input hardware 36 may comprise CD-ROM drives or disk drives 24. In conjunction with display terminal 26, keyboard 28 may also be used as an input device.
Output hardware 46, coupled to computer 11 by output lines 40, may similarly be implemented by conventional devices. By way of example, output hardware 46 may include CRT display terminal 26 for displaying a graphical representation of a binding pocket of this invention using a program such as QUANTA® [Molecular Simulations, Inc., San Diego, Calif.; now part of Accelrys, San Diego, Calif.] as described herein. Output hardware might also include a printer 42, so that hard copy output may be produced, or a disk drive 24, to store system output for later use. Output hardware may also include a display terminal, a CD or DVD recorder, ZIP™ or JAZ™ drive, or other machine-readable data storage device.
In operation, CPU 20 coordinates the use of the various input and output devices 36, 46, coordinates data accesses from mass storage 24 and accesses to and from working memory 22, and determines the sequence of data processing steps. A number of programs may be used to process the machine-readable data of this invention. Such programs are discussed in reference to the computational methods of drug discovery as described herein. Specific references to components of the hardware system 10 are included as appropriate throughout the following description of the data storage medium.
The magnetic domains of coating 102 of medium 100 are polarized or oriented so as to encode in a manner that may be conventional, machine readable data such as that described herein, for execution by a system such as system 10 of
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 data defines the above-mentioned binding pockets by comprising the structure coordinates of said amino acid residues according to
To use the structure coordinates generated for PLK3 or PLK3 homologue, one of its binding pockets, motifs, domains, or portion thereof, it is at times necessary to convert them into a three-dimensional shape or to generate 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 of molecules or portions thereof from a set of structure coordinates. In one embodiment, the three-dimensional structure may be displayed as a graphical representation.
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 of any of the molecule or molecular complexes, or binding pockets thereof, that are described herein.
In certain embodiment, this invention also provides a computer for producing a three-dimensional structure of:
According to other embodiments, the computer is also for producing the three-dimensional structure of the aforementioned molecules and molecular complexes and comprises the corresponding machine-readable data storage mediums. In one embodiment, the three-dimensional structure is displayed as a graphical representation.
In one embodiment, the structure coordinates of said molecules or molecular complexes are produced by homology modeling of at least a portion of the structure coordinates of
For example, since the protein sequence of the catalytic domains of PLK3 and homologues thereof can be aligned relative to each other, it is possible to construct models of the structures of PLK3 homologues, particularly in the regions of the active site, using the PLK3 structure. Software programs that are useful in homology modeling include XALIGN [Wishart, D. S. et al., Comput. Appl. Biosci., 10, pp. 687-88 (1994)] and CLUSTAL W Alignment Tool [Higgins D. G. et al., Methods Enzymol, 266, pp. 383-402 (1996)]. 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 D. G. et al., Methods Enzymol, 266, pp. 383-402 (1996), which is incorporated by reference] can be used. To model the amino acid side chains of homologous PLK3 proteins, the amino acid residues in PLK3 can be replaced, using a computer graphics program such as “O” [Jones et al, (1991) Acta Crest. 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.
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 of the above molecules or molecular complexes (e.g., PLK3, homologues and portions thereof), or binding pockets thereof, can be stored in a machine-readable storage medium, which is capable of displaying three-dimensional structural information or a graphical three-dimensional representation of the structure.
Rational Drug Design
The PLK3 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. In certain embodiments, the computer is programmed with software to translate those coordinates into the three-dimensional structure of PLK3.
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 PLK3 may inhibit or activate PLK3 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.
Thus, according to another embodiment, the invention provides a method for designing, selecting and/or optimizing a chemical entity that binds to the molecule or molecular complex comprising the steps of:
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 structure or graphical representation; subtract distances between atoms; calculate chemical energies for a PLK3 molecule, molecular complex or homologues thereof; or calculate or minimize energies of an association of PLK3 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.; now part of Accelrys, San Diego, Calif.], O [Jones et al., Acta Crystallogr. A47, pp. 110-119 (1991)] and RIBBONS [Carson, J. Appl. Crystallogr., 24, pp. 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.
Another embodiment of the invention provides a method for evaluating the potential of a chemical entity to associate with the molecule or molecular complex as described previously.
This method comprises the steps of: a) employing computational means to perform a fitting operation between the chemical entity and the molecule or molecular complex described before; b) analyzing the results of said fitting operation to quantify the association between the chemical entity and the molecule or molecular complex; and, optionally, c) outputting said quantified association to a suitable output hardware, such as a CRT display terminal, a printer, a CD or DVD recorder, ZIP™ or JAZ™ 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 molecule or molecular complex prior to step a). In one embodiment, the method is for evaluating the ability of a chemical entity to associate with the binding pocket of a molecule or molecular complex.
In another embodiment, the method comprises the steps of:
In another embodiment, the invention provides a method of using a computer for evaluating the ability of a chemical entity to associate with the molecule or molecular complex, wherein said computer comprises a machine-readable data storage medium comprising a data storage material encoded with said structure coordinates defining said binding pocket and means for generating a three-dimensional graphical representation of the binding pocket, and wherein said method comprises the steps of:
The above method may further comprise the steps of:
Alternatively, the structure coordinates of the PLK binding pockets may be utilized in a method for identifying an agonist or antagonist of a molecule comprising a binding pocket of PLK. In certain embodiments, the method comprises steps of:
In one embodiment, step a) is performed using a graphical representation of the binding pocket or portion thereof of the molecule or molecular complex.
In one embodiment, the three-dimensional structure is displayed as a graphical representation.
In another embodiment, the method comprises the steps of:
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 PLK3 or PLK3-like binding pockets, motifs and domains.
Applicants' elucidation of binding pockets on PLK3 provides the necessary information for designing new chemical entities and compounds that may interact with PLK3 or PLK3-like substrate or ATP-binding pockets, in whole or in part.
Assays to determine if a compound binds to PLK3 are well known in the art and are exemplified below.
The design of chemical entities that bind to or inhibit PLK3 binding pockets according to this invention generally involves consideration of two factors. First, the entity must be capable of physically and structurally associating with parts or all of the PLK3 binding pockets. Non-covalent molecular interactions important in this association include hydrogen bonding, van der Waals interactions, hydrophobic interactions and electrostatic interactions.
Second, the entity must be able to assume a conformation that allows it to associate with the PLK3 binding pockets directly. Although certain portions of the entity will not directly participate in these associations, those portions of the entity may still influence the overall conformation of the molecule. This, in turn, may have a significant impact on potency. Such conformational requirements include the overall three-dimensional structure and orientation of the chemical entity in relation to all or a portion of the binding pocket, or the spacing between functional groups of an entity comprising several chemical entities that directly interact with the PLK3 or PLK3-like binding pockets.
The potential inhibitory or binding effect of a chemical entity on PLK3 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 PLK3 binding pockets, testing of the entity is obviated. However, if computer modeling indicates a strong interaction, the compound may then be synthesized and tested for its ability to bind to a PLK3 binding pocket. This may be achieved by testing the ability of the molecule to inhibit PLK3 using the assays described in Example 9. In this manner, synthesis of inoperative compounds may be avoided.
A potential inhibitor of a PLK3 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 PLK3 binding pockets.
One skilled in the art may use one of several methods to screen chemical entities or fragments for their ability to associate with a PLK3 binding pocket. This process may begin by visual inspection of, for example, a PLK3 binding pocket on the computer screen based on the PLK3 structure coordinates in
Specialized computer programs may also assist in the process of selecting fragments or chemical entities. These include:
Once suitable chemical entities or fragments have been selected, they can be assembled into a single compound or complex of compounds. 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 PLK3. This would be followed by manual model building using software such as QUANTA® and DISCOVERY STUDIO® [Molecular Simulations, Inc., San Diego, Calif.; now part of Accelrys, San Diego, Calif.], SYBYL® [Tripos Associates, St. Louis, Mo.] or MAESTRO [Schrödinger, Portland, Oreg. 97204, USA], or OPENEYE [Copyright © 1997-2006, OpenEye Scientific Software, SanteFe, N. Mex. 87508, USA].
Useful programs to aid one of skill in the art in building an inhibitor of a PLK3 binding pocket in a step-wise fashion, including one fragment or chemical entity at a time, include:
Instead of proceeding to build an inhibitor of a PLK3 binding pocket in a step-wise fashion one fragment or chemical entity at a time as described above, inhibitory or other PLK 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:
Other molecular modeling techniques may also be employed in accordance with this invention [see, e.g., N. C. Cohen et al., “Molecular Modeling Software and Methods for Medicinal Chemistry, J. Med. Chem., 33, pp. 883-894 (1990); see also, M. A. Navia and M. A. Murcko, “The Use of Structural Information in Drug Design”, Current Opinions in Structural Biology, 2, pp. 202-210 (1992); L. M. Balbes et al., “A Perspective of Modern Methods in Computer-Aided Drug Design”, Reviews in Computational Chemistry, Vol. 5, K. B. Lipkowitz and D. B. Boyd, Eds., VCH, New York, pp. 337-380 (1994); see also, W. C. Guida, “Software For Structure-Based Drug Design”, Curr. Opin. Struct. Biology, 4, pp. 777-781 (1994)].
Once a chemical entity has been designed or selected by the above methods, the efficiency with which that chemical entity may bind to a PLK3 binding pocket may be tested and optimized by computational evaluation. For example, an effective PLK3 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 PLK3 binding pocket inhibitors should preferably be designed with a deformation energy of binding of not greater than about 10 kcal/mole, more preferably, not greater than 7 kcal/mole. PLK3 binding pocket inhibitors may interact with the binding pocket in more than one conformation that is similar in overall binding energy. In those cases, the deformation energy of binding is taken to be the difference between the energy of the free entity and the average energy of the conformations observed when the inhibitor binds to the protein.
An entity designed or selected as binding to a PLK3 binding pocket may be further computationally optimized so that in its bound state it would preferably lack repulsive electrostatic interaction with the target enzyme and with the surrounding water molecules. Such non-complementary electrostatic interactions include repulsive charge-charge, dipole-dipole and charge-dipole interactions.
Specific computer software is available in the art to evaluate compound deformation energy and electrostatic interactions. Examples of programs designed for such uses include: Gaussian 94, revision C [M. J. Frisch, Gaussian, Inc., Pittsburgh, Pa. ©1995]; AMBER, version 4.1 [P. A. Kollman, University of California at San Francisco, ©1995]; QUANTA®/CHARMm® [Accelrys, San Diego, Calif.]; Insight II®/Discovery Studio® [Accelrys, San Diego, Calif. ©2001, 2002]; DelPhi [Accelrys, San Diego, Calif. ©2001, 2002]; and AMSOL [Quantum Chemistry Program Exchange, Indiana University]. These programs may be implemented, for instance, using a Silicon Graphics® workstation such as an INDIGO2 with “IMPACT™” graphics. Other hardware systems and software packages will be known to those skilled in the art.
Another approach enabled by this invention, is the computational screening of small molecule databases for chemical entities or compounds that can bind in whole, or in part, to a PLK 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 [E. C. Meng et al., J. Comp. Chem., 13, pp. 505-524 (1992)].
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 compound by determining and evaluating the three-dimensional structures of successive sets of protein/compound complexes.
According to another embodiment, the invention provides compounds which associate with a PLK3 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 compound by determining and evaluating the three-dimensional structures of successive sets of protein/compound 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. Alternatively, a pre-formed protein crystal is soaked in the presence of an inhibitor, thereby forming a protein/compound complex and obviating the need to crystallize each individual protein/compound complex.
In one embodiment, this invention provides for a method of designing a compound which binds to a catalytic domain of a Polo-like Kinase 3 (PLK3) protein comprising a binding pocket, wherein said catalytic domain is characterized by:
In another embodiment, the invention provides for a method of identifying a potential inhibitor of a Polo-like Kinase 3 (PLK3) protein comprising a binding pocket, wherein said method comprising the steps of:
In another embodiment, this invention provides for a method of identifying a candidate inhibitor of a molecule or molecular complex comprising a binding pocket or domain selected from the group consisting of:
In one embodiment, this invention provides for a method for identifying a candidate inhibitor that interacts with a binding site of a human Polo-Like Kinase 3 (PLK3) protein or a homologue thereof, comprising the steps of:
In another embodiment, this invention provides the methods of identifying above, wherein the binding site of said human Polo-Like Kinase 3 (PLK3) protein or said homologue thereof determined in step (d) comprises the structure coordinates according to
In another embodiment, this invention provides the methods of identifying above, wherein the binding site of said human Polo-Like Kinase 3 (PLK3) protein or said homologue thereof determined in step (d) comprises the structure coordinates according to
In another embodiment, this invention provides the methods of identifying above, wherein the binding site of said human Polo-Like Kinase 3 (PLK3) protein or said homologue thereof determined in step (d) comprises the structure coordinates according to
In another embodiment, this invention provides the methods of identifying above, wherein the binding site of said human Polo-Like Kinase 3 (PLK3) protein or said homologue thereof determined in step (d) comprises the structure coordinates according to
One embodiment of this invention provides for a method of identifying compounds that bind Polo-like Kinase 3 (PLK3) protein comprising:
One embodiment of this invention provides for a method of identifying compounds that bind Polo-like Kinase 3 (PLK3) protein selectively comprising:
Another embodiment of this invention provides for a method of identifying compounds that bind Polo-like Kinase 3 (PLK3) protein comprising:
Another embodiment of this invention provides for a method of identifying compounds that bind Polo-like Kinase 3 (PLK3) protein selectively comprising:
Structure Determination of Other Molecules
The structure coordinates set forth in
According to an alternate embodiment, the machine-readable data storage medium comprises a data storage material encoded with a first set of machine readable data which comprises the Fourier transform of at least a portion of the structure coordinates set forth in
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, wherein said computer comprises:
For example, the Fourier transform of at least a portion of the structure coordinates set forth in
Therefore, in another embodiment this invention provides a method of utilizing molecular replacement to obtain structural information about a molecule or molecular complex whose structure is unknown comprising the steps of:
In one embodiment, the method is performed using a computer. In another embodiment, the molecule is selected from the group consisting of PLK3 and PLK3 homologues. In another embodiment, the molecule is a PLK3 molecular complex or homologue thereof.
By using molecular replacement, all or part of the structure coordinates of the PLK3 as provided by this invention (and set forth in
Molecular replacement provides an accurate estimation of the phases for an unknown structure. Phases are a factor in equations used to solve crystal structures that can not be determined directly. Obtaining accurate values for the phases, by methods other than molecular replacement, is a time-consuming process that involves iterative cycles of approximations and refinements and greatly hinders the solution of crystal structures. However, when the crystal structure of a protein containing at least a homologous portion has been solved, the phases from the known structure provide a satisfactory estimate of the phases for the unknown structure.
Thus, this method involves generating a preliminary model of a molecule or molecular complex whose structure coordinates are unknown, by orienting and positioning the relevant portion of the PLK3 according to
The structure of any portion of any crystallized molecule or molecular complex that is sufficiently homologous to any portion of the PLK3 can be resolved by this method.
In a preferred embodiment, the method of molecular replacement is utilized to obtain structural information about a PLK3 homologue. The structure coordinates of PLK3 as provided by this invention are particularly useful in solving the structure of PLK3 complexes that are bound by ligands, substrates and inhibitors.
Furthermore, the structure coordinates of PLK3 as provided by this invention are useful in solving the structure of PLK3 proteins that have amino acid substitutions, additions and/or deletions (referred to collectively as “PLK3 mutants”, as compared to naturally occurring PLK3). These PLK3 mutants may optionally be crystallized in co-complex with a chemical entity, such as a non-hydrolyzable ATP analog or a suicide substrate. The crystal structures of a series of such complexes may then be solved by molecular replacement and compared with that of wild-type PLK3. 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 PLK3 and a chemical entity or compound.
The structure coordinates are also particularly useful in solving the structure of crystals of PLK3 or PLK3 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 PLK3 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 PLK3 inhibition activity.
All of the complexes referred to above may be studied using well-known X-ray diffraction techniques and may be refined versus 1.5-3.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)), CNS (Brunger et al., Acta Crystallogr. D. Biol. Crystallogr., 54, pp. 905-921, (1998)) or CNX (Accelrys, ©2000, 2001). This information may thus be used to optimize known PLK inhibitors, and more importantly, to design new PLK inhibitors.
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.
Human cDNA encoding the PLK3 kinase domain (amino acids 48-332) was obtained by PCR from cDNA library using the following gene specific oligonucleotides; 5′CAT ATG GAC CCC GGG CGC CTC ATC ACG G 3′ (pBev10-TOPO Forward primer D48; SEQ ID NO:2) and 5′CTA TTG GGA CTG TCA CGC AGC TGC TG 3′ (pBev10-TOPO Reverse Primer P332; SEQ ID NO:3). The PCR product was gel purified on a 0.8% agarose gel, excised from the gel and solubilized using the QIAGEN® QIAEX® II Gel Extraction Kit following the manufacturer's instructions. The PCR product was resuspended in TE pH8.0 buffer (30 ul). This fusion product was subcloned downstream of the polyhedrin promoter in the baculovirus donor vector, pBEV10TOPO, using the BamHI and EcoRI sites. The vector incorporated an N-terminal hexa-histidine purification tag and thrombin cleavage site.
The expression of PLK3 was carried out using standard procedures known in the art.
pBEV10TOPO is a Bac-to-Bac compatible pBEV10 vector TOPO® adapted by Invitrogen Corporation via Invitrogen's Custom TOPO® Cloning Adaptation Service. Recombinant virus was generated according to the manufacturer's recommendations. These initially transfected Spodoptera frugiperda (Sf9) cells were tested for the expression of PLK3 (D48-P332) protein by loading a crude extract of the transfected insect cells onto an SDS-PAGE (sodium dodecyl sulfate polyacrylamide gel electrophoresis) gel and immunoblot analysis using an anti-His (Sigma) antibody. Upon confirmation of the expression of the PLK3 (D48-P332) protein the virus was further amplified two times to obtain high titer stocks and used for optimisation of expression studies in Hi5 and Sf9 insect cells (1.5×106 cells/ml) at 27° C. with shaking at 100 rpm, using defined volumes of virus. After infection, cells were harvested at regular intervals of 24, 48 and 72 hours and optimum expression was determined by SDS-PAGE gel and immunoblot analysis using an anti-His (Sigma) antibody. Large-scale cultivation/expression was conducted with Hi5 insect cells (1.5×106 cells/ml) using a multiplicity of infection (M.O.I.) of between 2 and 5 of recombinant PLK3 (D48-P332) virus particles/cell, incubated at 27° C., 100 rpm and harvested 48 hours post-infection.
The PLK3 (D48-P332) cell pellets were resuspended into Lysis buffer (50 mM Hepes, pH 8.0, 300 mM NaCl, 5 mM β-mercaptoethanol, 10% Glycerol (v/v), 0.25% Tween, 5 mM Imidazole, 0.5 mM Sodium orthovanadate, 50 mM sodium fluoride, 10 mM β-glycerophosphate, Protease inhibitor Cocktail I (Roche Diagnostics GmbH, Germany) and Phosphatase Inhibitor) and disrupted by dounce homogenization, on ice. Resuspended cells were further mechanically disrupted using a MICROFLUIDIZER® (MICROFLUIDICS™, Newton, Mass.). The cell debris was removed by centrifugation (21,000 rpm, 15 min at 4° C.) and the supernatant incubated for 2.5 hours at 4° C. with pre-equilibrated Nickel-NTA metal affinity resin. The Ni-NTA resin was collected by centrifugation (1000 g, 4 min) and the non-specifically bound protein was removed by washing with 30× bead volume of Lysis buffer. The PLK3 (D48-P332) protein was eluted 3 times with Lysis buffer containing 200 mM Imidazole and a final elution was carried out using Lysis buffer containing 500 mM Imidazole. The N-terminal hexa-histidine tagged protein was then purified by size-exclusion on a Superdex 200(26/60) column (Amersham Biotech, Sweden) pre-equilibrated in Gel Filtration Buffer (50 mM Hepes, pH 8.0, 300 mM NaCl, 2 mM DTT and 10% Glycerol). The N-terminal hexa-histidine tag was removed by a 4° C. overnight incubation using 30 Units thrombin (Sigma) per mg of eluted protein and successful cleavage analysed by SDS-PAGE gel. The cleaved protein was then isolated by further purification using cation exchange chromatography. Protein post size exclusion was diluted with Buffer A (20 mM Hepes, pH 7.3, 10% Glycerol and 2 mM DTT) to reach a final salt concentration of 50 mM. Protein was then applied to a 10 ml MUSTANG™ S module (Pall Corporation) pre-equilibrated with Buffer A at a flow rate of 5 ml/min. Unbound protein was washed extensively with 10CV (column volumes) of Buffer. Bound protein was eluted using a salt gradient of 0-100% in 25 minutes at a flow rate of 3 ml/min. The resultant protein fractions were analysed by SDS-PAGE gel and fractions containing purified PLK3 (D48-P332) protein (amino acids 48-322 of SEQ ID NO:1) were pooled accordingly. The purified protein sample was dialysed against 25 mM Tris pH8.0 containing 200 mM NaCl, 10% glycerol and 2 mM DTT at 4° C. and concentrated to 10 mg/ml for crystallization. The crystallisable PLK3 (D48-P332) protein (amino acids 48-322 of SEQ ID NO:1) was found to elute between 49-54% NaCl.
Activity of the PLK3 constructs was determined using a standard coupled enzyme system (Fox et al., Protein Sci., 7, pp. 2249 (1998)). Reactions were carried out in a solution containing 100 mM HEPES (pH 7.5), 10 mM MgCl2, 25 mM NaCl, 2 mM DTT and 3% DMSO. Final peptidic substrate (full length Myelin Basic Protein (MBP), Vertex Pharmaceuticals Inc., Cambridge, Mass.) concentration in the assay was 15 μM. Reactions were carried out at 30° C. in the presence of 500 nM PLK3 construct and a titration of ATP (Sigma Chemicals, St Louis, Mo.) at final assay concentrations spanning 0 to 500 μM. Final concentrations of the components of the coupled enzyme system were 2.5 mM phosphoenolpyruvate, 300 μM NADH, 60 μg/ml pyruvate kinase and 20 μg/ml lactate dehydrogenase. An assay stock buffer solution was prepared containing all of the reagents listed above with the exception of ATP and DMSO. The assay stock buffer solution (60 μl) was incubated in a 96 well plate with 2 μl DMSO. The reaction was initiated by the addition of 5 μl of ATP (final assay concentrations spanning 0 to 500 μM). Rates of reaction were obtained using a Molecular Devices Spectramax plate reader (Sunnyvale, Calif.) over 10 min at 30° C. The ATP Km and Vmax values were determined from the rate data as a function of ATP concentration using computerized nonlinear regression (PRISM® 4.0a, Graphpad Software, San Diego, Calif.).
Crystals of PLK3-inhibitor complex crystals were formed by co-crystallizing the protein with adenosine. Adenosine was added to the PLK3 protein solution immediately after the final protein concentration step (Example 1) and immediately prior to setting up the crystallization drop. The PLK3-adenosine co-complex was prepared by incubating protein on ice with 3 mM adenosine for 30 minutes. Specifically, 5.72 μl protein stock solution and 1.2 μl of adenosine at a concentration of 50 mM was added to 14.28 μl of PLK3 at a concentration of 11.7 mg/ml and left on ice for 30 minutes prior to setting up the hanging drops.
Crystallization of PLK3 was carried out using the hanging drop vapor diffusion technique. The PLK3 formed rod-like crystals over a reservoir containing 900-1100 mM sodium di-hydrogen phosphate/100-900 mM di-potassium hydrogen phosphate (i.e., the concentration of sodium di-hydrogen phosphate plus the concentration of di-potassium hydrogen phosphate is 2M), 100 mM MES pH6.6 and 10 mM DTT. The crystallization droplet contained 0.50 μl of 8 mg ml-1 protein solution and 0.50 μl of reservoir solution. Crystals formed in approximately 72 hours.
The formed crystals were transferred to a reservoir solution containing 25% glycerol. After soaking the crystals in 25% glycerol for less than 2 minutes, the crystals were scooped up with a cryo-loop, frozen in liquid nitrogen and stored for data collection.
The PLK3—adenosine structure was solved by molecular replacement as describe below using X-ray diffraction data collected at SRS Daresbury, station 14.1. A summary of data collection is given in Table 2. The diffraction images were processed with the program MOSFLM [A. G. Leslie, Acta Cryst. D, 55, pp. 1696-1702 (1999)] and the data was scaled using SCALA [Collaborative Computational Project, N., Acta Cryst. D, 50, pp. 760-763 (1994)].
The data statistics, unit cell parameters and spacegroup of the PLK3—adenosine crystal structure is given in Table 2. The starting phases for the PLK3 complexes were obtained by molecular replacement using as a search model the coordinates of a PLK3 homology model based on Wee1 (Protein Data Bank (PDB) code 1X8B) in the program AMoRe [J. Navaza, Acta. Cryst. A, 50, pp. 157-163 (1994)]. The asymmetric unit contained a single PLK3 complex. Multiple rounds of rebuilding with QUANTA® [Molecular Simulations, Inc., San Diego, Calif.; now part of Accelrys Inc.] and initial refinement with CNX [Accelrys Inc., San Diego, Calif. ©2000] resulted in a model that included residues 52 to 327 of PLK3. The model after initial refinement had a crystallographic R-factor of 30.2% and R-free of 33.2% (see initial refinement in Table 3 for details for data 20-2.0 Å). Further refinement was performed with BUSTER-TNT (Global Phasing Ltd., Cambridge, UK) on data from 13-2.1 Å with I/σ(I) of 10.2/1.9, and Rmerge*(%) (overall/outer shell) of 5.0/40.9 and resulted in a final model that included residues 52 to 332 of PLK3. Information on the further refinement is also given in Table 3.
PLK3 has the typical bi-lobal catalytic kinase fold or structural domain [S. K. Hanks, et al., Science, 241, pp. 42-52 (1988); Hanks, S. K. and A. M. Quinn, Meth. Enzymol., 200, pp. 38-62 (1991)] with a β-strand sub-domain (residues 52-137) at the N-terminal end and an α-helical sub-domain at the C-terminal end (residues 146-327). The ATP-binding pocket is at the interface of the α-helical and β-strand sub-domains, and is bordered by the glycine rich loop and the hinge.
Comparison with other kinases such as Itk, p38 and cdk3 revealed that the structure of PLK3 resembles closely the substrate-bound, activated, form of a kinase. The overall topology of the kinase domain is similar to other serine/threonine kinases, particularly Wee1 and PDK1, and distinct from tyrosine kinases such as Itk, MKK2 (Table 1).
The coordinates of
The atomic coordinates in
One skilled in the art would therefore be able to use the X-ray crystallography data to determine possible antigenic sites in the PLK3 kinase domain. Possible antigenic sites are exposed, small and mobile regions on the kinase surface which have atomic B-factors of greater than 80 Å2 in
This process may be aided by using a computer comprising a machine-readable data storage medium encoded with a set of machine-executable instructions, wherein the recorded instructions are capable of displaying a three-dimensional representation of the PLK3 or a portion thereof.
Compounds were screened for their ability to inhibit PLK3 using a radioactive-phosphate incorporation assay. Assays were carried out in a mixture of 25 mM HEPES (pH 7.5), 10 mM MgCl2, and 1 mM DTT. Final substrate concentrations were 75 μM [γ-33P]ATP (60 mCi 33P ATP/mmol ATP, Amersham Pharmacia Biotech/Sigma Chemicals) and 110 μM peptide (SAM68 protein Δ332-443). Assays were carried out at 25° C. in the presence of 5 nM PLK3 (S38-A340). An assay stock buffer solution was prepared containing all of the reagents listed above, with the exception of ATP and the test compound of interest. 30 μL of the stock solution was placed in a 96 well plate followed by addition of 2 μL of DMSO stock containing serial dilutions of the test compound (typically starting from a final concentration of 10μM with 2-fold serial dilutions) in duplicate (final DMSO concentration 5%). The plate was pre-incubated for 10 minutes at 25° C. and the reaction initiated by addition of 8 μL [γ-33P]ATP (final concentration 75 μM).
The reaction was stopped after 60 minutes by the addition of 100 μL 0.14M phosphoric acid. A multiscreen phosphocellulose filter 96-well plate (Millipore, Cat no. MAPHN0B50) was pretreated with 100 μL 0.2M phosphoric acid prior to the addition of 125 μL of the stopped assay mixture. The plate was washed with 4×200 μL 0.2M phosphoric acid. After drying, 100 μL Optiphase ‘SuperMix’ liquid scintillation cocktail (Perkin Elmer) was added to the well prior to scintillation counting (1450 Microbeta Liquid Scintillation Counter, Wallac).
After removing mean background values for all of the data points, Ki(app) data were calculated from non-linear regression analysis of the initial rate data using the Prism software package (GraphPad PRISM® version 3.0cx for Macintosh, GraphPad Software, San Diego Calif., USA).
Compounds were screened for their ability to inhibit PLK1 using a radioactive-phosphate incorporation assay. Assays were carried out in a mixture of 25 mM HEPES (pH 7.5), 10 mM MgCl2, and 1 mM DTT. Final substrate concentrations were 50 μM [γ-33P]ATP (136 mCi 33P ATP/mmol ATP, Amersham Pharmacia Biotech/Sigma Chemicals) and 10 μM peptide (SAM68 protein Δ332-443). Assays were carried out at 25° C. in the presence of 15 nM PLK1 (A20-K338). An assay stock buffer solution was prepared containing all of the reagents listed above, with the exception of ATP and the test compound of interest. 30 μL of the stock solution was placed in a 96 well plate followed by addition of 2 μL of DMSO stock containing serial dilutions of the test compound (typically starting from a final concentration of 10 μM with 2-fold serial dilutions) in duplicate (final DMSO concentration 5%). The plate was pre-incubated for 10 minutes at 25° C. and the reaction initiated by addition of 8 μL [γ-33P]ATP (final concentration 50 μM).
The reaction was stopped after 60 minutes by the addition of 100 μL 0.14M phosphoric acid. A multiscreen phosphocellulose filter 96-well plate (Millipore, Cat no. MAPHN0B50) was pretreated with 100 μL 0.2M phosphoric acid prior to the addition of 125 μL of the stopped assay mixture. The plate was washed with 4×200 μL 0.2M phosphoric acid. After drying, 100 μL Optiphase ‘SuperMix’ liquid scintillation cocktail (Perkin Elmer) was added to the well prior to scintillation counting (1450 Microbeta Liquid Scintillation Counter, Wallac).
After removing mean background values for all of the data points, Ki(app) data were calculated from non-linear regression analysis of the initial rate data using the Prism software package (GraphPad PRISM® version 3.0cx for Macintosh, GraphPad Software, San Diego Calif., USA).
Compounds were screened for their ability to inhibit PLK1 using a radioactive-phosphate incorporation assay. Assays were carried out in a mixture of 25 mM HEPES (pH 7.5), 10 mM MgCl2, 0.1% BSA, and 2 mM DTT. Final substrate concentrations were 150 μM [γ-33P]ATP (115 mCi 33P ATP/mmol ATP, Amersham Pharmacia Biotech/Sigma Chemicals) and 300 μM peptide (KKKISDELMDATFADQEAK) (SEQ ID NO:4. Assays were carried out at 25° C. in the presence of 4 nM PLK1. An assay stock buffer solution was prepared containing all of the reagents listed above, with the exception of ATP and the test compound of interest. 30 μL of the stock solution was placed in a 96 well plate followed by addition of 2 μL of DMSO stock containing serial dilutions of the test compound (typically starting from a final concentration of 10 μM with 2-fold serial dilutions) in duplicate (final DMSO concentration 5%). The plate was pre-incubated for 10 minutes at 25° C. and the reaction initiated by addition of 8 μL [γ-33P]ATP (final concentration 150 μM).
The reaction was stopped after 90 minutes by the addition of 100 μL 0.14M phosphoric acid. A multiscreen phosphocellulose filter 96-well plate (Millipore, Cat no. MAPHN0B50) was pretreated with 100 μL 0.2M phosphoric acid prior to the addition of 125 μL of the stopped assay mixture. The plate was washed with 4×200 μL 0.2M phosphoric acid. After drying, 100 μL Optiphase ‘SuperMix’ liquid scintillation cocktail (Perkin Elmer) was added to the well prior to scintillation counting (1450 Microbeta Liquid Scintillation Counter, Wallac).
After removing mean background values for all of the data points, Ki(app) data were calculated from non-linear regression analysis of the initial rate data using the Prism software package (GraphPad PRISM® version 3.0cx for Macintosh, GraphPad Software, San Diego Calif., USA).
Compounds were screened for their ability to inhibit PLK2 using a radioactive-phosphate incorporation assay. Assays were carried out in a mixture of 25 mM HEPES (pH 7.5), 10 mM MgCl2, 0.1% BSA, and 2 mM DTT. Final substrate concentrations were 200 μM [γ-33P]ATP (57 mCi 33P ATP/mmol ATP, Amersham Pharmacia Biotech/Sigma Chemicals) and 300 μM peptide (KKKISDELMDATFADQEAK) (SEQ ID NO:4). Assays were carried out at 25° C. in the presence of 25 nM PLK2. An assay stock buffer solution was prepared containing all of the reagents listed above, with the exception of ATP and the test compound of interest. 30 μL of the stock solution was placed in a 96 well plate followed by addition of 2 μL of DMSO stock containing serial dilutions of the test compound (typically starting from a final concentration of 10 μM with 2-fold serial dilutions) in duplicate (final DMSO concentration 5%). The plate was pre-incubated for 10 minutes at 25° C. and the reaction initiated by addition of 8 μL [γ-33P]ATP (final concentration 200 μM).
The reaction was stopped after 90 minutes by the addition of 100 μL 0.14M phosphoric acid. A multiscreen phosphocellulose filter 96-well plate (Millipore, Cat no. MAPHN0B50) was pretreated with 100 μL 0.2M phosphoric acid prior to the addition of 125 μL of the stopped assay mixture. The plate was washed with 4×200 μL 0.2M phosphoric acid. After drying, 100 μL Optiphase ‘SuperMix’ liquid scintillation cocktail (Perkin Elmer) was added to the well prior to scintillation counting (1450 Microbeta Liquid Scintillation Counter, Wallac).
After removing mean background values for all of the data points, Ki(app) data were calculated from non-linear regression analysis of the initial rate data using the Prism software package (GraphPad PRISM® version 3.0cx for Macintosh, GraphPad Software, San Diego Calif., USA).
Compounds are screened for their ability to inhibit PLK4 using a radioactive-phosphate incorporation assay. Assays are carried out in a mixture of 8 mM MOPS (pH 7.5), 10 mM MgCl2, 0.1% BSA and 2 mM DTT. Final substrate concentrations are 15 μM [γ-33P]ATP (227 mCi 33P ATP/mmol ATP, Amersham Pharmacia Biotech/Sigma Chemicals) and 300 μM peptide (KKKMDATFADQ) (SEQ ID NO:5). Assays are carried out at 25° C. in the presence of 25 nM PLK4. An assay stock buffer solution is prepared containing all of the reagents listed above, with the exception of ATP and the test compound of interest. 30 μL of the stock solution is placed in a 96 well plate followed by addition of 2 μL of DMSO stock containing serial dilutions of the test compound (typically starting from a final concentration of 10 μM with 2-fold serial dilutions) in duplicate (final DMSO concentration 5%). The plate is pre-incubated for 10 minutes at 25° C. and the reaction initiated by addition of 8 μL [γ-33P]ATP (final concentration 15 μM).
The reaction is stopped after 180 minutes by the addition of 100 μL 0.14M phosphoric acid. A multiscreen phosphocellulose filter 96-well plate (Millipore, Cat no. MAPHN0B50) is pretreated with 100 μL 0.2M phosphoric acid prior to the addition of 125 μL of the stopped assay mixture. The plate is washed with 4×200 μL 0.2M phosphoric acid. After drying, 100 μL Optiphase ‘SuperMix’ liquid scintillation cocktail (Perkin Elmer) is added to the well prior to scintillation counting (1450 Microbeta Liquid Scintillation Counter, Wallac).
After removing mean background values for all of the data points, Ki(app) data are calculated from non-linear regression analysis of the initial rate data using the Prism software package (GraphPad PRISM® version 3.0cx for Macintosh, GraphPad Software, San Diego Calif., USA).
Structure Refinement
The present application claims the benefit under 35 U.S.C. §119 of U.S. Provisional patent application No. 60/813,427, filed Jun. 14, 2006, the entire contents of which are incorporated herein by reference.
Number | Name | Date | Kind |
---|---|---|---|
4886646 | Carter et al. | Dec 1989 | A |
5096676 | McPherson et al. | Mar 1992 | A |
5130105 | Carter et al. | Jul 1992 | A |
5221410 | Kushner et al. | Jun 1993 | A |
5400741 | DeTitta et al. | Mar 1995 | A |
5884230 | Srinivasan et al. | Mar 1999 | A |
7432260 | Wang et al. | Oct 2008 | B2 |
20030077681 | Cogswell | Apr 2003 | A1 |
20050261836 | Meng et al. | Nov 2005 | A1 |
20080201123 | Cosgrove | Aug 2008 | A1 |
Number | Date | Country |
---|---|---|
WO 03048340 | Jun 2003 | WO |
WO 2004024759 | Mar 2004 | WO |
WO 2005047526 | May 2005 | WO |
WO 2006053125 | May 2006 | WO |
Number | Date | Country | |
---|---|---|---|
20080293579 A1 | Nov 2008 | US |
Number | Date | Country | |
---|---|---|---|
60813427 | Jun 2006 | US |