Molecular Docking Study Characterization of Rare Flavonoids at the Nac-Binding Site of the First Bromodomain of BRD4 (BRD4 BD1)

ε-N-Acetylation of lysine residues (Kac) is one of the most frequently occurring posttranslational modifications (PTMs) which control gene transcription and a vast array of diverse cellular functions. Bromodomains are epigenetic regulators involved in posttranslationalmodification. In silico docking studies were carried out to evaluate the binding potential of selected rare flavonoids on to Nac binding site of BD1 domain of BRD4 BET family proteins. Rare flavonoids like 3-O-acetylpinobanksin, naringenin triacetate, and kaempferol tetraacetate were found to occupy the WPF shelf and at the same time they exhibited a better binding affinity with multiple crystal structures of first bromodomain BRD4 (BRD4 BD1) when compared with the known inhibitors.


Introduction
The recent discovery of promising small molecule inhibitors for a class of nonenzymatic chromatin regulators, the BET bromodomains, suggests that future drug discovery for epigenetic therapy will involve the modulation of protein-protein interactions and multiprotein complexes.Also, it is expected that target-based discovery strategies will be increasingly complemented by approaches based on chemical probes generated by phenotypic or mechanistic cell based screening [1].
The bromodomain (BD) and extra terminal (BET) proteins, which comprise four members in human viz., BRD2, BRD3, BRD4, and the testis-specific BRDT [2] with each containing two N-terminal Bromodomains (BD).BRD4 and BRD2 are key mediators of transcriptional elongation by recruiting the positive transcription elongation factor complex (P-TEFb).The P-TEFb core complex is composed of cyclin-dependent kinase-9 (CDK9) and its activator cyclin T. CDK9 phosphorylates the RNA polymerase II (RNAPII) Cterminal domain, a region that contains 52 heptad repeats.RNAPII undergoes sequential phosphorylation at Ser5 during promoter clearance and at Ser2 by P-TEFb at the start of elongation.It has been shown that BRD4 couples P-TEFb to acetylated chromatin through its BRDs.Interestingly, in contrast to other BRD containing proteins and transcription factors, BET proteins remain associated with condensed and hypoacetylated mitotic chromosomes [3] suggesting a role in epigenetic memory [4,5].
Acetylation of lysine residues is a posttranslational modification with broad relevance to cellular signaling and disease biology.The principal readers of -N-acetyl-lysine (Kac) marks are bromodomains (BRDs), which are a diverse family of evolutionary conserved protein-interaction modules.The conserved BRD fold contains a deep, largely hydrophobic acetyl-lysine binding site, which represents an attractive pocket for the development of small pharmaceutically active molecules [6].
The bromodomain-containing protein-4 (BRD4) is a class of transcriptional regulators whose members are present in animals, plants, and fungi [7].The BET proteins typically have two tandem N-terminal bromodomains followed by an ET domain.As predicted by the presence of bromodomains, these proteins have been found to be associated with acetylated chromatin.BET proteins are involved in diverse cellular phenomena such as meiosis, cell-cycle control, and homeosis and have been suggested to modulate chromatin structure The emphasized bromodomain-1 (BD1) shows N-terminus, alpha-Z, ZA loop, alpha-A, AB loop, alpha-B, BC loop, and alpha C. The schematic diagram is based on as per the information from [49].
and affect transcription via a sequence-independent mechanism [8]. Figure 1 shows the schematic diagram of BRD4 with BD1, BD2 domains, and an extra terminal domain (ET).
Every first bromodomain BRD4 (BRD4 BD1) of BET has an N-terminus, alpha-Z helix, alpha-A helix, ZA loop, AB loop, alpha-B helix, BC loop, and alpha C helix.In these, the four identified distinct BET bromodomains have 80% similar amino acid sequences in BD1 domain and at the same time they conserve some amino acid residues.
-N-Acetylation of lysine residues (Kac) is one of the most frequently occurring posttranslational modifications (PTMs) in proteins [9].The secondary structure with helices and loops is shown in Figure 2 along with hydrophobic site for histone deacetylation (Kac site) which is projected.
Bromodomain inhibitors introduced so far are mostly BRD4 BD1 inhibitors.Among them, JQ1 is a BRD4 bromodomain inhibitor and it is highly active against human leukemia OCI-AML3 mutation lines such as nucleophosmin (NPM1) and DNA methyltransferase 3 (DNMT3A).JQ1 causes caspase 3/7-mediated apoptosis and DNA damage response in these cells.JQ1 prevented BRD4-mediated recruitment of p53 to chromatin targets following its activation in OCI-AML3 cells resulting in cell cycle arrest and apoptosis in a c-MYC-independent manner [10].GSK525762 binds to the acetylated lysine recognition motifs on the bromodomain of BET proteins, thereby preventing the interaction between the BET proteins and acetylated histone peptides.This disrupts chromatin remodeling and gene expression.Prevention of the expression of certain growth-promoting genes may lead to an inhibition of tumor cell growth.GW841819X displayed activity in vivo against NUTmidline carcinoma, multiple myeloma, mixed-lineage leukemia, and acute myeloid leukemia [11].Benzodiazepines and I-BET762 are BET bromodomain inhibitors now in clinical trial for nuclear protein in testis (NUT) midline carcinoma and other cancers [12][13][14].MS435 and MS436 effectively inhibit BRD4 activity in NF-B-directed production of nitric oxide and proinflammatory cytokine interleukin-6 in murine macrophages.MS435 and MS436 represent a new class of bromodomain inhibitors where MS436 is the best inhibitor [15].
X-ray crystallography shows that the chromen-4-one scaffold represents a new bromodomain pharmacophore and establishes LY294002 as a dual kinase and BET bromodomain inhibitor, whereas LY303511 exhibits anti-inflammatory and antiproliferative effects similar to the recently discovered BET inhibitors.LY294002 is a 2-morpholin-4-yl-7-phe-nyl-4Hchromen-4-one.Proteins BRD2, BRD3, and BRD4 comprise a family of targets structurally unrelated to PI3K [16].Figure 3 shows the structures of some known inhibitors of the first bromodomain of BRD4 (BRD4 BD1).
On this basis, some acetyl forms of flavonoids were selected and their binding affinity was analyzed.Selection of rare flavonoids is based on the structural aspects.Figure 4 shows the selective ligands under investigation.To support Lipinski's rule of five, Table 1 was listed out with chemical properties of ligands (selective rare flavonoids) like molecular formula, molar mass, Log P, and so forth.This study was aimed at understanding the differences in the presence of acetyl groups at various positions around the 1, 4-benzopyrone nucleus.In such perspective, acetyl form of naringenin, jaceidin, pinobanksin, kaempferol, artemetin, and hydroxy flavone were identified and carried over for in silico screening against bromodomains.The ligands were examined for Lipinski's rule of five.Rule of five is used to determine if a chemical compound, with a certain pharmacological or biological activity, has some properties that would make it a likely orally active drug in humans.It explains a molecule's pharmacokinetic properties, namely, absorption, distribution, metabolism, and excretion, but this rule lacks to explain whether the molecule would be pharmacologically active.In the discovery setting "the rule of 5" predicts that poor absorption or permeation is more likely when there are more than 5 hydrogen bond donors, 10 hydrogen bond acceptors, the molecular weight is greater than 500, and the calculated Log P (CLogP) is greater than 5 [17].According to the drug likeness properties, the ligands showed zero violation of the Lipinski rule of five, except kaempferol tetraacetate and jacedin triacetate.
In this study, the specific bromodomains BD1 domain of BRD4 BET family was selected.For validation, the specific pdb crystal structure with its cocrystallized ligand was redocked.Then a set of crystal structures with their cocrystallized ligands (known inhibitors) were evaluated for their effectiveness in binding.Multiple crystal structures of first bromodomain (BD1) of BRD4 BET family with different resolutions was downloaded from RCSB protein data bank and investigating ligands (rare flavonoids) with lowest binding energy was re-evaluated on different crystal structures.

Materials and Methods
The crystal structures of bromodomains used for evaluation in this study were obtained from RCSB protein data bank [18].Ligands were designed using Chemsketch [19,20] and   The chemical properties mentioned above for the ligands are available from the link: http://www.rsc.org.
designed ligands were optimized for their geometry with MOPAC [21] in Chem3D Ultra 9.0 platform [22], molecular docking was performed using LeDock [23], and overlapping analysis was done with PyMOL [24].Nonbond interactions of ligands with bromodomains were viewed through Accelrys' discovery studio 4.0 [25] and ligand map was generated using Molegro Molecular viewer 2.5 [26,27], Ligplot is generated using Ligplot plus 1.4.5 [28].MOPAC uses semiempirical theory to calculate the electronic part, in order to obtain molecular orbitals, the heat of formation, and its derivative with respect to molecular geometry [21].MOPAC in Chem3D Ultra 9.0 employs AM1 (Austin Model 1) method and closed shell wave function to minimize energy to minimum RMS gradient of 0.100.AM1, belong to the family of NDDO (Neglect of Diatomic Differential Overlap) methods.In these methods, all terms arising from the overlap of two atomic orbitals which are on different centers or atoms are set to zero.As this is not the forum for developing the ideas of Hartree-Fock theory [29], knowledge about the derivation of the Roothaan-Hall equations [30] will be assumed and description of the methods will start with the final equations.Properties like heat of formation, gradient norm, charges, cosmos ovation in water, electrostatic potential, molecular surfaces, spin density, hyperfine coupling constants, and polarizabilities were solved and, additionally, Mulliken charges were calculated.Mulliken charges [31] provide means of estimating partial atomic charges and are routinely used as variables in linear multiple regression QSAR procedures.3D structures obtained after MOPAC calculations were visualized through Accelrys' discovery Studio 4.0.

Docking Software.
LeDock (http://www.lephar.com/) is developed for fast, flexible, and accurate docking of small molecule ligands into a protein target with extremely simplified interface.LeDock utilizes a combination of simulated annealing and genetic algorithm for the optimization of position, orientation, and rotatable bonds of the ligand [23].Its built-in scoring function, adapted from a transferable and accurate energy function [32], mainly consists of Vdw interaction energy and a unique empirical-based hydrogen bonding energy.

First Bromodomain of BRD4 (BRD4 BD1
) Crystal Structure Energy Optimization.Different crystal structures of first bromodomain of BRD4 (BRD4 BD1) of BET family proteins were optimized for their geometry with Open Babel 2.3.2 [33,34].It uses a general Amber force field (GAFF) for energy minimization.GAFF is designed to be compatible with existing Amber force fields for proteins and nucleic acids and has parameters for most organic and pharmaceutical molecules that are composed of H, C, N, O, S, P, and halogens.It uses a simple functional form and a limited number of atom types but incorporates both empirical and heuristic models to estimate force constants and partial atomic charges [35].Optimized crystal structures were utilized for docking evaluation.

Validation of Computational Molecular Docking.
To know the accuracy of molecular docking, the cocrystallized ligand has to be retrieved and redocked onto the active site.Root-mean-square deviation (RMSD) between the docked structure and the original conformation of the inhibitor in each complex was calculated.Crystal structure (4CFK) with its cocrystallized ligand was selected and validation was performed.

Docking Methodology 2.5.1. Preparation of Protein Molecule Using Lepro.
Lepro is an advanced program in refining the protein molecule.It has a dual advantage in detecting the binding site and as well as removal of water molecules.Even though the binding site (the space search (grid) for the ligand) is recognized in the presence of a cocrystallized ligand, it conserves precision.The protein molecule is loaded and hydrogen is added prior to docking.Lepro writes refined protein molecule as a ".pdb" file (e.g., pro.pdb), but docking parameters as ".in" file (dock.in).
Pro.pdb and dock.in are the files needed in order for LeDock to carry out its docking evaluation.After adding hydrogen, the binding site was determined.

Docking of the Ligands Using LeDock.
Refined protein molecule (pro.pdb) and docking input files (dock.in)were uploaded before selecting a ligand.The ligand in a SYBYL MOL2 (ligand.mol2)format is an only ideal file that can be recognized by the LeDock.After constructing all necessary files, the docking is started using start "docking icon." The number of poses generated and their binding affinity (docking score) will be generated.Maximum of 20 poses will be generated depending upon torsions in the ligands.

Results and Discussion
As quoted above, LeDock is a built-in scoring function adapted from a transferable and accurate energy function mainly consisting of Vdw interaction energy and a unique empirical-based hydrogen bonding energy.The docking score is calculated as given in (1).Consider where Δ is binding energy; Δ  is interaction energy;  HB is hydrogen bonding penalty; , , and  are fitting parameters.
LeDock generates 20 different poses of each ligand depending upon torsions.There were no less than 6 poses on the above said criteria at the vicinity of the binding site for the investigated ligands 3.1.Validation of Molecular Docking.Validation of docking methodology was performed in order to measure the reliability of molecular docking.Root-mean-square deviation  (RMSD) between the docked structure and the original conformation of the inhibitor in each complex was calculated.
If the RMSD of the docked pose is less than or equal to 1.0 Å from the experimentally observed conformation, the prediction is regarded to be successful.LY294002 is a chromen-4-one compound exhibiting a docking score of −5.30 kcal/mol.Among the selective ligands, a prediction has been made to show their binding affinity towards gorge area.Figure 5 shows the docked conformation and original conformation of LY294002 with their respective docked energy.
In LY294002, the major functional group involved in hydrogen bond interaction was oxygen atom of ketonic (-C=O) group at C5 position and the rest of the aromatic rings interaction dealt with hydrophobic interaction.There was only a slight variation in occupying the WPF shelf, which can be evidenced from the ring shift and torsional shift of the redocked ligand from native ligand retrieved from the crystal structure.This reliability holds strong to perform the docking analysis with investigational ligands.
Asn140 donates its hydrogen to make a bond with oxygen atom of LY294002 at C5 position containing ketonic group.The former interaction of redocked ligand caused a shift of 0.15 Å when compared with the cocrystallized ligand.Table 2 shows the distance of H-bond formed by the cocrystallized ligand and the redocked ligand.The ring shift has maximum influence in deciding the position of redocked ligand.In such perspective ring shift was an appropriate consideration in the case of hydrophobic interaction.
Cocrystallized and redocked ligand conformation has its hydrophobic interaction in the same fashion and has limited ring shifts with slight torsion at the morpholinyl linkage.The hydrophobic interaction generated by the rings of LY294002 at the WPF shelf was directed at amino acids like Val87, Cys136, Leu92, Ile146, and Phe83.Table 3 shows the position (XYZ) of the ligand at the WPF shelf along with the type of hydrophobic bond.This explains the ring shift from the centre of gorge area with respect to the cocrystallized ligand.

Binding Affinity of Known Inhibitors and Investigational
Ligands.Except JQ1, the binding affinity of known inhibitors was found to be less in comparison with flavonoids under study.
The resolution of crystal structures used and their cocrystallized ligand has been listed out in Table 2 along with docked energy (Δ), amino acid residue interacted, and its distance in formation of hydrogen bond ( Å).
Repeated docking analysis showed the reliability in docking and the known inhibitors occupying the WPF shelf with maximum conformity.Among them JQ1 exhibited a docked energy of −6.52 kcal/mol.In this conformation, there is a well oriented H-bond interaction with Asn140.The Ligplot (Figure 6) shows the occupancy of the known inhibitors at WPF shelf.
Whereas, di-acetyl group substituted naringenin and 3-O-acetyl pinobanksin substituted with an extra acetyl group at C7 positionlag behind in scoring as −5.95 and −5.82 kcal/mol, respectively.
Native flavone moiety with differing hydroxy, methoxy, and acetyl substitution on C5, C7, and C8, namely, 5-hydroxy-7-acetoxy-8-methoxy flavone, exhibits a satisfactory binding  as the ZA and BC loops form the acetyl-lysine binding pocket located at one end of the helix bundle.The amino acid residues in the acetyl-lysine binding pocket are highly conserved with over 90% sequence identity between the two bromodomains in each BET protein.Of these is a highly conserved Asn residue that is essential for lysine-acetylated histone recognition by forming a hydrogen bond to the acetyl amide group of the acetylated lysine [15].BRD4 BD1 binds to the acetylated H3 tail polypeptide at residues 12-19.Specifically, the conserved Asn140 in BD1 binds to the N zeta -acetylated Lys14 in H3.Intermolecular hydrogen bonds are formed between the (O) of the acetylated Lys14 and Asn140 of the BC loop.Water molecules surrounding the acetyl moiety orientate it so that there are hydrogen bonds between the hydroxy group of the Asn side chain and the backbone amide of the Ala residue in the histone peptide.Additional contacts with the acetylated lysine were mediated by Tyr139 and Ile46 and with Cys139 in some BRD4 BD1 [36].
The crystal structure with pdb code 4CFK was selected for docking analysis with rare flavonoids.In those selected pdb's, there was a conserved intermolecular hydrogen bonding between the amino acid residues in protein structure that was observed.Important hydrogen bond interactions within the amino acid residues of bromodomain at the Kac site were Asn140 and Cys136.Focusing on crystal structure of 4CFK, it showed an N-terminus commencing with Ser42 and continue as alpha-A helix at Asn61.Helical structure extends up to His77, where it has a slight coil and lengthens up to Asp106.In between this helix, residues form an AB loop (Gln78-Met105) and continue as alpha-B helix to descend at Asn116.The successive residues Asn117-Ala122 form a BC loop and then from Gln123 extend as alpha C helix up to Asn162.
Validation of docking methodology was performed in order to measure the reliability of molecular docking.Rootmean-squared deviation (RMSD) between the docked structure and the original conformation of the inhibitor in each complex was calculated.The RMSD of the docked pose is less than 1.0 Å from the experimentally observed conformation; this exhibits a successful prediction.From the binding poses and superimposition of docked conformation of the ligand, the prediction was said to be accurate.
Selected flavonoids contain methoxy, hydroxyl, and acetyl groups.Binding characteristics at the acetylation gorge of BRD4 BD1 were evaluated in comparison with known inhibitors Selective rare flavonoids showed their interaction with the amino acids residues of predicted binding site of BC loop.Ligand map (Figures 7,11,and 14) shows the hydrogen bond and steric interactions of the ligand with the bromodomain.It is a planar ligand map where only one set of interactions are visible and it has been shown in following discussion.
The binding free energy was mainly accounted for by hydrophobic interactions, while electrostatic forces might contribute to the specificity of the interaction.Thus affinity constant is independent of the ionic strength [37].Hence hydrophobic interactions have been visualized and their importance is discussed throughout the binding interaction analysis.
C2 substituted (2R, 3R)-2-cyclohexyl moiety has its unidirectional hydrophobic interactions from Pro82 and a bidirectional force of attraction with Leu92 and Val87.Figure 9 shows benzopyrone nucleus under the influence of Val87 holding a hydrophobic interaction with Ile146 and Cys136.
There was a predicted binding mode of acetyl mimics in BRD4.It is characterized by a lipophilic sandwich of their bicyclic core between residues Val87, Leu92, Leu94, and Tyr97 on one side and Phe83 and Ile146 on the other side of the binding pocket [38].The carbonyl oxygen of both compounds is engaged in hydrogen bonding interactions with the highly conserved Asn140, a typical feature among bromodomain inhibitors as most of them are acetyl-lysine mimics [39][40][41][42][43][44][45].
An investigation on compounds like acetyl mimics was predicted to further occupy the hydrophobic WPF shelf, which is an important region for ligand design to gain potency [38].
Kaempferol tetraacetate is a 3,5-bis(acetyloxy)-2-[4-(acetyloxy) cyclohexyl]-4-oxo-4a,5,6,7,8,8a-hexahydro-4Hchromen-7-yl acetate.Acetyl group of 4-(acetyloxy) cyclohexyl moiety interacts at Tyr97 (2.50 Å) through hydrogen bonding and as well as with a steric interaction.In case of jaceidin triacetate, at the same position there is a void of interaction with expected Tyr97; this might be due to the presence of methoxy group at the adjacent carbon (C5  ).This screening effect would prohibit strain along the ring and make both the oxygen atoms at carbonyl oxygen (O-C=O) of acetyl group at C7 induce a force of nonbond hydrogen bond interaction at Asn140 (2.02 Å)and Steric interaction with Val87.Hydrophobic moments along the interaction field were found to have the same symmetry like 3-Oacetylpinobanksin but due to the position of acetyl groups their nonbond interaction varies considerably and elicits free energy of binding.This explains the strain created along the bonds due to methoxy groups and acetyl groups in the ligands at the gorge environment.
RVX-208 binds BRD4 BD1 in a mostly hydrophobic pocket where it is flanked by Val87, Leu92, and Leu94 side chains on one side of the compound and Pro82 and Leu146 side chain on the other.It is also within van der Waals distances to the four aromatic residues Tyr97, Tyr139, Phe83, and Trp81.RVX-208 is a quinazolin-4-one it forms a bidentate set of two hydrogen bond interactions with the conserved Asn140 amide side chain and the amino acid recognized by acetylated lysine of the natural substrate, as well as an extensive water network.The quinazolin-4-one scaffold ketone and methoxy oxygens both form a watermediated interaction with the Tyr97 hydroxyl group [47].
Diazobenzene inhibitors bind across the acetyl-lysine binding pocket in the ZA and BC loops in the BRD4 BD1.The hydroxyl group in the A ring forms a hydrogen bond to the amide of the conserved Asn140 and another water (W1) mediated hydrogen bond to the phenoxyl group for Tyr97 [15].Some set of structurally optimized flavonoids bearing a common flavone nucleus with acetyl, methoxy, and hydroxyl groups specifically at C5, C7, and C8 positions like 5-hydroxy-7-acetoxy-8-methoxy flavone, 5-acetoxy-7-hyroxy flavone, and 5,7-diacetoxy flavone were evaluated through docking for their affinity in accessing the lysine binding site of bromodomains.
Being a best scored ligand, 3-O-acetylpinobanksin was examined for redocking on multiple crystal structures with different resolution, namely, 3U5J, 3MXF, 4CFK, and 4C66 to know the reliability in their access towards -N-acetyllysine binding site.Table 5 shows their binding energy with its interacting amino acids and Figure 15 shows the ligand map of 3-O-acetylpinobanksin with different crystal structures.
The critical amino acid residue of importance is Asn140 AND Tyr97.On analysing the interaction of investigational ligands with 4CFK and specific interaction of 3-O-acetylpinobanksin with different crystal structures, there was a slight difference in the locations of the interacting (H-bonding) amino acid residues with the ligands.This was evidenced by the distance of hydrogen bonds formed between the ligand and the residues.3-O-Acetylpinobanksin and 5-hydroxyl-7-acetoxy-8-methoxy flavone shared a common distance of 1.94 Å whereas naringenin diacetate, 3,7-Oacetylpinobanksin made interaction at a distance of 1.84 Å.
There was a 0.10 Å difference in every set of interactions among the investigational ligands.On comparison, the specific interaction of 3-O-acetylpinobanksin with different crystal structures as shown in Table 5 also exhibited a 0.10 Å distance variation.
Overlapping analysis of investigational ligands with a reported known inhibitor, LY294002, with a benzopyran nucleus is done in order to know about the structural efficiency of rare flavonoids in occupying the binding site area.Then analysis is made on the necessary structural importance of rare flavonoids in bromodomain binding.Figure 16 shows the evidence for superimposition of 3-O-acetylpinobanksin and LY294002.
Flavonoids occupy the WPF shelf and nonbonding interaction was also found at these areas with the investigational It is concluded that selective rare flavonoids like 3-O-acetylpinobanksin, naringenin diacetate, and kaempferol tetraacetate were found to occupy the WPF shelf with maximum fidelity and with satisfactory binding energy (Δ).Substitution of acetyl, methoxy, and hydroxyl group on benzopyrone structures shows influence in binding interaction at Kac site.Most of the investigational ligands interacted at Asn140, Tyr97, Met132, and Lys91 through nonbonding interactions like steric interactions and hydrophobic interaction.In that perspective, these compounds would act as a better BRD4 BD1 inhibitor.Further in vitro and in vivo studies can be carried out in order to produce supporting information and can target the inhabitation of bromodomains, the epigenetic readers of lysine acetylation.

Figure 1 :
Figure 1: Schematic diagram of human BRD4 with two N-terminal bromodomains (BD1 and BD2) with an extra terminal (ET) domain.The emphasized bromodomain-1 (BD1) shows N-terminus, alpha-Z, ZA loop, alpha-A, AB loop, alpha-B, BC loop, and alpha C. The schematic diagram is based on as per the information from [49].

Figure 2 :Figure 3 :
Figure 2: Secondary structure of bromodomain BD1-BRD4 with water molecules showing N-terminus, -Z helix, -A helix, ZA loop, AB loop, -B helix, BC loop, and -C helix.Adjacent image with white arrows below shows a secondary structure and surface diagram which in turn emphasizes a circled area of hydrophobic sites.Diagram of secondary structure generated using Accelrys' discovery studio 4.0 [25].

2. 1 .
Ligand Preparation Using MOPAC.MOPAC (Molecular Orbital Package) is a semiempirical quantum chemistry program based on Dewar and Thiel's NDDO approximation.

Figure 5 :
Figure 5: Validation of molecular docking (LeDock).The magenta shows the native ligand (LY294002) pose retrieved from the crystal structure 4CFK and red colour shows the docked pose of the ligand (LY294002) using LeDock.(the docked files have been displayed using PyMOL and the retrieved ligand from cocrystallized structure in its conformation has been loaded to superimpose on each other).

Figure 6 :
Figure 6: Ligplot image showing the occupancy and interaction of known inhibitors with its crystal structures.The image is generated using Ligplot plus version 1.4.5.

Figure 8 :Figure 9 : 3 -
Figure 8: Emphasized image shows the binding interaction of 3-O-acetylpinobanksin at N-acetyl-lysine binding site of BD1 domain of BRD4 family protein.Image generated using Accelrys' discovery studio 4.0 [25].(a) Secondary structure showing the ligand at the WPF shelf with surface around cavity.(b) The cavity area is shown as a transparent surface with ligand.(c) The emphasized image showing the H-bond interaction (dark blue lines) and hydrophobic interaction (sky blue cylindrical lines) of amino acid residues [shown as sticks (0.20 Å) and lines, resp.] with the ligands.

Figure 12 :Figure 13 : 5 -
Figure 12: The hydrophobic moments showing the multiple nonbond interaction with artemetin acetate at the gorge area interactions.

Figure 16 :
Figure 16: Image displayed is showing the crystal structure of 4CFK containing its cocrystallized ligand (LY294002) overlapped onto the energy minimized 4CFK with docked 3-O-acetylpinobanksin.Red colour is LY294002 and red salmon 3-O-acetylpinobanksin; those colors mentioned are the same for protein structure and cocrystallized ligand [image is generated using PyMOL viewer].

Table 1 :
Table showing the list of ligands (with CAS number) used for investigation and their chemical properties.

Table 2 :
Validation set showing the docking score of known inhibitors onto the crystal structure of multiple BRD4 (BD1).

Table 3 :
Redocked ligand and cocrystallized ligand position at grid centre of crystal structure of BRD4 BD1.(The amino acid involved in hydrophobic interaction with LY294002 and its positioning with respect to the free energy of binding is shown according to its -axis.)

Table 4 :
Docking score with the hydrogen bond interaction of rare flavonoids onto the detected binding site in crystal structure of 4CFK.
affinity but the values of 5, 7-diacetoxy flavone, 5-acetoxy-7-hyroxy flavone show a reduction in docking score.It is evident from Table4which shows their docking score along with the dissociation constant.The top scored ligand (3-O-acetylpinobanksin) among 10 ligands was again subjected to redocking onto different crystal structures like 4C66, 4NUC, 3U5J, and 3MXF.Three of the crystal structures exhibited docking score (Table5) ≥ 6.40 kcal/mol.All bromodomains share a conserved left-handed helix bundle that is made out of four  helices, named Z, A, B, and C, respectively.The interhelical loop regions known