Advances in Bioinformatics

Advances in Bioinformatics / 2014 / Article

Research Article | Open Access

Volume 2014 |Article ID 104823 |

Ridla Bakri, Arli Aditya Parikesit, Cipta Prio Satriyanto, Djati Kerami, Usman Sumo Friend Tambunan, "Utilization of Boron Compounds for the Modification of Suberoyl Anilide Hydroxamic Acid as Inhibitor of Histone Deacetylase Class II Homo sapiens", Advances in Bioinformatics, vol. 2014, Article ID 104823, 10 pages, 2014.

Utilization of Boron Compounds for the Modification of Suberoyl Anilide Hydroxamic Acid as Inhibitor of Histone Deacetylase Class II Homo sapiens

Academic Editor: Antoine van Kampen
Received03 Jun 2014
Revised12 Aug 2014
Accepted12 Aug 2014
Published24 Aug 2014


Histone deacetylase (HDAC) has a critical function in regulating gene expression. The inhibition of HDAC has developed as an interesting anticancer research area that targets biological processes such as cell cycle, apoptosis, and cell differentiation. In this study, an HDAC inhibitor that is available commercially, suberoyl anilide hydroxamic acid (SAHA), has been modified to improve its efficacy and reduce the side effects of the compound. Hydrophobic cap and zinc-binding group of these compounds were substituted with boron-based compounds, whereas the linker region was substituted with p-aminobenzoic acid. The molecular docking analysis resulted in 8 ligands with Δ value more negative than the standards, SAHA and trichostatin A (TSA). That ligands were analyzed based on the nature of QSAR, pharmacological properties, and ADME-Tox. It is conducted to obtain a potent inhibitor of HDAC class II Homo sapiens. The screening process result gave one best ligand, Nova2 (513246-99-6), which was then further studied by molecular dynamics simulations.

1. Introduction

Cervical cancer is cause by human papillomavirus (HPV) and in the second rank as a cause of cancer death in women worldwide [1]. Cervical cancer occurs in the cervical region, which is located in the hollow area between the vagina and the uterus or commonly called cervix. Cervical cancer can be contagious among all women; a ratio of 1 out of every 4 women is likely to suffer from it [2].

Based on data from the World Health Organization, in 2008, it is estimated 530,232 cases of cervical cancer in the world, with 275,008 mortality cases [3]. Through these data, the estimated global mortality rate from cervical cancer is 50% [2].

HPV is a virus of the family Papillomaviridae and has a nonenveloped, icosahedral-shaped capsid and the double stranded circular DNA as its genetic material [46]. It is 7,800–7,900 base pairs long with a 55 nm diameter [7, 8]. HPV has more than 100 different genotypes, and over 40 types of it can infect any part of the epithelial and mucosal lining of the anogenital tissue [9]. The HPV virus is divided into two classes, namely, low-risk HPV (e.g., HPV-6 and HPV-11) and high-risk HPV (e.g., HPV 16 and HPV 18) [10]. Low-risk HPV usually causes a bulge impact on disease areas such as anogenital condylomata (wart) that usually grows on the cervix and vulva [11].

HPV genome is divided into 3 regions, namely, upstream regulatory (URR, noncoding), early gene, and late gene regions [12]. Proteins E6 and E7 oncogenes can make HPV-infected cells to become immortal [13]. E6 protein is associated with ubiquitin (protein ligase), which in turn interacts with p53. It results in the degradation process in the proteosome. E6 also increases the activity of telomerase and induces the creation of immortal cells [14]. E7 protein interacts with the retinoblastoma protein (Rb) and releases the E2F transcription factor that induces expression of genes involved in the process of cell proliferation [15]. E7 oncoprotein can interact directly with the interferon regulatory factor- (IRF-) 1 tumor suppressor proteins that inhibit the performance of the release of E2F and E7, thus increasing the transcriptional activity of cells containing the HPV genome [16]. E6 and E7 activities generally cause epigenetic changes that interfere with the process of cell regulation, apoptosis, DNA-repair processes, hormonal response, and cell differentiation processes, which lead to cervical cancer [17, 18].

Oncoprotein of HPV E6 and E7 in particular has a correlation with the enzyme activity of histone deacetylase (HDAC) [19]. HDAC is a medium for binding with oncogene transcription of genes with the aim of transforming the processes of cells into the media of the viral proliferation [20].

HDAC is an enzyme with EC number 3.5.1 which acts as a catalyst for histone deacetylase [21]. In eukaryotic cells, it is useful for removing acetyl groups from lysine amino acid on a histone tail and wrapping the histones around DNA, thus interfering with the process of gene transcription by binding with transcription factor [22, 23]. In general, there are two regulation processes of gene expression and DNA replication by regulation of chromatin structure [23]. The process of protein acetylation of histone and nonhistone was carried out by the histone acetyl transferases (HATs) and histone deacetylase by histone deacetylase (HDACs) enzyme [24]. These two enzymes are working as opposites because HATs will cause chromatin structure to stretch into euchromatin [25]. It provides space for the specific enzyme or other protein complexes involved in gene expression that serves to increase the activity of transcription and DNA repair [26]. While HDAC causes the release of an acetyl group on the N-acetyl lysine that is available on the histone tail, it causes the DNA to form loops on the histone called heterochromatin [27]. Hence, the transcription of DNA is obstructed and gene expression does not occur properly, thus causing the transformation of normal cells into cancer cells [17, 28]. HDAC inhibition can inhibit the proliferation of epigenetic gene transcription of HPV that causes cancer cells broke down to apoptosis [29].

Suberoyl anilide hydroxamic acid (SAHA) has been through the stages of clinical trials and approved by the U.S. Food and Drug Administration (FDA) as cancer drug [30]. SAHA has the ability to inhibit HDAC, and it could interact with the HDAC metalloenzyme site [31]. Hence, the Zn2+ ion lies at the basis of the metalloenzyme site of HDAC [32]. The following is an explanation of each unit in the design of HDAC inhibitors.

(1) Zinc Binding Group (ZBG). It is a site where a ligand (inhibitor) interacts with Zn2+ cofactor contained in the HDAC formed charge relays system with amino acid residues [33]. In general, compounds that can interact with the Zn2+ cofactor are nucleophilic compounds, for example, hydroxyl, carbonyl, thiol, carboxylic, and sulfonyl [34].

(2) Linker. It is the liaison between CAP and ZBG that forms a short-chain hydrocarbons, long-chain hydrocarbons or aromatic such as butane, fatty acids, γ-aminobutyric acid (GABA), p-aminobenzoic acid (PABA), furans, and others [35]. The linker is able to interact with amino acid residues found in the cylinder pocket of HDAC enzymes [36].

(3) Hydrophobic Cap (CAP). It is a group of compounds that are used to design a cap which is generally a hydrophobic compound that has properties of high lipophilicity [37]. It easily reacts with the surface of the active site and closes the entry point to the enzyme substrate. In general, hydrophobic cap is composed of phenyl, benzyl, furans, polycyclic, and so forth [37, 38].

The reason boron compounds are selected to be substituted at the ZBG and CAP is because the clinical trials have shown that consumption of boron may prevent cervical cancer caused by HPV [37, 39]. By consuming the boron content of 84.1 mg per day, it could prevent cervical cancer [40]. The forms of a functional group of boron compounds that have been proven to have therapeutic effects till date are diazoborin, boronic acid, boronic ester, and benzoxaborole [41]. Carborane has been found to be useful as a good inhibitor and has high lipophilicity properties which are useful for binding with the receptor binding site on the hydrophobic active site of the enzyme [42]. Carborane in the closed form can increase receptor affinity and activity of the enzyme with hydrophobic ligand binding cavity. Therefore, it can inhibit the enzyme activity that contributes to a disruption in the disease [43].

2. Materials and Methods

This research method was developed based on established pipeline of our group [4448]. This is an in silico research that involves the use of computerized system. Each query was done using the online and offline software. In this study, MOE 2012.10, ACDLabs, ChemSketch, Toxtree, and VegaZZ were used. Multiple sequence alignment was done in HDAC class II sequences of Homo sapiens to obtain the conserved region. Then, the HDAC class II Homo sapiens enzyme was computed for its homology modeling with SWISS-MODEL server. As a result, ligand inhibitors for HDAC class II Homo sapiens were produced. After both the ligand and the enzyme were ready, the molecular docking simulations were performed. The result of molecular docking simulation was forwarded to the analysis of the existing parameters, namely, pharmacological analysis, ADMET testing, and bioavailability. Further test was carried out to examine the thermodynamic stability of ligand in the presence of solvents with molecular dynamics simulations.

3. Results and Discussion

3.1. Results

The ZBG would be substituted with boronic acid and carborane, part of the linker would be substituted with p-aminobenzoic acid compounds, and hydrophobic parts of cap would be substituted with organoboron compounds from the website of organoborons database ( (Figure 1).

The sequences of HDAC class II Homo sapiens were searched in a protein sequence database. They could be accessed via the National Center for Biotechnology Information (NCBI). HDAC class II consists of six types of enzymes, namely, HDAC4, HDAC5, HDAC6, HDAC7, HDAC9, and HDAC10. The whole isoform of protein sequences has been encoded in the NCBI Reference Sequence (NCBI RefSeq), GenBank, and UniProt Knowledge Base (UniProtKB)/SWISS-PROT.

After conducting multiple sequence alignment, conserved region sequences were obtained. The obtained sequences of HDAC enzyme code are seen in Table 1. Furthermore, the sequences were piped into the Basic Local similarity Alignment Search Tools (BLAST) which could be accessed through the NCBI website ( BLAST is useful for comparing sequences derived from the conserved region of the existing protein database. The BLAST protein code is written in Table 1.

EnzymeSequence code


Furthermore, 3D structure modeling was conducted by using SWISS-MODEL server. The predicted 3D image crystal structure of HDAC class II Homo sapiens results can be seen in Figures 2(a) to 2(f).

The modeling of the data was also obtained in determining the active site of the enzyme from each HDAC class II Homo sapiens enzyme. The active site of each enzyme is listed in Table 2.

EnzymeActive site

HDAC4Zn2+ cofactor with Asp196, His198, and Asp290 residue
HDAC5Zn2+ cofactor with Asp870, His872, and Asp964
HDAC6Zn2+ cofactor with Cys5, His7, and Cys78
HDAC7Zn2+ cofactor with Asp707, His709, and Asp801
HDAC9Zn2+ cofactor with Asp820, His822, and Asp914
HDAC10Zn2+ cofactor with Asp172, His174, and Asp265

Furthermore, molecular docking simulation has been conducted and produced 8 best ligands. The obtained ligands were having Δ value lower than SAHA and trichostatin A (TSA) as standards; they are listed in Table 3. The low Δ value would facilitate the spontaneous ligand binding reaction with the enzyme to form stable complexes. Deriving from Δ value, the inhibition constants of complex formation are in Table 4. The lower the Δ value the greater the p value. Docking results also included visualization of the interaction between the ligand with the target enzyme. Ligands would interact with amino acid residues that were owned by the enzyme and also with the active site of the enzyme. Table 5 presents the interaction between multiple ligands with the enzyme.

LigandGibbs free energy (Δ) towards HDAC class II Homo sapiens

Nova2 (9058064-6)−25,83−28,94−29,43−31,33−25,24−24,01
Nova2 (95752-88-8)−18,25−31,67−28,87−28,14−26,74−25,96
Nova2 (88765-82-6)−23,63−29,22−28,20−25,41−24,20−23,32
Nova2 (Unique10)−23,79−29,14−27,58−26,66−24,39−20,90
Nova2 (16876-27-0)−16,96−27,72−25,88−24,31−24,80−22,84
Nova2 (513246-99-6)−22,69−28,57−27,01−26,72−24,92−21,75
Nova2 (Unique80)−22,14−25,00−26,80−24,51−23,43−20,83
Nova2 (279262-23-6)−22,13−27,46−26,04−27,09−24,25−20,68

Note: numbers in bold show the lowest binding energy.

LigandInhibition constant (p)

Nova2 (9058064-6)18,8221,0921,4422,8218,3917,49
Nova2 (95752-88-8)13,3023,0721,0320,5019,4818,92
Nova2 (88765-82-6)17,2121,2920,5418,5117,6316,99
Nova2 (Unique10)17,3321,2320,0919,4217,7715,23
Nova2 (16876-27-0)12,3620,2018,8617,7118,0716,64
Nova2 (513246-99-6)16,5320,8119,6819,4718,1515,85
Nova2 (Unique80)16,1318,2119,5317,8617,0715,18
Nova2 (279262-23-6)16,1320,0118,9719,7317,6715,07

LigandDocking interactions

Nova2 (9058064-6)Zn2+, Arg37, Trp195, Asp220Zn2+, His833, His872, Trp869, Leu892, Tyr937, Phe983Zn2+, His7, Leu8, Ala10Zn2+, Arg547, His843, Asp801, His709, Arg731, Trp706Zn2+, Arg606, Asp820, Tyr887, Ala911, Gly912, Asp914Zn2+, Arg30, Trp171, His174, Arg196, Phe254, Asp265

Nova2 (95752-88-8)Zn2+, Arg37Zn2+, His833, Asp870, His872, Tyr937, Trp869, Asp964Zn2+, His7, Val9Zn2+, Arg547, Tyr774, Gly799, Asp801, His843Zn2+, Arg660, Phe833, Gly912, Asp914Zn2+, Arg30, Trp171, Gly263, Phe282, Gly303, Gly304

Nova2 (88765-82-6)Zn2+, Asp290, Glu329, Trp195, Phe289Zn2+, His872, Arg710, Asp870, Asp964, Gly1005Zn2+, Arg 47, Ala80Zn2+, Arg547, Ala798, Phe800, Gly841, His843, Zn2+, Asp820, Leu842, Glu853, Asp914Zn2+, Asp172, Asp265

Nova2 (Unique10)Zn2+, Arg37, trp195, Ser287, Gly288Zn2+, His833, His872, Asp964, Leu973Zn2+Zn2+, Arg547, Asp707, Trp706, Asp801, Gly799, His843, Zn2+, Arg660, Gly791, Asp914Zn2+, Asp172, Tyr238, Asp265, Glu302

Nova2 (16876-27-0)Zn2+, Gly288Zn2+, Asp870, His872, Arg710Zn2+Zn2+, Arg547Zn2+, Arg660, Asp914Zn2+, Asp172, Asp265

Nova2 (513246-99-6)Zn2+, Ser287, Gly288, Asp290, Leu304Zn2+, His833, Asp870, His872, Gly962Zn2+, Ser72, Trp74, Tyr77Zn2+, His709, Gly799, Asp801, His843Zn2+, Arg660, Arg778, His822, Asp914, His956Zn2+, Trp171, Asp172, Tyr238, ASp265, Ser266

Nova2 (Unique80)Zn2+, Gly288, Asp290Zn2+, Arg710, Trp869, Ala961, Phe963Zn2+, Leu64, Typ74Zn2+, Val708, Gly812, His843Zn2+, His822, Asp914, Gly 954Zn2+, Trp171, Asp172, Trp228, Tyr238, Gly303

Nova2 (279262-23-6)Zn2+, Ser287, Asp290, Phe209, Ala327Zn2+, Arg710, Asp870, Gly962, Gly1005Zn2+, Val9Zn2+, Arg547, Trp706, Asp801, His843Zn2+, Arg660, Asp820, Asp914Zn2+, Asp172, Tyr238, Phe264, Asp265

SAHAZn2+, Asp196, His198, Asp290Zn2+, Asp870, His872, Asp964, Gly1005Zn2+, Cys5, His7, Cys75, Cys78Zn2+, Asp707, His709, Asp801, Tyr813Zn2+, Arg660, Gly954, Gly955Zn2+, Trp171, Asp172, Val173, His174, Asp265, Phe282

TSAZn2+, Asp196, His198, Asp290, Ser287Zn2+, Asp870, His872, Asp964, Arg710Zn2+, Cys75, Cys78Zn2+, Asp707, His709, Asp801Zn2+, Asp820, His822, Asp914, Gly955Zn2+, Asp172, Val173, His174, Asp265, Phe282

Furthermore, the previous modified ligands were forwarded for pharmacological analysis. The analysis was carried out by using Lipinski’s rule of five Egan’s, and Veber’s rules to determine the best drug candidates in its stability and oral bioavailability. According to these rules, the drugs should have a molecular weight of less than 500 Dalton (Da), Log values ​​of less than 5, the number of hydrogen bond donors of less than 5, the number of hydrogen bond acceptors of less than 10, polar surface area of less than 140 Å2, and rotation of the ligand binding compound of less than 10 [49, 50]. The test of pharmacological analysis was conducted using FAF-Drugs2 online software. Table 6 shows the results of pharmacological analysis of each ligand.

LigandMwTPSARot. bondHBDHBAViolation

Nova2 (9058064-6)425,63−1,52141,7510681
Nova2 (95752-88-8)465,742,4088,699370
Nova2 (88765-82-6)423,661,26102,689470
Nova2 (Unique10)423,702,4490,4612460
Nova2 (16876-27-0)409,631,1281,798270
Nova2 (513246-99-6)425,701,2893,708470
Nova2 (Unique80)471,792,5079,909270
Nova2 (279262-23-6)437,691,2890,909370

Note: Mw: molecular weight; : octanol-water partition coefficient; TPSA: topological polar surface area; Rot. bond: the amount of rotatable bond; HBD: the amount of hydrogen donor atom; HBA: the amount of hydrogen atom acceptor; violation: the amount of the violation to Lipinski’s rule.

As seen in Table 7, the data show that the best 8 ligands have good oral bioavailability with parameters based on the existing rules. Furthermore, an analysis of health impact of ligands has been completed by observing its absorption, distribution, metabolism, excretion, and toxicology (ADMET) properties. The analysis was carried out using Toxtree 2.6.0 software and ACD/I-Lab with the parameter of Benigni-Bossa rule. The method involved analyzing the groups of ligands that have fragments containing substances that cause mutagenic and carcinogenic effects in the cells of the body. Examples of the compounds that could cause these effects are acyl halides, benzyl, esters, epoxides, aliphatic halogen, alkyl nitrites, quinones, hydrazine, polycyclic aromatic hydrocarbons, tiocarbamate, aromatic amines, hydroxylamine, and so forth. Analysis of Toxtree software assessment was generated from data obtained by test using Salmonella typhimurium organism or commonly referred to as the AMES test. Table 8 shows the results of ligands that have been analyzed with the Toxtree software. The result is the best 8 ligands did not have mutagenic and carcinogenic effects.

LigandEgan’s ruleVeber’s ruleBioavability
% (Oral) > 30%% (Oral) > 70%

Nova2 (9058064-6)GoodLow0.030.01
Nova2 (95752-88-8)GoodGood0.850.48
Nova2 (88765-82-6)GoodGood0.850.27
Nova2 (Unique10)GoodGood0.760.27
Nova2 (16876-27-0)GoodGood0.850.48
Nova2 (513246-99-6)GoodGood0.760.27
Nova2 (Unique80)GoodGood0.760.23
Nova2 (279262-23-6)GoodGood0.590.27

LigandPotential carcinogen based QSARGenotoxic carcinogenicityNongenotoxic carcinogenicityPotential S. typhimurium mutagenicity

Nova2 (9058064-6)NoNegativeNegativeNo
Nova2 (95752-88-8)NoNegativeNegativeNo
Nova2 (88765-82-6)NoNegativeNegativeNo
Nova2 (Unique10)NoNegativeNegativeNo
Nova2 (16876-27-0)NoNegativeNegativeNo
Nova2 (513246-99-6)NoNegativeNegativeNo
Nova2 (Unique80)NoNegativeNegativeNo
Nova2 (279262-23-6)NoNegativeNegativeNo

Furthermore, a probability analysis of ligand’s side effects on the human health was conducted. The computational analysis was carried out by using the ACD/I-Lab software and the results for the safest ligands candidate are shown in bold, at Table 9.

LigandThe probability of side effects ()
BloodCardiovascular systemGastrointestinal systemKidneyLiverLung

Nova2 (9058064-6)
Nova2 (95752-88-8)
Nova2 (88765-82-6)
Nova2 (Unique10)
Nova2 (16876-27-0)
Nova2  (513246-99-6)
Nova2 (Unique80)
Nova2 (279262-23-6)

Table 9 shows that the entire modified ligand had adverse effects on the gastrointestinal tract, but it was not a problem because in order to distribute the drug, utilization of drug delivery technology could be in place to target the receptor. After passing the test, it was determined that the best ligand is Nova2 (513246-99-6). This is because of lower Δ value than the standard. Its value is almost the same in all of the HDAC class II enzymes. The best ligands were tested using molecular dynamics simulations to look at their stability due to the changes of solvent as well as temperature. The process was divided into three phases with the first phase of initialization temperature of 300 K, equilibration stage and heating temperature of 310 K, and the last stage of production. The simulated stage happened when the drug met the solvent and the occurrence of temperature changed and when the drug was distributed and reached its intended target. The dynamic simulation result is shown in Figure 3 as RMSD versus time (ps) curve at the molecular dynamic stage of the best drug candidate, Nova2 (513246-996), against the target enzyme.

The graph shows that Nova2 ligand (513246-99-6) was stable at the time of 5000 ps. However, progression in HDAC6 shows an increase in the curve due to the shallowness of the binding pocket; the ligand was somewhat less stable due to the influence of the solvent, which was outside the surface of the enzyme.

3.2. Discussion

Due to the versatility of organoboron compounds, they have proved themselves to be useful at the field or chemical science [51]. The synthesizing of organoboron compounds has been proven to be physically possible [52]. The flexibility of organoboron as electrophile and nucleophile compound adds a tunable property in the drug design [51]. The efficacy of organoboron compounds has been proven in the fungicide and bactericide-based in vitro experiment [53]. More specifically, the closo-carborane compounds have indeed become a logical choice for drug design due to their biological activity [43]. The hydrophobic property of closo-carborane is making chemical modifications on various compounds feasible, as a measure to observe its pharmacochemistry properties [54]. Hence, Velcade © is the only boron-based therapeutic in the market that is useful to threat multiple myolema, and there are several more that undergo clinical trials, including Talabostat © as lung cancer drug candidate [41, 55]. In this end, it is expected that there will be more boron based drugs in the market.

The computational measure for docking of simple boronic acid based compounds was already utilized and has paved the way for serious organoboron based rational drug design [56]. The starting point of the development of our pipeline was the optimization of docking method, which would be improved later on [57]. Our experience in HPV drug design was based on design of organic compounds as drug candidates [46, 58]. However, as carbon and boron based compounds have some similarity of physicochemical properties, the utilization of boron as carbon substitute for drug candidate becomes more feasible [59]. Now, closo-carborane and boronic acid based lead compounds were successfully designed based on our established methodology. The low toxicity and degradability of boronic acid into safer boric acid made it an environmental friendly compound [60]. It is also observed that closo-carborane could improve the hydrophobic interaction with enzyme [43]. Due to the robustness of organoboron compounds, the existing pipeline could be applied “as it is” without any major modification, for working with organometalloid compounds. Moreover, the designed organoboron compounds are still below the threshold of Lipinski’s rule molecular weight barrier, so it will be much simpler to develop them [59]. Thus, the Nova2 (513246-99-6) that was a combined compound from SAHA and (4-piperazin-1-ylphenyl) boronic acid was chosen as the most feasible drug candidate [61]. The slightly acidic properties of the boronic acid functional groups, combined with the electronegativity tendency from its nitrogen atom, may contribute to the inhibition activity to the Zn2+ metalloenzyme pocket [62]. However, the acidic property of Nova2 (513246-99-6) should be taken into account for its oral delivery measure, as it could inflict certain hazard for heartburn patient. In this end, prodrug construction should be considered in its synthesis strategy [63].

The wet laboratory experiment that is working with the interaction of protein and organoboron compound was already proved to be feasible [64]. Moreover, organoboron compound is just starting to be applied as radiotherapy agent [65]. The synthesis pathway for both boronic acid and closo-carborane derivatives are already applied by some research group [66, 67]. Thus, in order to soften the complexity of the synthesis, a prediction method will be utilized to evaluate the synthesis accessibility [6870]. To this end, by applying the information from in silico results, it is expected that the laboratory synthesis and bioassay experimentation for organoboron compound should be straightforward and not difficult.

4. Conclusion

Modification of ZBG and CAP at SAHA with a boron compound and the p-aminobenzoic (PABA) as linker turns out to yield results as expected, as they show better HDAC inhibition than SAHA. This modification resulted in a total of 1,100 ligands. After going through molecular docking simulations, the top 8 ligands were obtained; they have much lower Δ than standard SAHA and TSA. The ligands were Nova2 (9058064-6), Nova2 (95752-88-8), Nova2 (88765-82-6), Nova2 (unique10), Nova2 (16876-27-0), Nova2 (513246-99-6), Nova2 (unique80), and Nova2 (279262-23-6). After undergoing the screening process of QSAR, pharmacological properties, and ADME-Tox, 4 best ligands were collected that can be used as drug candidates. They were Nova2 (95752-88-8), Nova2 (16876-27-0), Nova2 (513246-99-6), and Nova2 (279262-23-6). After a thorough analysis, Nova2 (513246-99-6) was concluded as the best drug candidate, and it was piped to the molecular dynamics simulation process.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.


The authors would like to thank The Directorate of Higher Education, Ministry of Culture and Education, for giving us Decentralization Grant no. 1087/H2.R12/HKP.05.00/2014. Usman Sumo Friend Tambunan and Arli Aditya Parikesit supervised the research, Ridla Bakri added more relevant bioinorganic chemistry materials and proof read the paper with Niken Widiyanti, Cipta Prio Satriyanto worked on the technical and experimental details, Arli Aditya Parikesit was responsible for writing the paper, and Djati Kerami was giving critical suggestion for improving our pipeline. The funding source had no involvement in study design, in the collection, analysis and interpretation of data, in the writing of the report and in the decision to submit the paper for publication.


  1. WHO, Cervical Cancer Summary Report, 2010.
  2. U. S. F. Tambunan and A. A. Parikesit, “HPV bioinformatics: in silico detection, drug design and prevention agent development,” in Topics on Cervical Cancer With an Advocacy for Prevention, R. Rajkumar, Ed., pp. 237–252, Intech Publishing, Rijeka, Croatia, 2012, View at: Google Scholar
  3. Cancer Research UK, “CancerStats—Cancer Worldwide,” 2013. View at: Google Scholar
  4. H. zur Hausen, “Papillomaviruses in the causation of human cancers—a brief historical account,” Virology, vol. 384, no. 2, pp. 260–265, 2009. View at: Publisher Site | Google Scholar
  5. H. zur Hausen, “Human papillomavirus & cervical cancer,” Indian Journal of Medical Research, vol. 130, no. 3, p. 209, 2009. View at: Google Scholar
  6. H. Zur Hausen, “Papillomaviruses and cancer: from basic studies to clinical application,” Nature Reviews Cancer, vol. 2, no. 5, pp. 342–350, 2002. View at: Publisher Site | Google Scholar
  7. P. T. Soliman, B. M. Slomovitz, and J. K. Wolf, “Mechanisms of cervical cancer,” Drug Discovery Today: Disease Mechanisms, vol. 1, pp. 253–258, 2004. View at: Publisher Site | Google Scholar
  8. J. Paavonen, “Human papillomavirus infection and the development of cervical cancer and related genital neoplasias,” International Journal of Infectious Diseases, vol. 11, supplement 2, pp. S3–S9, 2007. View at: Publisher Site | Google Scholar
  9. H. Trottier and A. N. Burchell, “Epidemiology of mucosal human papillomavirus infection and associated diseases,” Public Health Genomics, vol. 12, no. 5-6, pp. 291–307, 2009. View at: Publisher Site | Google Scholar
  10. S. R. Pagliusi, J. Dillner, M. Pawlita, W. G. V. Quint, C. M. Wheeler, and M. Ferguson, “Chapter 23: international Standard reagents for harmonization of HPV serology and DNA assays-an update,” Vaccine, vol. 24, no. 3, pp. S193–S200, 2006. View at: Publisher Site | Google Scholar
  11. M. Steben and E. Duarte-Franco, “Human papillomavirus infection: epidemiology and pathophysiology,” Gynecologic Oncology, vol. 107, no. 2, pp. S2–S5, 2007. View at: Publisher Site | Google Scholar
  12. E. Boccardo, A. P. Lepique, and L. L. Villa, “The role of inflammation in HPV carcinogenesis,” Carcinogenesis, vol. 31, no. 11, pp. 1905–1912, 2010. View at: Publisher Site | Google Scholar
  13. H. Zur Hausen, “Papillomaviruses causing cancer: evasion from host-cell control in early events in carcinogenesis,” Journal of the National Cancer Institute, vol. 92, no. 9, pp. 690–698, 2000. View at: Publisher Site | Google Scholar
  14. L. Gewin, H. Myers, T. Kiyono, and D. A. Galloway, “Identification of a novel telomerase repressor that interacts with the human papillomavirus type-16 E6/E6-AP complex,” Genes and Development, vol. 18, no. 18, pp. 2269–2282, 2004. View at: Publisher Site | Google Scholar
  15. M. W. Łuczak and P. P. Jagodzinski, “Apicidin down-regulates human papillomavirus type 16 E6 and E7 transcripts and proteins in SiHa cervical cancer cells,” Cancer Letters, vol. 272, no. 1, pp. 53–60, 2008. View at: Publisher Site | Google Scholar
  16. K. Münger, A. Baldwin, K. M. Edwards et al., “Mechanisms of human papillomavirus-induced oncogenesis,” Journal of Virology, vol. 78, no. 21, pp. 11451–11460, 2004. View at: Publisher Site | Google Scholar
  17. G. Herbein and D. Wendling, “Histone deacetylases in viral infections,” Clinical Epigenetics, vol. 1, no. 1-2, pp. 13–24, 2010. View at: Publisher Site | Google Scholar
  18. A. Szalmás and J. Kónya, “Epigenetic alterations in cervical carcinogenesis,” Seminars in Cancer Biology, vol. 19, no. 3, pp. 144–152, 2009. View at: Publisher Site | Google Scholar
  19. P. Finzer, A. Aguilar-Lemarroy, and F. Rösl, “The role of human papillomavirus oncoproteins E6 and E7 in apoptosis,” Cancer Letters, vol. 188, no. 1-2, pp. 15–24, 2002. View at: Publisher Site | Google Scholar
  20. J. S. Park, E. J. Kim, H. J. Kwon, E. S. Hwang, S. E. Namkoong, and S. J. Um, “Inactivation of interferon regulatory factor-1 tumor suppressor protein by HPV E7 oncoprotein. Implication for the E7-mediated immune evasion mechanism in cervical carcinogenesis,” The Journal of Biological Chemistry, vol. 275, no. 10, pp. 6764–6769, 2000. View at: Publisher Site | Google Scholar
  21. G. Bouchain and D. Delorme, “Novel hydroxamate and anilide derivatives as potent histone deacetylase inhibitors: synthesis and antiproliferative evaluation,” Current Medicinal Chemistry, vol. 10, no. 22, pp. 2359–2372, 2003. View at: Publisher Site | Google Scholar
  22. G. P. Delcuve, D. H. Khan, and J. R. Davie, “Roles of histone deacetylases in epigenetic regulation: emerging paradigms from studies with inhibitors,” Clinical Epigenetics, vol. 4, no. 1, article 5, 2012. View at: Publisher Site | Google Scholar
  23. T. Kouzarides, “Chromatin modifications and their function,” Cell, vol. 128, no. 4, pp. 693–705, 2007. View at: Publisher Site | Google Scholar
  24. M. H. Kuo and C. D. Allis, “Roles of histone acetyltransferases and deacetylases in gene regulation,” Bioessays, vol. 20, pp. 615–626, 1998. View at: Publisher Site | Google Scholar
  25. C. L. Woodcock and R. P. Ghosh, “Chromatin higher-order structure and dynamics,” Cold Spring Harbor Perspectives in Biology, vol. 2, Article ID a000596, 2010. View at: Publisher Site | Google Scholar
  26. G. Legube and D. Trouche, “Regulating histone acetyltransferases and deacetylases,” EMBO Reports, vol. 4, no. 10, pp. 944–947, 2003. View at: Publisher Site | Google Scholar
  27. B. R. Keppler and T. K. Archer, “Chromatin-modifying enzymes as therapeutic targets: part 2,” Expert Opinion on Therapeutic Targets, vol. 12, no. 11, pp. 1457–1467, 2008. View at: Publisher Site | Google Scholar
  28. C. Dell'Aversana, I. Lepore, and L. Altucci, “HDAC modulation and cell death in the clinic,” Experimental Cell Research, vol. 318, pp. 1229–1244, 2012. View at: Publisher Site | Google Scholar
  29. K. Darvas, S. Rosenberger, D. Brenner et al., “Histone deacetylase inhibitor-induced sensitization to TNFα/TRAIL- mediated apoptosis in cervical carcinoma cells is dependent on HPV oncogene expression,” International Journal of Cancer, vol. 127, no. 6, pp. 1384–1392, 2010. View at: Publisher Site | Google Scholar
  30. M. S. Finnin, J. R. Donigian, A. Cohen et al., “Structures of a histone deacetylase homologue bound to the TSA and SAHA inhibitors,” Nature, vol. 401, no. 6749, pp. 188–193, 1999. View at: Publisher Site | Google Scholar
  31. K. Chen, L. Xu, and O. Wiest, “Computational exploration of zinc binding groups for HDAC inhibition,” Journal of Organic Chemistry, vol. 78, pp. 5051–5055, 2013. View at: Publisher Site | Google Scholar
  32. S.-K. Anandan, J. S. Ward, R. D. Brokx, M. R. Bray, D. V. Patel, and X. Xiao, “Mercaptoamide-based non-hydroxamic acid type histone deacetylase inhibitors,” Bioorganic and Medicinal Chemistry Letters, vol. 15, no. 8, pp. 1969–1972, 2005. View at: Publisher Site | Google Scholar
  33. S. Thangapandian, S. John, Y. Lee, V. Arulalapperumal, and K. W. Lee, “Molecular modeling study on tunnel behavior in different histone deacetylase isoforms,” PLoS ONE, vol. 7, Article ID e49327, 2012. View at: Publisher Site | Google Scholar
  34. P. Bertrand, “Inside HDAC with HDAC inhibitors,” European Journal of Medicinal Chemistry, vol. 45, no. 6, pp. 2095–2116, 2010. View at: Publisher Site | Google Scholar
  35. L. Marek, A. Hamacher, F. K. Hansen et al., “Histone deacetylase (HDAC) inhibitors with a novel connecting unit linker region reveal a selectivity profile for HDAC4 and HDAC5 with improved activity against chemoresistant cancer cells,” Journal of Medicinal Chemistry, vol. 56, no. 2, pp. 427–436, 2013. View at: Publisher Site | Google Scholar
  36. H. Rajak, P. Kumar, P. Parmar et al., “Appraisal of GABA and PABA as linker: design and synthesis of novel benzamide based histone deacetylase inhibitors,” European Journal of Medicinal Chemistry, vol. 53, pp. 390–397, 2012. View at: Publisher Site | Google Scholar
  37. P. C. Chen, V. Patil, W. Guerrant, P. Green, and A. K. Oyelere, “Synthesis and structure-activity relationship of histone deacetylase (HDAC) inhibitors with triazole-linked cap group,” Bioorganic and Medicinal Chemistry, vol. 16, no. 9, pp. 4839–4853, 2008. View at: Publisher Site | Google Scholar
  38. D. Huang, X. Li, Y. Wei, and Z. Xiu, “A novel series of l-2-benzyloxycarbonylamino-8-(2-pyridyl)- disulfidyloctanoic acid derivatives as histone deacetylase inhibitors: design, synthesis and molecular modeling study,” European Journal of Medicinal Chemistry, vol. 52, pp. 111–122, 2012. View at: Publisher Site | Google Scholar
  39. M. Calvaresi and F. Zerbetto, “In silico carborane docking to proteins and potential drug targets,” Journal of Chemical Information and Modeling, vol. 51, no. 8, pp. 1882–1896, 2011. View at: Publisher Site | Google Scholar
  40. M. Korkmaz, E. Uzgören, S. Bakirdere, F. Aydin, and O. Y. Ataman, “Effects of dietary boron on cervical cytopathology and on micronucleus frequency in exfoliated buccal cells,” Environmental Toxicology, vol. 22, no. 1, pp. 17–25, 2007. View at: Publisher Site | Google Scholar
  41. S. J. Baker, C. Z. Ding, T. Akama, Y.-K. Zhang, V. Hernandez, and Y. Xia, “Therapeutic potential of boron-containing compounds,” Future Medicinal Chemistry, vol. 1, pp. 1275–1288, 2009. View at: Publisher Site | Google Scholar
  42. J. F. Valliant, K. J. Guenther, A. S. King et al., “The medicinal chemistry of carboranes,” Coordination Chemistry Reviews, vol. 232, pp. 173–230, 2002. View at: Publisher Site | Google Scholar
  43. F. Issa, M. Kassiou, and L. M. Rendina, “Boron in drug discovery: carboranes as unique pharmacophores in biologically active compounds,” Chemical Reviews, vol. 111, no. 9, pp. 5701–5722, 2011. View at: Publisher Site | Google Scholar
  44. U. S. F. Tambunan, I. K. Limans, and A. A. Parikesit, “Molecular dynamic simulation of complex NS2B-NS3 DENV2 protease with potential inhibitors of cyclic tripeptide,” African Journal of Biotechnology, vol. 12, no. 28, pp. 4419–4431, 2013. View at: Publisher Site | Google Scholar
  45. A. A. Parikesit, K. Kinanty, and U. S. F. Tambunan, “Screening of commercial cyclic peptides as inhibitor envelope protein dengue virus (DENV) through molecular docking and molecular dynamics,” Pakistan Journal of Biological Sciences, vol. 16, pp. 1836–1848, 2013. View at: Publisher Site | Google Scholar
  46. U. S. F. Tambunan, A. A. Parikesit, T. Prasetia, and D. Kerami, “In silico molecular interaction studies of suberoylanilide hydroxamic acid and its modified compounds with histones deacetylase class II Homo sapiens as curative measure towards cervical cancer,” Engineering, vol. 5, pp. 203–206, 2013. View at: Publisher Site | Google Scholar
  47. U. S. F. Tambunan, A. Randy, and A. A. Parikesit, “Design of Candida antarctica lipase B thermostability improvement by introducing extra disulfide bond into the enzyme,” OnLine Journal of Biological Sciences, vol. 14, no. 2, pp. 108–118, 2014. View at: Publisher Site | Google Scholar
  48. U. S. F. Tambunan, N. Apriyanti, A. A. Parikesit, W. Chua, and K. Wuryani, “Computational design of disulfide cyclic peptide as potential inhibitor of complex NS2B-NS3 dengue virus protease,” African Journal of Biotechnology, vol. 10, no. 57, pp. 12281–12290, 2011. View at: Google Scholar
  49. C. A. Lipinski, F. Lombardo, B. W. Dominy, and P. J. Feeney, “Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings,” Advanced Drug Delivery Reviews, vol. 23, no. 1–3, pp. 3–25, 1997. View at: Publisher Site | Google Scholar
  50. D. Lagorce, O. Sperandio, H. Galons, M. A. Miteva, and B. O. Villoutreix, “FAF-Drugs2: Free ADME/tox filtering tool to assist drug discovery and chemical biology projects,” BMC Bioinformatics, vol. 9, article 396, 2008. View at: Publisher Site | Google Scholar
  51. N. A. Petasis, “Expanding roles for organoboron compounds—versatile and valuable molecules for synthetic, biological and medicinal chemistry,” Australian Journal of Chemistry, vol. 60, no. 11, pp. 795–798, 2007. View at: Publisher Site | Google Scholar
  52. M. Penner, V. Rauniyar, L. T. Kaspar, and D. G. Hall, “Catalytic asymmetric synthesis of palmerolide a via organoboron methodology,” Journal of the American Chemical Society, vol. 131, no. 40, pp. 14216–14217, 2009. View at: Publisher Site | Google Scholar
  53. V. P. Singh, R. V. Singh, and J. P. Tandon, “Stereochemical and biochemical aspects of some organoboron complexes of sulphur donor ligands,” Journal of Inorganic Biochemistry, vol. 39, pp. 237–245, 1990. View at: Publisher Site | Google Scholar
  54. R. Tiwari, K. Mahasenan, R. Pavlovicz, C. Li, and W. Tjarks, “Carborane clusters in computational drug design: a comparative docking evaluation using AutoDock, FlexX, Glide, and Surflex,” Journal of Chemical Information and Modeling, vol. 49, no. 6, pp. 1581–1589, 2009. View at: Publisher Site | Google Scholar
  55. C. C. Cunningham, “Talabostat,” Expert Opinion on Investigational Drugs, vol. 16, no. 9, pp. 1459–1465, 2007. View at: Publisher Site | Google Scholar
  56. A. Tafi, M. Agamennone, P. Tortorella, S. Alcaro, C. Gallina, and M. Botta, “AMBER force field implementation of the boronate function to simulate the inhibition of β-lactamases by alkyl and aryl boronic acids,” European Journal of Medicinal Chemistry, vol. 40, no. 11, pp. 1134–1142, 2005. View at: Publisher Site | Google Scholar
  57. U. S. F. Tambunan, F. Fadilah, and A. A. Parikesit, “Bioactive compounds screening from Zingiberaceae family as influenza a/swine flu virus neuraminidase inhibitor through docking approach,” OnLine Journal of Biological Sciences, vol. 10, no. 4, pp. 151–156, 2010. View at: Google Scholar
  58. U. S. F. Tambunan, R. Bakri, T. Prasetia, A. A. Parikesit, and D. Kerami, “Molecular dynamics simulation of complex Histones Deacetylase (HDAC) Class II Homo Sapiens with suberoylanilide hydroxamic acid (SAHA) and its derivatives as inhibitors of cervical cancer,” Bioinformation, vol. 9, pp. 696–700, 2013, View at: Google Scholar
  59. P. Hunter, “Not boring at all. Boron is the new carbon in the quest for novel drug candidates,” EMBO Reports, vol. 10, no. 2, pp. 125–128, 2009. View at: Publisher Site | Google Scholar
  60. D. Hall, “Structure, properties, and preparation of boronic acid derivatives. Overview of their reactions and applications,” in Boronic Acids: Preparation and Applications in Organic Synthesis, Medicine, chapter 1, pp. 1–99, 2006. View at: Publisher Site | Google Scholar
  61., “4-piperazin-1-ylphenylboronic acid,” synthesis and properties,, 2003, View at: Google Scholar
  62. C. Chen, P. Dagneau, E. J. J. Grabowski et al., “Practical asymmetric synthesis of a potent Cathepsin K inhibitor. Efficient palladium removal following Suzuki coupling,” Journal of Organic Chemistry, vol. 68, no. 7, pp. 2633–2638, 2003. View at: Publisher Site | Google Scholar
  63. A. T. de Albuquerque Silva, A. Silva, M. C. Chung et al., “Advances in prodrug design,” Mini-Reviews in Medicinal Chemistry, vol. 5, no. 10, pp. 893–914, 2005. View at: Publisher Site | Google Scholar
  64. Y. Endo, T. Yoshimi, K. Kimura, and A. Itai, “Protein kinase C modulators bearing dicarba-closo-dodecaborane as a hydrophobic pharmacophore,” Bioorganic & Medicinal Chemistry Letters, vol. 9, no. 17, pp. 2561–2564, 1999. View at: Publisher Site | Google Scholar
  65. R. Tiwari, A. Toppino, H. K. Agarwal et al., “Synthesis, biological evaluation, and radioiodination of halogenated closo -carboranylthymidine analogues,” Inorganic Chemistry, vol. 51, no. 1, pp. 629–639, 2012. View at: Publisher Site | Google Scholar
  66. D. S. Wilbur, D. K. Hamlin, R. R. Srivastava, and M. Chyan, “Synthesis, radioiodination, and biodistribution of some nido- and closo-monocarbon carborane derivatives,” Nuclear Medicine and Biology, vol. 31, no. 4, pp. 523–530, 2004. View at: Publisher Site | Google Scholar
  67. G. A. Molander, S. L. J. Trice, and S. D. Dreher, “Palladium-catalyzed, direct boronic acid synthesis from aryl chlorides: a simplified route to diverse boronate ester derivatives,” Journal of the American Chemical Society, vol. 132, no. 50, pp. 17701–17703, 2010. View at: Publisher Site | Google Scholar
  68. P. Ertl and A. Schuffenhauer, “Estimation of synthetic accessibility score of drug-like molecules based on molecular complexity and fragment contributions,” Journal of Cheminformatics, vol. 1, no. 1, article no. 8, 2009. View at: Publisher Site | Google Scholar
  69. K. Boda, T. Seidel, and J. Gasteiger, “Structure and reaction based evaluation of synthetic accessibility,” Journal of Computer-Aided Molecular Design, vol. 21, no. 6, pp. 311–325, 2007. View at: Publisher Site | Google Scholar
  70. J. C. Baber and M. Feher, “Predicting synthetic accessibility: application in drug discovery and development,” Mini-Reviews in Medicinal Chemistry, vol. 4, no. 6, pp. 681–692, 2004. View at: Google Scholar

Copyright © 2014 Ridla Bakri et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2020, as selected by our Chief Editors. Read the winning articles.