{"id":39171,"date":"2021-06-30T11:28:35","date_gmt":"2021-06-30T11:28:35","guid":{"rendered":"https:\/\/biomedpharmajournal.org\/?p=39171"},"modified":"2021-07-13T09:41:17","modified_gmt":"2021-07-13T09:41:17","slug":"molecular-docking-studies-of-ros-agent-from-quinone-family-to-reductase-enzymesimplication-in-finding-anticancer-drug-candidate","status":"publish","type":"post","link":"https:\/\/biomedpharmajournal.org\/staging\/vol14no2\/molecular-docking-studies-of-ros-agent-from-quinone-family-to-reductase-enzymesimplication-in-finding-anticancer-drug-candidate\/","title":{"rendered":"Molecular Docking Studies of ROS Agent from Quinone Family to Reductase Enzymes:Implication in Finding Anticancer Drug Candidate"},"content":{"rendered":"<p style=\"text-align: justify;\"><strong>Introduction<\/strong><\/p>\n<p style=\"text-align: justify;\">Drugs belonging to the quinone family have been successfully explored for their anti-tumor activity. Understanding the cellular metabolism of their cytotoxic compounds is importance in the field of oncology for finding new therapeutic strategy. Generally, cytotoxicity of quinones is related to its rapid redox cycling which associated with Reactive Oxygen Species (ROS) generation to provoke oxidative stress. Quinones undergo either one-electron reduction to form radical semiquinones or two-electron reduction to form less toxic compound hydroquinones.<sup>1,2<\/sup>Those formations may react with oxygen to produce superoxide radicals O<sub>2<\/sub><sup>.-<\/sup>. Those formations have a short lifetime in the cells and reactive.<sup>3<\/sup>The ROS radical consequences on the occurrence a series of chain reactions of neutral species in the cells or transformation into other types of ROS such as hydrogen peroxide formation. ROS is highly damaging for cells such as vandalized an important macromolecules which consequences on the disruption of their normal functions in the cell. Excessive ROS production will be detrimental to cells by damaging lipid membranes, proteins, and DNA.<sup>3\u20135<\/sup>Normal and cancer cells have different ROS limit level. The researchers reported that cancer cells have an endogenous ROS concentrations greater than normal cells.<sup>6,7<\/sup> Increasing the amount of ROS in cancer cells does not necessarily affect to cell death, due to differences limit of oxidative stress in cancercells. Therefore, it becomes important to study the strong stress agents to kill cancer cells.<\/p>\n<p style=\"text-align: justify;\">The use of compounds that produce ROS as their main mechanism for cell death, are often selected as a drug clinically, although their mechanism to produce ROS has not clearly understood yet. Some quinones are used as chemotherapy, such as mitoxantrone and menadione. Their reduction ability is believed to be responsible in ROS production. Mitoxantrone is anticancer agents that used clinically for the breast cancer, leukemia, lymphoma, and a variety of other malignant tumors treatment.<sup>8,9<\/sup> Menadione has two capacities considered in doing lethal action as an anticancer agent; involved in the redox cycle related to ROS production<sup>10,11<\/sup>and menadione has arylation capacity to form conjugation product with thiols proteins in the cell caused a decline in the levels of antioxidants in the cells.<sup>11\u201314<\/sup>Two capacities of menadione has been confirmed by in silico study that showed menadione can interact with NADH domain (correlated with menadione induced ROS production) or FAD domain of AIF (correlated with arylation capacity of menadione).<sup>15<\/sup><\/p>\n<p style=\"text-align: justify;\">The difference of quinone types and its environment affect the oxygen radicals rate, which consequences to cell damage.<sup>16<\/sup> Cytotoxic level of quinones is also affected by enzymatic reducing systems.<sup>17<\/sup><sup>,<\/sup><sup>18<\/sup>For instance, P450 reductase is an enzyme that play a major role in one-electron reduction of anticancer compounds. Other enzyme involved in one-electron reduction is ubiquinone oxidoreductase contained in the mitochondria (mitochondrial NADH-dependent ubiquinone oxidoreductase)<sup>18<\/sup>or the AIF.<sup>19<\/sup> While, two-electron reduction of quinones generally are catalyzed by NAD(P)H:quinone oxidoreductase 1 (NQO1, DT-diaphorase). This enzyme catalyzes the formation of a neutral compound_hydroquinone, which is less reactive than radical semiquinone.<sup>20<\/sup> The molecular mechanism of quinones to those enzymes has not yet understood. The aim of this study is to understand molecular interaction of quinones to reductase enzymes\u00a0 such as cytochrome P450 reductase, ubiquinone reductase (NQO1), and apoptosis inducing factor (AIF) which is also reported as NADH:quinones reductase.<sup>19<\/sup><\/p>\n<p style=\"text-align: justify;\"><strong>Material and Methods<\/strong><\/p>\n<p style=\"text-align: justify;\">The crystal structure of\u00a0 NQO1, P450 reductase, and Apoptosis Inducing Factor (AIF) with Protein Data Bank (PDB) code subsequently 2F1O, 1AMO, and 3GD4 in complex with its ligands were retrieved from the PDB (http:\/\/www.rscb.org\/pdb)<sup>21,22<\/sup>Three original ligands of each enzyme used in this study; Bishydroxy [2H-1-benzopyran-2-one,1,2-benzopyrone] (C<sub>19<\/sub>H<sub>12<\/sub>O<sub>6<\/sub>, DTC), Nicotinamide Adenine Dinucleotide Phosphate (C<sub>21<\/sub>H<sub>28<\/sub>N<sub>7<\/sub>O<sub>17<\/sub>P<sub>3<\/sub>, NAP), and Nicotinamide Adenine Dinucleotide (NADH) (NADH).<\/p>\n<p style=\"text-align: justify;\"><strong>General Procedure<\/strong><\/p>\n<p style=\"text-align: justify;\"><strong>Protein and ligand preparation<\/strong><\/p>\n<p style=\"text-align: justify;\">The complexes, heteroatoms and water molecules bound to the receptor molecules were removed from each protein structure. Finally hydrogen atoms were merged to the target receptor molecules. Protein preparation have done using Autodock Vina software.<sup>23<\/sup>The ligands were taken from each crystal structure of the enzymes. DTC is ligand of NQO1, NAP is a ligand of P450, and NADH is a ligand of AIF. Quinone compounds were used as ligands are 1,4-benzoquinone, 2-methyl-1,4-benzoquinone, 1,4-naphthoquinone, 2-methyl-1,4-naphthoquinone (menadione), antraquinone, 1,4-dihydroxy-5,8-bis[2-(2-hydroxyethylamino)ethylamino]-anthracene-9,10-dione (mitoxantrone). Ligand structures were sketched using marvin sketch software. Dreiding energy lowest was built to optimize all ligands using marvin sketch software. The optimized ligands were prepared for docking studies.<\/p>\n<p style=\"text-align: justify;\"><strong>Docking<\/strong><strong> analysis<\/strong><\/p>\n<p style=\"text-align: justify;\">The docking analysis of AIF with ligands was carried out by auto dock vina docking software which is most commonly available software. Grid resolution was set to 1 \u00c5, located to NADH binding domain.Docking analysis of AIF-menadione and dreiding energy lowest of menadione were taken from previous report that showed -6.3 kcal\/mol and 34.05 kcal\/mol respectively.<sup>15<\/sup><\/p>\n<p style=\"text-align: justify;\"><strong>Superposition analysis<\/strong><\/p>\n<p style=\"text-align: justify;\">Superposition of 3D structure of receptors were studied in order to know the folding similarity of both protein, especially on transferase binding domain. This analysis was performed by FAT CAT rigid\/flexible-structure alignment.<sup>24<\/sup><\/p>\n<p style=\"text-align: justify;\"><strong>Results and discussion<\/strong><\/p>\n<p style=\"text-align: justify;\"><strong>Validation of docking system<\/strong><\/p>\n<p style=\"text-align: justify;\">This report begins by showing the validation result of docking program settingin order to confirm the feasibility of system used. Ligand-receptor interaction on crystal structure was used as reference, such as a number of contact residues, the type of contact residues, the number of hydrogen bonds, and Root Mean Standard Deviation (RMSD). This information was retrieved from PDB. All parameters were compared to docking experiment of each ligand-receptor interaction. When the experiment has high similarity with the reference and the value of RMSD \u2264 2, the experimental system will be approved. The result showed that the setup of docking program used in this experiment is good enough in predicting the receptor-ligand interactions. (Table 1, Figure 1)<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig1.jpg\"><img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-39174\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig1-150x150.jpg\" alt=\"Vol14No2_Mol_Hes_fig1\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig1-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig1-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig1.jpg 693w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/a><\/td>\n<td><strong>Figure 1: Ligand position relative to its receptor. P450 reductase and NAP (experiment or refference) (A); NQO1 and DTC (experiment or\u00a0 refference) (B);\u00a0 AIF and NADH (experiment\u00a0\u00a0or refference) (C). <\/strong><strong>RMSD was employed base on the position of ligand on crystal structure\u00a0(as reference) compare to docking result using autodock vina software<\/strong><strong>.<\/strong><\/p>\n<p><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig1.jpg\" target=\"_blank\">Click here to view figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\"><strong>Table 1: Validation parameters of docking program<\/strong><\/p>\n<table style=\"width: 95%;\" border=\"1\" cellspacing=\"0\" cellpadding=\"4\">\n<tbody>\n<tr>\n<td style=\"text-align: center;\" width=\"234\"><strong>Receptor-Ligand<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"98\"><strong>P450-NAP<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"98\"><strong>NQO1-DTC<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"98\"><strong>AIF-NADH<\/strong><\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"234\"><strong>% similarity of contact residues <\/strong><\/td>\n<td style=\"text-align: center;\" width=\"98\">86%<\/td>\n<td style=\"text-align: center;\" width=\"98\">100%<\/td>\n<td style=\"text-align: center;\" width=\"98\">73%<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"234\"><strong>RMSD<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"98\">0.97<\/td>\n<td style=\"text-align: center;\" width=\"98\">1.68<\/td>\n<td style=\"text-align: center;\" width=\"98\">1.51<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\">Caption: % similarity of contact residues was measured by identify the type of contact residues of receptor on experiment compare to contact residues of receptor on crystal structure (as reference) RMSD was employed base on the position of ligand on crystal structure (as reference) compare to docking result using autodockvina software.<\/p>\n<p style=\"text-align: justify;\">After validating the docking system, three reductase enzymes as receptor and quinones as ligands were prepared in order to get the lowest energy for interaction. Three dimensional structure of P450 reductase (1AMO), NQO1 (2J1O) and AIF (3GD4) were prepared by removing their original ligand and water molecules from the complexes. (Figure 2) Quinone compounds were prepared as ligand by optimizing the structure using dreiding energy force field. (Figure 2)<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig2.jpg\"><img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-39175\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig2-150x150.jpg\" alt=\"Vol14No2_Mol_Hes_fig2\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig2-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig2-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig2.jpg 637w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/a><\/td>\n<td><strong>Figure 2: Three dimension (3D) structure of receptors and optimized ligands. The crystal structure of\u00a0 receptors (NQO1, P450 reductase, and AIF) were retrieved from the PDB.<\/strong><\/p>\n<p><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig2.jpg\" target=\"_blank\">Click here to view figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\"><strong>The number of ring affect to quinone\u2019s affinity to enzymes. <\/strong><\/p>\n<p style=\"text-align: justify;\">Docking results between reductase enzymes and quinones are tabulated in Table 2. Docking analysis showed that anthraquinone have higher affinity of all enzymes compared toother quinones. Affinity energy of anthraquinone&lt;naphtaquinone&lt;benzoquinone. The number of benzene rings seems affect to the affinity, affinity will increase with increasing number of rings. These results seem to be in line with previous studies which revealed that AIF was able to perform the most efficient enzymatic activity against benzoquninone, followed by naphthoquinone and anthraquinone[19]. The lower the affinity of benzoquinone to AIF, with the same rate constant of the enzyme, it will make AIF perform more efficiently.<\/p>\n<p style=\"text-align: justify;\"><strong>Table 2: Affinity energy of quinones to P450 reductase, NQO1, and AIF<\/strong><\/p>\n<table style=\"width: 95%;\" border=\"1\" cellspacing=\"0\" cellpadding=\"4\">\n<tbody>\n<tr>\n<td style=\"text-align: center;\" width=\"208\"><strong>Receptor\/Ligand<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"158\"><strong>P450 reductase<\/strong><\/p>\n<p><strong>(kcal\/mol)<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"102\"><strong>NQO1<\/strong><\/p>\n<p><strong>(kcal\/mol)<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"102\"><strong>AIF<\/strong><\/p>\n<p><strong>(kcal\/mol)<\/strong><\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"208\">1,4-benzoquinone<\/td>\n<td style=\"text-align: center;\" width=\"158\">-5.2<\/td>\n<td style=\"text-align: center;\" width=\"102\">-4.3<\/td>\n<td style=\"text-align: center;\" width=\"102\">-4.8<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"208\">2-methyl-1,4-benzoquinone<\/td>\n<td style=\"text-align: center;\" width=\"158\">-5.1<\/td>\n<td style=\"text-align: center;\" width=\"102\">-4.3<\/td>\n<td style=\"text-align: center;\" width=\"102\">-4.8<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"208\">1,4-naphthoquinone<\/td>\n<td style=\"text-align: center;\" width=\"158\">-7.1<\/td>\n<td style=\"text-align: center;\" width=\"102\">-5.9<\/td>\n<td style=\"text-align: center;\" width=\"102\">-6.5<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"208\">2-methyl-1,4-naphthoquinone<\/td>\n<td style=\"text-align: center;\" width=\"158\">-6.7<\/td>\n<td style=\"text-align: center;\" width=\"102\">-5.5<\/td>\n<td style=\"text-align: center;\" width=\"102\">-6.9<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"208\">Antraquinone<\/td>\n<td style=\"text-align: center;\" width=\"158\">-8.4<\/td>\n<td style=\"text-align: center;\" width=\"102\">-6.4<\/td>\n<td style=\"text-align: center;\" width=\"102\">-8.1<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"208\">Mitoxantrone<\/td>\n<td style=\"text-align: center;\" width=\"158\">-7.8<\/td>\n<td style=\"text-align: center;\" width=\"102\">-5.2<\/td>\n<td style=\"text-align: center;\" width=\"102\">-7.8<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"208\">NADH<\/td>\n<td style=\"text-align: center;\" width=\"158\">n.d<\/td>\n<td style=\"text-align: center;\" width=\"102\">-6.6<\/td>\n<td style=\"text-align: center;\" width=\"102\">-9.3<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"208\">NAP<\/td>\n<td style=\"text-align: center;\" width=\"158\">-8.8<\/td>\n<td style=\"text-align: center;\" width=\"102\">n.d<\/td>\n<td style=\"text-align: center;\" width=\"102\">n.d<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"208\">DTC<\/td>\n<td style=\"text-align: center;\" width=\"158\">n.d<\/td>\n<td style=\"text-align: center;\" width=\"102\">-7.2<\/td>\n<td style=\"text-align: center;\" width=\"102\">n.d<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\">Caption: n.d: not determined<\/p>\n<p style=\"text-align: justify;\">P450 reductase-ligand interaction analysis showed that naphthoquinones and anthraquinones were held by 2 or 3 contact residues more than benzoquinones. (Figure 3) Even, other anthraquinone, mitoxantrone, interacts with 12 contact residues. The number of contact residues may effect to the bigger affinity of mitoxantrone to P450 reductase. While antraquinone, the number of benzene ring appears have a greater contribution (compared to Van der Waals interaction) to its affinity energy value. (Table 3)<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig3.jpg\"><img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-39176\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig3-150x150.jpg\" alt=\"Vol14No2_Mol_Hes_fig3\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig3-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig3-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig3.jpg 740w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/a><\/td>\n<td><strong>Figure 3: Binding pocket of P450 reductase to each quinones<\/strong><\/p>\n<p><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig3.jpg\" target=\"_blank\">Click here to view figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\"><strong>Table 3: Ligands properties and its contact residues on P450 reductase<\/strong><\/p>\n<table style=\"width: 95%;\" border=\"1\" cellspacing=\"0\" cellpadding=\"4\">\n<tbody>\n<tr>\n<td style=\"text-align: center;\" width=\"143\"><strong>ligands<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"216\"><strong>Contact residues on P450 reductase<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"132\"><strong>Water Solubility, pH 7.4<\/strong><\/td>\n<td style=\"text-align: center;\" width=\"113\"><strong>Polar Surface Area (PSA)<\/strong><\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"143\">1,4-benzoquinone<\/td>\n<td style=\"text-align: center;\" width=\"216\">ser596,\u00a0 val605, tyr604<\/td>\n<td style=\"text-align: center;\" width=\"132\">High (higher than 0.06 mg\/ml)<\/td>\n<td style=\"text-align: center;\" width=\"113\">34.14<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"143\">2-methyl-1,4-benzoquinone<\/td>\n<td style=\"text-align: center;\" width=\"216\">val605, tyr604<\/td>\n<td style=\"text-align: center;\" width=\"132\">High (higher than 0.06 mg\/ml)<\/td>\n<td style=\"text-align: center;\" width=\"113\">34.14<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"143\">1,4-naphthaquinone<\/td>\n<td style=\"text-align: center;\" width=\"216\">gly565, ser596, val605, tyr604, met636<\/td>\n<td style=\"text-align: center;\" width=\"132\">High (higher than 0.06 mg\/ml)<\/td>\n<td style=\"text-align: center;\" width=\"113\">34.14<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"143\">2-methyl-1,4-naphthaquinone<\/td>\n<td style=\"text-align: center;\" width=\"216\">val605, tyr604, ser596, gly565, met636<\/td>\n<td style=\"text-align: center;\" width=\"132\">High (lower than 0.01 mg\/ml)<\/td>\n<td style=\"text-align: center;\" width=\"113\">34.14<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"143\">anthraquinone<\/td>\n<td style=\"text-align: center;\" width=\"216\">cys566, ser596, tyr604, gln606, met636<\/td>\n<td style=\"text-align: center;\" width=\"132\">Low<\/td>\n<td style=\"text-align: center;\" width=\"113\">34.14<\/td>\n<\/tr>\n<tr>\n<td style=\"text-align: center;\" width=\"143\">Mitoxantrone<\/td>\n<td style=\"text-align: center;\" width=\"216\">arg298, leu300, pro533, gly534, tyr564, cys566, asp572,ser596, tyr604, val605, met636, asn635<\/td>\n<td style=\"text-align: center;\" width=\"132\">High (higher than 0.06 mg\/ml)<\/td>\n<td style=\"text-align: center;\" width=\"113\">152.11<\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\">Caption: Contact residues of enzyme were performed as results of docking analysis. Ligand properties parameters (water solubility and polar surface area) were measured using marvin sketch software.<\/p>\n<p style=\"text-align: justify;\">The presence of methyl group on each quinone, seems not affect to the affinity (Table 2, Figure 2). Other factors such as water solubility and Polar Surface area (PSA) do not appear to contribute to the affinity of ligand to P450. (Table 3) Docking analysis of quinones to other reductase enzymes (NQO1 and AIF) indicate the similar result.<\/p>\n<p style=\"text-align: justify;\"><strong>Binding energy analysis showed that generally, quinones prefere to bind to cytochrome P450 reductase rather than to NQO1 or AIF.<\/strong><\/p>\n<p style=\"text-align: justify;\">In general, quinones prefer bound to P450 reductase rather than to AIF or NQO1 (Table 2). It may be affected by reduction potential value of each quinones. Benzoquinone has reduction potential higher followed by naphthoquinone then antraquinone.<sup>19<\/sup>Reduction process is often achieved stepwise via one-electron reduction. The intermediate product (semiquinone molecules) is a free radical. The reduction potentials of such one-electron couples (quinone\/semiquinone) are value in predicting the direction of many free-radical reactions. Meaning, quinone which has higher reduction potential\u2019s value tends to interact with enzymes that facilitate one-electron reduction. It can explain why all ligands (benzoquinone or naphthaquinone or anthraquinone) have better affinity to P450 reductase then followed by AIF, rather than NQO1. This result confirmed the ability of AIF to catalyze one-electron reduction.<sup>19<\/sup>Even, 2-methyl-1,4-naphthoquinone showed different phenomenon. Here, 2-methyl-1,4-naphthoquinone interacts to AIF better than interacts to P450 reductase.Analysis of ligand-receptor interaction revealed also that there is a residue that consistent interact to each quinones. Tyr604 support quinone interaction using Van der Waals interaction. It appears that Tyr604 pose 90<sup>o<\/sup> to central ring on each quinone. Thus, Tyr604 residues may acts as recognizer residue the target compound of P450 reductase. (Table 3)Here also showed that the energy values of quinones-AIF close to quinones-P450 reductase (rather than to quinones-NQO1). The structural alignment of three enzymes was applied to explain this situation. FAT CAT program revealed that the folding of AIF close to P450 reductase rather than NQO1. (Figure 4)<\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig4.jpg\"><img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-39177\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig4-150x150.jpg\" alt=\"Vol14No2_Mol_Hes_fig4\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig4-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig4-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig4.jpg 738w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/a><\/td>\n<td><strong>Figure 4: Folding alignment between AIF (3GD4) to P450 reductase (1AMO) or NQO1 (2F1O). AIF has 65% folding similarity to P450 reductase (A). AIF has 35% folding similarity to NQO1 (B).<\/strong><\/p>\n<p><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig4.jpg\" target=\"_blank\">Click here to view figure<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\"><strong>Residues analysis confirmed that reduction of quinone is NAD(P)H dependent<\/strong><\/p>\n<p style=\"text-align: justify;\">Several contact residues to the binding of NAD(P)H as redox couple of quinones, were used to do analysis NAD(P)H-dependent on quinones-enzymes interaction. (Table 2) Here, 2-methyl-1,4-benzoquinone was used as a ligand model in order to understand the relationship between ligand to its redox couple on each reductase enzymes. Receptor-ligand interaction analysis showed that contact residues of target ligand are fit with contact residues of its couple redox. (Figure 5)This result revealed that the target ligand has high possibility to compete with redox couples, although it requires more support factor to overcome the obstacle threshold (ex. concentration), given the affinity energy of 2-methyl-1,4-benzoquinone is under their redox partner.These results confirmed that AIF can act as NADH:quinone reductase.<sup>19<\/sup><\/p>\n<table style=\"width: 70%;\" border=\"1\" cellpadding=\"5\">\n<tbody>\n<tr>\n<td><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig5.jpg\"><img decoding=\"async\" class=\"alignnone size-thumbnail wp-image-39178\" src=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig5-150x150.jpg\" alt=\"Vol14No2_Mol_Hes_fig5\" width=\"150\" height=\"150\" srcset=\"https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig5-150x150.jpg 150w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig5-256x256.jpg 256w, https:\/\/biomedpharmajournal.org\/staging\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig5.jpg 784w\" sizes=\"(max-width: 150px) 100vw, 150px\" \/><\/a><\/td>\n<td><strong>Figure 5: Docking interaction analysis of of 2-methyl-1,4-benzoquinone to P450 reductase, NQO1, and AIF<\/strong><\/p>\n<p><a href=\"https:\/\/biomedpharmajournal.org\/wp-content\/uploads\/2021\/06\/Vol14No2_Mol_Hes_fig5.jpg\" target=\"_blank\">Click here to view figure\u00a0<\/a><\/td>\n<\/tr>\n<\/tbody>\n<\/table>\n<p style=\"text-align: justify;\"><strong>Conclusion<\/strong><\/p>\n<p style=\"text-align: justify;\">Exploration of compounds that has capacity to kill cancer cells is challenge. One is to find potential anti-cancer agent which is associated with its ability to stimulate oxidative stress through generate ROS production. Given that cancer cells have higher threshold of ROS than normal cells, finding the ROS agent is absolutely necessary. One type of compound that is widely explored for its anticancer abilities are quinone family. Quinone have reduction capacity to form radical semiquinone or the less-reactive compound, hydroquinone. These reductions is associated with reducing enzymes including P450 reductase, NQO1, and AIF. This report revealed that in general quinones prefer bound to cytochrome P450 reductase rather than NQO1 or AIF. Several factors indicated influence to this process; the amount of quinone ring and reduction potential value. Functional groups of the ligand and Polar Surface area (PSA) do not appear to contribute to the affinity. The lower affinity were found for quinones-NQO1 interaction compared to quinones-AIF or Quinones-P450 reductase interaction. Residues analysis confirmed that reduction of quinone is NAD(P)H:dependent. Further analysis revealed that all quinones have high possibility to compete with their redox couples, although it requires more support factor to overcome the obstacle threshold. Finally, all describes that all quinones on this study have the potential as anticancer.<\/p>\n<p style=\"text-align: justify;\"><strong>Acknowledgment<\/strong><\/p>\n<p style=\"text-align: justify;\">The authors would like to humbly acknowledge The Indonesian Ministry of Research, Technology and Higher Education (Ristekdikti) for financial support during the research.<\/p>\n<p style=\"text-align: justify;\"><strong>Conflict of Interest<\/strong><\/p>\n<p style=\"text-align: justify;\">No potential conflict of interest was reported by the authors.<\/p>\n<p style=\"text-align: justify;\"><strong>Funding Source<\/strong><\/p>\n<p style=\"text-align: justify;\">The present study was funded by The Indonesian Ministry of Research, Technology and Higher Education (Ristekdikti).<\/p>\n<p style=\"text-align: justify;\"><strong>References<\/strong><\/p>\n<ol>\n<li style=\"text-align: justify;\">O\u2019Brien, P. J. Molecular mechanisms of quinone cytotoxicity. <em>Chem. Interact.<\/em>80, 1\u201341 (1991).<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/0009-2797(91)90029-7\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Morgan WA1, Hartley JA, C. G. Quinone-induced DNA single strand breaks in rat hepatocytes and human chronic myelogenous leukaemic K562 cells. <em>Biochem Pharmacol.<\/em>44 (2), 215\u201321 (1992).<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/0006-2952(92)90003-2\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Pelicano Helena, Carney Dennis, H. P. ROS stress in cancer cells and therapeutic implications. <em>Drug Resist. Updat.<\/em>7, 97\u2013110 (2004).<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/j.drup.2004.01.004\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">R., S. E. Metal ion-catalyzed oxidation of proteins: biochemical mechanism and biological consequences. <em>Free Radic Biol Med<\/em>9, 315\u201325 (1990).<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/0891-5849(90)90006-5\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Hensley Kenneth, Robinson Kent A., Gabbita S. Prasad, Salsman S., F. R. A. Reactive oxygen species, cell signaling, and cell injury. <em>Free Radic Biol Med<\/em>28, 1956\u201362 (2000).<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/S0891-5849(00)00252-5\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Hileman E.O., Liu J., Albitar M., Keating M.J., H. P. Intrinsic oxidative stress in cancer cells: a biochemical basis for therapeutic selectivity. <em>Cancer Chemother Pharmacol<\/em>53, 209\u201319 (2004).<br \/>\n<a href=\"https:\/\/doi.org\/10.1007\/s00280-003-0726-5\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Behrend L., Henderson G., Z. R. M. Reactive oxygen species in oncogenic transformation. <em>Biochem. Soc. Trans.<\/em>31, 1441\u20134 (2003).<br \/>\n<a href=\"https:\/\/doi.org\/10.1042\/bst0311441\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Ehninger Gerhard, Schuler Ulrich, Proksch Barbara, Zeller Klaus-Peter, B. J. Pharmacokinetics and Metabolism of Mitoxantrone A Review. <em>Clin. Pharmacokinet.<\/em>18, 365\u201380 (1990).<br \/>\n<a href=\"https:\/\/doi.org\/10.2165\/00003088-199018050-00003\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Faulds Diana, Balfour Julia A., Chrisp Paul, L. H. D. Mitoxantrone A Review of its Pharmacodynamic and Pharmacokinetic Properties, and Therapeutic Potential in the Chemotherapy of Cancer. <em>Drugs<\/em>41, 400\u2013449 (1991).<br \/>\n<a href=\"https:\/\/doi.org\/10.2165\/00003495-199141030-00007\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Daugas, E. <em>et al.<\/em> Apoptosis-inducing factor (AIF): A ubiquitous mitochondrial oxidoreductase involved in apoptosis. <em>FEBS Lett.<\/em>476, 118\u2013123 (2000).<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/S0014-5793(00)01731-2\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Wiraswati, H. L. <em>et al.<\/em> Apoptosis inducing factor (AIF) mediates lethal redox stress induced by menadione. <em>Oncotarget<\/em>7, 76496\u201376507 (2016).<br \/>\n<a href=\"https:\/\/doi.org\/10.18632\/oncotarget.12562\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Di Monte Donato, Ross David, Bellomo Giorgio, Eklow Lena, O. S. Alterations in intracellular thiol homeostasis during the metabolism of menadione by isolated rat hepatocytes. <em>Biochem Biophys<\/em>235, 334\u2013342 (1984).<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/0003-9861(84)90206-6\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Wang, X. <em>et al.<\/em> Mechanism of arylating quinone toxicity involving Michael adduct formation and induction of endoplasmic reticulum stress. <em>Proc. Natl. Acad. Sci.<\/em>103, 3604\u20133609 (2006).<br \/>\n<a href=\"https:\/\/doi.org\/10.1073\/pnas.0510962103\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Scott, G. Vitamin K3 (menadione) induced oncosis associated with keratin 8 phosphorylation and histone H3 arylation. <em>Mol. Pharmacol.<\/em>68, 606\u2013615 (2005).<br \/>\n<a href=\"https:\/\/doi.org\/10.1124\/mol.105.013474\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Wiraswati, H. L., Martoprawiro, M. A. &amp; Warganegara, F. M. Apoptosis Inducing Factor ( AIF ) Stabilizes Menadione- Conjugate Product in Programmed Cell Death. 10, 237\u2013245 (2017).<br \/>\n<a href=\"https:\/\/doi.org\/10.20902\/IJPTR.2017.10430\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Rockwell S., Sartorelli A. C., Tomasz M., K. K. A. Cellular pharmacology of quinone bioreductive alkylating agents. <em>Cancer Metastasis Rev.<\/em>12, 165\u201376 (1993).<br \/>\n<a href=\"https:\/\/doi.org\/10.1007\/BF00689808\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Gibson N.W., Hartley J.A, Siegel D., R. D. Relationship between DT-diaphorase-mediated metabolism of a series of aziridinylbenzoquinones and DNA damage and cytotoxicity. <em>Mol. Pharmacol.<\/em>42, 531\u201336 (1992).<\/li>\n<li style=\"text-align: justify;\">Mustofa, H. N., Wiraswati, H. L., Ekawardhani, S., Fianza, P. I. &amp; Windria, S. Anticancer Activities of Saponins and Quinones Group through Oxidative Stress and Glycolysis Inhibition via in Silico Studies. <em>Biomed. Pharmacol. J.<\/em>13, 999\u20131010 (2020).<br \/>\n<a href=\"https:\/\/doi.org\/10.13005\/bpj\/1969\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Iyanagi Takashi, Y. I. One-electron-transfer reactions in biochemical systems V. Difference in the mechanism of quinone reduction by the NADH dehydrogenase and the NAD(P)H dehydrogenase (DT-diaphorase). <em>Biochim. Biophys. Acta<\/em>216, 282\u201394 (1970).<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/0005-2728(70)90220-3\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Misevi\u010dien, L., Anusevi\u010dius, \u017d., \u0160arlauskas, J., Sevrioukova, I. F. &amp; \u010cnas, N. Redox reactions of the FAD-containing apoptosis-inducing factor (AIF) with quinoidal xenobiotics: A mechanistic study. <em>Arch. Biochem. Biophys.<\/em>512, 183\u2013189 (2011).<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/j.abb.2011.05.015\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">L., G. P. The role of NAD(P)H oxidoreductase (DT-Diaphorase) in the bioactivation of quinone-containing antitumor agents. <em>Free Radic. Biol. Med.<\/em>29, 263\u201375 (2000).<br \/>\n<a href=\"https:\/\/doi.org\/10.1016\/S0891-5849(00)00314-2\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Wang, M. <em>et al.<\/em> Three-dimensional structure of NADPH-cytochrome P450 reductase: Prototype for FMN- and FAD-containing enzymes. <em>Proc. Natl. Acad. Sci.<\/em>94, 8411\u20138416 (1997).<br \/>\n<a href=\"https:\/\/doi.org\/10.1073\/pnas.94.16.8411\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Asher Gad, Dym Orlt, Tsvetkov Peter, Adler Julia, S. Y. The crystal structure of NAD(P)H quinone oxidoreductase 1 in complex with its potent inhibitor dicoumarol. <em>Biochemistry<\/em>45, 6372\u201378 (2006).<br \/>\n<a href=\"https:\/\/doi.org\/10.1021\/bi0600087\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Trott Oleg, O. A. J. autodock vina:improving the speed and accuracy of docking with a new scoring function, efficient optimization and multithreading. <em>J Comput Chem.<\/em>31, 455\u201361 (2010).<br \/>\n<a href=\"https:\/\/doi.org\/10.1002\/jcc.21334\" target=\"_blank\">CrossRef<\/a><\/li>\n<li style=\"text-align: justify;\">Ye Y., G. A. Flexible structure alignment by chaining aligned fragment pairs allowing twists. <em>Bioinformatics<\/em>2, 246\u201355 (2003).<br \/>\n<a href=\"https:\/\/doi.org\/10.1093\/bioinformatics\/btg1086\" target=\"_blank\">CrossRef<\/a><\/li>\n<\/ol>\n","protected":false},"excerpt":{"rendered":"<p>Introduction Drugs belonging to the quinone family have been successfully  [&#8230;]<\/p>\n","protected":false},"author":14,"featured_media":0,"comment_status":"closed","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[92],"tags":[],"class_list":["post-39171","post","type-post","status-publish","format-standard","hentry","category-vol14no2"],"_links":{"self":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/39171","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/users\/14"}],"replies":[{"embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/comments?post=39171"}],"version-history":[{"count":5,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/39171\/revisions"}],"predecessor-version":[{"id":39848,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/posts\/39171\/revisions\/39848"}],"wp:attachment":[{"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/media?parent=39171"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/categories?post=39171"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/biomedpharmajournal.org\/staging\/wp-json\/wp\/v2\/tags?post=39171"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}