Combinatorial Chemistry Computational System and Enhanced Selection Method
20210134398 ยท 2021-05-06
Inventors
Cpc classification
G16C20/10
PHYSICS
International classification
Abstract
A method for identifying a potentially useful molecular combination includes applying a selection procedure to a compound to identify a first set of candidate molecules, the procedure including providing a chemical synthesis scheme, a virtual scaffold molecule of the compound, and a virtual reactant fragment to react with the scaffold molecule according to the scheme; preparing the reactant fragment and the scaffold molecule for analyzing combinations of them; designating a remaining scaffold subset and a remaining fragment subset if a product molecule can be formed from them; rotating the fragment subset about an axis connecting the scaffold subset and the fragment subset incrementally through 360 degrees; and identifying potentially useful combinations of the reactant fragment and the scaffold molecule; identifying a set of combinatorial fragments from the first set of candidates; and applying the selection procedure to the set of combinatorial fragments to identify a second set of candidate molecules.
Claims
1. A method for identifying one or more potentially useful molecular combinations comprising: applying a selection procedure to a compound of interest to identify a first set of one or more candidate molecules, the selection procedure comprising: providing a chemical synthesis scheme for a compound of interest, a virtual scaffold molecule of the compound of interest, and a virtual reactant fragment to react with the virtual scaffold molecule according to the chemical synthesis scheme; preparing the virtual reactant fragment and the virtual scaffold molecule for analyzing combinations of the virtual reactant fragment and the virtual scaffold molecule; designating a remaining scaffold subset and a remaining fragment subset if a product molecule can be formed from the virtual scaffold molecule and the virtual reactant fragment; rotating the remaining fragment subset about an axis connecting the remaining scaffold subset and the remaining fragment subset through 360 degrees in increments of less than or equal to 5 degrees; and identifying potentially useful combinations of the virtual reactant fragment and the virtual scaffold molecule, by: recording as a potential product increment each increment at which a steric collision is not detected; and recording a separation distance between the remaining fragment subset and the remaining scaffold subset at each increment and identifying a set of product increments for which the separation distances are less than or equal to a predetermined criterion distance to identify the one or more potentially useful molecular combinations; identifying a set of combinatorial fragments from the first set of one or more candidates; and applying the selection procedure to the set of combinatorial fragments to identify a second set of one or more candidate molecules that are the one or more potentially useful molecular combinations.
2. The method of claim 1, wherein: the preparing the virtual reactant fragment and the virtual scaffold molecule comprises providing a three-dimensional coordinate system for the virtual reactive fragment.
3. The method of claim 1, wherein: the preparing the virtual reactant fragment and the virtual scaffold molecule comprises identifying a fragment alignment atom and a fragment root atom in the virtual reactant fragment; and the preparing the virtual reactant fragment and the virtual scaffold molecule comprises: identifying a scaffold alignment atom and a scaffold root atom in the virtual scaffold molecule; and providing a three-dimensional coordinate system for the virtual scaffold molecule and aligning the scaffold root atom with an origin and the scaffold alignment atom with an x-axis.
4. The method of claim 3, wherein the preparing the virtual reactant fragment and the virtual scaffold molecule comprises: aligning the fragment alignment atom with the scaffold root atom; and aligning the fragment root atom with the scaffold alignment atom.
5. The method of claim 4, wherein the axis connecting the remaining scaffold subset and the remaining fragment subset is defined by the scaffold root atom and the virtual root atom.
6. The method of claim 1, wherein: the identifying potentially useful combinations further comprises creating a product file for a configuration of the remaining fragment subset and the remaining scaffold subset at each increment of the set of product increments.
7. A non-transitory computer-readable medium encoded with a computer program for execution by a processor for identifying one or more potentially useful molecular combinations, the computer program comprising instructions for: applying a selection procedure to a compound of interest to identify a first set of one or more candidate molecules, the selection procedure comprising: receiving a chemical synthesis scheme for a compound of interest, a virtual scaffold molecule of the compound of interest, and a virtual reactant fragment to react with the virtual scaffold molecule according to the chemical synthesis scheme; receiving input to prepare the virtual reactant fragment and the virtual scaffold molecule for analyzing combinations of the virtual reactant fragment and the virtual scaffold molecule; designating a remaining scaffold subset and a remaining fragment subset if a product molecule can be formed from the virtual scaffold molecule and the virtual reactant fragment; rotating the remaining fragment subset about an axis connecting the remaining scaffold subset and the remaining fragment subset through 360 degrees in increments of less than or equal to 5 degrees; identifying potentially useful combinations of the virtual reactant fragment and the virtual scaffold molecule by: recording as a potential product increment each increment at which a steric collision is not detected; and recording a separation distance between the remaining fragment subset and the remaining scaffold subset at each increment and identifying a set of product increments for which the separation distances are less than or equal to a predetermined criterion distance, to identify the first set of one or more candidate molecules; and identifying a set of combinatorial fragments from the first set of one or more candidates; and applying the selection procedure to the set of combinatorial fragments to identify a second set of one or more candidate molecules that are the one or more potentially useful molecular combinations.
8. The medium of claim 7, wherein: the preparing the virtual reactant fragment and the virtual scaffold molecule comprises providing a three-dimensional coordinate system for the virtual reactive fragment.
9. The medium of claim 7, wherein: the preparing the virtual reactant fragment and the virtual scaffold molecule comprises identifying a fragment alignment atom and a fragment root atom in the virtual reactant fragment; and the preparing the virtual reactant fragment and the virtual scaffold molecule comprises: identifying a scaffold alignment atom and a scaffold root atom in the virtual scaffold molecule; and providing a three-dimensional coordinate system for the virtual scaffold molecule and aligning the scaffold root atom with an origin and the scaffold alignment atom with an x-axis.
10. The medium of claim 9, wherein the preparing the virtual reactant fragment and the virtual scaffold molecule comprises: aligning the fragment alignment atom with the scaffold root atom; and aligning the fragment root atom with the scaffold alignment atom.
11. The medium of claim 10, wherein the axis connecting the remaining scaffold subset and the remaining fragment subset is defined by the scaffold root atom and the virtual root atom.
12. The medium of claim 7, wherein: the identifying potentially useful combinations further comprises creating a product file for a configuration of the remaining fragment subset and the remaining scaffold subset at each increment of the set of product increments.
13. An apparatus for identifying one or more potentially useful molecular combinations comprising: a processor; a memory communicably coupled to the processor; an output device communicably coupled to the processor; and a non-transitory computer-readable medium encoded with a computer program for execution by the processor that causes the processor to: apply a selection procedure to a compound of interest to identify a first set of one or more candidate molecules, the selection procedure comprising: receiving a chemical synthesis scheme for a compound of interest, a virtual scaffold molecule of the compound of interest, and a virtual reactant fragment to react with the virtual scaffold molecule according to the chemical synthesis scheme; receiving input to prepare the virtual reactant fragment and the virtual scaffold molecule for analyzing combinations of the virtual reactant fragment and the virtual scaffold molecule; designating a remaining scaffold subset and a remaining fragment subset if a product molecule can be formed from the virtual scaffold molecule and the virtual reactant fragment; rotating the remaining fragment subset about an axis connecting the remaining scaffold subset and the remaining fragment subset through 360 degrees in increments of less than or equal to 5 degrees; and identifying potentially useful combinations of the virtual reactant fragment and the virtual scaffold molecule, by: recording as a potential product increment each increment at which a steric collision is not detected; and recording a separation distance between the remaining fragment subset and the remaining scaffold subset at each increment and identifying a set of product increments for which the separation distances are less than or equal to a predetermined criterion distance, to identify the first set of one or more candidate molecules; and identify a set of combinatorial fragments from the first set of one or more candidates; and apply the selection procedure to the set of combinatorial fragments to identify a second set of one or more candidate molecules that are the one or more potentially useful molecular combinations.
14. The apparatus of claim 13, wherein: the preparing the virtual reactant fragment and the virtual scaffold molecule comprises providing a three-dimensional coordinate system for the virtual reactive fragment.
15. The apparatus of claim 13, wherein: the preparing the virtual reactant fragment and the virtual scaffold molecule comprises identifying a fragment alignment atom and a fragment root atom in the virtual reactant fragment; and the preparing the virtual reactant fragment and the virtual scaffold molecule comprises: identifying a scaffold alignment atom and a scaffold root atom in the virtual scaffold molecule; and providing a three-dimensional coordinate system for the virtual scaffold molecule and aligning the scaffold root atom with an origin and the scaffold alignment atom with an x-axis.
16. The apparatus of claim 15, wherein the preparing the virtual reactant fragment and the virtual scaffold molecule comprises: aligning the fragment alignment atom with the scaffold root atom; and aligning the fragment root atom with the scaffold alignment atom.
17. The method of claim 16, wherein the axis connecting the remaining scaffold subset and the remaining fragment subset is defined by the scaffold root atom and the virtual root atom.
18. The apparatus of claim 13, wherein: the identifying potentially useful combinations further comprises creating a product file for a configuration of the remaining fragment subset and the remaining scaffold subset at each increment of the set of product increments.
Description
BRIEF DESCRIPTION OF THE DRAWINGS
[0014] For a more complete understanding of the features and advantages of the present invention, reference is now made to the detailed description of the invention along with the accompanying figures, in which:
[0015]
[0016]
[0017]
[0018]
[0019]
[0020]
[0021]
[0022]
[0023]
[0024]
[0025]
[0026]
[0027]
[0028]
[0029]
[0030]
[0031]
[0032]
[0033]
[0034]
[0035]
[0036]
[0037]
[0038]
[0039]
[0040]
[0041]
[0042]
[0043]
[0044]
[0045]
[0046]
[0047]
[0048]
[0049]
[0050]
[0051]
[0052]
[0053]
[0054]
[0055]
[0056]
[0057]
[0058]
[0059]
[0060]
DETAILED DESCRIPTION OF THE INVENTION
[0061] Illustrative embodiments of the system of the present application are described below. In the interest of clarity, not all features of an actual implementation are described in this specification. It will of course be appreciated that in the development of any such actual embodiment, numerous implementation-specific decisions must be made to achieve the developer's specific goals, such as compliance with system-related and business-related constraints, which will vary from one implementation to another. Moreover, it will be appreciated that such a development effort might be complex and time-consuming but would nevertheless be a routine undertaking for those of ordinary skill in the art having the benefit of this disclosure.
[0062] In the specification, reference may be made to the spatial relationships between various components and to the spatial orientation of various aspects of components as the devices are depicted in the attached drawings. However, as will be recognized by those skilled in the art after a complete reading of the present application, the devices, members, apparatuses, etc. described herein may be positioned in any desired orientation. Thus, the use of terms such as above, below, upper, lower, or other like terms to describe a spatial relationship between various components or to describe the spatial orientation of aspects of such components should be understood to describe a relative relationship between the components or a spatial orientation of aspects of such components, respectively, as the device described herein may be oriented in any desired direction.
[0063] An objective of the present invention is an increase in the efficiency of drug development programs by performing iterative molecular syntheses using efficient computational approaches, instead of synthesizing large sets of compounds related to a molecule of interest using organic chemistry synthesis methods. The present invention allows the synthesis of large numbers of variants in silico to inform choices of which of the vast numbers of possible compounds to synthesize for experimentation, saving expense and time in the drug development process. After in silico production, the virtual compound variants are computationally assessed for any predicted improvements in pharmacological characteristics, including any characteristic that can be calculated, for example, physicochemical data, such as total polar surface area, log P values, molecular weight, etc. and more complex indicators of improved drug-like characteristics, for example, predicted toxicities, mutagenicity, likelihood of inducing potential drug-drug interactions (cytochrome P450 isozyme substrate character, etc.), increased binding affinities to targeted proteins, decreased bonding affinities to undesired protein targets, and more.
[0064] A prior art method for combinatorial in silico drug-lead optimization exists.
[0065] In the course of a typical virtual-computational or nonvirtual-wet-laboratory-based drug discovery program, initial candidate hit molecules are identified as molecules that may interact with and possibly inhibit or otherwise affect a given drug target (usually a cellular protein or a protein from a pathogen). The virtually identified molecules or those identified through conventional wet-lab procedures such as high throughput screens are then normally tested in wet-lab experimentation to verify their functional properties in a relevant assay that when successful, provides support for their interaction with, and effects on, the target. It is these initial hit molecules that are molecularly varied for potential identification of variant molecules that have more optimal molecular, biochemical and pharmacological properties using the prior art method 100.
[0066] The prior art method 100 (also, herein, chemical variation procedure 100) is used as a virtual chemical synthesis procedure in embodiments of the present invention a first time to produce thousands to millions of virtual variant candidate molecules in short periods of time, e.g., one to three days. From analyses performed by those skilled in organic chemical synthesis and/or medicinal chemistry, synthetic methods to synthesize the previously identified candidate molecule(s) are devised or elucidated from the literature. From the so-elucidated chemical synthetic pathways of the candidate molecules, a set of varied precursor molecules (combinatorial fragments) can be identified that would likely be capable of reaction to generate variants of the originally identified candidate hit molecule. By applying blocks 105-135 of prior art method 100, these fragments can be combined to identify a large set of thousands or millions of variant candidate molecules that are derivatives of the originally identified candidates that may possess one or more potentially useful molecular combinations.
[0067] The second set of one or more variant candidate molecules with predicted improvements are then synthesized by organic chemistry synthetic methods and tested using actual biochemical, biophysical, pharmacological, cell biological, and/or animal experimentation, including the chemical synthesis schemes identified originally identified for a compound of interest.
[0068] Another flow diagram for the prior art method 100 is shown in
[0069] The subsequent testing of these new variants of the originally identified hit molecule using computational molecular docking programs to identify improved ligand interactions (blocks 8-11 of prior art
[0070]
[0071] Prior art for one type of intelligent pre-selection of compounds in
[0072] The type of intelligent pre-selection of compounds of interest as described in the prior art of Brewer et al. 2014, Follit et al. 2015 and Nanayakkara et al. 2018 as well as other selection schemes can be used with the prior art described in
[0073] The present invention, illustrated in a method embodiment 200 in
[0074] In the selection procedure 100 used in embodiments of the present invention, preparing the virtual reactant fragment and the virtual scaffold molecule for analyzing combinations of the virtual reactant fragment and the virtual scaffold molecule includes providing a three-dimensional coordinate system for the virtual reactive fragment. This may be done, for example, by identifying a three-dimensional coordinate system from a database of commercially available compounds such as zinc.docking.com or www.emolecules.com. SMILES-formatted or SDF-formatted 3D coordinate files for a virtual reactive fragment can be converted to PDB-type files by a program such as the openbabel program (O'Boyle et al. 2011) or the cactvs program (Ihlenfeldt et al. 2002; Ihlenfeldt et al. 1994).
[0075] In the selection procedure 100 used in embodiments of the present invention, preparing the virtual reactant fragment and the virtual scaffold molecule for analyzing combinations of the virtual reactant fragment and the virtual scaffold molecule also includes marking, via one or more marking programs, the virtual reactant fragment and the virtual scaffold molecule: identifying an alignment atom and a root atom in each of the virtual reactant fragment and the virtual scaffold molecule. In each case, a root atom is an atom that remains after the reaction and an alignment atom is an atom that does not remain after the reaction. A combination of bash (main_new10.bsh) and tcl (various generic_*.tcl) scripts running a program such as Visual Molecular Dynamics (Humphrey et al. 1996) may be used to mark the virtual reactant fragment and the virtual scaffold molecule. Using these tools, the alignment and the root atoms are marked in the beta column of respective 3D coordinate files. Each chemistry step which is inherent in marking is specifically programmed. This programming allows the selection procedure and embodiments of the present invention to go beyond the pre-specified click chemistries that limit some other prior art methods, but that do not limit the present invention.
[0076] In the selection procedure 100 used in embodiments of the present invention, preparing the virtual reactant fragment and the virtual scaffold molecule for analyzing combinations of the virtual reactant fragment and the virtual scaffold molecule further includes moving the root atom of the virtual scaffold molecule to the origin of the x,y,z coordinate system assigned to that molecule and orienting the alignment atom of the virtual scaffold molecule to the x-axis of that coordinate system. This preparation also includes aligning the alignment atom of the virtual reactant fragment with the root atom of the virtual scaffold molecule and aligning the root atom of the virtual reactant fragment with the alignment atom of the virtual scaffold molecule. A combination of bash and tcl scripts running a program such as Visual Molecular Dynamics may be used to prepare the virtual reactant fragment and the virtual scaffold molecule as described herein.
[0077] A virtual scaffold molecule of a compound of interest is a computer model of a molecule such as SMU 29 or SMU 45. A virtual reactant fragment is a computer model of a fragment that can be reacted with a compound of interest. The prior art method 100 and the prior art method 150 make use of modeling a reaction of the virtual reactant fragment with the virtual scaffold molecule according to the chemical synthesis scheme.
[0078] In the selection procedure 100 used in embodiments of the present invention, at each increment, steric hindrances of atoms may be assessed and recorded. If rotation through 360 degrees is completed and no solution without a steric hindrance is detected, a product of the remaining scaffold subset and the remaining fragment subset and is identified as failed. If increments of rotation are found that result in no steric hindrances, a separation distance between the remaining fragment subset and the remaining scaffold subset at each such increment is recorded and a set of product increments for which the separation distances satisfy one or more predetermined criteria is identified. For example, the increments with the most widely separated subset groups may be determined as preferred and identified as a preferred set of product increments.
[0079] Alternatively, in the selection procedure 100 used in embodiments of the present invention, at each increment, steric collisions of atoms may be assessed and recorded. If rotation through 360 degrees is completed and no solution without a steric collision is detected, a product of the remaining scaffold subset and the remaining fragment subset and is identified as failed. If increments of rotation are found that result in no steric collisions, a separation distance between the remaining fragment subset and the remaining scaffold subset at each such increment is recorded and a set of product increments for which the separation distances satisfy one or more predetermined criteria is identified. For example, the increments with the most widely separated subset groups may be determined as preferred and identified as a preferred set of product increments.
[0080] A virtual scaffold molecule of a compound of interest is a computer model of a molecule such as SMU 29 or SMU 45. A virtual reactant fragment is a computer model of a fragment that can be reacted with a compound of interest. The selection procedure used in embodiments of the present invention makes use of modeling a reaction of the virtual reactant fragment with the virtual scaffold molecule according to the chemical synthesis scheme.
[0081] In the selection procedure 100 used in embodiments of the present invention, analyzing combinations of the virtual reactant fragment and the virtual scaffold molecule includes identifying a remaining scaffold subset and a remaining fragment subset if a product molecule can be formed from the virtual scaffold molecule and the virtual reactant fragment. Then, the virtual scaffold molecule and the virtual reactant fragment are rotated about each other in increments around an axis defined by the root atoms of each virtual structure.
[0082] Further, the selection procedure 100 used in embodiments of the present invention uses the predicted binding energies and predicted fold improvements (defined as the ratio of the estimated K.sub.d of the compound of interest to the predicted K.sub.d of candidates for synthesis and testing) as factors in assessing the suitability of candidates for synthesis and testing.
[0083]
[0084]
[0085]
[0086] Where a chemical synthesis scheme has one or more intermediate steps, as, for example the scheme for SMU 29 has, the method of the present invention can be repeated as necessary. Three-dimensional coordinates for a virtual reactant fragment for an intermediate step may be obtained as described herein. This intermediate virtual reactant fragment and the last product from the last analysis, as the intermediate virtual scaffold molecule of the intermediate compound of interest, can be prepared and analyzed as described herein. Potentially useful combinations of the intermediate virtual reactant fragment and the intermediate virtual scaffold molecule can be identified as described herein.
[0087] The skilled artisan will recognize that methods and systems of the present invention allows the synthesis of large numbers of variants in silico to inform choices of which of the vast numbers of possible compounds to synthesize for experimentation, saving expense and time in the drug development process.
[0088] Example of Use of an Embodiment of the Present Invention.
[0089] Resistances of cancers to chemically unrelated anti-cancer drugs are frequently caused by the expression of members of the ABC transporter superfamily, including ABCB1 (P-glycoprotein or P-gp).sup.1-3 and/or ABCG2 (the breast cancer resistance protein or BCRP).sup.3, 4. The phenomenon of multidrug resistance (MDR) remains a major obstacle in the treatment of both adult and pediatric cancers.sup.5-7. Despite close to 40 years of intense research, no inhibitor of these proteins has yet been approved for clinical use.sup.8-13. The reasons for failure in clinical trials are multifaceted: some could be attributed to flaws in trial design, others reported tumor penetration problems, and still others failed because of drug-drug interactions and associated toxicities. A significant number of failures of the clinical trials may be attributed to the fact that many of the assessed MDR-inhibitors were also good transport substrates of the pumps. This latter characteristic likely required elevated systemic inhibitor concentrations for efficacy that may have resulted in significant off-target toxicities. The fact that a Phase III trial using the immunosuppressant, cyclosporine A, led to improved patient outcomes in poor-risk acute myeloid lymphoma patients.sup.14 suggests, however, that despite the limited success in finding clinically successful inhibitors of MDR pumps, these proteins are important targets for drug discovery and development.
[0090] One of the biggest obstacles to developing effective inhibitors of membrane proteins like P-gp and BCRP has been the dearth of detailed structural and mechanistic knowledge, the lack of which made targeting these pumps ineffective. Over the last several years, however, significant advances in knowledge of the structure.sup.18-20 and mechanism.sup.21-25 of P-gp, BCRP and other related pumps have emerged that will enable the design of potent inhibitors of the pumps that may prove to be more successful in clinical applications.
[0091] Using the evolutionary relationship of different ABC transporters and the structural knowledge of both prokaryotic and eukaryotic ABC transporters, dynamic models of human P-gp 2 were created and a putative catalytic cycle.sup.22 was simulated that correlated well with published biochemical and biophysical studies as well as with the recently elucidated outward facing structure of the human P-gp.sup.16. These conformationally dynamic models of human P-gp in ultrahigh throughput in silico screenings were previously used to identify and characterize inhibitors of P-glycoprotein.sup.26. One desirable characteristic of the inhibitors identified in these screens was that these potential P-gp inhibitors should not be transport substrates of the pump.sup.26. For this reason, drug-like molecules that were computationally predicted to interact well with the nucleotide binding domains were computationally counter-screened for interactions with the drug binding parts of the protein. Compounds were discarded from further evaluation, if significant binding to the drug binding sites was predicted by the in silico docking calculations.sup.26. Using this approach, three compounds were initially identified and characterized that reversed the multidrug resistance phenotype of various cancer cell lines in both conventional and microtumor-spheroid cell cultures.sup.26-28. The efficacy of the compounds was assessed using P-gp overexpressing, multidrug resistant prostate and ovarian cancer cells in culture. All three compounds were observed to reverse multidrug resistance by increasing lethality of various chemotherapeutics. The increased lethality was correlated with increased cellular retention of the chemotherapeutics when inhibitor was present.sup.28. Importantly, studies also showed that these three P-gp inhibitors were not significantly transported by P-gp.sup.28, supporting the premise that these inhibitors would not bind effectively to the drug binding domains of the pump.sup.26.
[0092] One of these compounds that reversed MDR phenotypes in cancer cells was predicted to be an allosteric inhibitor of P-glycoprotein (SMU-29, compound 29, or 29 herein, 2-[(5-cyclopropyl-4H-1,2,4-triazol-3-yl)sulfanyl]-N-[2-phenyl-5-(2,4,5-trimethylphenyl)-pyrazol-3-yl,
[0093] Inspection of the computational docking of 29 at the putative allosteric site on P-gp (
[0094] A number of computational approaches exist that seek to analyze the potential chemical space of a hit compound by creating virtual libraries of variants of a given hit molecule.sup.29-32. These approaches mostly vary in the chemistry that is applied to perform the synthesis reactions. Presented herein are initial efforts at creating variants of compound 29 with optimized binding affinity to P-glycoprotein using a novel computer-aided and structure-based approach that was applied to the Western end of compound 29. Results of the virtual synthesis of a moderate number of variants of 29 and the virtual screening of these variants with structural models of P-gp are reported. A small portion of the nearly infinite chemical space around hit compound 29 was synthesized and assessed for reversal of the MDR phenotype in a multidrug resistant prostate cancer cell line that over-expresses P-gp (DU145TXR.sup.33). Using the same cell line, the inhibition of P-gp-catalyzed pumping of a P-gp substrate by these 29 variants was also assessed. In addition, biochemical analysis of the mode of P-gp inhibition was performed for all 29-variants using ATP hydrolysis and ATP binding assays as in.sup.26. After initial evaluation of the computationally predicted inhibitor variants in these assays as well as of the physicochemical properties of these variants, we developed a new, structure-based rational design to synthesize and analyze a small number of 29-derivatives with different structural and physicochemical characteristics. These compounds were not initially computationally evaluated using the subtractive binding routines as described in.sup.26, but were chosen mostly for the shape and size of the Western half of the molecule as well as for their physicochemical characteristics like polar surface area and solubility. All of the novel 29-variants were experimentally assessed for their potential of being transported by P-gp. The work led to the discovery of several variants of P-gp inhibitor hit compound 29 with improved efficacy in reversing MDR in P-glycoprotein over-expressing cancer cells by inhibiting P-gp catalyzed substrate pumping.
[0095] Results. Virtual synthesis of novel variants of the Western half of SMU-29 using the ChemGen computational suite. Evaluation of the fit of compound 29 into a putative allosteric binding site on P-gp as visualized from the results of docking studies the inventors recognized that if variants of inhibitor 29 were made larger and more hydrophobic they would likely fill the relatively large hydrophobic pocket in the protein where the cyclopropyl group of 29 interacts (
[0096] Some of the potential chemical space of these 647 ChemGen-produced Western variants of 29 (Group 1 compounds) is visualized in
[0097] Docking to P-glycoprotein and chemical synthesis of compound 29 variants. The 647 Group 1 molecules created by ChemGen were used in docking studies to the same structural model of P-gp that was employed previously and led to the identification of compound 29 as a potential inhibitor.sup.26 (see steps 8 through 10 of
TABLE-US-00001 TABLE 1 ChemGen/docking routine identified variants of P-gp inhibitor 29. Topological ratio polar estimated estimated K.sub.d 29/ Molecular surface Synthesized G.sub.binding K.sub.d K.sub.d weight area Consensus Variant name variant (kcal/mol) (nM) variant (Da) (.sup.2) logP opt_0009_19 12.4 0.8 50 599 173.8 3.4 opt_0005_53 29-541 12.0 1.6 25 522 132.5 3.9 opt_0009_29 29-216 11.8 2.3 17 558 89.3 6.6 opt_0010_33 11.8 2.3 17 576 151.5 3.8 opt_0005_47 11.8 2.3 17 591 138.6 4.3 opt_0000_30 11.7 2.7 15 541 157 3.4 opt_0002_34 11.7 2.7 15 579 124 4.7 opt_0005_27 11.7 2.7 15 584 131.4 4.4 opt_0010_49 29-551 11.6 3.2 13 627 114.2 5.9 opt_0009_33 29-231 11.5 3.8 11 535 101.3 5.7 opt_0004_12 29-227 11.5 3.8 11 587 92.5 6.4 opt_0006_13 11.5 3.8 11 596 130.1 4.6 ZINC08767731 29 10.1 40 1 461 113.8 4.0
[0098] Estimated K.sub.d values for the ligand interactions with P-gp were calculated from the lowest estimated binding G values from the AutoDock calculations. The ratio of K.sub.d values is shown as a relative value for increased affinity exhibited by the respective variants over the parent P-gp inhibitor compound 29. Molecular weights, topological polar surface areas, and consensus log P values were calculated at the SwissADME website (http://www.swissadme.c/) as described herein.
[0099] One additional compound (29-551) which contains abromo-substituent and has a molecular weight of 627 Da was added for consideration.
[0100] Syntheses were performed using the scheme shown in
[0101]
[0102] The Group 1 ChemGen/docking routine produced variants of P-gp inhibitor 29 resensitize a multidrug resistant, P-gp overexpressing prostate cancer cell line to paclitaxel. The mitochondrial reduction potential of cells is often used as an indicator for cell viability using MTT assays.sup.4. Using these assays it was observed that the five Group 1 derivatives of 29 predicted through the ChemGen/docking routine and synthesized here, 216, 227, 231, 541 and 551, re-sensitized the P-gp overexpressing prostate cancer cell line, DU145TXR.sup.33, to the chemotherapeutic, paclitaxel (
TABLE-US-00002 TABLE 2 Increased toxicity of paclitaxel to DU145TXR in the presence of P-gp inhibitors identified by the ChemGen/docking routine. Resensitization to paclitaxel with the indicated treatment and fold increased sensitivity in the presence of inhibitors PTX PTX + 29 PTX + 216 alone Fold Fold PTX + 227 IC.sub.50 IC.sub.50 Fold vs IC.sub.50 Fold vs IC.sub.50 Fold Inhibitor PTX PTX vs PTX + PTX vs PTX + PTX vs concentration (nM) (nM) PTX 29 (nM) PTX 29 (nM) PTX 3 M 2120 629 3 1 266 8 2.4 164 13 5 M 2120 194 11 1 154 14 1.3 52 41 7 M 2120 59 36 1 62 34 1 6 353 10 M 2120 21 101 1 57 37 0.4 4 530 Resensitization to paclitaxel with the indicated treatment and fold increased sensitivity in the presence of inhibitors PTX + 227 PTX + 231 PTX + 541 PTX + 551 Fold Fold Fold Fold vs IC.sub.50 Fold vs IC.sub.50 Fold vs IC.sub.50 Fold vs Inhibitor PTX + PTX vs PTX + PTX vs PTX + PTX vs PTX + concentration 29 (nM) PTX 29 (nM) PTX 29 (nM) PTX 29 3 M 3.8 153 14 4.1 426 5 1.5 56 38 11 5 M 3.7 46 46 4.2 21 101 9.2 20 106 9.7 7 M 10 7 303 8.4 17 125 3.4 9 236 6.6 10 M 5 2 1060 11 17 125 3.4 6 353 3.5
[0103] Cytotoxicity of the chemotherapeutic, paclitaxel (PTX) to P-gp overexpressing prostate cancer cells, DU145TXR, was determined in the absence and presence of the ChemGen designed 29-variants, 216, 227, 231, 541 and 551. For each experimental compound IC.sub.50 values of PTX alone or in the presence of inhibitors, fold improvement of PTX sensitivity in the presence of inhibitor, and fold improvement of PTX sensitivity by variants compared to parental compound 29 are given.
[0104] At increasing concentrations, the efficacy of variant 216 in increasing paclitaxel 10 toxicity decreased when compared to the parental compound 29. At the highest concentration tested (10 M), the 216 variant was observed to be somewhat less effective than parental compound 29 (0.4 fold compared to the 1 fold of paclitaxel+29). Unlike 216, variants 227, 231, 541 and 551 were more effective than 29 at all concentrations tested. At 3 M, the presence of variants 227, 231, 541 and 551 resulted in 4 to 11-fold decreased paclitaxel IC.sub.50 when compared to parental compound 29 and up to 38-fold overall sensitization to paclitaxel when compared to paclitaxel alone (compound 551). At higher concentrations (5 to 10 M), addition of these variants resulted in increased paclitaxel toxicity and decreased paclitaxel IC.sub.50 of up to 500-1000-fold at 10 M, as compared to 100-fold sensitization caused by the parental compound 29 at 10 M. This data indicated that variants 227, 231, 541 and 551 were better re-sensitizers of the multidrug resistant cells to paclitaxel than the original compound 29 at all concentrations tested, while variant 216 appeared to be marginally better than 29 at lower concentrations. These data demonstrate that the ChemGen generated and docking analyses selected Group 1 variants of compound 29 had increased affinity for P-gp resulting in improved efficacy for reversing chemotherapy resistance in a P-gp overexpressing cancer cell line than did the parental compound.
[0105] Accumulation and cellular retention of calcein AM in P-gp overexpressing prostate cancer cells upon incubation with Group 1 SMU-29 variants. Calcein AM accumulation assays have been used by us previously to evaluate P-gp-substrate accumulation in real time in the presence or absence of P-gp inhibitors.sup.28. For these assays, P-gp overexpressing DU145TXR cells were incubated with the respective inhibitors in the presence of the P-gp substrate, calcein AM. Inhibition of P-gp leads to cellular accumulation of calcein AM and to cleavage of its acetoxymethyl ester groups, resulting in the generation of the highly fluorescent compound, calcein. The anionic calcein is not a substrate of P-gp and remains in the cells. In these assays, the relative fluorescence of cellular calcein was measured over time and the results of these assays are shown in
[0106] To test whether the lower accumulation of calcein in the presence of the Group 1 29 variants was the result of retention of the compounds in the cellular membrane due to their mostly increased log P values relative to 29 (Table 1), similar calcein accumulation assays were performed after a 6-hour pre-incubation with the 29-variants and parental compound 29,
[0107] Assessing the roles of polarity and size of 29 variants in improving efficacy of compound 29 variants. To assess the contributions of overall hydrophobicity and size of the western halves of the 29-variants on efficacy in inhibiting P-gp and reversing multidrug resistance in cancer cells, five structural derivatives of 29 (Group 2 variants,
TABLE-US-00003 TABLE 3 29-variants differing in overall shape, size and polarity that did not undergo docking routine. Topological ratio polar estimated estimated K.sub.d 29/ Molecular surface Synthesized G.sub.binding K.sub.d K.sub.d weight area Consensus variant name (kcal/mol) (nM) variant (Da) (.sup.2) logP 29-238 9.6 92 0.4 575 129.0 4.8 29-255 11.0 9 4.4 542 142.5 5.3 29-278 9.9 55 0.7 543 101.3 5.7 29-280 9.6 92 0.4 575 101.3 6.0 29-286 9.9 55 0.7 517 101.3 5.2 ZINC08767731, 10.1 40 1 461 108.6 3.8 29
[0108] Estimated K.sub.d values for the ligand interactions with P-gp were calculated from the lowest estimated binding G values from the AutoDock calculations. The ratio of K.sub.d values is given as a relative value for potentially changed affinities exhibited by the respective variants over the parent P-gp inhibitor compound 29. Molecular weights, topological polar surface areas, and consensus log P values were calculated at the SwissADME website (http://www.swissadme.ch/) as described herein.
[0109] The variations of structure in these molecules were made in the western half of the molecule (
[0110] Table 3 shows that all five of these Group 2 variants are somewhat larger than 29, but the calculated log P values for these Group 2 compounds (Table 3) are closer to that of 29 than the log P values of the Group 1 ChemGen generated variants with the exception of 541 (Table 1). The consensus log P values.sup.43, 44 of 238 and 255 were calculated to be 4.8 and 5.3, while those of 278, 280 and 286 were calculated to be 5.7, 6.0 and 5.2, respectively. The topological polar surface areas (TPSA) of 238 and 255 were calculated to be higher than those of 278, 280 and 286. These TPSA values were also higher than those of the Group 1 variants, 216, 227 and 231 (compare Tables 1 and 3). Of the five structural variants in Table 3, two had increased calculated topological polar surface areas and three had reduced calculated topological polar surface areas when compared with 29.
[0111] Subsequent docking of the Group 2 variants to the putative allosteric inhibitor binding site of P-gp showed that the western portions of compounds 238, 255 and 286 could penetrate deeper into the hydrophobic void than does compound 29, but that compounds 278 and 280 did not seem to penetrate into the hydrophobic void as well as does compound 29 (compare
[0112] Effects of Group 2 variants on the paclitaxel sensitivity of P-glycoprotein overexpressing prostate cancer cells, DU145TXR. MTT assays were again used to assess the efficacy of the new structural 29-variants on sensitizing the chemotherapy-resistant prostate cancer cell line, DU145TXR, to paclitaxel, at concentrations of 3 M, 5 M, 7 M, and 10 M as shown in
TABLE-US-00004 TABLE 4 Topological ratio polar estimated estimated K.sub.d 29/ Molecular surface Synthesized G.sub.binding K.sub.d K.sub.d weight area Consensus variant name (kcal/mol) (nM) variant (Da) (.sup.2) logP 29-238 9.6 92 0.4 575 129.0 4.8 29-255 11.0 9 4.4 542 142.5 5.3 29-278 9.9 55 0.7 543 101.3 5.7 29-280 9.6 92 0.4 575 101.3 6.0 29-286 9.9 55 0.7 517 101.3 5.2 ZINC08767731, 10.1 40 1 461 108.6 3.8 29
[0113] The data suggest that the presence of all of the variants except for 286 increased the paclitaxel toxicities to these P-gp overexpressing cells. Compounds 238 and 255 increased paclitaxel toxicities to the greatest extents: At 5 M concentration, compound 238 decreased the paclitaxel IC.sub.50 by about a thousand-fold, similar to compound 255 at 7 M. Comparison with parental compound 29 showed that 238 improved sensitization of DU145TXR to paclitaxel by 3.4-fold at 3 M, and by around 100-fold at 5 and 7 M. At 10 M, compound 238 resulted in a 300-fold decreased paclitaxel IC.sub.50 of DU145TXR compared to parental compound 29 at the same concentration. Compound 255 also showed substantial improvement in sensitization of DU145TXR to paclitaxel that was comparable to 238 at 3 M, but the effect was not as pronounced at higher concentrations. The other three compounds, 278, 280 and 286, were similar or less effective than the Group 1 variants, 216, 227 and 231. In addition, compound 280 seemed to exhibit some toxicity to the multidrug resistant cancer cells as judged by the lowered cell viability at very low concentrations of paclitaxel in the presence of 280 (
[0114] Accumulation and cellular retention of calcein AM in DU145TXR in the presence of 29-variants from Group 2.
[0115] Evaluation of mode of inhibition of P-glycoprotein by Group 1 and Group 2 variants of compound 29. To assess the mode of inhibition of P-gp by the novel variants of P-gp inhibitor 29, ATP hydrolysis by P-gp was evaluated in the presence or absence of the variants. Both basal ATP hydrolysis (assayed in the absence of added transport substrate) and stimulated ATP hydrolysis (assayed in the presence of the P-gp transport substrate, verapamil) were assessed as described in reference 2.sup.66. Murine P-gp (MDR3) expressed in Pichia pastoris that had all naturally occurring cysteines replaced with alanine, was used.sup.45, 46. It is widely assumed that the ATP hydrolytic rate of P-gp is stimulated in the presence of transport substrates when compared to ATP hydrolysis in the absence of transport substrates. Assays comparing these rates can therefore be useful not only in identifying inhibitors of P-gp catalyzed ATP hydrolysis, but also to potentially infer whether a compound might be a transport substrate if basal ATPase rates are stimulated by the addition of a compound.
[0116] Effects of compound 29 variants on verapamil-stimulated ATP hydrolysis by P-gp. The effects of compound 29 variants on P-gp ATP hydrolysis rates assayed in the presence of verapamil (a good substrate for transport by P-gp) are presented in Table 5 A (Stimulated ATPase).
TABLE-US-00005 TABLE 5 Mode of Inhibition of Cysteineless Mouse MDR3 P-glycoprotein by Group 1 and Group 2 compound 29 variants. Cellular accumulation:ratio Stimulated of plus ATPase Effect on Basal ATPase Effect on Tariquidar (% of DMSO | stimulated (% of DMSO | basal over no Compound significance) ATPase significance) ATPase Tariquidar DMSO 100 8 100 7 SMU29 49 2 ** inhibitor 95 11 NS none 1.0 Group 1 - ChemGen and docking selected SMU29-216 108 2 NS none 105 7 NS none 1.1 SMU29-227 50 2 ** inhibitor 88 4 NS none 1.0 SMU26-231 141 18 * stimulator 70 0 * inhibitor 0.9 SMU29-541 68 7 * inhibitor 101 12 NS none 1.0 SMU29-551 62 5 ** inhibitor 150 6 ** stimulator 1.1 Group 2 - Rationally designed/no docking selection SMU29-238 194 22 ** stimulator 1143 46 ** stimulator 1.9 SMU29-255 123 15 NS none 355 7 *** stimulator 1.2 SMU29-278 41 7 ** inhibitor 78 4 * inhibitor 1.0 SMU29-280 116 4 * stimulator 148 4 * stimulator 1.2 SMU29-286 98 2 NS none 143 32 NS none 1.3 Cellular SL-ANP accumulation:ratio Maximum binding to of plus ATP binding P-gp Tariquidar Transport (mol SL-ANP Apparant Effect on over no substrate bound/mol Kd SL-ANP Compound Tariquidar for P-gp P-gp) (M) binding DMSO 1.8 0.1 36.5 3.6 SMU29 NS no 71.0 12.6 no Group 1 - ChemGen and docking selected SMU29-216 NS no 1.7 0.1 23.1 3.9 no SMU29-227 NS no 1.8 0.1 36.9 4.1 no SMU26-231 NS no 1.8 0.1 22.2 3.8 marginally SMU29-541 NS no 1.8 0.1 24.1 4.0 no SMU29-551 * no 1.7 0.1 22.8 3.6 no Group 2 - Rationally designed/no docking selection SMU29-238 *** yes 1.2 0.1 20.1 4.0 yes SMU29-255 * yes 1.9 0.1 25.6 4.7 no SMU29-278 NS no 1.3 0.1 22.9 4.5 yes SMU29-280 NS yes 1.6 0.1 21.6 4.6 marginally SMU29-286 * yes 1.9 0.1 25.7 4.9 no
indicates data missing or illegible when filed
[0117] ATP hydrolysis assays using purified P-glycoprotein were performed without added transport substrate (basal ATPase) or in the presence of verapamil (Stimulated ATPase). Results are presented compared to DMSO control standard deviation (three independent experiments with duplicate samples). Basal activity of P-gp was 20 to 30 nmol/min mg, verapamil-stimulated rates were 200 to 400 nmol/min mg P-gp. Stimulation of basal ATPase by 29-variants was used as an indicator that a compound may be a P-gp transport substrate. Effects on stimulated P-gp ATPase activity indicated whether a compound directly interfered with ATP usage by the protein (***, p<0.001; **, p<0.01; *, p<0.1; NS, not significant). Quantitative cellular accumulation of 29-variants was performed using LC-MS/MS and is presented as a ratio of the cellular amounts of 29-variants in the presence of P-gp inhibitor, tariquidar, divided by amounts accumulated in its absence. A ratio >1 indicates that the compound likely is a transport substrate of P-gp (***, very significant; *, significant; NS, not significant). Binding of an ATP analog, SL-ATP, to P-gp was used to determine whether ATP binding to P-gp was affected by the 29-variants. Values +/standard deviations are shown for at least three different P-gp preparations and three independent SL-ATP titration experiments. The values for SL-ATP binding in the presence of 29 were taken directly from Brewer et al. (2014).
[0118] The respective percent ATPase activity is shown, normalized to the ATPase in the presence of DMSO carrier/no added experimental compound. Interestingly, the Group 1 compounds differed in their effects on stimulated ATPase: 216 did not affect stimulated ATP hydrolysis activities, while compounds 227, 541 and 551 inhibited activity similar to parental compound 29. Compound 231 slightly stimulated ATP hydrolysis rates in the presence of verapamil. For Group 2 compounds, 238 stimulated the stimulated ATPase rates by about two-fold, while variant 280 showed only a slight stimulation of hydrolysis rates and compounds 255 and 286 had no significant effect. Only compound 278 of the Group 2 variants inhibited verapamil-stimulated ATP hydrolysis by P-gp similar to the parental compound 29.
[0119] Effects on basal ATP hydrolysis rates of compound 29 variants. Group 1 compounds 216, 227 and 541 did not significantly affect basal ATP hydrolysis by P-gp, while compound 231 inhibited the basal ATPase rates of P-gp. Only 551 of the Group 1 molecules stimulated basal ATPase activities of P-gp. Of the Group 2 compounds, 238, 255 stimulated basal ATPase by 10 and 3 fold respectively. Compounds 280 and 286 stimulated basal ATPase only marginally or with no statistical significance. Only compound 278 inhibited basal ATPase of P-gp. The relatively strong activation of basal ATPase by compounds 238 and 255 was suggestive that these two compounds and potentially to a lesser extent, compound 280, may be transport substrates of the pump. Compound 278 was not indicated to be a good transport substrate for P-gp since it inhibited basal ATPase by P-gp.
[0120] Intracellular accumulation of compound 29 variants. Cell accumulation assays for each of the 29 variants were performed as in reference.sup.28 to more directly assess whether the compounds were indeed transport substrates for P-gp. These assays measured the intracellular accumulation of the experimental compounds using LC-MS/MS methods after incubation with the P-gp over-expressing cell line, DU145TXR, in the absence and presence of the strong P-gp inhibitor, tariquidar.sup.47 (TQR). Low levels of cellular accumulation of a compound in the absence of tariquidar accompanied by much higher levels of accumulation in the presence of tariquidar suggests that the compound in question may be a transport substrate of P-gp. In other words, if a compound is an effective transport substrate for P-gp, active P-glycoprotein in these cells would limit intracellular accumulation, while inhibited P-gp would result in higher intracellular concentrations. Daunorubicin (DNR) is an example of a good transport substrate for P-gp and showed very strong cellular accumulation in these assays when P-gp was inhibited by tariquidar, but much less accumulation in the cells when P-gp was not inhibited (see
[0121] None of the Group 1 molecules tested resulted in intracellular accumulations that were considerably different in the absence versus presence of TQR, similar to the parental compound 29 (
[0122] Of the Group 2 compounds, variant 238 showed a very large and significant increase in intracellular accumulation in the presence of TQR (
[0123] To assess whether the observed discrepancies of compounds stimulating basal P-gp ATPase activity but not being transport substrates of the human pump in the cell culture assessments were due to the fact that these biochemical assays used a cysteineless variant of the mouse MDR3 P-glycoprotein, the experiments were repeated using normal human MDR1 P-gp. In order to stabilize the human protein for the activity assays, the protein was assembled into membrane nanodiscs as described herein. The results of the experiments are shown in Table 6.
TABLE-US-00006 TABLE 6 Effects of Group 1 and Group 2 compound 29 variants on ATP Hydrolysis by Normal Human MDR1 P-glycoprotein. Stimulated ATPase Effect on Basal ATPase Effect on (% of DMSO | stimulated (% of DMSO | basal Compound significance) ATPase significance) ATPase DMSO 100 6 100 5 SMU29 62 3 ** inhibitor 88 9 NS none Group 1 - ChemGen and docking selected SMU29-216 64 6 ** inhibitor 69 8 ** inhibitor SMU29-227 62 5 ** inhibitor 83 7 NS none SMU29-231 40 3 **** inhibitor 70 4 *** inhibitor SMU29-541 66 6 ** inhibitor 87 9 NS none SMU29-551 27 1 *** inhibitor 72 9 * inhibitor Group 2 - Rationally designed/no docking selection SMU29-238 61 7 ** inhibitor 89 7 NS none SMU29-255 59 6 ** inhibitor 84 9 NS none SMU29-278 63 6 ** inhibitor 74 6 ** inhibitor SMU29-280 54 7 *** inhibitor 76 8 * inhibitor SMU29-286 59 4 *** inhibitor 95 3 NS none
[0124] ATP hydrolysis assays using purified P-glycoprotein were performed without added transport substrate (basal ATPase) or in the presence of verapamil (Stimulated ATPase). Results are presented compared to DMSO control standard deviation (three independent experiments with duplicate samples). The specific basal activity of normal human MDR1 P-gp was between 123 and 193 nmol min.sup.1g.sup.1, and transport substrate (verapamil) stimulated activity was between 193-263 nmol min.sup.1mg.sup.1. Effects on stimulated P-gp ATPase activity indicated whether a compound directly interfered with ATP usage by the protein (***, p<0.001; **, p<0.01; *, p<0.1; NS, not significant).
[0125] Interestingly, neither compound 29 nor any of is variants had a stimulatory effect on basal ATPase activity of the normal human protein reconstituted into membrane nanodiscs, while all of them significantly inhibited transport substrate (verapamil) stimulated activity. The results clearly indicate that the source (human vs. mouse) and potentially also the membrane environment of P-glycoprotein strongly affects the overall behavior of potential biochemical inhibitors.
[0126] Effects of 29 variants on binding of an ATP-analog to purified P-glycoprotein. ATP binding in the presence of the 29 variants was assessed in titration assays using a spin-labeled analog of ATP, 2,3-SL-ATP (2,3-(2,2,5,5,-tetramethyl-3-pyrroline-1-oxyl-3-carboxylic acid ester) ATP; (2,3 indicates a rapid equilibrium of the ester bond between the C2 and C3 of the ribose moiety)).sup.48-50, and electron spin resonance spectroscopy as described in.sup.26. Due to the lower stability of the human P-glycoprotein in the extended times needed for these experiments, the cysteineless mouse protein was used here. The goal was to assess whether binding of the 29 variants to P-gp affected nucleotide binding to the protein. Results of these assays are presented in Table 5A. Except for compounds 238 and 278, neither of which initially underwent the selective docking routines used for Group 1 compounds, none of the novel inhibitors affected maximal binding of the ATP analog or the apparent K.sub.d, showing that the inhibitors were indeed targeted to the putative allosteric binding site on P-gp. Compounds 238 and 278 reduced SL-ATP binding to about 1 mol SL-ANP (adenine nucleotide with an undefined number of phosphoryl groups) bound per mol of enzyme, suggesting that these inhibitors may also interact with the nucleotide binding sites or may indirectly induce changes in nucleotide binding to P-gp that affect ATP binding.
[0127] In addition to evaluating the effects of the decreased hydrophobicity of the Group 2 variants on the reversal of MDR in cell-based assays, it was also of interest to evaluate whether or not the docking routines employed for the Group 1ChemGen derived variants of 29 were better able to predict compounds that were not transport substrates of P-gp when compared to the Group 2 compounds. All five of the Group 1 molecules that were chosen though the subtractive docking routine and that were predicted to not interact well with the drug binding domains were observed to not be transport substrates. When Group 2 molecules were docked to a structural model of P-gp that was essentially identical to the one initially used to identify parental compound 29 as a P-gp inhibitor.sup.26 as well as to choose Group 1 variants, four of the five molecules were predicted to interact well with the drug binding domains and one was not predicted to bind well (data not shown). Three of these five predictions were confirmed by the LC MS/MS experiments described above. All in all, the experimental data for Groups 1 and 2 suggested that the subtractive docking method was predicting the correct outcome (transport substrate vs. no transport substrate) in 4 out of 5 cases or at 80%.
[0128] While optimization efforts of hit compound 29 did not include aspects of compound toxicity, it seems of interest to note that only one of the 29-variants (compound 541) showed some toxicity in cell viability assays in the P-gp overexpressing DU45TXR.sup.13 cancer cells when administered in the absence of chemotherapeutic. No significant toxicity of the compounds was observed in non-cancerous human lung fibroblast cells, HFL-1.sup.51 (data not shown). In addition, toxicity of the chemotherapeutic, paclitaxel, was not increased in the presence of 29 or 29-variants in cells that do not overexpress P-gp, i.e. HFL-1 and the not chemotherapy resistant, not P-gp overexpressing prostate cancer line, DU145.sup.52 (data not shown). The overall results indicate that increased lethality of paclitaxel to the P-gp overexpressing cells was due to the inhibition of the pump and increased accumulation of paclitaxel to therapeutic levels within the cells. It should be noted that compounds 541 and 551 were not assessed in these latter experiments.
[0129] Using computational approaches to create novel variants of hit molecules from drug discovery programs. In drug development, often large sets of molecules that are related to a molecule of interest are synthesized to identify derivative molecules with better drug-like characteristics than those of the originally identified molecule. The required organic chemistry syntheses are expensive, costly in time, can be very laborious, and require the skills and time of highly qualified chemists. Often hundreds or thousands of compounds are synthesized with only a few variants showing desired improvements in pharmacological characteristics. Even in these early steps of a medicinal chemistry project, these efforts add to the already significant costs of drug development.sup.53. While some drug-like characteristics can be estimated from computational approaches, most require biochemical, biophysical, pharmacological, cell biological and/or animal experimentation to assess potential improvement over the parental compound, which again increases the time and expenditures required for each potential lead compound.
[0130] A number of virtual chemical synthesis computer programs have been previously described. Some use fragments annotated with reaction rules.sup.29 or compound scaffolds with chemically reactive linkers.sup.30, and still others use popular click chemistries that can easily translate into the laboratory.sup.31, 32 just to mention only a few. To make more informed choices about which of the vast numbers of possible compound variants to synthesize for subsequent testing, a set of computational routines (collectively called ChemGen) have been written and developed to synthesize in silico what can be very large numbers of variant compounds. The methods differ from predecessor methods in that retrosynthetic approaches to the discovered hit molecule synthetic routes are mimicked in the computations. This results in advantageous translation to actual chemical syntheses of identified variants of interest, is not constrained to one or a few chemical reaction types, but can theoretically encompass any chemical reaction. A disadvantage is that each reaction type must be programmed ahead of its implementation, but the ChemGen platform may be adapted to new chemistries relatively easily.
[0131] Next, the inventors determined whether an increase in the efficiency of drug development can be achieved by performing iterative virtual molecular syntheses using efficient computational approaches instead of physically synthesizing a large set of compounds related to a molecule of interest. After production of the virtual compound variant library, the new molecules were computationally assessed for predicted improvements in any pharmacological characteristics that can be calculated, including relatively simple physiochemical data (topological polar surface area or TPSA, log P values, molecular weight, etc.) as well as more complex indicators of improved drug-like characteristics such as predicted toxicities, mutagenicity, likelihood of inducing potential drug-drug interactions (cytochrome P450 isozyme substrate character, etc.). Other important factors that can be calculated and that are valuable for decision making are predicted increased binding affinities to targeted proteins as well as potentially decreased binding affinities to undesired protein targets of the drug lead compounds. In recent years, machine learning machine methods have been employed for predicting toxicities, various ADME characteristics and even protein-ligand binding affinities of molecules of potential interest.sup.54-59.
[0132] By creating hit variants computationally and then assaying themagain computationallyfor improved characteristics, variants that do not possess the desired improved characteristics can be eliminated from consideration before any actual organic synthetic chemistry is performed. This path can lead to expedited and much more cost-effective syntheses of a relatively small number of potentially improved hit-variants. The latter part of this approach, namely computational counter-selection against compounds with characteristics that are undesirable, was used by us previously to identify molecules that inhibited P-gp catalysis, but that were not transport substrates of P-gp.sup.26-28. Counter selections such as these, used to eliminate from consideration compounds with undesirable target interactions, can be extended to any property of a molecule that is calculable. When coupled with virtual synthesis of hit variants, increasing the efficiency and cost-effectiveness of synthesis programs is practically assured.
[0133] Ligand docking methods as described in.sup.26 have previously led us to identify the P-glycoprotein inhibitor, compound 29, that served as the initial hit for further drug development. Compound 29 and several other hits discovered were evaluated in biochemical and biophysical studies for their mechanism of inhibition of P-gp action.sup.26, as well as for their potential to reverse multidrug resistance in different cancer cell lines in culture.sup.27, 28 Compound 29 was chosen here as an initial compound for further development mostly for the fact that binding of compound 29 was predicted to be at an allosteric site, away from the nucleotide binding sites of P-gp.sup.26. Biophysical assessment using electron spin resonance spectroscopy and a spin-labeled ATP analog suggested that ATP binding was not affected in the presence of the inhibitor, while ATP hydrolysis assays showed inhibition of ATPase activity.sup.26. A putative mechanism for P-gp inhibition by 29 can be envisioned when comparing the position to which 29 docked with high affinity.sup.26 in
[0134] A closer evaluation of the high affinity allosteric docking site of compound 29 to P-gp revealed a relatively large hydrophobic pocket where the cyclopropyl moiety of the western half of the molecule interacted with the protein (
[0135] Five of the 29 variants from Table 1 (216, 227, 231, 541 and 551) were chosen for actual chemical synthesis mostly based on the perceived ease of synthesis and expense of precursor fragments. All three variants added more volume to the western half of the molecules and all but 541 had lower TPSA and higher log P values than the original compound 29. Closer inspection of the docking poses of the three variants (
[0136] Cell viability assays using a P-gp overexpressing prostate cancer cell line indicated that all five of the 29 variants had improved characteristics over the parental 29 for causing cell mortality (summarized in Table 2). This remarkable result, that five out of five Group 1ChemGen produced variants showed between 2.4-fold and 11-fold improvement over the performance of 29 in reversing P-glycoprotein-conveyed multidrug resistance phenotype (summarized in Table 2), underscores the utility of the ChemGen virtual synthesis approach and the employed docking selections for increasing affinity to P-gp. It is a reasonable conclusion that the larger Group 1 variants were able to interact more strongly with the protein as depicted in
[0137] In assays designed to allow quantification of the accumulation of the P-gp transport substrate, calcein AM, in cells that over-express P-gp, however, the larger, more hydrophobic Group 1 variants did not perform better than 29. Comparison of the calculated consensus log P values for variants 216, 227, 231 and 551 (6.6, 6.4, 5.7, and 5.9 respectively, Table 1) with the log P of the parental compound 29 (4.0, Table 1) led us to ask whether the lack of efficacy in inhibiting P-gp-catalyzed export of calcein AM may have been due to 29 variants being too hydrophobic to efficiently transfer across the cellular membrane to the cytosol-located nucleotide binding domains of P-gp for efficacious inhibition in the nucleotide binding domain of the protein to occur. The relatively short incubation times used in the calcein AM assays as shown in
[0138] Since the ChemGen/docking routine computational approach for selecting variants at the putative allosteric site resulted in larger but mostly also more hydrophobic western fragments in variants 216, 227, 231, 541 and 551 and since increased hydrophobicity may have limited utility of the compounds, the goal was to rationally design other variants of the western portion of compound 29 that might result in improved interactions with P-gp as well as more favorable physicochemical characteristics. In this light, five additional western derivatives of 29 were synthesized with varying shape, size and physicochemical properties, compounds 238, 255, 278, 280 and 286 (Group 2 variants). Again, the compounds were chosen for relative ease of synthesis as well as availability and expense of precursors. Unlike the previous set of ChemGen generated 29 variants (Group 1 variants), this latter set did not initially undergo the docking routines that selected for high affinity binding to the nucleotide binding domains and low affinity to the drug binding domains that were used to identify the parental compound 29 and the Group 1 variants. Instead these variants were chosen by visual inspection of the putative binding site, as well as by the calculated physicochemical properties of the resulting variants.
[0139] Efficacy of rationally designed Group 2 variants in reversing MDR caused by P-gp over-expression. If compounds 216, 227, 231 and 551 were too hydrophobic for favorable passage through the cellular membrane, thereby causing lowered efficacy of P-gp inhibition and decreased calcein AM accumulation in the P-gp overexpressing cells, then the more hydrophilic Group 2 derivatives, 238 and 255, as judged by their TPSA values and possibly 286 as judged by its log P value, might show increased re-sensitization of MDR cells to paclitaxel as well as improved calcein AM retention as compared to the other derivatives. These predictions were also based on the observed docking poses shown in
[0140] Effectiveness of the selective docking routines used on Group 1 variants to ensure targeting to the nucleotide binding domains and avoid the drug transport domains of P-gp. It was also of interest to investigate whether the subtractive docking routine performed on Group 1 variants added value to the overall process in selecting variants for synthesis and subsequent testing. Specifically, the docking selections employed were aimed at identifying compound 29 variants that preferentially bound to the nucleotide binding domains and were therefore not good transport substrates of P-gp. Comparison of the results from Group 1 compounds to the more traditional rationally designed Group 2 variants that were chosen without additional input on docking derived binding affinities to different substructures of the protein was therefore made. As can be seen in Table 5A, only compound 551 of the Group 1 variants selected through the ChemGen/docking routines stimulated basal ATP hydrolysis rates, and compound 231 inhibited basal ATPase activities. Intracellular accumulation assays which more directly measure whether a compound is a transport substrate of P-gp (see reference.sup.28) confirmed the predictions (Table 5,
[0141] On the other hand, intracellular accumulation data (summarized in Table 5) and stimulation of basal ATP hydrolysis activity of the mouse cysteine-less P-gp that was induced by Group 2 compounds 238, 255 and 286 (also summarized in Table 5) correlated quite well. Although the intracellular accumulation for compound 280 with or without addition of tariquidar were not observed to be significantly different, compound 280 did stimulate basal ATP hydrolysis in the absence of any other added transport substrate, which may indicate that this compound may also be a transport substrate of P-gp. However, the lack of stimulation of the human P-gp lets this conclusion be more doubtful. Compound 278 of Group 2 did not show differences in intracellular accumulation with or without tariquidar, but showed strong inhibition of both basal and verapamil-stimulated ATP hydrolysis of both orthologous enzymes. Compound 278 is therefore the only member of Group 2 that one can conclude is definitely not a transport substrate of P-gp, while 4 of the 5 Group 2 compounds either were transport substrates of P-gp or were potentially transport substrates of P-gp. In contrast to these conclusions, none of the five Group 1 molecules were transported by P-gp. Although the general case cannot be statistically proven since the number of synthesized and tested variants was too small, these studies show that the ChemGen produced variants that were selected against interactions at drug transporting structures on P-gp were much more likely not to be transport substrates of P-gp (five of five Group 1 variants were not observed to be transport substrates) than the Group 2 molecules that did not undergo the docking counter-selection procedures (where 4 out of 5 molecules were likely or potentially transport substrates for P-gp). It is very clear therefore that the docking and counter-selections were very effective in identifying 29 variants that were not transported by P-gp and that these procedures out-performed the rationally designed molecules that were not subjected to these selections.
[0142] Mechanism of action of the 29 variants. The effects of the 29 variants on verapamil-stimulated ATP hydrolysis rates catalyzed by P-gp (Table 5A and B) are consistent with the effects observed for the compounds in reversing MDR caused by P-gp (Tables 2 and 4). Group 1 compound 216 and Group 2 compound 286 had little effect in either the stimulated ATPase or the MDR reversal assays, likely indicating that these compounds do not strongly interact with P-gp under the conditions tested in these assays.
[0143] To assess whether the 29 variants also interacted with the nucleotide binding sites of P-glycoprotein, titration experiments using an electron spin resonance (ESR) active ATP analog, SL-ATP were performed, where the amount of P-gp bound SL-ANP was determined in the presence of the 29 variants (Table 5A). The results showed that none of the Group 1 variants significantly affected ATP binding while two of the Group 2 variants reduced ATP binding to about 1 mol per mol of protein. This again indicates that the selection process through ChemGen/selective docking was much more predictive of the effects the potential inhibitors have on the enzyme.
[0144] Among all compounds tested, the strongest reversal of MDR and the strongest stimulator of both basal and verapamil-stimulated ATPase activities was variant 238 of the Group 2 molecules. Compound 238, while strongly reversing MDR, was also the best transport substrate of all the variants tested, a characteristic that is deemed undesirable for any clinically relevant P-gp modulator lead as discussed above.
[0145] Group 2 compound 255 reversed MDR and moderately accelerated verapamil-stimulated ATPase rates, an effect that may be related to its role as a transport substrate of P-gp. Compound 280 affected MDR, ATP hydrolysis, as well as SL-ANP binding, but was not strong enough in any of these assays to warrant further study. Finally, Group 2 compound 278 affected MDR relatively weakly, but did inhibit both ATPase hydrolysis and SL-ANP binding to P-gp and therefore may warrant further study. Of the 10 new variants of the P-glycoprotein inhibitor compound 29 that were experimentally tested, compounds 227, 278, 541 and 551 inhibited verapamil-stimulated ATPase activity of the mouse cysteine-less P-gp similarly to that observed with the parental compound 29. However, the results using normal human P-gp reconstituted in a more native-like nanodisc environment indicate that relying on ATP hydrolysis assays, especially when not the same isoform of the enzyme is used and when the enzymes are in different environments (mixed micelles vs. nanodiscs), interpretations of results may not be as clear cut as often assumed.
[0146] It was found that the virtual synthesis of hit variants with a program suite like ChemGen combined with a novel selection of characteristics predicted from docking experiments, i.e., increased affinity to targeted structures and decreased affinity against sub-structures that should be avoided as employed here, resulted in efficient and cost-effective identification of five out of five variants assessed that met these goals, was demonstrated in every example tested so date. The comparison of physicochemical characteristics of the resulting new variants and the rational alteration of structure to investigate changed solubility properties (TPSA and log P) without the aid of computational prediction of desired properties, i.e. avoidance of transport structures on P-gp, led to the synthesis of molecules of which a majority was not targeted as desired, while every molecule that underwent the previous selection process inhibited as predicted. This is therefore a good example of how computational evaluation and selection of potential inhibitors before their actual synthesis adds to the speed and overall success rates in identifying hit-to-lead variants that possess desired characteristics.
[0147] Materials and methods.
[0148] Virtual synthesis of compound 29 derivatives using ChemGen. Because of the nature of the putative allosteric site shown for compound 29 and in light of rational design considerations, only a moderate set of variants were virtually synthesized using the ChemGen program. The ChemGen program and its use are described in detail herein below. In the virtual syntheses performed here, a scaffold molecule, 2-chloro-N-[1-phenyl-3-(2,4,5-trimethylphenyl)-1H-pyrazol-5-yl]acetamide, which is equivalent to the chloroacetamide that retains the Eastern substituent group from compound 29 was reacted with approximately 650 thiol compounds obtained by simple structural searches for thiols from the clean drug-like commercially available molecule set at the ZINC database.sup.34. The 2-chloro-N-[1-phenyl-3-(2,4,5-trimethylphenyl)-1H-pyrazol-5-yl]acetamide scaffold and the thiol precursor molecules were marked for reaction in ChemGen as described for the second reaction shown in
[0149] In silico docking of compound 29 variants to a model of human P-gp. AutoDock Vina.sup.60 and AutoDock 4.2.sup.61-63 were used with a model of P-glycoprotein with drug binding domains open to the outside and nucleotide binding sites fully formed that was extracted from targeted molecular dynamics trajectories as described in.sup.21, 22. This conformation of P-gp is one that is very similar to the homologous Sav1866 crystal structure reported by Dawson and Locher.sup.64 and was the conformation with which compound 29 was originally identified.sup.26. Ligand docking was limited to a volume equivalent to 202426 .sup.3 centered on the putative allosteric site of P-gp (see
[0150] Calculation of physicochemical properties. The SWISS-ADME server at http://www.swissadme.ch/ was used for the calculation of the physicochemical properties of the compounds as discussed in.sup.44.
[0151] Imaging of P-glycoprotein and ligands. The VMD (Visual Molecular Dynamics) program suite was used extensively in this work for the analysis of structural data and for the presentation of visual images.sup.65 and included the SURF surface representation program.sup.66 as well as the pdb2pgr.sup.67, 68 and APBS.sup.69, 70 programs for electrostatic/solvation calculations.
[0152] Synthetic procedures: All synthetic procedures and analyses of products are provided herein below.
[0153] Cell lines and cell culture. The chemotherapeutic sensitive DU145 human prostate cancer cells.sup.52 as well as the multidrug resistant sub-line, DU145TXR.sup.33 were generous gifts from Dr. Evan Keller (University of Michigan, Ann Arbor, Mich.). The multidrug resistant DU145TXR was maintained under positive selection pressure by supplementing complete medium with 10 nM paclitaxel (Acros Organics, NJ). Both cell lines were maintained in complete media consisting of RPMI-1640 with L-glutamine, 10% fetal bovine serum (FBS; BioWest, Logan, Utah), 100 U/mL penicillin and 100 g/mL streptomycin in a humidified incubator at 37 C. and 5% CO.sub.2. The noncancerous human fetal lung cell line, HFL1.sup.51, was kindly provided by Dr. Robert Harrod (Southern Methodist University, Dallas, Tex.) and maintained in complete media consisting of F-12K with L-glutamine, 10% FBS (BioWest, Logan, Utah), 100 U/mL penicillin, and 100 g/mL streptomycin in a humidified incubator at 37 C. and 5% CO.sub.2. To promote attachment of HFL1 cells, growth surfaces were treated with 0.1 mg/mL rat tail collagen (BD Biosciences, Palo Alto, Calif.) in 0.02 N acetic acid for 10 min and rinsed with PBS prior to use. Cell culture materials were purchased from Corning Inc. (Corning, N.Y.) unless otherwise stated.
[0154] MTT cell viability assay. Cells were trypsinized from monolayers and seeded with 3000 cells in 150 L of complete medium in a 96 well plate. After 24 hours, cells were treated for 48 hours with paclitaxel (Acros Organics, NJ) and/or P-gp inhibitory compounds dissolved in DMSO, or DMSO controls. All additions were diluted into complete medium. After 48 hours of treatment, MTT assays were performed as described.sup.41 using 5 mg/mL of MTT (Acros Organics, NJ) solution prepared in PBS (137 mM NaCl, 2.7 mM KCl, 10 mM Na.sub.2HPO.sub.4, 1.8 mM KH.sub.2PO.sub.4, pH 7.4). After 4 hours of incubation with MTT, the media was removed and the formazan crystals were dissolved in 100 L of DMSO. The absorbance at 570 nm was then measured using a BioTek Cytation 5 imaging multi-mode reader (Bio-Tek, Winooski, Vt.). Percent viability was calculated using DMSO treated cells as representative for 100% viability, according to Equation 1. Background absorbance was determined using MTT and complete medium without cells and subtracted from all the test values. Equation 1:
[0155] The results were plotted as the mean with standard deviation (SD) of eight replicates per concentration from at least two independent experiments with n=8. The graphical representations and IC.sub.50 values were determined using four parameter variable slope non-linear regression, using the following equation: Y=bottom+(top-bottom)/(1+10{circumflex over ()}((log IC50X)*Hill Slope) (GraphPad Prism, La Jolla Calif., USA, Version 6.05). The reported fold sensitization was calculated as follows, per Equation 2:
[0156] Calcein AM assay. To assess inhibition of P-gp-catalyzed transport of the P-gp pump substrate, calcein AM, DU145TXR cells were seeded in 96 wells plates and allowed to grow in complete medium until confluency was reached. Medium was removed, and cells were treated without 2 sM P-gp inhibitory compounds and 1 g/mL calcein AM (Life Technologies, OR) and diluted into phenol red free RPMI 1640 media. To study the effect of pre-incubation of compounds, cells were treated with just P-gp inhibitory compounds and incubated at 37 C. for six hours before adding calcein AM. Fluorescence excitation at 485 nm with a 20 nm gate and at emission at 535 nm with a 20 nm gate was measured using a BioTek Cytation 5 imaging multi-mode reader (Bio-Tek, Winooski, Vt.) over 60 minutes in 20 minute intervals. DMSO was used as vehicle. Results were plotted as the mean with standard deviation (SD) of six replicates per concentration and are representative of at least two independent experiments.
[0157] Cellular Accumulation Assays for Experimental P-gp Inhibitors. Cells used (DU145TXR), cell culturing, cell exposure to compounds, cellular handling and extractions were performed as described in.sup.28. LC-MS/MS methods were performed as described in.sup.71 and as modified in.sup.28.
[0158] P-glycoprotein Purification. Cysteine-less MDR3 and the human normal MDR1 P-glycoprotein was recombinantly expressed in the yeast Pichia pastoris essentially as in.sup.45, 46 and used for assaying ATP hydrolysis and ATP binding to the protein in the presence of the 29-variants. Purification of the protein was performed as described.sup.49 with small modifications resulting in highly enriched P-gp in mixed micelles containing dodecyl maltoside (DDM) and lysophosphatidyl choline.sup.26.
[0159] Nanodisc assembly. Human P-gp in mixed detergent micelles obtained during protein purification, was reconstituted into nanodiscsas described in references.sup.73, 74 with small modifications. P-gp was assembled with membrane scaffold protein, MSP1E3D1 (Sigma-Aldrich) expressed in BL21 (DE3) and L-alpha-phosphatidylcholine (Sigma-Aldrich), at a ratio of 1:10:500 (P-gp:MSP:PC) in 50 mM Tris-CL (pH 8). The detergent was removed with Bio-Beads SM-2 Adsorbent Media (BioRad). Ni-NTA Agarose (Qiagen) chromatography was used to purify P-gp reconstituted nanodiscs using 6 bed volumes of start buffer (20% (v/v) glycerol, 50 mM Tris-CL pH 7.5 at 4 C., 50 mM NaCl), and 5 bed volumes of elution buffer (20% (v/v) glycerol, 50 mM Tris-Cl pH 7.5 at 4 C., 50 mM NaCl, 300 mM imidazole).
[0160] ATPase Activity Assays. ATP hydrolysis activity was measured using a coupled enzyme assay.sup.72 as modified in reference.sup.26. The specific basal activity of the mouse MDR3 cysteineless P-gp was between 20 and 30 nmol min.sup.1mg.sup.1, and transport substrate (verapamil) stimulated activity was 200-400 nmol min.sup.1mg.sup.1. The specific basal activity of normal human MDR1 P-gp was between 123 and 193 nmol min.sup.1 mg.sup.1, and transport substrate (verapamil) stimulated activity was between 193-263 nmol min.sup.1mg.sup.1.
[0161] ESR Measurements. ESR measurements were as described in.sup.26. The amount of protein-bound spin-labeled (SL) adenine nucleotide was determined as the difference between the known total concentration of SL-ATP (2,3-(2,2,5,5,-tetramethyl-3-pyrroline-1-oxyl-3-carboxylic acid ester) ATP.sup.48) added and the free spin-labeled nucleotide (SL-ANP) observed in the experiment. Hyperbolic curve fitting of the results was performed using GraphPad Prism 7 to determine maximum binding and apparent affinity for the spin-labeled nucleotide. The equation used for the fitting the curves was y=P*x/(P2+x), where P1 corresponds to the maximum binding of SL-ANP (moles of SL-ANP bound per mole P-gp), and P2 equals the apparent dissociation constant for SL-ANP. To quantify the amount of free SL-ANP, standard curves were established where the signal amplitude of the high field signal of the ESR spectra of free SL-ANP in the absence of protein was correlated to the concentration of SL-ANP added. All ESR measurements were performed using a Bruker EMX 6/1 ESR spectrophotometer operating in X-band mode and equipped with a high sensitivity cavity. Spectra were acquired at a microwave frequency of 9.33 GHz, microwave power of 12.63 mW, 100 kHz modulation frequency and a resolution of 1024 points. The centerfield of the scan was at 3325 G and an area of 100 G was scanned. The peak to peak modulation amplitude was 1G and the time constant was set to 10.240 ms. The conversion time was 163.84 ms, resulting in a total time sweep of 167.772 s. The signal gain was adjusted for the SL-ATP concentrations in the different experiments.
[0162] Synthetic Procedures. General Materials and Methods.
[0163] The reactions were performed under nitrogen and dried glassware. Reagents were purchased from Sigma-Aldrich (St. Louis, Mo.), Alfa Aesar (Ward Hill, Mass.), EMD Millipore (Billerica, Mass.), Oakwood Chemical (West Columbia, S.C.), and Cayman Chemical (Ann Arbor, Mich.). Silica gel P60 (SiliCycle) was used for column chromatography and Analytical Chromatography TLC Silica gel 60 F.sub.254 (Merck Millipore, Darmstadt, Germany) was used for analytical thin layer chromatography. .sup.1H NMR and .sup.13C NMR spectra were used for analyzed the compounds by using CDCl.sub.3 (Cambridge Isotope Laboratories, Cambridge, Mass.) on a JEOL 500 MHz and BRUKER 400 MHz spectrometer in the Department of Chemistry at Southern Methodist University. Chemical abbreviations are used as follows: CH.sub.2Cl.sub.2, dichloromethane; EtOAc, ethyl acetate; THF, tetrahydrofuran; DMF, dimethylformamide; H.sub.2O, water; HBTU, O-benzotriazole-N,N,N,N-tetramethyl-uronium-hexafluoro-phosphate; DIPEA, N,N-diisopropylethylamine; KOH potassium hydroxide; DMSO, dimethylsulfoxide N.sub.2, nitrogen. High resolution mass spectroscopy was performed on a Shimadzu IT-TOF (ESI source) and low resolution mass spectroscopy was performed on a Shimadzu LCMS-8050 Triple Quadrupole LCMS (ESI source) or a Shimadzu Matrix Assisted Laser Desorption/Ionization MS (MALDI) at the Shimadzu Center for Advanced Analytical Chemistry at the University of Texas, Arlington.
[0164]
[0165]
[0166]
[0167]
[0168]
[0169]
General Procedure for the Synthesis of 2-Chloro-Acetamide Derivatives
[0170] Each substituted amine (1.0 equiv) was dissolved in 3 mL of anhydrous THF, followed by addition of one equivalent of Et.sub.3N. After placing the reaction mixture in an ice bath, 2-chloroacetyl chloride (1.2 equiv) was added dropwise for one hour and the reaction was stirred overnight at room temperature. After being concentrated, CH.sub.2Cl.sub.2 and water were added and the organic compounds were extracted three times with CH.sub.2Cl.sub.2. The organic layers were washed with brine, dried over Na.sub.2SO.sub.4, filtered, and concentrated.
[0171]
[0172]
[0173]
[0174]
[0175]
[0176]
[0177]
[0178]
[0179] General Synthesis for S.sub.N2 Coupling of Alkyl Thiols.
[0180] The reaction was performed by dissolving the thiol (1 equiv) in 3 mL of DMF (deoxygenated by bubbling N.sub.2) and adding K.sub.2CO.sub.3 (2 equiv). The chloroacetamide (1.2 equiv) was then added to the reaction mixture and stirred overnight at room temperature. The reaction was diluted in EtOAc and washed with water. The water layer was extracted three times with EtOAc, washed with brine, dried over Na.sub.2SO.sub.4, filtered and concentrated to give the crude product, which was purified as indicated.
[0181]
[0182]
[0183]
[0184]
[0185]
[0186]
[0187]
[0188]
[0189] Synthetic Procedures for Aromatic Sulfide Derivatives.
[0190]
[0191]
[0192]
[0193]
[0194]
[0195] It will be understood that particular embodiments described herein are shown by way of illustration and not as limitations of the invention. The principal features of this invention can be employed in various embodiments without departing from the scope of the invention. Those skilled in the art will recognize, or be able to ascertain using no more than routine experimentation, numerous equivalents to the specific procedures described herein. Such equivalents are considered to be within the scope of this invention and are covered by the claims.
[0196] All publications and patent applications mentioned in the specification are indicative of the level of skill of those skilled in the art to which this invention pertains. All publications and patent applications are herein incorporated by reference to the same extent as if each individual publication or patent application was specifically and individually indicated to be incorporated by reference.
[0197] The use of the word a or an when used in conjunction with the term comprising in the claims and/or the specification may mean one, but it is also consistent with the meaning of one or more, at least one, and one or more than one. The use of the term or in the claims is used to mean and/or unless explicitly indicated to refer to alternatives only or the alternatives are mutually exclusive, although the disclosure supports a definition that refers to only alternatives and and/or. Throughout this application, the term about is used to indicate that a value includes the inherent variation of error for the device, the method being employed to determine the value, or the variation that exists among the study subjects.
[0198] As used in this specification and claim(s), the words comprising (and any form of comprising, such as comprise and comprises), having (and any form of having, such as have and has), including (and any form of including, such as includes and include) or containing (and any form of containing, such as contains and contain) are inclusive or open-ended and do not exclude additional, unrecited elements or method steps. In embodiments of any of the compositions and methods provided herein, comprising may be replaced with consisting essentially of or consisting of. As used herein, the phrase consisting essentially of requires the specified integer(s) or steps as well as those that do not materially affect the character or function of the claimed invention. As used herein, the term consisting is used to indicate the presence of the recited integer (e.g., a feature, an element, a characteristic, a property, a method/process step, or a limitation) or group of integers (e.g., feature(s), element(s), characteristic(s), property(ies), method/process(s) steps, or limitation(s)) only.
[0199] The term or combinations thereof as used herein refers to all permutations and combinations of the listed items preceding the term. For example, A, B, C, or combinations thereof is intended to include at least one of: A, B, C, AB, AC, BC, or ABC, and if order is important in a particular context, also BA, CA, CB, CBA, BCA, ACB, BAC, or CAB. Continuing with this example, expressly included are combinations that contain repeats of one or more item or term, such as BB, AAA, AB, BBC, AAABCCCC, CBBAAA, CABABB, and so forth. The skilled artisan will understand that typically there is no limit on the number of items or terms in any combination, unless otherwise apparent from the context.
[0200] As used herein, words of approximation such as, without limitation, about, substantial or substantially refers to a condition that when so modified is understood to not necessarily be absolute or perfect but would be considered close enough to those of ordinary skill in the art to warrant designating the condition as being present. The extent to which the description may vary will depend on how great a change can be instituted and still have one of ordinary skill in the art recognize the modified feature as still having the required characteristics and capabilities of the unmodified feature. In general, but subject to the preceding discussion, a numerical value herein that is modified by a word of approximation such as about may vary from the stated value by at least 1, 2, 3, 4, 5, 6, 7, 10, 12 or 15%.
[0201] All of the devices and/or methods disclosed and claimed herein can be made and executed without undue experimentation in light of the present disclosure. While the devices and/or methods of this invention have been described in terms of particular embodiments, it will be apparent to those of skill in the art that variations may be applied to the compositions and/or methods and in the steps or in the sequence of steps of the method described herein without departing from the concept, spirit and scope of the invention. All such similar substitutes and modifications apparent to those skilled in the art are deemed to be within the spirit, scope, and concept of the invention as defined by the appended claims.
[0202] Furthermore, no limitations are intended to the details of construction or design herein shown, other than as described in the claims below. It is therefore evident that the particular embodiments disclosed above may be altered or modified and all such variations are considered within the scope and spirit of the disclosure. Accordingly, the protection sought herein is as set forth in the claims below.
[0203] Modifications, additions, or omissions may be made to the systems and apparatuses described herein without departing from the scope of the invention. The components of the systems and apparatuses may be integrated or separated. Moreover, the operations of the systems and apparatuses may be performed by more, fewer, or other components. The methods may include more, fewer, or other steps. Additionally, steps may be performed in any suitable order.
[0204] To aid the Patent Office, and any readers of any patent issued on this application in interpreting the claims appended hereto, applicants wish to note that they do not intend any of the appended claims to invoke 35 U.S.C. 112(f) as it exists on the date of filing hereof unless the words means for or step for are explicitly used in the particular claim.
REFERENCES
[0205] Brewer, F. K., Follit, C. A., Vogel, P. D., & Wise, J. G. (2014). In silico Screening for Inhibitors of P-Glycoprotein that Target the Nucleotide Binding Domains. Molecular pharmacology, 86(6), 716-726. [0206] Cramer, R. D., Patterson, D. E., Clark, R. D., Soltanshahi, F. & Lawless, M. S. Virtual compound libraries: a new approach to decision making in molecular discovery research. Journal of chemical information and computer sciences 38, 1010-1023 (1998). [0207] Durrant, J. D. & McCammon, J. A. AutoClickChem: click chemistry in silico. PLoS Comput Biol 8, e1002397 (2012). [0208] Follit, C. A., Brewer, F. K., Wise, J. G., & Vogel, P. D. (2015). In silico identified targeted inhibitors of P-glycoprotein overcome multidrug resistance in human cancer cells in culture. Pharmacol Res Perspect, 3(5), e00170. doi:10.1002/prp2.170 [0209] Humphrey, W., Dalke, A. & Schulten, K. VMD: visual molecular dynamics. J Mol Graph 14, 33-38 (1996). [0210] Ihlenfeldt, W. D., Voigt, J. H., Bienfait, B., Oellien, F. & Nicklaus, M. C. Enhanced CACTVS browser of the Open NCI Database. J Chem Inf Comput Sci 42, 46-57 (2002). [0211] Ihlenfeldt, W. D., Takahashi, Y., Abe, H. & Sasaki, S.-i. Computation and management of chemical properties in CACTVS: An extensible networked approach toward modularity and compatibility. Journal of chemical information and computer sciences 34, 109-116 (1994). [0212] Krier, M. et al. Design of small-sized libraries by combinatorial assembly of linkers and functional groups to a given scaffold: application to the structure-based optimization of a phosphodiesterase 4 inhibitor. Journal of medicinal chemistry 48, 3816-3822 (2005). [0213] Melnikov, A. A., Palyulin, V. A. & Zefirov, N. S. Generation of molecular graphs for QSAR studies: an approach based on supergraphs. Journal of chemical information and modeling 47, 2077-2088 (2007). [0214] Morris, G. M., R. Huey, W. Lindstrom, M. F. Sanner, R. K. Belew, D. S. Goodsell and A. J. Olson (2009). AutoDock4 and AutoDockTools4: Automated docking with selective receptor flexibility. J Comput Chem 30(16): 2785-2791. [0215] Nanayakkara, A. K., Follit, C. A., Chen, G., Williams, N. S., Vogel, P. D., & Wise, J. G. (2018). Targeted inhibitors of P-glycoprotein increase chemotherapeutic-induced mortality of multidrug resistant tumor cells. Sci Rep, 8(1), 967. [0216] O'Boyle, N. M. et al. Open Babel: An open chemical toolbox. Journal of cheminformatics 3, 33 (2011). [0217] Ousterhout, J. K., and Jones, K. (2010) Tcl and the Tk toolkit, 2nd ed., Addison Wesley, Upper Saddle River, N.J. [0218] Pearlman, R. S. & Smith, K. Novel software tools for chemical diversity. Perspectives in Drug Discovery and Design 9, 339-353 (1998). [0219] Srinivasan, R., Uttamchandani, M. & Yao, S. Q. Rapid assembly and in situ screening of bidentate inhibitors of protein tyrosine phosphatases. Organic letters 8, 713-716 (2006). [0220] Trott, O. and A. J. Olson (2010). AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem 31(2): 455-461. [0221] Wang, Z., H. Sun, X. Yao, D. Li, L. Xu, Y. Li, S. Tian and T. Hou (2016). Comprehensive evaluation of ten docking programs on a diverse set of protein-ligand complexes: the prediction accuracy of sampling power and scoring power. Physical Chemistry Chemical Physics 18(18):12964-12975.
EXAMPLE REFERENCES
[0222] [1] Ferreira, R. J., dos Santos, D. J., and Ferreira, M. J. (2015) P-glycoprotein and membrane roles in multidrug resistance, Future Med Chem 7, 929-946. [0223] [2] Sharom, F. J. (2008) ABC multidrug transporters: structure, function and role in chemoresistance, Pharmacogenomics 9, 105-127. [0224] [3] Jaramillo, A. C., Al Saig, F., Cloos, J., Jansen, G., and Peters, G. J. (2018) How to overcome ATP-binding cassette drug efflux transporter-mediated drug resistance?, Cancer Drug Resistance 1, 6-29. [0225] [4] Mao, Q., and Unadkat, J. D. (2015) Role of the breast cancer resistance protein (BCRP/ABCG2) in drug transportan update, The AAPS journal 17, 65-82. [0226] [5] Zhao, J. (2016) Cancer stem cells and chemoresistance: The smartest survives the raid, Pharmacol Ther 160, 145-158. [0227] [6] Wijdeven, R. H., Pang, B., Assaraf, Y. G., and Neefjes, J. (2016) Old drugs, novel ways out: Drug resistance toward cytotoxic chemotherapeutics, Drug Resist Updat 28, 65-81. [0228] [7] Fruci, D., Cho, W. C., Nobili, V., Locatelli, F., and Alisi, A. (2016) Drug Transporters and Multiple Drug Resistance in Pediatric Solid Tumors, Current drug metabolism 17, 308-316. [0229] [8] Tamaki, A., Ierano, C., Szakacs, G., Robey, R. W., and Bates, S. E. (2011) The controversial role of ABC transporters in clinical oncology, Essays in biochemistry 50, 209-232. [0230] [9] Palmeira, A., Sousa, E., Vasconcelos, M. H., and Pinto, M. M. (2012) Three decades of P-gp inhibitors: skimming through several generations and scaffolds, Current medicinal chemistry 19, 1946-2025. [0231] [10] Binkhathlan, Z., and Lavasanifar, A. (2013) P-glycoprotein inhibition as a therapeutic approach for overcoming multidrug resistance in cancer: current status and future perspectives, Current cancer drug targets 13, 326-346. [0232] [11] Kathawala, R. J., Gupta, P., Ashby, C. R., and Chen, Z.-S. (2015) The modulation of ABC transporter-mediated multidrug resistance in cancer: A review of the past decade, Drug Resistance Updates 18, 1-17. [0233] [12] Gottesman, M. M., and Pastan, I. H. (2015) The Role of Multidrug Resistance Efflux Pumps in Cancer: Revisiting a JNCI Publication Exploring Expression of the MDR1 (P-glycoprotein) Gene, J Natl Cancer Inst 107. [0234] [13] Marin, J. J., Briz, O., Rodriguez-Macias, G., Diez-Martin, J. L., and Macias, R. I. (2016) Role of drug transport and metabolism in the chemoresistance of acute myeloid leukemia, Blood Rev 30, 55-64. [0235] [14] List, A. F., Kopecky, K. J., Willman, C. L., Head, D. R., Persons, D. L., Slovak, M. L., Dorr, R., Karanes, C., Hynes, H. E., Doroshow, J. H., Shurafa, M., and Appelbaum, F. R. (2001) Benefit of cyclosporine modulation of drug resistance in patients with poor-risk acute myeloid leukemia: a Southwest Oncology Group study, Blood 98, 3212-3220. [0236] [15] Li, J., Jaimes, K. F., and Aller, S. G. (2014) Refined structures of mouse P-glycoprotein, Protein Sci 23, 34-46. [0237] [16] Kim, Y., and Chen, J. (2018) Molecular structure of human P-glycoprotein in the ATP-bound, outward-facing conformation, Science, eaar7389. [0238] [17] Pan, L., and Aller, S. G. (2015) Equilibrated atomic models of outward-facing P-glycoprotein and effect of ATP binding on structural dynamics, Scientific reports 5, 7880. [0239] [18] Jin, M. S., Oldham, M. L., Zhang, Q., and Chen, J. (2012) Crystal structure of the multidrug transporter P-glycoprotein from Caenorhabditis elegans, Nature 490, 566-569. [0240] [19] Taylor, N. M. I., Manolaridis, I., Jackson, S. M., Kowal, J., Stahlberg, H., and Locher, K. P. (2017) Structure of the human multidrug transporter ABCG2, Nature 546, 504-509. [0241] [20] Johnson, Z. L., and Chen, J. (2017) Structural Basis of Substrate Recognition by the Multidrug Resistance Protein MRP1, Cell 168, 1075-1085 e1079. [0242] [21] Wise, J. G. (2012) Catalytic transitions in the human MDR1 P-glycoprotein drug binding sites, Biochemistry 51, 5125-5141. [0243] [22] McCormick, J. W., Vogel, P. D., and Wise, J. G. (2015) Multiple Drug Transport Pathways through Human P-Glycoprotein, Biochemistry 54, 4374-4390. [0244] [23] Frank, G. A., Shukla, S., Rao, P., Borgnia, M. J., Bartesaghi, A., Merk, A., Mobin, A., Esser, L., Earl, L. A., and Gottesman, M. M. (2016) Cryo-EM analysis of the conformational landscape of human P-glycoprotein (ABCB1) during its catalytic cycle, Molecular pharmacology 90, 35-41. [0245] [24] Khunweeraphong, N., Stockner, T., and Kuchler, K. (2017) The structure of the human ABC transporter ABCG2 reveals a novel mechanism for drug extrusion, Scientfic reports 7, 13767. [0246] [25] Zoghbi, M. E., Mok, L., Swartz, D. J., Singh, A., Fendley, G. A., Urbatsch, I. L., and Altenberg, G. A. (2017) Substrate-induced conformational changes in the nucleotide-binding domains of lipid bilayer-associated P-glycoprotein during ATP hydrolysis, J Biol Chem 292, 20412-20424. [0247] [26] Brewer, F. K., Follit, C. A., Vogel, P. D., and Wise, J. G. (2014) In silico Screening for Inhibitors of P-Glycoprotein that Target the Nucleotide Binding Domains, Molecular pharmacology 86, 716-726. [0248] [27] Follit, C. A., Brewer, F. K., Wise, J. G., and Vogel, P. D. (2015) In silico identified targeted inhibitors of P-glycoprotein overcome multidrug resistance in human cancer cells in culture, Pharmacol Res Perspect 3, e00170. [0249] [28] Nanayakkara, A. K., Follit, C. A., Chen, G., Williams, N. S., Vogel, P. D., and Wise, J. G. (2018) Targeted inhibitors of P-glycoprotein increase chemotherapeutic-induced mortality of multidrug resistant tumor cells, Scientific Reports 8, 967. [0250] [29] Cramer, R. D., Patterson, D. E., Clark, R. D., Soltanshahi, F., and Lawless, M. S. (1998) Virtual compound libraries: a new approach to decision making in molecular discovery research, Journal of chemical information and computer sciences 38, 1010-1023. [0251] [30] Krier, M., de Arajo-Jnior, J. X., Schmitt, M., Duranton, J., Justiano-Basaran, H., Lugnier, C., Bourguignon, J.-J., and Rognan, D. (2005) Design of small-sized libraries by combinatorial assembly of linkers and functional groups to a given scaffold: application to the structure-based optimization of a phosphodiesterase 4 inhibitor, Journal of medicinal chemistry 48, 3816-3822. [0252] [31] Srinivasan, R., Uttamchandani, M., and Yao, S. Q. (2006) Rapid assembly and in situ screening of bidentate inhibitors of protein tyrosine phosphatases, Organic letters 8, 713-716. [0253] [32] Durrant, J. D., and McCammon, J. A. (2012) AutoClickChem: click chemistry in silico, PLoS computational biology 8, e1002397. [0254] [33] Takeda, M., Mizokami, A., Mamiya, K., Li, Y. Q., Zhang, J., Keller, E. T., and Namiki, M. (2007) The establishment of two paclitaxel-resistant prostate cancer cell lines and the mechanisms of paclitaxel resistance with two cell lines, Prostate 67, 955-967. [0255] [34] Irwin, J. J., Sterling, T., Mysinger, M. M., Bolstad, E. S., and Coleman, R. G. (2012) ZINC: A Free Tool to Discover Chemistry for Biology, Journal of chemical information and modeling 52, 1757-1768. [0256] [35] Moriarty, N. W., Grosse-Kunstleve, R. W., and Adams, P. D. (2009) electronic Ligand Builder and Optimization Workbench (eLBOW): a tool for ligand coordinate and restraint generation, Acta crystallographica. Section D, Biological crystallography 65, 1074-1080. [0257] [36] Yang, Z. H., Son, K. I., Li, S. Q., Zhou, B. N., and Xu, J. X. (2014) Specific 1,2-Hydride Shift in the Boron Trifluoride Catalyzed Reactions of Aromatic Aldehydes with Diazoacetonitrile: Simple Synthesis of beta-Ketonitriles, Eur J Org Chem, 6380-6384. [0258] [37] Su, W. N., Lin, T. P., Cheng, K. M., Sung, K. C., Lin, S. K., and Wong, F. F. (2010) An Efficient One-Pot Synthesis of N-(1,3-Diphenyl-H-Pyrazol-5-yl)amides, J Heterocyclic Chem 47, 831-837. [0259] [38] Hilton, S. T., Motherwell, W. B., Potier, P., Pradet, C., and Selwood, D. L. (2005) Observations on the reactivity of thiyl radicals derived from 3,6-epidithiodiketopiperazine-2,5-diones and related congeners, Bioorganic & Medicinal Chemistry Letters 15, 2239-2242. [0260] [39] Robertson, D. W., Krushinski, J. H., Beedle, E. E., Leander, J. D., Wong, D. T., and Rathbun, R. C. (1986) Structure-Activity-Relationships of (Arylalkyl)Imidazole AnticonvulsantsComparison of the (Fluorenylalkyl)Imidazoles with Nafimidone and Denzimol, Journal of Medicinal Chemistry 29, 1577-1586. [0261] [40] van Meerloo, J., Kaspers, G. J., and Cloos, J. (2011) Cell sensitivity assays: the MTT assay, Methods Mol Biol 731, 237-245. [0262] [41] Berridge, M. V., and Tan, A. S. (1993) Characterization of the cellular reduction of 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT): subcellular localization, substrate dependence, and involvement of mitochondrial electron transport in MTT reduction, Arch Biochem Biophys 303, 474-482. [0263] [42] Kupcsik, L. (2011) Estimation of cell number based on metabolic activity: the MTT reduction assay, Methods Mol Biol 740, 13-19. [0264] [43] Mannhold, R., Poda, G. I., Ostermann, C., and Tetko, I. V. (2009) Calculation of Molecular Lipophilicity: State-of-the-Art and Comparison of Log P Methods on More Than 96,000 Compounds, Journal of pharmaceutical sciences 98, 861-893. [0265] [44] Daina, A., Michielin, O., and Zoete, V. (2017) SwissADME: a free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules, Scientfic reports 7, 42717. [0266] [45] Lerner-Marmarosh, N., Gimi, K., Urbatsch, I. L., Gros, P., and Senior, A. E. (1999) Large scale purification of detergent-soluble P-glycoprotein from Pichia pastoris cells and characterization of nucleotide binding properties of wild-type, Walker A, and Walker B mutant proteins, J Biol Chem 274, 34711-34718. [0267] [46] Tombline, G., Urbatsch, I. L., Virk, N., Muharemagic, A., White, L. B., and Senior, A. E. (2006) Expression, purification, and characterization of cysteine-free mouse P-glycoprotein, Arch Biochem Biophys 445, 124-128. [0268] [47] Martin, C., Berridge, G., Mistry, P., Higgins, C., Charlton, P., and Callaghan, R. (1999) The molecular interaction of the high affinity reversal agent XR9576 with P-glycoprotein, Br J Pharmacol 128, 403-411. [0269] [48] Streckenbach, B., Schwarz, D., and Repke, K. R. (1980) Analysis of phosphoryl transfer mechanism and catalytic centre geometries of transport ATPase by means of spin-labelled ATP, Biochim Biophys Acta 601, 34-46. [0270] [49] Delannoy, S., Urbatsch, I. L., Tombline, G., Senior, A. E., and Vogel, P. D. (2005) Nucleotide binding to the multidrug resistance P-glycoprotein as studied by ESR spectroscopy, Biochemistry 44, 14010-14019. [0271] [50] Vogel, P. D., Nett, J. H., Sauer, H. E., Schmadel, K., Cross, R. L., and Trommer, W. E. (1992) Nucleotide binding sites on mitochondrial F1-ATPase. Electron spin resonance spectroscopy and photolabeling by azido-spin-labeled adenine nucleotides support an adenylate kinase-like orientation, J Biol Chem 267, 11982-11986. [0272] [51] Breul, S. D., Bradley, K. H., Hance, A. J., Schafer, M. P., Berg, R. A., and Crystal, R. G. (1980) Control of collagen production by human diploid lung fibroblasts, J Biol Chem 255, 5250-5260. [0273] [52] Stone, K. R., Mickey, D. D., Wunderli, H., Mickey, G. H., and Paulson, D. F. (1978) Isolation of a human prostate carcinoma cell line (DU 145), Int J Cancer 21, 274-281. [0274] [53] Morgan, S., Grootendorst, P., Lexchin, J., Cunningham, C., and Greyson, D. (2011) The cost of drug development: A systematic review, Health Policy 100, 4-17. [0275] [54] Merk, D., Friedrich, L., Grisoni, F., and Schneider, G. (2018) De novo design of bioactive small molecules by artificial intelligence, Molecular informatics 37, 1700153. [0276] [55] Chen, H., Engkvist, O., Wang, Y., Olivecrona, M., and Blaschke, T. (2018) The rise of deep learning in drug discovery, Drug discovery today. [0277] [56] Gomes, J., Ramsundar, B., Feinberg, E. N., and Pande, V. S. (2017) Atomic convolutional networks for predicting protein-ligand binding affinity, arXiv preprint arXiv:1703.10603. [0278] [57] Goh, G. B., Hodas, N. O., and Vishnu, A. (2017) Deep learning for computational chemistry, Journal of computational chemistry 38, 1291-1307. [0279] [58] Altae-Tran, H., Ramsundar, B., Pappu, A. S., and Pande, V. (2017) Low data drug discovery with one-shot learning, ACS central science 3, 283-293. [0280] [59] Subramanian, G., Ramsundar, B., Pande, V., and Denny, R. A. (2016) Computational modeling of p-secretase 1 (BACE-1) inhibitors using ligand based approaches, Journal of chemical information and modeling 56, 1936-1949. [0281] [60] Trott, O., and Olson, A. J. (2010) AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading, J Comput Chem 31, 455-461. [0282] [61] Goodsell, D. S., Morris, G. M., and Olson, A. J. (1996) Automated docking of flexible ligands: applications of AutoDock, J Mol Recognit 9, 1-5. [0283] [62] Morris, G. M., Goodsell, D. S., Huey, R., and Olson, A. J. (1996) Distributed automated docking of flexible ligands to proteins: parallel applications of AutoDock 2.4, J Comput Aided Mol Des 10, 293-304. [0284] [63] Osterberg, F., Morris, G. M., Sanner, M. F., Olson, A. J., and Goodsell, D. S. (2002) Automated docking to multiple target structures: incorporation of protein mobility and structural water heterogeneity in AutoDock, Proteins 46, 34-40. [0285] [64] Dawson, R. J. P., and Locher, K. P. (2006) Structure of a bacterial multidrug ABC transporter, Nature 443, 180-185. [0286] [65] Humphrey, W., Dalke, A., and Schulten, K. (1996) VMD: visual molecular dynamics, J Mol Graph 14, 33-38. [0287] [66] Varshney, A., Brooks, F. P., and Wright, W. V. (1994) Computing smooth molecular surfaces, IEEE Computer Graphics and Applications 14, 19-25. [0288] [67] Dolinsky, T. J., Nielsen, J. E., McCammon, J. A., and Baker, N. A. (2004) PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations, Nucleic acids research 32, W665-W667. [0289] [68] Dolinsky, T. J., Czodrowski, P., Li, H., Nielsen, J. E., Jensen, J. H., Klebe, G., and Baker, N. A. (2007) PDB2PQR: expanding and upgrading automated preparation of biomolecular structures for molecular simulations, Nucleic acids research 35, W522-W525. [0290] [69] Baker, N. A., Sept, D., Joseph, S., Holst, M. J., and McCammon, J. A. (2001) Electrostatics of nanosystems: application to microtubules and the ribosome, Proceedings of the National Academy of Sciences 98, 10037-10041. [0291] [70] Jurrus, E., Engel, D., Star, K., Monson, K., Brandi, J., Felberg, L. E., Brookes, D. H., Wilson, L., Chen, J., and Liles, K. (2018) Improvements to the APBS biomolecular solvation software suite, Protein Science 27, 112-128. [0292] [71] Wu, C.-Y., Feng, Y., Cardenas, E. R., Williams, N., Floreancig, P. E., De Brabander, J. K., and Roth, M. G. (2012) Studies toward the unique pederin family member psymberin: structure-activity relationships, biochemical studies, and genetics identify the mode-of-action of psymberin, Journal of the American Chemical Society 134, 18998-19003. [0293] [72] Vogel, G., and Steinhart, R. (1976) ATPase of Escherichia coli: purification, dissociation, and reconstitution of the active complex from the isolated subunits, Biochemistry 15, 208-216. [0294] [73] Denisov, I. G., Grinkova, Y. V., Lazarides, A. A., and Sligar, S. G. (2004) Directed self-assembly of monodisperse phospholipid bilayer Nanodiscs with controlled size, J Am Chem Soc 126, 3477-3487. [0295] [74] Ritchie, T. K., Grinkova, Y. V., Bayburt, T. H., Denisov, I. G., Zolnerciks, J. K., Atkins, W. M., and Sligar, S. G. (2009) Chapter 11Reconstitution of membrane proteins in phospholipid bilayer nanodiscs, Methods Enzymol 464, 211-231.
SUPPORTING INFORMATION REFERENCES
[0296] [1] Humphrey, W., Dalke, A., and Schulten, K. (1996) VMD: visual molecular dynamics, J Mol Graph 14, 33-38. [0297] [2] O'Boyle, N. M., Banck, M., James, C. A., Morley, C., Vandermeersch, T., and Hutchison, G. R. (2011) Open Babel: An open chemical toolbox, Journal of cheminformatics 3, 33. [0298] [3] Ihlenfeldt, W. D., Takahashi, Y., Abe, H., and Sasaki, S.-i. (1994) Computation and management of chemical properties in CACTVS: An extensible networked approach toward modularity and compatibility, Journal of chemical information and computer sciences 34, 109-116. [0299] [4] Ihlenfeldt, W.-D., Voigt, J. H., Bienfait, B., Oellien, F., and Nicklaus, M. C. (2002) Enhanced CACTVS browser of the Open NCI Database, J Chem Inf Comput Sci 42, 46-57. [0300] [5] Ousterhout, J. K., and Jones, K. (2010) Tcl and the Tk toolkit, 2nd ed., Addison Wesley, Upper Saddle River, N.J. [0301] [6] Moriarty, N. W., Grosse-Kunstleve, R. W., and Adams, P. D. (2009) electronic Ligand Builder and Optimization Workbench (eLBOW): a tool for ligand coordinate and restraint generation, Acta crystallographica. Section D, Biological crystallography 65, 1074-1080. [0302] [7] Wang, J., Wang, W., Kollman, P. A., and Case, D. A. (2001) Antechamber: an accessory software package for molecular mechanical calculations, J. Am. Chem. Soc 222, U403. [0303] [8] Schmidt, M. W., Baldridge, K. K., Boatz, J. A., Elbert, S. T., Gordon, M. S., Jensen, J. H., Koseki, S., Matsunaga, N., Nguyen, K. A., Su, S. J., Windus, T. L., Dupuis, M., and Montgomery, J. A. (1993) General Atomic and Molecular Electronic-Structure System, Journal of Computational Chemistry 14, 1347-1363. [0304] [9] Nanayakkara, A. K., Follit, C. A., Chen, G., Williams, N. S., Vogel, P. D., and Wise, J. G. (2018) Targeted inhibitors of P-glycoprotein increase chemotherapeutic-induced mortality of multidrug resistant tumor cells, Scientific reports 8, 967. [0305] [10] Takeda, M., Mizokami, A., Mamiya, K., Li, Y. Q., Zhang, J., Keller, E. T., and Namiki, M. (2007) The establishment of two paclitaxel-resistant prostate cancer cell lines and the mechanisms of paclitaxel resistance with two cell lines, Prostate 67, 955-967. [0306] [11] Yang, Z. H., Son, K. I., Li, S. Q., Zhou, B. N., and Xu, J. X. (2014) Specific 1,2-Hydride Shift in the Boron Trifluoride Catalyzed Reactions of Aromatic Aldehydes with Diazoacetonitrile: Simple Synthesis of beta-Ketonitriles, Eur J Org Chem, 6380-6384. [0307] [12] Su, W. N., Lin, T. P., Cheng, K. M., Sung, K. C., Lin, S. K., and Wong, F. F. (2010) An Efficient One-Pot Synthesis of N-(1,3-Diphenyl-H-Pyrazol-5-yl)amides, J Heterocyclic Chem 47, 831-837. [0308] [13] Saravanan, G., Selvaraju, R., and Nagarajan, S. (2012) Synthesis of Novel 2-Iminothiazolidin-4-Ones, Synthetic Commun 42, 3361-3367. [0309] [14] Marcu, A., Schurigt, U., Muller, K., Moll, H., Krauth-Siegel, R. L., and Prinz, H. (2016) Inhibitory effect of phenothiazine- and phenoxazine-derived chloroacetamides on Leishmania major growth and Trypanosoma brucei trypanothione reductase, European Journal of Medicinal Chemistry 108, 436-443. [0310] [15] Zhao, P.-L., Chen, P., Li, Q., Hu, M.-J., Diao, P.-C., Pan, E.-S., and You, W.-W. (2016) Design, synthesis and biological evaluation of novel 3-alkylsulfanyl-4-amino-1,2,4-triazole derivatives, Bioorganic & Medicinal Chemistry Letters 26, 3679-3683. [0311] [16] Heil, M. L., Cosford, N. D., Ardecky, R., and Zou, J. (2017) Quinazolinone analogs and use of quinazolinone analogs for treating or preventing certain viral infections, (Office, U. S. P. a. T., Ed.), Sanford-Burnham Prebys Medical Discovery Institute, Southern Research Institute, US. [0312] [17] Toscan, C. E., Rahimi, M., Bhadbhade, M., Pickford, R., McAlpine, S. R., and Lock, R. B. (2015) Thioimidazoline based compounds reverse glucocorticoid resistance in human acute lymphoblastic leukemia xenografts, Organic & biomolecular chemistry 13, 6299-6312. [0313] [18] Beisel, T., and Manolikakes, G. (2013) A Lewis Acid Palladium(II)-Catalyzed Three-Component Synthesis of alpha-Substituted Amides, Organic Letters 15, 6046-6049. [0314] [19] Li, Y.-l., Qi, X.-y., Jiang, H., Deng, X.-d., Dong, Y.-p., Ding, T.-b., Zhou, L., Men, P., Chu, Y., Wang, R.-x., Jiang, X.-c., and Ye, D.-y. (2015) Discovery, synthesis and biological evaluation of 2-(4-(N-phenethylsulfamoyl)phenoxy)acetamides (SAPAs) as novel sphingomyelin synthase 1 inhibitors, Bioorganic & Medicinal Chemistry 23, 6173-6184. [0315] [20] Robertson, D. W., Krushinski, J. H., Beedle, E. E., Leander, J. D., Wong, D. T., and Rathbun, R. C. (1986) Structure-Activity-Relationships of (Arylalkyl)Imidazole AnticonvulsantsComparison of the (Fluorenylalkyl)Imidazoles with Nafimidone and Denzimol, Journal of Medicinal Chemistry 29, 1577-1586. [0316] [21] Gianella-Borradori, M., Christou, I., Bataille, C. J., Cross, R. L., Wynne, G. M., Greaves, D. R., Russell, A. J. (2016) Ligand-based virtual screening identifies a family of selective cannabinoid receptor 2 agonists. Biorganic & Medicinal Chemistry 23, 241-263. [0317] [22] Tao, C., Wang, Q., Nallan, L., Polat, T., Koroniak, L. and Desai, N., California Capital Equity LLC, 2012. Triazine Derivatives and their Therapeutical Applications. U.S. patent application Ser. No. 13/376,964.