This invention relates to methods of identifying and designing cell surface receptor inhibitors, particularly integrin inhibitors.
This research was supported by grant numbers DK48549 and HL70219 awarded by the National Institutes of Health (NIDDK and NHLBI).
Integrins are adhesion receptors that mediate vital bi-directional signals during morphogenesis, tissue remodeling and repair (reviewed in 2). Integrins are heterodimers formed by noncovalent association of an α and a β subunit, both type I membrane proteins with large extracellular segments. In mammals, eighteen α and eight β subunits assemble into 24 different receptors. Integrins depend on divalent-cations to bind their extracellular ligands. Although these ligands are structurally diverse, they all employ an acidic residue during integrin recognition. Specificity for a particular ligand is then determined by additional contacts with the integrin. High affinity binding of integrins to ligands is usually not constitutive, but is elicited in response to cell “activation” signals (so-called “inside-out” signaling) that alter the tertiary and quaternary structure of the extracellular region, making the integrin ligand-competent. Ligand binding in turn induces structural rearrangements in integrins that trigger “outside-in” signaling (reviewed in 4).
Integrins can be grouped into two classes based on the presence or absence of α˜180 amino acid A-type domain (αA or I domain; see 18). In the nine αA-containing integrins (αA-integrins), αA is the major ligand binding site. Thus, isolated αA binds directly and in a divalent-cation-dependent manner to physiologic ligands with affinity equal to that of the respective ligand-competent heterodimer.18 The structures of isolated αA domains in “liganded” (high-affinity) and “unliganded” (low-affinity) conformations have revealed how this domain interacts with ligands.19,30 A metal ion is coordinated at the ligand-binding interface of αA through a conserved five amino acid motif, the metal-ion-dependent-adhesion-site (MIDAS), and the metal coordination is completed by a glutamate from the ligand or, in its absence, by a water molecule.19,20,29 In αA-lacking integrins, we proposed that ligand recognition is supported by an αA-like domain (αA) present in all integrin β subunits.19 The a subunit is also believed to participate in ligand recognition; we have recently defined its binding site to the prototypical ligand RGD.23
The invention features methods for identifying compounds which bind and modulate integrins. The compounds identified by the invention can inhibit integrins in a manner referred to herein as “deadbolt inhibition.” The region involved in deadbolt inhibition includes the βA domain strand-F/α7 loop contacted by the CD loop of the β tail domain (βTD) of an unliganded integrin (e.g., αVβ3). The contact region between the βA domain and the βTD acts as a regulatable deadbolt to lock the βA domain in an inactive state by preventing strand-F/α7 loop movement associated with the activation-initiated inward movement of the α1 helix. The deadbolt contact is stabilized in place by an additional contact on the same side between the βTD and the hybrid domain and by an ionic bond through an ADMIDAS cation, linking the strand-F/α7 loop to the α1 helix of βA. The compounds (e.g., mimetics of this loop, of the ADMIDAS cation, or of the βTD/βA contact) identified by the methods of the invention act to stabilize the CD loop in locking the integrin in an inactivated state.
The invention features methods for evaluating potential of a compound to associate with a molecule or molecular complex comprising a non-ligand binding site of an integrin βA domain, the method comprising: (a) employing computational means to perform a fitting operation between the compound and β tail domain (βTD) contact region on the strand-F/α7 loop of an unliganded integrin (e.g., αVβ3); and (b) analyzing the results of the fitting operation to quantify said association potential. In certain embodiments, the compound mimics the interaction of a peptide, comprising the amino acid sequence C663VVRFQYYE671D672S673S674G675KSILYVVEEPEC687 or a fragment thereof, or K618KFDREPYMTENTCNR633YCRD or a fragment thereof, with the strand-F/α7 loop of a βA domain of the integrin.
In another aspect, the invention features a method for identifying a candidate selective modulator of the activity of an integrin, the method comprising: (a) modeling test compounds that fit spatially and preferentially into a βA domain non-ligand binding site of an integrin of interest using an atomic structural model of the integrin βA domain, wherein the atomic structural model is generated using amino acid sequence comprising C663VVRFQYYE671D672S673S674G675KSILYVVEEPEC687 or a fragment thereof, or K618KFDREPYMTENTCNR633YCRD or fragment thereof; (b) screening test compounds in a biological assay for integrin activation characterized by binding of a test compound to the βA domain non-ligand binding site of the integrin; and (c) identifying a test compound that selectively modulates the activity of the integrin, and optionally, (d) screening an identified test compound in a biological assay for its ability to prevent interaction of the βA domain and the βTD domain by binding of the identified test compound to the βA domain non-ligand binding site of the integrin, and (e) identifying the screened test compound as a compound capable of selectively modulating the activity of an integrin from the pool of test compounds. In various embodiments: the modulating comprises inhibiting the integrin activation and the modulating comprises inhibiting ligand binding to the integrin.
In another aspect, the invention features a method of identifying a candidate inhibitor of the activity of an integrin, the method comprising: (a) introducing into a suitable computer program information defining a non-ligand binding site of an integrin βA domain, the information comprising a conformation defined by the coordinate atoms as in Table 1 and Table 2, wherein the program displays the three-dimensional structure thereof; (b) creating a three dimensional structure of a test compound in the computer program; (c) superimposing the model of the test compound on the model of the non-ligand binding site of the integrin βA domain; and (c) assessing whether the test compound model fits spatially into the non-ligand binding site. In one embodiment, the atoms include those in Table 3 and Table 4.
In another embodiment, the invention features a method for identifying a candidate integrin modulating compound, the method comprising: (a) generating a three-dimensional structure of the βTD contact region of strand-F/α7 loop of an βA domain of a non-ligand bound integrin; (b) employing the three dimensional structure to design or select the candidate integrin modulating compound, and (c) identifying said candidate integrin modulating compound by the data obtained by steps (a) and (b). In one embodiment, the method also includes: (d) synthesizing the integrin modulating compound; and (e) determining the ability of the integrin inhibitor to bind to the integrin by contacting the modulating compound with the integrin.
In another aspect, the invention features a method for identifying a candidate integrin modulating compound, the method comprising (a) expressing recombinant integrin fragments containing βA domain and βTD domains and (b) employing the protein-protein interaction of these two domains to in screening assays to identify modulators of the βA domain and βTD domain interaction. The method can further include: (c) synthesizing the integrin modulating compound; (d) determining the ability of the integrin modulator to bind to the integrin by measuring the interaction of the modulating compound with the integrin; and (e) using the modulating compound as a basis for drug design.
The invention also features a method of inhibiting activation of an integrin the method comprising contacting a compound to with an integrin thereby locking the βA domain structure of the integrin into a non-activatable form. In various embodiments, the compound mimics an intrachain ligand in its interaction with the integrin, the intrachain ligand comprises the sequence of SEQ ID No:1, or a fragment thereof, or SEQ ID No. 2, or a fragment thereof; and the intrachain ligand is a member of the statin family.
In another aspect, the invention features a method of identifying an integrin modulator comprising: (a) selecting a potential inhibitor by performing rational drug design with the three-dimensional structural coordinates of Table 1 and Table 2, wherein selecting is performed in conjunction with computer modeling; (b) contacting the potential inhibitor with an integrin domain; and (c) detecting the ability of the potential inhibitor for inhibiting the integrin. In various embodiments, detecting the ability of the potential inhibitor for inhibiting the integrin in step (c) is performed using a ligand binding assay; detecting the ability of the potential inhibitor for inhibiting the integrin in step (c) is performed using a cellular-based assay; and the method further includes: (d) growing a supplemental crystal comprising a complex formed between the integrin domain and a first potential inhibitor from step (a), the supplemental crystal effectively diffracts X-rays for the atomic coordinates of the complex a resolution of greater than 4.0 Å; (e) determining the three-dimensional structure of the supplemental crystal; (f) selecting a second potential inhibitor by performing rational drug design with the three-dimensional structure determined for the supplemental crystal, wherein selecting is performed in conjunction with computer modeling;
(g) contacting the second potential inhibitor with the integrin domain; and (h) detecting the ability of the second potential inhibitor for inhibiting the integrin.
Also within the invention are compounds identified by the methods of the invention, with the proviso that the compound is not lovastatin, and pharmaceutical compositions comprising such compounds and a pharmaceutically acceptable carrier.
The invention also features a method for modulating, inhibiting or stimulating binding of ligands or associated proteins to integrins by modifying the interaction of integrin beta-A domain (βA) with the beta-tail domain (βTD). In certain embodiments the integrin is selected from the group consisting of: αVβ1, αVβ3, αVβ5, αVβ6, αVβ8, α3β1, α4β1, α5β1, α6β1, α6β4, α7β1, α9β1, α4β7, gp3b3a, α1β1, α2β1, α10β1, α11β1, LFA-1, MAC-1, or α150β95; the interaction of βA with βTD is investigated using either computational or biochemical or biophysical techniques; and either the βA, the βTD or both serve as structures on which the computational or biochemical or biophysical techniques are based.
The invention also features a method for discovering pharmacologically relevant substances including antibodies, small molecules, polypeptides, peptides, and peptide mimetics, which can perturb or stimulate the interaction of integrin beta-A domain (βA) with the beta-tail domain (βTD). In certain embodiments, the pharmacologically relevant substances modulate (e.g., inhibit) the interaction of the native integrin with its ligands or cofactors.
Definitions
As used herein, “integrin” and “integrin receptor” are used interchangeably. “Integrin” or “integrin receptor” refers to any of the many cell surface receptor proteins, also referred to as adhesion receptors which bind to extracellular matrix ligands or other cell adhesion protein ligands thereby mediating cell-cell and cell-matrix adhesion processes. The integrins are encoded by genes belonging to a gene superfamily and are typically composed of heterodimeric transmembrane glycoproteins containing α- and β-subunits. Integrin subfamilies contain a β-subunit combined with different α-subunits to form adhesion protein receptors with different specificities. The integrins are grouped into two classes, those containing the αA domain and those that do not contain the αA domain. Both classes have a βA domain.
“Compounds” refer to a chemical entity that can comprise a peptide or polypeptide, including antibodies, phage display antibodies, and their biologically active fragments, a small molecules, (e.g., chemically synthesized or of natural origin), or synthetic peptides or polypeptides (e.g., non-naturally occurring polypeptides, e.g., peptoids or peptidomimetics). The compound has the ability to bind the integrin at a non-ligand binding site of the βA domain. The “non-ligand binding site of an integrin βA domain” includes the site with which the CD loop of the β tail domain (βTD) of the β subunit that contacts the strand-F/α7 loop in the βA domain of the β subunit of an unliganded integrin. As such, ligand refers to the naturally occurring ligand that binds the integrin so that it can perform its physiological function or functions. This site is distinguishable from the ligand binding site. Compounds can include synthetic and naturally occurring mimetics of the CD loop or fragments of the CD loop of the β tail domain (βTD) that contacts the strand-F/α7 loop in the βA domain of an unliganded integrin, the α1/strand-A loop or fragments of that loop that contact the hybrid domain of an unliganded integrin, and the ADMIDAS coordination site in the βA domain of an unliganded integrin. A “mimetic” has structural similarity or has similar binding properties as the entity it is mimicking. Compounds can comprise or consist of the mimetic. According to the invention, “compounds” preferably refers to chemical entities as described above, selected from the group consisting of
“Modulation” refers to regulating or changing the activatability or ligand binding of an integrin. For example, a modulator can inhibit or promote the activation of the integrin or it can prevent or promote ligand binding to the integrin.
A “peptidomimetic” refers to a chemical variant of a polypeptide or a peptide in which the side chains of the polypeptide or peptide are substantially maintained in the variant, yet the chemical backbone of the peptidomimetic is altered relative to the polypeptide or peptide in at least one peptide bond.
A “peptoid” is an oligomer of N-substituted glycines. A peptoid can be synthesized from a variety of different N-alkylglycines that have side chains similar to amino acid side chains, (e.g., as described in Simon et al., (1992) PNAS 89:9367-9371). It can serve as a motif for the generation of chemically diverse libraries of novel molecules. As an alternative to natural polymers, it is a modular system that allows one to synthesize monomers in large amounts. The monomers have a wide variety of functional groups presented as side chains off of an oligomeric backbone, the linking chemistry is high yielding and amenable to automation. The linkage in a peptoid is resistant to hydrolytic enzymes such as proteases. Another advantage is that the monomers are achiral.
A “small molecule” is a molecule of less than 32 kDa, e.g., 0.5 kDa, 1 kDa, 5 kDa, 10, kDa, 15 kDa, 20 kDa, 25 kDa, 30 kDa, or 32 kDa.
Unless otherwise defined, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Although methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present invention, suitable methods and materials are described below. All publications, patent applications, patents, and other references mentioned herein are incorporated by reference in their entirety. In case of conflict, the present specification, including definitions, will control. In addition, the materials, methods, and examples are illustrative only and not intended to be limiting.
The details of one or more embodiments of the invention are set forth in the accompanying drawings and the description below. Other features, objects, and advantages of the invention will be apparent from the description and drawings, and from the claims.
The invention is based, in part, on the identification of a mechanism for “deadbolt inhibition” of integrins and features methods for identifying compounds which bind to integrins (e.g., αVβ3), particularly compounds that inhibit activation of the integrin and inhibit ligand binding to the integrin thus providing deadbolt inhibition. The structural information in
The methods of the invention entail identification of compounds (e.g., small molecules naturally occurring polypeptides, and non-naturally occurring peptides or polypeptides (e.g., peptoids, peptidomimetics, antibodies, and the like) having a particular structure. The methods rely on the use of precise structural information derived from x-ray crystallographic studies which include the βA domain of an integrin (e.g., αVβ3). These crystallographic data permit the identification of atoms in the compound (e.g., naturally occurring and non-naturally occurring polypeptides, small molecules, and the like) that are important for integrin binding and integrin modulation (e.g. inhibition). More importantly, these data define a three-dimensional array of the important contact atoms. Other molecules which include a portion in which the atoms have a three-dimensional arrangement similar to some or all of these contact atoms are likely to be capable of binding an integrin and acting as an integrin modulator (e.g., inhibitor). The methods can also employ isolated integrin βA domains defined on the basis of their X-ray crystallographic or other high-resolution molecular structural data (e.g., via solution NMR), or on the basis of their homology to the structure of integrin αVβ3 as isolated domain or as fusion partner (e.g., with the Fc domain of immunoglobulin), and use conventional inhibitor screens based on the inhibition or stimulation of interaction of the βTD domain with the βA domain.
A recent study, based on analysis of negatively-stained EM (electron micrograph) images of extracellular integrins in solution, found that Mn2+ and RGD (arginine-glycine-aspartate) elicit global changes in the integrin, causing it to switch from the bent to the knee-straightened form, regardless of whether its legs are chemically tied24. Straightening at the knees was proposed to be caused by an outward swing of the hybrid domain relative to βA, driven by a downward slide of the βA's C-terminal α7 helix. Since the bent form of extracellular or membrane-bound αVβ3 (both produced by locking the head to the legs through a disulfide bond introduced based on the αVβ3 atomic structure) was found inactive but could be activated by chemical reduction of disulfides24, it was proposed that this “switchblade” model underlies integrin activation. While appealing, this model does not explain a number of observations. First, inside-out activation is rapid (<1 sec) and reversible. The switchblade model does not explain how the extended structure refolds within seconds, especially within the crowded confines of the extracellular matrix. Second, movement at the knees alone is insufficient to extend the 20 Å cryoEM structure of the native inactive αIIbβ345. Third, analysis of mAb epitope maps of active an inactive native αIIbβ3 suggest that both can be in the compact (bent) form82. Fourth, activation and extension can be dissociated: extracellular α5β1 tied at its legs can assume the straight conformation when unlocked, yet it is inactive72. Furthermore, fixing the structure in its bent unliganded state by artificial disulphides is expected to keep it unliganded, as tertiary and quaternary changes of the integrin are needed and occur even in the protein crystal lattice to allow RGD ligand binding (see above). Thus it may be that flexion/extension at the knees is an outcome, rather than a cause, of activation which may contribute to post-ligand binding event(s).
The flexible CD loop of the βTD contacts strand-F/α7 loop (Ser674 of the former contacts Val332 in βA) in unliganded αVβ3 (Table 1)(
Significant progress has recently been made in elucidating the crystal structure of the integrin ectodomain in its unliganded and RGD-liganded states and the NMR structure of its cytoplasmic tails “feet”.22,23,44 The liganded structure helped define the basis of cation-dependent ligand binding specificity and helped in the development of models of the interaction of an integrin with more complex ligands. Comparison of the unliganded and liganded structures elucidated the tertiary changes leading to ligand binding and provided a minimalist view of the associated quaternary changes. This information resulted in the new “deadbolt inhibition” model for inside-out activation. Given the broad spectrum of diseases caused or exacerbated by dysregulated integrin function, this structural information provides a rich source for the rational development of specific therapeutic drugs targeting these receptors.
The αVβ3 structure has a curious analogy to the lovastatin binding site. Due to severe folding of the beta-3 chain, a membrane proximal loop of the βTD contacts the membrane distal F7 fold of the βA domain, at a position near where lovastatin contacts CD11a. On binding RGD-ligand, this contact is broken in a movement of 0.3 nm allowing movement within the βA domain. RGD-ligands can activate integrins. It is likely that this intrachain contact thus solves the problem of integrin activation for non-α-A domain integrins, explaining how an intracellular signal can communicate via 5 domains to the ligand-binding site, that is 23 nm away along the chain.
This mechanism suggests that the slot formed by the F1-F7 helices and the F7 pleat of the non-α A domains can be a general novel target, and will allow a new category of drugs to be rationally designed. These drugs will prevent integrin activation by locking the βA domain structure into a non-activatable form. The drugs would necessarily be completely different, and function by a completely different mechanism to those acting at the ligand binding site, and also, because they are mimics of an intrachain ligand, are also different from the lovastatin inhibition, that acts on extracellular activators of the integrin alpha chain. As the A domains have been successfully crystallized as independent domains, a rapid elucidation of the structure of these domains from as yet uncharacterized integrins can be made (i.e., by production of the βA domains as recombinant protein, and determination of their crystal structure), to allow rational drug design.
Starting points for rational drug design include: 1) the amino acid sequence of the βTD loops in β3 (C663VVRFQYYE671D672S673S674G675 KSILYVVEEPEC687: SEQ ID No. 1 especially the central amino acids 671-675; and K618KFDREPYMTENTCNR633YCRD: SEQ ID NO.2 especially arginine 633), 2) the homologous loops in the other relevant integrin β chains, 3) the high resolution NMR or X-ray crystallographic structures of the βA domains intact in the full-length integrin or isolated as a chain fragment containing the βA domain and variable fragments of the chain or of other appropriate fusion partners to allow rapid crystallization, and 4) structural homologues of the statin family of drugs. Screening for drug candidates in high-throughput assays can utilize the following: 1) assays based on ligand binding to integrins or recombinant ligand-binding integrin fragments (e.g., vitronectin, fibrinogen, or thrombospondin binding assays and β3A domains/βTD domain binding assays), 2) perturbation of the NMR signal, 3) perturbation at the isolated integrin receptor or βA domain-containing receptor fragment, or 4) via attachment at the cellular level.
Computer Modeling
The methods of the invention employ computer-based methods for identifying compounds having a desired structure. These computer-based methods fall into two broad classes: database methods and de novo design methods. Database methods fall in two main classes, those based on a compound (i.e., a ligand of a binding site alone) or those based on the three dimensional structure of the binding site. In the former approach, the compound of interest is compared to all compounds present in a database of chemical structures and compounds whose structure is in some way similar to the compound of interest are identified. In the latter approach, all compounds in a database are docked by appropriate computer software into the binding site, and their degree of fit is evaluated and ranked. The structures in the database are based on either experimental data, generated by NMR or x-ray crystallography, or modeled three-dimensional structures based on two-dimensional protein or DNA sequence data. In de novo design methods, models of compounds whose structure is in some way similar to the compound of interest are generated by a computer program using information derived from known structures (e.g., data generated by x-ray crystallography and/or theoretical rules). Such design methods can build a compound having a desired structure in either an atom-by-atom manner or by assembling stored small molecular fragments.
The success of both database and de novo methods in identifying compounds with activities similar to the compound of interest depends on the identification of the functionally relevant portion of the compound of interest. For drugs, the functionally relevant portion is referred to a pharmacophore. A pharmacophore then is an arrangement of structural features and functional groups important for biological activity.
Not all identified compounds having the desired pharmacophore will act as an integrin modulator. The actual activity can be finally determined only by measuring the activity of the compound in relevant biological assays. However, the methods of the invention are extremely valuable because it can be used to greatly reduce the number of compounds which must be tested to identify an actual mimetic.
Programs suitable for generating predicted three-dimensional structures from two-dimensional data include: Concord (Tripos Associated, St. Louis, Mo.), 3-D Builder (Chemical Design Ltd., Oxford, U.K.), Catalyst (Bio-CAD Corp., Mountain View, Calif.), Daylight (Abbott Laboratories, Abbott Park, Ill.).
Programs suitable for searching three-dimensional databases to identify molecules bearing a desired pharmacophore include: MACCS-3D and ISIS/3D (Molecular Design Ltd., San Leandro, Calif.), ChemDBS-3D (Chemical Design Ltd., Oxford, U.K.), Sybyl/3DB Unity (Tripos Associates, St. Louis, Mo.). Programs suitable for pharmacophore selection and design include: DISCO (Abbott Laboratories, Abbott Park, Ill.), Catalyst (Bio-CAD Corp., Mountain View, Calif.), and ChemDBS-3D (Chemical Design Ltd., Oxford, U.K.).
Databases of chemical structures are available from Cambridge Crystallographic Data Centre (Cambridge, U.K.) and Chemical Abstracts Service (Columbus, Ohio).
De novo design programs include Ludi (Biosym Technologies Inc., San Diego, Calif.) and Aladdin (Daylight Chemical Information Systems, Irvine Calif.), LEGEND (Nishibata, Y., Itai, A., Tetrahedron, 47, 8985 (1991))(Molecular Simultations, Burlington, Mass.), and LeapFrog (available from Tripos associates, St. Louis, Mo.).
Upon determination of the three-dimensional structure of an integrin, a potential modulator can be evaluated by any of several methods, alone or in combination. Such evaluation can utilize visual inspection of a three-dimensional representation of the binding site on the integrin, based on the x-ray coordinates of a crystal described herein, on a computer screen. Evaluation, or modeling, can be accomplished through the use of computer modeling techniques, hardware, and software known to those of ordinary skill in the art. This can additionally involve model building, model docking, or other analysis of protein-ligand interactions using software including, for example, QSC, GOLD (Jones et al., J. Mol. Biol., 245, 43-53, 1995), FlexX (Lengauer, Rarey, 1996), Autodock (Morris et.al., 1998), GLIDE, Modeler, or Sybyl, followed by energy minimization and molecular dynamics with standard molecular mechanics forcefields including, for example, CHARMM and AMBER. The three-dimensional structural information of an unliganded integrin, (e.g., the CD loop of the βTD contacting the strand-F/α7 loop of the βA domain in an unliganded integrin) can also be utilized in conjunction with computer modeling to generate computer models of other unliganded integrins. Computer models of unliganded integrin structures can be created using standard methods and techniques known to those of ordinary skill in the art, including software packages described herein.
Once the three-dimensional structure of a crystal, or solution structure via NMR, comprising an unliganded integrins determined, a potential non-ligand site binder (e.g., a binder that mimics the βTD binding of the strand-F/α7 loop of the βA domain) is examined through the use of computer modeling using a docking program such as QSC, GOLD, FlexX, or Autodock to identify potential non-ligand binding site binders to ascertain how well the shape and the chemical structure of the potential ligand will interact with the binding site. Computer programs can also be employed to estimate the attraction, repulsion, and steric hindrance of the two binding partners (i.e., the non-ligand binding site and the modulating binder). Generally tighter fit, lower steric hindrances, and greater attractive force between the potential ligand and the allosteric binding site are consistent with a tighter binding constant between the two. Furthermore, the more specificity in the design of a potential drug, the more likely that the drug will also not interact with other proteins. This will minimize potential side-effects due to unwanted interactions with other proteins.
A variety of methods are available to one skilled in the art for evaluating and virtually screening molecules or chemical fragments appropriate for associating with a protein, particularly an integrin. Such association can be in a variety of forms including, for example, steric interactions, van der Waals interactions, electrostatic interactions, solvation interactions, charge interactions, covalent bonding interactions, non-covalent bonding interactions (e.g., hydrogen-bonding interactions), entropically or enthalpically favorable interactions, and the like.
Numerous computer programs are available and suitable for rational drug design and the processes of computer modeling, model building, and computationally identifying, selecting and evaluating potential modulating compounds in the methods described herein. These include, for example, QSC (WO 01/98457), FlexX, Autodock, Glide, Accelrys' Discovery Studio, or Sybyl. Potential inhibitors can also be computationally designed “de novo” using such software packages as QSC (WO 01/98457), Accelrys' Discovery Studio, Sybyl, ISIS, ChemDraw, or Daylight. Compound deformation energy and electrostatic repulsion, can be evaluated using programs such as GAUSSIAN 92, AMBER, QUANTA/CHARMM, AND INSIGHT II/DISCOVER.
There are a number of ways to select moieties to fill individual binding pockets. These include QUANTA [Molecular Simulations, Inc., Burlington, Mass., 1992], SYBYL [Molecular Modeling Software, Tripos Associates, Inc., St. Louis, Mo., 1992], AMBER [S. J. Weiner, P. A. Kolhman, D. A. Case, U. C. Singh, C. Ghio, G. Alagona, and P. Weiner, J. Am. Chem. Soc., vol. 106, pp. 765-784 (1984)], or CHARMM [B. R. Brooks, R. E. Bruccoleri, B. D. Olafson, D. J. States, S Swaminathan, and M. Karplus, J. Comp. Chem. vol. 4, pp. 187-217 (1983)]. This modelling step may be followed by energy minimization with standard molecular mechanics forcefields such as CHARMM and AMBER. In addition, there are a number of more specialized computer programs to assist in the process of selecting the binding moieties of this invention.
These include:
1.) GRID (Goodford, P. J. A Computational Procedure for Determining Energetically Favorable Binding Sites on Biologically Important Macromolecules. J. Med. Chem., 28, pp. 849-857 (1985)). GRID is available from Oxford University, Oxford, UK.
2.) MCSS (Miranker, A.; Karplus, M. Functionality Maps of Binding Sites: A Multiple Copy Simultaneous Search Method. Proteins: Structure, Function and Genetics, 11, pp. 29-34 (1991)). MCSS is available from Molecular Simulations, Burlington, Mass.
3.) AUTODOCK (Goodsell, D. S.; Olsen, A. J. Automated Docking of Substrates to Proteins by Simmulated Annealing. PROTEINS: Structure, Function and Genetics, 8, pp. 195-202 (1990)). AUTODOCK is available from the Scripps Research Institute, La Jolla, Calif.
4.) DOCK (Kuntz, I. D.; Blaney, J. M.; Oatley, S. J.; Langridge, R.; Ferrin, T. E. A Geometric Approach to Macromolecule-Ligand Interactions. J. Mol. Biol., 161, pp. 269-288 (1982)). DOCK is available from the University of California, San Francisco, Calif.
5.) GOLD (Jones et al., J. Mol. Biol., 245, 43-53, 1995). GOLD is available from the Cambridge Crystallography Data Centre, Camdridge, UK.
6.) FlexX (T. Lengauer and M. Rarey, Computational Methods for Biomolecular Docking, Current Opinion in Structural Biology, Vol. 6, pp. 402-406, 1996). FlexX is available through Tripos Associated, St. Louis, Mo.
These computer evaluation and modeling techniques can be performed on any suitable hardware including for example, workstations available from Silicon Graphics, Sun Microsystems, and the like. These techniques, methods, hardware and software packages are representative and are not intended to be comprehensive listing. Other modeling techniques known in the art can also be employed in accordance with this invention. See for example, QSC (WO 01/98457), FlexX, Autodock, Glide, Accelrys' Discovery Studio, or Sybyl and software identified at various internet sites (e.g., netsci.org/Resources/Software/Modeling/CADD/
ch.cam.ac.uk/SGTL/software.html
cmm.info.nih.gov/modeling/universal_software.html
dasher.wustl.edu/tinker/
zeus.polsl.gliwice.pl/˜nikodem//linux4chemistry.html
nyu.edu/pages/mathmol/software.html
msi.umn.edu/user_support/software/MolecularModeling.html
us.expasy.org/
sisweb.com/software/model.htm).
A potential integrin modulator is selected by performing rational drug design with the three-dimensional structure (or structures) determined for the site on the βA domain, at which the βTD is in contact, described herein, in conjunction with or solely by computer modeling and methods described above. The potential modulator is then obtained from commercial sources or is synthesized from readily available starting materials using standard synthetic techniques and methodologies known to those of ordinary skill in the art. The potential inhibitor is then assayed to determine its ability to modulate the target (e.g., integrin, e.g., αVβ3) and/or integrin pathway.
A potential inhibitor can also be selected by screening a library of compounds (e.g., a combinatorial library, e.g., a mass-coded combinatorial library). The library of compounds can be screened by affinity screening in which members with the greatest affinity to a particular integrin at the new non-ligand binding site can be selected.
Once suitable binding moieties have been selected, they can be assembled into a single modulating binder. This assembly may be accomplished by connecting the various moieties to a central scaffold. The assembly process may, for example, be done by visual inspection followed by manual model building, again using software such as Quanta or Sybyl. A number of other programs may also be used to help select ways to connect the various moieties. These include: CAVEAT (Bartlett, P. A.; Shea, G. T.; Telfer, S. J.; Waterman, S. CAVEAT: A Program to Facilitate the Structure-Derived Design of Biologically Active Molecules. In “Molecular Recognition in Chemical and Biological Problems,” Special Pub., Royal Chem. Soc., 78, pp. 182-196 (1989))(available from the University of California, Berkeley, Calif.), 3D Database systems such as MACCS-3D (MDL Information Systems, San Leandro, Calif. (reviewed by Martin (Martin, Y. C. 3D Database Searching in Drug Design. J. Med. Chem., 35, pp. 2145-2154 (1992))), and HOOK (available from Molecular Simulations, Burlington, Mass.).
A number of techniques commonly used for modeling drugs may be employed (for a review, see: Cohen, N.C.; Blaney, J. M.; Humblet, C.; Gund, P.; Barry, D.C., “Molecular Modeling Software and Methods for Medicinal Chemistry”, J. Med. Chem., 33, pp. 883-894 (1990)). There are likewise a number of examples in the chemical literature of techniques that can be applied to specific drug design projects (for a review, see: Navia, M. A. and Murcko, M. A., “The Use of Structural Information in Drug Design”, Current Opinions in Structural Biology, 2, pp. 202-210 (1992)). Some examples of these specific applications include: Baldwin, J. J. et al., “Thienothiopyran-2-sulfonamides: Novel Topically Active Carbonic Anhydrase Inhibitors for the Treatment of Glaucoma”, J. Med. Chem., 32, pp. 2510-2513 (1989); Appelt, K. et al., “Design of Enzyme Inhibitors Using Iterative Protein Crystallographic Analysis”, J. Med. Chem., 34, pp. 1925-1934 (1991); and Ealick, S. E. et al., “Application of Crystallographic and Modeling Methods in the Design of Purine Nucleotide Phosphorylase Inhibitors” Proc. Nat. Acad. Sci USA, 88, pp. 11540-11544 (1991).
A variety of conventional techniques may be used to carry out each of the above evaluations as well as the evaluations necessary in screening a candidate compound in modulation (e.g., inhibition) of an integrin. Generally, these techniques involve determining the location and binding proximity of a given moiety, the occupied space of a bound modulator (e.g., inhibitor), the deformation energy of binding of a given compound and electrostatic interaction energies. Examples of conventional techniques useful in the above evaluations include: quantum mechanics, molecular mechanics, molecular dynamics, Monte Carlo sampling, systematic searches and distance geometry methods (G. R. Marshall, Ann. Ref. Pharmacol. Toxicol., 27, p. 193 (1987)). Specific computer software has been developed for use in carrying out these methods. Examples of programs designed for such uses include: Gaussian 92, revision E.2 (M. J. Frisch, Gaussian, Inc., Pittsburgh, Pa. © 1993); AMBER, version 4.0 (P. A. Kollman, University of California at San Francisco, © 1993); QUANTA/CHARMM [Molecular Simulations, Inc., Burlington, Mass. © 1992]; and Insight IL/Discover (Biosysm Technologies Inc., San Diego, Calif. © 1992). These programs may be implemented, for instance, using a Silicon Graphics Indigo 2 workstation or IBM RISC/6000 workstation model 550. Other hardware systems and software packages will be known and of evident applicability to those skilled in the art.
Conventional Screening
In addition to the computer-based technologies, the invention shows that the βTD binding interaction with the βA domain can be used in conventional drug library screening, to directly identify compounds capable of modulating the interaction between these two domains.
These assays may be based on binding and interaction assays where one partner is marked (e.g., by biotin or fluorescent labelling) and the other partner immobilized (e.g., on 96-well ELISA plates). Compound libraries are screened for their ability to enhance or block the interaction between the immobilized and the added biotinylated or fluorescent partner. Binding interaction is measured by anti-biotin antibodies, or fluorescence spectrometry analogous to the method described in detail for the αVβ3-vitronectin binding assay, below. Many other labelling technologies are usable in this method (e.g., radioactive marking, proximity assay).
Alternatively the interaction between the domains can be generated in a yeast two-hybrid system, using the βA domain (or larger protein fragment containing that domain) as bait and the βTD domain (or larger protein fragment containing that domain) as prey. Compound libraries can be tested for their ability to perturb the transcription of a suitable marker gene on a Gal4 promoter.
Compound Synthesis
Synthetic chemistry transformations and protecting group methodologies (protection and deprotection) useful in synthesizing the modulating compounds described herein are known in the art and include, for example, those such as described in R. Larock, Comprehensive Organic Transformations, VCH Publishers (1989); T. W. Greene and P. G. M. Wuts, Protective Groups in Organic Synthesis, 2nd ed., Jolm Wiley and Sons (1991); L. Fieser and M. Fieser, Fieser and Fieser's Reagents for Organic Synthesis, John Wiley and Sons (1994); and L. Paquette, ed., Encyclopedia of Reagents for Organic Synthesis, John Wiley and Sons (1995), and subsequent editions thereof.
The modulating compounds described herein can contain one or more asymmetric centers and thus occur as racemates and racemic mixtures, single enantiomers, individual diastereomers and diastereomeric mixtures. All such isomeric forms of these compounds are expressly included in the present invention. The modulating compounds described herein can also be represented in multiple tautomeric forms, all of which are included herein. The modulating compounds can also occur in cis- or trans- or E- or Z-double bond isomeric forms. All such isomeric forms of such modulating compounds are expressly included in the present invention.
Conservative Amino Acid Substitution
Polypeptide mimetic compounds can have a different amino acid content as the polypeptide of SEQ ID No. 1 or SEQ ID No.2 and serve as a useful mimetic. Substitution mutants can include amino acid residues that represent either a conservative or non-conservative change (or, where more than one residue is varied, possibly both). A “conservative” substitution is one in which one amino acid residue is replaced with another having a similar side chain. Families of amino acid residues having similar side chains have been defined in the art. These families include amino acids with basic side chains (e.g., lysine, arginine, histidine), acidic side chains (e.g., aspartic acid, glutamic acid), uncharged polar side chains (e.g., glycine, asparagine, glutamine, serine, threonine, tyrosine, cysteine), nonpolar side chains (e.g., alanine, valine, leucine, isoleucine, proline, phenylalanine, methionine, tryptophan), beta-branched side chains (e.g., threonine, valine, isoleucine) and aromatic side chains (e.g., tyrosine, phenylalanine, tryptophan, histidine). The invention includes polypeptides that include one, two, three, five, or more conservative amino acid substitutions, where the resulting mutant polypeptide binds a non-ligand binding site of the integrin βA domain (e.g., strand-F/α7 loop in the βA domain of unliganded αVβ3).
Fragments or other mutant nucleic acids can be made by mutagenesis techniques well known in the art, including those applied to polynucleotides, cells, or organisms (e.g., mutations can be introduced randomly along all or part of the nucleic acid encoding the polypeptide of SEQ ID No. 1 by saturation mutagenesis), and the resultant proteins can be screened for ability to inhibit integrin activation as seen in one or more of the following assays.
Integrin Inhibition Assays
The utility of the methods of identifying compounds and compounds of the present invention can be assessed by testing in one or more of the following assays as described in detail below and further described in U.S. Pat. No. 6,489,333: Purified αVβ3 (human placenta)—Vitronectin ELISA, αVβ3—Vitronectin Binding Assay, Human Aortic Smooth Muscle Cell Migration Assay, In Vivo Angiogenesis Model, Pig Restenosis Model, Mouse Retinopathy Model. The assays are assumed to be made appropriate for the integrin of interest and the following are not limiting and merely serve as examples. A compound identified by the present invention is considered to be active if it has an IC50 or Ki value of less than about 10 μM for the inhibition of αVβ3-Vitronectin Binding Assay, with compounds preferably having Ki values of less than about 0.1 βM. Tested compounds of the present invention are active in the αVβ3—Vitronectin Binding Assay as well as in cell-based assays of integrin adhesion mediated by the αVβ3-receptor. Generally, the assays can be adopted to more appropriately apply to the particular integrin of interest. For example, use of the appropriate ligand (e.g., RGD-containing, e.g., fibrinogen, vitronectin, fibronectin, thrombospondin, laminin, collagen, VCAM-1, ICAM-1, ICAM-2, Factor X, osteopontin, bone sialoprotein, or vWF), use of appropriate cell types, and use of appropriate conditions which is knowledge readily available to one skilled in the art. The assays that follow are useful for detecting inhibition or modulation of αVβ3, may also be directly applied to the testing of other integrins, or may be minimally modified for appropriate use for other integrins.
Purified αVβ3 (Human Placenta)—Vitronectin ELISA
The αVβ3 receptor can be isolated from human placental extracts prepared using octylglucoside. The extracts can be passed over an affinity column composed of anti-αVβ3 monoclonal antibody (LM609) to Affigel. The column can subsequently be washed extensively at pH 7 and pH 4.5 followed by elution at pH 3. The resulting sample can be concentrated by wheat germ agglutinin chromatography and can be identified by the presence of two bands on SDS gel and confirmed as αVβ3 by western blotting. The receptor can also be prepared in a soluble recombinant form using baculovirus expression as described.21
Affinity purified protein can be diluted at different levels and plated to 96 well plates. ELISA can be performed using fixed concentration of biotinylated vitronectin (approximately 80 nM/well). This receptor preparation is confirmed to contain the αVβ3 with no detectable levels of αVβ5 by gel (αVβ3) and by testing the effects of blocking antibodies for the αVβ3 or αVβ5 in the ELISA.
A submaximal concentration of biotinylated vitronectin can be selected based on a concentration response curve with a fixed concentration of receptor and variable concentrations of biotinylated vitronectin.
αVβ3—Vitronectin Binding Assay
Integrin-ligand binding interactions can be measured as detailed previously21. The purified receptor can be diluted with coating buffer (20 mM Tris HCl, 150 mM NaCl, 1.0 mM CaCl2, 1.0 mM MgCl26H2O, 10.0 μM MnCl2.4H2 O) and coated (100 μL/well) on Costar (3590) high capacity binding plates overnight at 4° C. The coating solution is discarded and the plates washed once with blocking/binding buffer (B/B buffer, 50 mM Tris HCl, 100 mM NaCl, 1.0 mM CaCl2, 1.0 mM MgCl2.6H2O, 10.0 μM MnCl2.4H2 O). Receptor is then blocked (200 μL/well) with 3.5% BSA in B/B buffer for 2 hours at room temperature. After washing once with 1.0% BSA in B/B buffer, biotinylated vitronectin (100 μL) and either inhibitor (11 μL) or B/B buffer w/1.0% BSA (11 μL) is added to each well. The plates are incubated 2 hours at room temperature. The plates are washed twice with B/B buffer and incubated 1 hour at room temperature with anti-biotin alkaline phosphatase (100 μL/well) in B/B buffer containing 1.0% BSA. The plates are washed twice with B/B buffer and alkaline phosphatase substrate (100 μL) is added. Color is developed at room temperature. Color development is stopped by addition of 2N NaOH (25 μL/well) and absorbance is read at 405 nm. The IC50 is the concentration of test substance needed to block 50% of the vitronectin binding to the receptor. A compound is considered to be active if it has an IC50 value of less than or equal to about 10 μM in the αVβ3—Vitronectin Binding Assay. Compounds with an IC50 less than 100 nM for the inhibition of vitronectin are desirable. Using the methodology described above, a number of compounds of the present invention can be found to exhibit an IC50 of less than or equal to about 10 μM, thereby confirming the utility of the compounds of the present invention as effective αVβ3 integrin inhibitors.
βTD-βA Domain Binding Assay
Purified βTD as Fc-fusion protein is immobilized and the interaction with the biotinylated βA domain is measured as described above for the integrin αVβ3—vitronectin binding interaction. The IC50 is the concentration of test substance needed to block 50% of the βTD binding to the receptor. A compound is considered to be active if it has an IC50 value of less than or equal to about 10 μM in the βTD-βA domain binding assay. Compounds with an IC50 less than 100 nM for the inhibition of βTD-βA interaction-are desirable.
Integrin Cell-Based Adhesion Assays
In the adhesion assays, a 96 well plate are coated with the appropriate ligand (e.g., fibrinogen, vitronectin, fibronectin, thrombospondin, laminin, collagen, VCAM-1, ICAM-1, ICAM-2, Factor X, osteopontin, bone sialoprotein, or vWF) for the integrin to be tested and incubated overnight at 4° C. The following day, the cells are harvested, washed, and loaded with a fluorescent dye. Compounds and cells are added together and then are immediately added to the coated plate. After incubation, loose cells are removed from the plate, and the plate (with adherent cells) is counted on a fluorometer. The ability of test compounds to inhibit cell adhesion by 50% is given by the IC50 value and represents a measure of potency of inhibition of integrin mediated binding. Compounds are tested for their ability to block cell adhesion using integrin interaction assays specific for the integrin of interest.
Platelet Aggregation Assay
Venous blood is obtained from anesthetized mongrel dogs or from healthy human donors who are drug- and aspirin-free for at least two weeks prior to blood collection. Blood is collected into citrated Vacutainer tubes. The blood is centrifuged for 15 minutes at 150×g (850 RPM in a Sorvall RT6000 Tabletop Centrifuge with H-1000 B rotor) at room temperature, and platelet-rich plasma (PRP) is removed. The remaining blood is centrifuged for 15 minutes at 1500×g (26, 780 RPM) at room temperature, and platelet-poor plasma (PPP) is removed. Samples are assayed on a PAP-4 Platelet Aggregation Profiler, using PPP as the blank (100% transmittance). 200 μL of PRP (5×108 platelets/mL) are added to each micro test tube, and transmittance is set to 0%. 20 μL of ADP (10 μM) is added to each tube, and the aggregation profiles are plotted (% transmittance versus time). Test agent (20 μL) is added at different concentrations prior to the addition of the platelet agonist. Results are expressed as % inhibition of agonist-induced platelet aggregation.
Human Aortic Smooth Muscle Cell Migration Assay
A method for assessing αVβ3-mediated smooth muscle cell migration and agents which inhibit αVβ3-mediated smooth muscle cell migration is described in Liaw et al., J. Clin. Invest. (1995) 95:713-724).
In Vivo Angiogenesis Model
A quantitative method for assessing angiogenesis and antiangiogenic agents is described in Passaniti et al., Laboratory Investigation (1992) 67:519-528.
Pig Restenosis Model
A method for assessing restenosis and agents which inhibit restenosis is described in Schwartz et al., J. Am. College of Cardiology (1992) 19:267-274.
Mouse Retinopathy Model
A method for assessing retinopathy and agents which inhibit retinopathy is described in Smith et al., Invest. Ophthal. & Visual Science (1994) 35: 101-111.
Contact Points of the “Deadbolt”
Table 1 below lists the contact points between the βTD and the βA domain. Table 2 below lists the more extensive contact points between the hybrid domain and the βTD. The amino acids listed in the source column of each table represent the contact points to which a modulator would be designed to bind. Modulators identified by the methods described herein would generally mimic the interactions of the peptides, c663 WRFQyyE671 D672s673S674G675KSILYVVEEPEC687 (SEQ ID No. 1), and K618KFDREPYMTENTCNRYCRD (SEQ ID NO. 2) specifically those listed in Table 1 and additionally those listed in Table 2, listed in the target column. The important amino acid contact to which an integrin modulator would bind include Ser673 as seen in Table 1. As seen in Table 2, the important amino acid contacts to which an integrin modulator would bind would include Arg633, Thr630, Glu628, and Arg636 (of the α1/strand A loop). The flexible CD loop of the βTD contacts strand-F/α7 loop (Ser674 of the former contacts Val332 in βA) in unliganded αVβ3 (Table 1)(
Interestingly, it was discovered that a D126D127/A126A127 mutant inactivates the ADMIDAS. Binding of the resulting membrane-bound integrin to physiologic ligand increases by 1.5-2-fold (
* weak bond
** strong bond
*** very strong bond
This application is a continuation (and claims the benefit of priority under 35 USC 120) of U.S. application Ser. No. 10/769,497 filed Jan. 30, 2004, which claims priority to U.S. application Ser. No. 60/444,609 filed Jan. 30, 2003. The disclosure of the prior applications are considered part of (and is incorporated by reference in) the disclosure of this application.
Number | Date | Country | |
---|---|---|---|
60444609 | Jan 2003 | US |
Number | Date | Country | |
---|---|---|---|
Parent | 10769497 | Jan 2004 | US |
Child | 11042964 | Jan 2005 | US |