The present invention relates generally to bioinformatics methods for predicting efficient protein inhibitors, compounds synthesized based on the methods and uses thereof. More specifically, the present invention relates to bioinformatics methods for predicting BRD2 inhibitors, BRD2 inhibitor compounds and uses thereof.
Small molecules that influence the epigenetic status of cells have attracted attention as potential anticancer therapeutic agents, and many of these compounds are in preclinical and clinical trials (Dawson and Kouzarides 2012). Among the most promising targets of epigenetic drugs are bromodomain extra terminal (BET) proteins (Belkina and Denis 2012; Prinjha, Witherington, and Lee 2012).
The mammalian BET protein family includes ubiquitously expressed Brd2, Brd3, Brd4, and testis-specific BrdT. BET proteins contain two N-terminal bromodomains (BDs; BD1 and BD2) and the extra terminal (ET) protein interaction domain. Tandem BDs recognize acetylated lysines at histones, and this recognition interaction might recruit BET proteins to chromatin.
A distinct feature of BET proteins is that their attachment to chromatin persists during mitosis, which might help to reactivate postmitotic transcription and suggests a role in epigenetic memory (Garcia-Gutierrez, Mundi, and Garcia-Dominguez 2012). BDs can selectively distinguish specific acetylated lysines; for example, BD2 of Brd2 interacts preferentially with H4K12ac or H4K5ac/K12ac, whereas BD1 of the same protein recognizes H4K5ac/K8ac (Umehara et al. 2010; Filippakopoulos and Knapp 2012).
In addition, Brd2 BDs have much higher affinity to H4K5ac/K8ac/K12ac/K16ac tetraacetylated peptides than to single acetylation marks (Filippakopoulos et al. 2012; LeRoy et al. 2012). In contrast to diverse BD targets, ET domains, which are highly conserved among individual BET proteins (>80% identity among Brd2, Brd3, and Brd4), interact with the same proteins, for example, NSD3 methyltransferase or JMJD6 lysyl hydroxylase (Webby et al. 2009; Rahman et al. 2011).
BET proteins are found in multiple protein complexes; for example, Brd2 associates with RNA polymerase II and TATA box binding protein (TBP)-associated factors (TAFs), histone acetylases (CBP/p300), histone deacetylases, and chromatin assembly and remodeling factors (Denis et al. 2006). The BET proteins are transcription regulators, and the role of Brd3 and Brd4 in transcription is quite well established. Brd3 specifically binds the GATA1 transcription factor and enhances expression of GATA1-dependent genes (Gamsjaeger et al. 2011; Lamonica et al. 2011). Brd4 stimulates transcription via recruitment of active P-TEFb and induces phosphorylation of the C-terminal domain of RNA polymerase II (Hargreaves, Horng, and Medzhitov 2009; Devaiah et al. 2012).
However, less is known about the transcription regulation by Brd2. Brd2 associates with E2F transcription factors and helps to recruit TBP to promoters (Denis et al. 2006; Peng et al. 2006). In vitro, Brd2 enhances the passage of RNA polymerase II through acetylated chromatin within the gene body (LeRoy, Rickards, and Flint 2008). The importance of Brd2 is underscored by several experiments showing that Brd2 is an essential gene. Brd2−/− mice display defects in the development of the neural tube and abnormal brain structures and die during embryogenesis (Gyuris et al. 2009; Shang et al. 2009). Heterozygous Brd2+/− mice are viable but have a low number of inhibitory neurons in brain and develop spontaneous seizures (Velíšek et al. 2011). In addition, several studies link the BRD2 gene to epilepsy in humans (Pal et al. 2003; Cavalleri et al. 2007). Of interest, mice with lowered Brd2 expression are extremely obese but are protected from type 2 diabetes (Wang et al. 2009).
WO2015035112A1 describes cancer therapies, which combine epigenetic modulating agent(s) with immune modulating agent(s) These agents were identified to provide an improved treatment regimen over single agent therapy. In particular embodiments, the invention provides for improved treatment of NSCLC in patients via administration of exemplary immune modulating agents anti-PD-1 antibody or anti-PD-L1 antibody, which were observed to show enhanced activity in combination with the exemplary epigenetic modulating agent 5-deoxyazacytidine. Further, expression markers of responsive neoplastic cells are also disclosed.
There still remains a need for improved methods for identifying candidate agents for treating diseases.
There further remains a need for improved methods for identifying candidate agents and combinations of agents for treating proliferative diseases.
Additionally, there remains a need for improved agents and combinations of agents for treating proliferative diseases.
It is an object of some aspects of the present invention to provide bioinformatics methods for predicting efficient protein inhibitors.
It is another object of some aspects of the present invention to provide bioinformatics methods for predicting efficient BRD2 protein inhibitors.
It is yet another object of some aspects of the present invention to provide efficient protein inhibitors.
It is yet a further object of some aspects of the present invention to provide efficient BRD2 protein inhibitors.
It is yet another object of some aspects of the present invention to synthesize efficient protein inhibitors.
It is yet a further object of some aspects of the present invention to synthesize efficient BRD2 protein inhibitors.
In some embodiments of the present invention, improved methods for identifying candidate agents for treating diseases are provided.
In some further embodiments of the present invention, improved methods for identifying candidate agents and combinations of agents for treating proliferative diseases are provided.
In some additional embodiments of the present invention, improved agents and combinations of agents for treating diseases are provided.
In some further additional embodiments of the present invention, improved agents and combinations of agents for treating proliferative diseases are provided.
In some embodiments of the present invention, improved bioinformatics methods are provided for predicting efficient protein inhibitors.
In some further embodiments of the present invention, improved bioinformatics methods are provided for predicting BRD2 protein inhibitors.
In other preferred embodiments of the present invention, a method and system is described for providing a BRD2 protein inhibitor.
In some embodiments of the present invention, a BRD2 protein inhibitor structure is predicted.
In some embodiments of the present invention, a BRD2 protein inhibitor structure is defined.
In additional embodiments of the present invention, a method of synthesis of a BRD2 protein inhibitor structure is provided.
In further embodiments of the present invention, methods of structure-based drug design are provided to find a new pharmaceutical compound that binds to a macromolecular receptor. Traditionally, design campaigns are carried out experimentally, a process that is both time-consuming and costly. Therefore, computational methods have been developed with the aim of reducing the time and cost to design new drugs. The choice of method depends on at which stage in the campaign it will be employed and there is typically a trade-off between accuracy and efficiency.
In the present invention, three main approaches for new drug searches are combined, namely:
a) a pharmacophore-based search; and
b) a docking method
The combination methods of the present invention are applied to rapid screening of relatively large collections of compounds.
The present invention provides bioinformatics methods for predicting efficient protein inhibitors, compounds synthesized based on the methods and uses thereof, the method including identifying a BRD2 inhibitor by performing a pharmacophore-based search to produce a first plurality of candidate agents, applying a docking method to the first plurality of candidate agents second plurality of candidate agents to provide at least three BRD2 inhibitors which are 2-acetyl-N-benzyl-4-(2-methanesulfonylphenyl)-1-methyl-1H-imidazole-5-carboxamide, 5-acetyl-N-benzyl-4-ethyl-2-(2-methanesulfonylphenyl)-1H-pyrrole-3-carboxamide, 1-[3-(2-hydroxyethyl)-5-(2-methanesulfonylphenyl)-4-[2-(pyridin-4-yl)ethynyl]-1H-pyrrol-2-yl]ethan-1-one configured to inhibit the BRD2 protein.
There is thus provided according to an embodiment of the present invention, a method for identifying a BRD2 inhibitor, the method including;
a. performing a pharmacophore-based search to produce a first plurality of candidate agents; and
b. applying a docking method to the first plurality of candidate agents to output a second plurality of candidate agents.
Notably, according to an embodiment of the present invention, the pharmacophore search was applied with at least one of the following rules;
Furthermore, according to an embodiment of the present invention, the pharmacophore search was applied with all of the following rules;
Additionally, according to an embodiment of the present invention, the at least one candidate agent is listed in Table 1.
Furthermore, according to an embodiment of the present invention, the at least one candidate agent is PDB 4UYG.
Moreover, according to an embodiment of the present invention, the at least one candidate agent is a compound of the formulae;
a derivative thereof, a salt thereof, a derivative salt thereof and combinations thereof.
Further, according to an embodiment of the present invention, the method identifies compounds including:
Furthermore, according to an embodiment of the present invention, the compound has at least one of the following properties;
Importantly, according to an embodiment of the present invention, the compound has all of the following properties;
There is thus provided according to an embodiment of the present invention a compound according to any of following formulae:
Furthermore, according to an embodiment of the present invention, the compound is a BRD2 inhibitor.
Moreover, according to an embodiment of the present invention, the compound includes
Furthermore, according to an embodiment of the present invention, the compound has at least one of the following properties;
Importantly, according to an embodiment of the present invention the compound has all of the following properties;
Additionally, according to an embodiment of the present invention, the compounds described herein are suitable for administration to a subject to treat any disease or condition for which BET bromodomain inhibition provides clinical benefits.
Furthermore, according to an embodiment of the present invention, the disease or condition is selected from cancer, autoimmune disorder, respiratory disorder, asthma, cardiovascular disorder, neurological disease, Alzheimer's disease, inflammatory condition, sepsis and infection.
Moreover, according to an embodiment of the present invention, the compound is applied in combination with other therapeutic agents.
Further, according to an embodiment of the present invention, the compound, as described herein is used as immunotherapy agent to boost the patient's immune response against tumors or infection.
Yet further, according to an embodiment of the present invention the compound, as described herein, is administered to a patient in combination with immunotherapy agent(s).
Additionally, according to an embodiment of the present invention the immunotherapy agent(s) include one or many compounds or antibodies targeting PD-1, PD-L1 or CTLA4.
The present invention will be more fully understood from the following detailed description of the preferred embodiments thereof, taken together with the drawings.
The invention will now be described in connection with certain preferred embodiments with reference to the following illustrative figures so that it may be more fully understood.
With specific reference now to the figures in detail, it is stressed that the particulars shown are by way of example and for purposes of illustrative discussion of the preferred embodiments of the present invention only and are presented in the cause of providing what is believed to be the most useful and readily understood description of the principles and conceptual aspects of the invention. In this regard, no attempt is made to show structural details of the invention in more detail than is necessary for a fundamental understanding of the invention, the description taken with the drawings making apparent to those skilled in the art how the several forms of the invention may be embodied in practice.
In the drawings:
In all the figures similar reference numerals identify similar parts.
In the detailed description, numerous specific details are set forth in order to provide a thorough understanding of the invention. However, it will be understood by those skilled in the art that these are specific embodiments and that the present invention may be practiced also in different ways that embody the characterizing features of the invention as described and claimed herein.
A structure-based design approach was used employing three-dimensional structural modeling of BRD2 protein active site and its binding relationship with inhibitor compounds to design the inventive compounds described herein. The crystal structures of the BRD2 protein were obtained from http://www.rcsb.org/pdb/home/home.do. The modeled protein molecules first were docked with their inhibitors with docking software.
The crystal structures of the BRD2 protein obtained from Protein Data Bank with natural ligand and inhibitors were used as templates for pharmacophore search via screening chemical databases. The active site of the BRD2 protein is well-known and conservative (Filippakopoulos et al., 2010), therefore all the structures for the search were chosen according to the its binding pocket and Lipinski rule of 5 (Al-Lazikani, 2004).
Docking into the site of binding allows calculating different energy contributions for ligand-receptor system. For example, on the basis of docking chemicals in the active site, it was determined that certain compound classes (4-[(2,2-disubstituted 2H-1,3-benzodioxol-5-yl)substituted]-6-ethane-1,2,5,7,8-pentasubstituted 1,2,3,4-tetrahydroquinoline) could replace acetylated lysine thereby allowing new hydrogen bonding and hydrophobic interactions within the pocket.
Docking procedure was provided with the usage of option of flexible chains (5 {acute over (Å)}) and the center for binding for each ligand was chosen according to the coordinates of active site of the protein. However, the possibility of allostery inhibition of the BRD2 protein was not taken into account.
Crystal structure of BRD2 protein was taken from open access Protein Data Bank database (PDB code 4UYG). A pharmacophore search was performed according to the following rules:
1. No more than 5 hydrogen bond donors (the total number of nitrogen-hydrogen and oxygen-hydrogen bonds).
2. No more than 10 hydrogen bond acceptors (all nitrogen or oxygen atoms).
3. A molecular mass less than 500 daltons.
4. An octanol-water partition coefficient log P not greater than 5.
5. A polar surface area no greater than 140 {acute over (Å)}2.
Pharmacophores were chosen as it's shown in
As a result, 16 hits were identified and 9 of them were chosen for the future analysis of the structures (Table 1). A docking procedure was provided as it was described before and molecules—potential inhibitors—as a template for the further design (highlighted in the table 2.1 with “+”). As the reference structure the one from the crystal structure was used and docked it into the binding site. The inhibitor has hydrogen bond with Asn429, Other interactions are mostly hydrophobic: with Met438, Pro371, Val435, Trp370, His433, Cys425 and Leu383. Also water molecule is coordinated inside the binding site forming hydrogen bonds with ligand and Tyr386.
As it can be seen from Table 1, the lowest energy has inhibitor from the crystal structure because protein structure is well optimized around it.
To identify new chemicals that satisfy these structural requirements, a de novo design approach was employed. Initially, a core of the molecule was chosen as structural analogs of PCID AC1NCAKO, PCID 01137466. Below, there is the list of non-limiting examples of possible derivatives of a candidate molecule.
Then, all the molecules obtained were filtered with the following parameters:
As a result of the methods of the present invention, more than 119000 structures (combinations) were obtained.
Representative compounds of the present invention include:
The following embodiments provide non-limiting examples of the present invention, which should not be deemed as limiting:
1. A method for identifying a BRD2 inhibitor, the method comprising:
a. performing a pharmacophore-based search to produce a first plurality of candidate agents; and
b. applying a docking method to said first plurality of candidate agents to output a second plurality of candidate agents.
2. A method according to embodiment 1, wherein said pharmacophore-based search comprises using a template selected from a crystal structures of a BRD2 protein, a natural ligand of a BRD2 protein and an inhibitor of a BRD2 protein.
3. A method according to embodiment 1, wherein said docking method comprises applying ligand-protein visual analysis software to said first plurality of candidate agents.
4. A method according to embodiment 1, wherein said docking method comprises applying docking software to said first plurality of candidate agents.
5. A method according to embodiment 2, wherein said pharmacophore search comprises at least one of the following rules:
a) No more than 5 hydrogen bond donors (the total number of nitrogen-hydrogen and oxygen-hydrogen bonds);
b) No more than 10 hydrogen bond acceptors (all nitrogen or oxygen atoms);
c) A molecular mass less than 500 daltons;
d) An octanol-water partition coefficient log P not greater than 5; and
e) A polar surface area no greater than 140 {acute over (Å)}2.
6. A method according to embodiment 5, wherein said pharmacophore search comprises 2, 3, 4 or all of the following rules:
a) No more than 5 hydrogen bond donors (the total number of nitrogen-hydrogen and oxygen-hydrogen bonds);
b) No more than 10 hydrogen bond acceptors (all nitrogen or oxygen atoms);
c) A molecular mass less than 500 daltons;
d) An octanol-water partition coefficient log P not greater than 5; and
e) A polar surface area no greater than 140 {acute over (Å)}2.
7. A method according to embodiment 5, wherein said at least one candidate agent is listed in Table 1.
8. A method according to embodiment 7, wherein said at least one candidate agent is PDB 4UYG.
9. A method according to embodiment 5, wherein said at least one candidate agent is a compound selected from one or more of the formulae:
a derivative thereof, a salt thereof, a derivative salt thereof and combinations thereof.
10. A method according to embodiment 9, wherein
b)—X: Aromatic, Heteroaromatic ring;
c)—W: alkyl, heteroalkyl, cyclic alkyl;
d)—A: C or N or H
e)—D: H, alkyl, heteroalkyl
f)—P, Q: H, alkyl, cycloalkyl, aryl, alkyl or heteroaryl
g)—J: alkyl aryl, alkyl heteroaryl, unsaturated alkyl and
h)—R: H, alkyl.
11. A method according to embodiment 9, wherein said compound has at least one of the following properties:
i. A molecular mass of less than 500;
i. has less than 5 H-bond donors; and
ii. has less than 10 H-bond acceptors.
12. A method according to embodiment 9, wherein said compound has two or all of the following properties:
ii. A molecular mass of less than 500;
iii. has less than 5 H-bond donors; and
iv. has less than 10 H-bond acceptors.
13. A compound of any one or more of the following formulae:
14. A compound according to embodiment 13, wherein said compound is a BRD2 inhibitor.
15. A compound according to embodiment 13, wherein
16. A compound according to embodiment 13, wherein said compound has at least one of the following properties:
i. A molecular mass of less than 500;
ii. has less than 5 H-bond donors; and
iii. has less than 10 H-bond acceptors.
17. A compound according to embodiment 13, wherein said compound has two or all of the following properties:
i. A molecular mass of less than 500;
ii. has less than 5 H-bond donors; and
iii. has less than 10 H-bond acceptors.
18. A method wherein any compound of embodiments 13-17 is administered to treat any disease or condition for which BET bromodomain inhibition provides clinical benefits.
19. A method of embodiment 18, wherein a disease or condition is cancer, autoimmune disorder, respiratory disorder, asthma, cardiovascular disorder, neurological disease, Alzheimer's disease, inflammatory condition, sepsis or infection.
20. A method of embodiment 19, wherein a compound of embodiments 13-19 is applied in combination with other therapeutic agents.
21. A method of embodiments 19 or 20, wherein a compound of any of embodiments 13-17 is used as immunotherapy agent to boost the patient's immune response against tumors or infection.
22. A method of embodiments 18-21, wherein a compound of embodiments 13-17 is administered to a patient in combination with immunotherapy agent(s).
23. A method of embodiment 22, wherein immunotherapy agent(s) include one or many compounds or antibodies targeting PD-1, PD-L1 or CTLA4.
The references cited herein teach many principles that are applicable to the present invention. Therefore the full contents of these publications are incorporated by reference herein where appropriate for teachings of additional or alternative details, features and/or technical background.
It is to be understood that the invention is not limited in its application to the details set forth in the description contained herein or illustrated in the drawings. The invention is capable of other embodiments and of being practiced and carried out in various ways. Those skilled in the art will readily appreciate that various modifications and changes can be applied to the embodiments of the invention as hereinbefore described without departing from its scope, defined in and by the appended claims.
Number | Date | Country | |
---|---|---|---|
62456877 | Feb 2017 | US |