Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptHHS Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Bioorg Med Chem. Author manuscript; available in PMC 2012 December 1.
Published in final edited form as:
PMCID: PMC3216315

Metabolic and Electrochemical Mechanisms of Dimeric Naphthoquinones Cytotoxicity in Breast Cancer Cells


Cancer cells reprogram their metabolism due to genetic alteration to compensate for increased energy demand and enhanced anabolism, cell proliferation, and protection from oxidative damage. Here, we assessed the cytotoxicity of three dimeric naphthoquinones against the glycolytic MCF-7 versus the oxidative MDA-453 breast carcinoma cell lines. Dimeric naphthoquinones 1 and 2 impaired MDA-453, but not MCF-7, cell growth at IC50 = 15 μM. Significant increase in reactive oxygen species, decrease in oxygen consumption and ATP production were observed in MDA-453 cells but not in MCF-7 cell. These findings suggest that oxidative stress and mitochondrial dysfunction are mechanisms by which these agents exert their cytotoxic effects. Cyclic voltammetry and semi-empirical molecular orbital calculations further characterized the electrochemical behavior of these compounds. These results also suggest that dimeric naphthoquinones may be used to selectively target cancer cells that depend on oxidative phosphorylation for energy production and macromolecular synthesis.

Keywords: Dimeric Naphthoquinones, Anticancer Agents, Cytotoxicity, Oxidative Stress, Tumor Metabolism

1. Introduction

Breast cancer is the most prevalent cancer in women in industrialized countries. In the United States, approximately 207,000 women were diagnosed with breast cancer and 40,000 died of it in 2010.1 Although advanced therapies have contributed to decreased mortality from breast cancer in recent years, the global burden of this heterogeneous disease highlights the need for new therapies against more aggressive and advanced forms of the disease.

In 2006, we prepared a series of novel dimeric naphthoquinones for the purpose of developing structurally less complex HIV-integrase inhibitors that mimicked the action of conocurvone, a naturally occurring trimeric naphthoquinone with potent anti-HIV activity.2 While the less complex dimeric naphthoquinones inhibited HIV-1 mediated cytopathogenicity, they also showed cytotoxicity in CEM-T4 lymphoblastic cells at varying concentrations, suggesting their possible use as anti-neoplastic agents.2

Recently, we reported that dimeric naphthoquinones possessed cytotoxic activity and synergistic effects with radiation against prostate cancer cells.3 Substitution patterns on the quinone double bond and aromatic ring of the quinone units were shown to influence cytotoxicity. Mechanistic studies using a chemical genetic screen in yeast revealed that their cytotoxicity was probably related to their ability to undergo redox cycling, produce free radicals and cause mitochondrial dysfunction.4

Targeting aerobic glycolysis (Warburg effect) has been the focus of numerous studies in cancer, but more recent data have identified a number of cancer cells that use mitochondria and oxidative phosphorylation as a main energy source.56 Our previous studies of the cytotoxicity of dimeric naphthoquinones in prostate cancer cells and yeast suggested that one mechanism by which these compounds affect cell survival may be through induction of mitochondrial dysfunction.34 We hypothesized that dimeric naphthoquinones would be selectively cytotoxic for cancer cells which utilize oxidative phosphorylation while cancers relying on glycolysis would not be affected by these agents.34

In the present study, we investigated how exposure to dimeric naphthoquinones affected MDA-453 and MCF-7 differently as two metabolically distinct breast cancer cell lines.7 We report that dimeric naphthoquinones 1–3 are selectively cytotoxic to a MDA-453, which depend on oxidative phosphorylation for energy production, with IC50’s in the micromolar range. We further show that exposure of MDA-453 breast cancer cells to dimeric naphthoquinones results in the accumulation of reactive oxygen species, depletion of intracellular oxygen and inhibition of ATP generation, while have no or minimal effect on MCF-7 breast cancer cells. The ability of the quinone dimers to undergo redox cycling was studied using cyclic voltammetry. Molecular orbital calculations were performed to predict lowest unoccupied molecular orbital (LUMO) energies of the dimers in an attempt to see if there is a correlation between their LUMO’s and cytotoxicity.

2. Results and Discussion

A high yielding one step synthesis of dimeric naphthoquinones 1–3 was described by us previously8 and has been subsequently utilized by other research groups910 to study the synthesis and biological properties of this unique class of compounds.

To determine whether dimeric naphthoquinones affected the viability of breast cancer cells, we performed a preliminary screen using three structurally related chloro-hydroxy dimeric naphthoquinones 1–3 in MDA-453 and MCF-7 breast cancer cell lines using an MTT assay. As mentioned previously, MDA-453 is a basal-type breast cancer cell line, which depends more on oxidative phosphorylation for survival, proliferation, invasion and energy production.7, 1112 MCF-7 is a luminal-type breast cancer cell line that depends more on glycolysis for survival.7, 1112 At concentrations of 40 μM or less dimeric naphthoquinones 1–3 significantly impaired conversion of MTT to formazan in MDA-453 breast cancer cells (p<0.05, compared to DMSO) but not in MCF-7 cells, suggesting selective, dose-dependent cytotoxicity by possibly affecting oxidative phosphorylation (Figure 2). Trypan blue exclusion assays for cell viability confirmed cell deaths. Dimeric quinones 1 and 2 showed cytotoxic effect against MDA-453 with IC50s of 15 and 17 μM, respectively, compared to the pyranylated dimeric quinone 3 with an IC50 of 34 μM. In contrast, approximately 80% of MCF-7 cells survived at 50 μM of dimeric naphthoquinones concentration (1, 80±4%; 2, 83±5%; 3, 85±7%).

Figure 2
Effect of dimeric naphthoquinones 1–3 on MDA-453 (A) and MCF-7 (B) breast cancer cell proliferation (compared to DMSO). All cells were grown at 1×105 cells/mL. Cell counts were performed in triplicate and shown as mean ± SD, and ...

To test if dimeric naphthoquinones selectively affect energetic metabolism in cancer cells by disrupting mitochondrial function, we performed a series of metabolic experiments in MDA-453 and MCF-7 breast cancer cells. Indeed, after treatment of MDA-453 breast cancer cells with dimeric naphthoquinones at a concentration of 10μM, quinones 1 and 2 induced a 53% and a 37% decrease in ATP in the remaining viable cells compared to DMSO, respectively (Figure 3). At 10 μM, quinone 3 did not decrease ATP production significantly. When the quinone concentration was increased to 20 μM, quinones 1 and 2 decreased ATP production by approximately 60%, while quinone 3 induced a 27% decrease in ATP compared to DMSO (Figure 3). As expected, treatment of glycolysis-dependent MCF-7 cells with dimeric naphthoquinones affected ATP production to a much lesser extent (Supplemental Figure 1). At concentrations greater than 20 μM, ATP production was not further decreased, suggesting the presence of other mechanisms involved in cytotoxic effects of these compounds.

Figure 3
Relative (compared to DMSO) decrease in ATP production in MDA-453 breast cancer cells by dimeric naphthoquinones 1–3.

To test if reactive oxygen species (ROS) increased in breast cancer cells treated with dimeric naphthoquinones, we preloaded MDA-453 breast cancer cells with 5-carboxy-2′,7′-dichlorodihydrofluorescein diacetate (DCFDA) followed by exposure to dimeric naphthoquinones. Intracellular oxidation of DCFDA by ROS yields a fluorescent compound which can be measured using flow cytometry. All three dimeric naphthoquinones increased production of ROS in MDA-453 cells as shown by a 50–63% (p<0.05) absolute increase in the percentage of fluorescent cells compared to DMSO (20%), suggesting that oxidative stress may contribute to cytotoxicity through free radical production in these cell lines (Figure 4).

Figure 4
Generation of free oxygen radicals by dimeric naphthoquinones 1–3 compared to DMSO in MDA-453 cells. ROS levels were determined by DCFDA fluorescence in the cells. Data are representative of triplicate samples of two separate experiments. For ...

Interestingly, when we performed similar experiments in MCF-7 cells, dimeric naphthoquinones did not increase ROS production significantly (Supplemental Figure 2). Furthermore, when we measured O2 consumption by MDA-453 cells in the presence and absence of dimeric naphthoquinones, we found that dimeric quinones 13 decreased O2 consumption to different degrees, suggesting interference with mitochondrial oxidative phosphorylation. As expected, dimeric naphthoquinones did not alter O2 consumption in MCF-7 cells which rely on non-oxygen non-mitochondrial glycolysis for survival (supplemental Figure 3).11 These results suggest that dimeric naphthoquinones may be used to selectively target cancer cells that depend on oxidative phosphorylation for energy production and macromolecular synthesis.

In view of the ability of dimeric naphthoquinones to selectively generate ROS in MDA-453 breast cancer cell lines, it was of interest to examine the electrochemical behavior of these compounds in order to determine whether they had the inherent ability to undergo reversible oxidation-reduction reactions. The correlation of electrochemical potentials of monomeric quinones and their inhibitory effects on Epstein–Barr virus activation,1315 and cytotoxic activities1617 have been reported. However, to the best of our knowledge, such correlations have not been investigated in dimeric naphthoquinones. Thus, 2mM stock solutions of the dimeric quinones were prepared in acetonitrile and analyzed by cyclic voltammetry using a standard Ag/AgCl reference electrode and tetrapropyl ammonium bromide as the supporting electrolyte.18 The cyclic voltammograms are shown in Figure 5.

Figure 5
Cyclic Voltammetry of 2 mM dimeric naphthoquinones in acetonitrile plus 0.1 M tetrapropyl ammonium bromide (TBAP), with a Ag/AgCl reference electrode (RE), and the scan rates between 0.1 volt per second to 2 volts per second, E = electrochemical potential, ...

The electrochemical behavior of 2-hydroxynaphthoquinones are complex.19 Quinones are known to undergo disproportionation reactions, which could account for the non-uniform peak shapes in the cyclic voltammograms of the hydroxyquinone dimers. In Figure 5, the first and second waves correspond to formation of the radical anion and dianion (presumably by reduction of the more electron deficient chloronaphthoquinone unit first), while the third wave corresponds to radial anion formation of the more electron rich hydroxyquinone unit at higher potential. The fourth wave was not definable. Table 1 shows cathodic (Ec) and anodic (Ea) peak potentials for each of the four redox steps (E1–E4) as well as reduction potentials (E1/2 = Ea + Ec/2).2021 Similar values for E11/2 reduction potentials for each of the dimeric quinones supports the idea that reduction of the chloroquinone occurs first because all three quinones dimers contain a structurally similar chloroquinone unit (i.e. with no substituents on the aromatic ring). The peaks selected for Table 1 are depicted in the supplemental Figure 4.

Table 1
Cathodic (Ec) and Anodic (Ea) Potentials

To study whether regioselective reduction of the chloroquinone unit could be predicted, we used semi-empirical molecular orbital calculations to determine the partial charges and molecular orbital energies of the quinone dimers (Figure 6).2224 The nearly identical partial charges of the carbon atoms of the four carbonyl groups at C-1, C-4, C-6 and C-8 (ranges between 0.34 to 0.39) revealed that oxygen n-donation to the carbonyl π-system was similar in all three dimeric naphthoquinones. Furthermore, the charges of the chlorinated carbons at C-2 (i.e. 0.17–0.18) were nearly identical in all three of the quinone dimers making it difficult to predict selective reduction or correlate cytotoxicity merely from partial charge distributions. Interestingly, calculation of the lowest unoccupied molecular orbital (LUMO) and the second LUMO (SLUMO) energies of the quinones (Figure 6) revealed a slightly narrower LUMO-SLUMO orbital energy gap [ΔE (LUMO − SLUMO)] in quinones 1 and 2 (5 kcal/mol) compared to quinone 3 (7.6 kcal/mol). This slightly narrow energy gap may make the SLUMO more accessible to the reacting cellular nucleophiles to form Michael adducts as additional mechanism for cytotoxicity.4 Images of the calculated LUMO and SLUMO for dimeric naphthoquinones are presented in supplemental Figure 5.

Figure 6
Semi-empirical partial charges and molecular orbital energy calculations for dimeric naphthoquinones using the PM3 Hamiltonian carried out with the Gaussian03 program package. All geometries were initially sketched and cleaned and molecular orbitals were ...

In conclusion, the potential mechanisms of cytotoxicity of three structurally related 2-chloro-2′-hydroxyl dimeric naphthoquinones in MDA-453 breast carcinoma cells was investigated. Dimeric naphthoquinone cytotoxicity is likely mediated by oxidative stress and mitochondrial dysfunction as evidenced by increased levels of ROS, and decreased ATP production and oxygen consumption. Cyclic voltammetry and semi-empirical molecular orbital calculations assisted to further characterize the electrochemical behavior of the quinone dimers. Our data suggests that participation in redox cycling by dimeric naphthoquinones generate significant levels of reactive oxygen radicals which may contribute to their cytotoxicities. Further in vitro and in vivo studies are warranted to continue to explore anti-neoplastic properties of this class of compounds in different cancer cell lines as single agents and in combination with other chemotherapeutics.

Material and Methods

2.1. Chemistry

Binaphthoquinones were synthesized and characterized using previously described methods, and dissolved in DMSO.8 We already performed and reported a comprehensive study of biological activities and physical properties of halohydroxy monomeric naphthoquinones. None of the monomeric quinones showed anti-viral or anti-neoplastic activities.2, 8, 2526

2.2. Cell lines

Breast cancer cell lines, MCF-7 and MDA-453 were grown in a humidified incubator at 37°C and 5% CO2 in DMEM media (Life Technology, Rockville, Maryland) containing 10% fetal bovine serum and 1% antibiotics.

2.3. Cytotoxic activity assays

Cell survival was assessed by MTT assays and trypan blue exclusion.27 For MTT assays, an equal numbers of cells were plated in three or more replicates and allowed to attach overnight. Cells were then treated with the vehicle (DMSO) or individual solutions of dimeric naphthoquinones at 5, 10, 20 and 50 μM in complete media for 4 hours. Following treatment, cells were incubated with 0.5 mg/ml MTT in media for 1 hour and then the media was removed and PBS was added. Absorbance was assayed at 570 nm (signal) and 690 nm (background).

2.4. ATP measurement

Cells were plated at 10,000 cells per well in a 96-well plate and treated in sextuplet with either DMSO or 5, 10, 20 or 50 μM of dimeric naphthoquinones for 4 hours. ATP production was measured using the Cell Titer-Glo Luminescent assay (Promega) and adjusted for viable cell number. The resultant luminescence was normalized to the average number of viable, trypan blue excluding cells per treatment condition.

2.5. Reactive Oxygen Species (ROS) measurement

Flow cytometric analysis of ROS production was performed using a FACScan instrument and CellQuest Pro software (Becton Dickinson) as previously described.11 To measure reactive oxygen species generation, breast cancer cells were preloaded with 5μM of 5-(and-6)-carboxy-2′,7′-dichlorodihydrofluorescein diacetate (DCFDA) for 30 minutes. Cells were then treated with DMSO or 20 μM of naphthoquinones for 2 hours, collected and analyzed using an excitation wavelength of 488 nm produced by an argon laser and measuring emission at 530nm.

2.6. Oxygen Consumption

Oxygen consumption was measured using a Clark-type oxygen electrode (Oxytherm System; Hansatech Instruments Ltd).11 Cells were treated with 20 μM dimeric naphthoquinones or DMSO for 3hours. Then the media were removed, the cells were detached using trypsin, washed with PBS and spun down. Then, 5×106 cells in 0.8 mL of PBS containing 10% FCS were placed in a chamber above a membrane that is permeable to oxygen. Oxygen diffuses through the membrane and is reduced at the cathode surface so that a current flows through the circuit, which is completed by a thin layer of KCl solution. The current that is generated bears a direct stoichiometric relation to the oxygen reduced and is converted to a digital signal. Determinations were done in duplicate.

2.7. Cyclic Voltammetry (CV)

Cyclic voltammetry was performed on a CV 130 system from Gamry Instruments (Warminster, PA, USA).28 The electrochemical cells were installed in a Faraday cage to reduce externally generated noise in the sensor circuits. The initial time delay before starting CV was 30 seconds. A three electrode cell (gold disk working electrode (WE), platinum wire counter electrode (CE), and Ag/AgCl reference electrode (RE)13, 29) was used for all electrochemical experiments. The WE was scanned between −0.5V and +0.5V or less versus RE. Resolution was set to collect a current data point every 2 mV. Under this condition, the scan rates ranged from 0.1 volt per second to 2 volts per second as indicated in the data. All solutions were deoxygenated for 10 minutes by bubbling of nitrogen gas. A blanket of nitrogen was maintained over the solution during all experiments. Test solutions contained 2 mM dimeric naphthoquinones and 0.1 M tetrapropyl ammonium bromide (TPAB) in anhydrous acetonitrile. The potentials of oxidation (anodic E, Epa) and reduction (cathodic E, Epc) values were measured. Eredox is calculated as Epa+Epc/2.2021

2.8. Semi-empirical molecular orbital calculations

Semi-empirical molecular orbital calculations using the PM3 Hamiltonian30 were carried out with the Gaussian03 program package.31 All geometries were initially sketched and cleaned and molecular orbitals were generated and viewed using GaussView 3.09. Geometry optimizations were performed in Gaussian03 using the very tight optimization convergence criteria with no symmetry constraints.

2.9. Statistical analysis

Statistical analyses were performed using Stata® Software version 10.1 (StataCorp, College Station, Tx). All p-values are two-sided and those < 0.05 were considered as statistically significant.

Figure 1
Dimeric naphthoquinones.

Supplementary Material


Supp. Figure 1:

Relative (compared to DMSO) decrease in ATP production in MCF-7 breast cancer cells by dimeric naphthoquinones 1–3.


Supp. Figure 2:

Generation of reactive oxygen species, which were determined by DCFDA fluorescence in MCF-7 cells untreated (0.02%) and treated with DMSO (0.3%), dimeric naphthoquinones 1 (19%), 2 (10%) and 3 (10%). Data are representative of triplicate samples of two separate experiments.


Supp. Figure 3:

Comparison of oxygen consumption of MDA-453 and MCF-7 cells was determined by a Clark-type oxygen electrode in the presence and the absence of dimeric naphthoquinones. Data are representative of duplicate experiments.


Supp. Figure 4:

Cyclic voltammetry of dimeric naphthoquinones 1–3 with demonstration of the peaks selected for Table 1.


Supp. Figure 5:

Images of the calculated LUMO and SLUMO for dimeric naphthoquinones 13. In quinones 1 and 2, the LUMO is localized on the quinone ring that contains hydroxyl group and the SLUMO is localized on the chloroquinone ring.


Funding for this study was provided in part by a Johns Hopkins Hospital Department of Pathology Start-up Fund and Richard Starr Ross Clinician Scientist Award (to MVR), and National Institute of Health (grant AI43687 and T-32 grant to the Division of Hematology, Department of Internal Medicine, Johns Hopkins University) to AE.


Publisher's Disclaimer: This is a PDF file of an unedited manuscript that has been accepted for publication. As a service to our customers we are providing this early version of the manuscript. The manuscript will undergo copyediting, typesetting, and review of the resulting proof before it is published in its final citable form. Please note that during the production process errors may be discovered which could affect the content, and all legal disclaimers that apply to the journal pertain.

References and notes

1. Jemal A, Siegel R, Xu J, Ward E. CA Cancer J Clin. 60:277. [PubMed]
2. Stagliano KW, Emadi A, Lu Z, Malinakova HC, Twenter B, Yu M, Holland LE, Rom AM, Harwood JS, Amin R, Johnson AA, Pommier Y. Bioorg Med Chem. 2006;14:5651. [PubMed]
3. Ross AE, Emadi A, Marchionni L, Hurley PJ, Simons BW, Schaeffer EM, Vuica-Ross M. BJU Int. 2010
4. Emadi A, Ross AE, Cowan KM, Fortenberry YM, Vuica-Ross M. PLoS One. 2010;5:e10846. [PMC free article] [PubMed]
5. Dang CV. Mol Cell Biol. 30:1300. [PMC free article] [PubMed]
6. Koppenol WH, Bounds PL, Dang CV. Nat Rev Cancer. 11:325. [PubMed]
7. Mazurek S, Michel A, Eigenbrodt E. J Biol Chem. 1997;272:4941. [PubMed]
8. Emadi A, Harwood JS, Kohanim S, Stagliano KW. Org Lett. 2002;4:521. [PubMed]
9. Crosby IT, Rose ML, Collis MP, de Bruyn PJ, Keep PLC, Robertson AD. Australian J Chem. 2008;61:768.
10. Crosby IT, Bourke DG, Jones ED, de Bruyn PJ, Rhodes D, Vandegraaff N, Cox S, Coates JA, Robertson AD. Bioorg Med Chem. 2010;18:6442. [PubMed]
11. Le A, Cooper CR, Gouw AM, Dinavahi R, Maitra A, Deck LM, Royer RE, Vander Jagt DL, Semenza GL, Dang CV. Proc Natl Acad Sci U S A. 2010;107:2037. [PubMed]
12. Kung HN, Marks JR, Chi JT. PLoS Genet. 2011;7:e1002229. [PMC free article] [PubMed]
13. Koyama J, Morita I, Kobayashi N, Osakai T, Hotta H, Takayasu J, Nishino H, Tokuda H. Cancer Lett. 2003;201:25. [PubMed]
14. Koyama J, Morita I, Kobayashi N, Osakai T, Hotta H, Takayasu J, Nishino H, Tokuda H. Cancer Lett. 2004;212:1. [PubMed]
15. Koyama J, Morita I, Kobayashi N, Osakai T, Nishino H, Tokuda H. Cancer Lett. 2005;225:193. [PubMed]
16. Pan SS, Gonzalez H. Mol Pharmacol. 1990;37:966. [PubMed]
17. Wu HQ, Huang ZS, Bu XZ, Shen YD, Zhang ZL, Xie BF, Liu ZC, Gu LQ, Chan AS. Eur J Med Chem. 2005;40:1341. [PubMed]
18. El Maali NA, Berg H. Electroanalysis. 1996;8:808.
19. Aguilar-Martinez M, Macias-Ruvalcaba NA, Bautista-Martinez JA, Gomez M, Gonzalez FJ, Gonzale I. Current Org Chem. 2004;8:1721.
20. Crawford PW, Carlos E, Ellegood JC, Cheng CC, Dong Q, Liu DF, Luo YL. Electrochimica. 1996;41:2399.
21. Rieger PH. Electrochemistry. New York: Chapman and Hall; 1994.
22. Stefanou V, Matiadis D, Melagraki G, Afantitis A, Athanasellis G, Igglessi-Markopoulou O, McKee V, Markopoulos J. Molecules. 2011;16:384. [PubMed]
23. Benigni R. Chem Rev. 2005;105:1767. [PubMed]
24. Cavelier G, Amzel LM. PROTEINS: Structure, Function, and Genetics. 2001;43:420. [PubMed]
25. Emadi A. PhD. Department of Biological, Chemical and Physical Sciences; 2004.
26. Ross AE, Emadi A, Marchionni L, Hurley PJ, Simons BW, Schaeffer EM, Vuica-Ross M. British Journal of Urology International. 2010 in press.
27. Mosmann T. J Immunol Methods. 1983;65:55. [PubMed]
28. Cunningham AJ, Underwood AL. Biochemistry. 1967;6:266. [PubMed]
29. Ferraz PAL, de Abreu FC, Pinto AV, Glezer V, Tonholo J, Goulart MOF. Journal of Electroanalytical Chemistry. 2001;507:275.
30. Stewart JJP. J Comp Chem. 1989;10:221.
31. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Montgomery JA, Jr, Vreven T, Kudin KN, Burant JC, Millam JM, Iyengar SS, Tomasi J, Barone V, Mennucci B, Cossi M, Scalmani G, Rega N, Petersson GA, Nakatsuji H, Hada M, Ehara M, Toyota K, Fukuda R, Hasegawa J, Ishida M, Nakajima T, Honda Y, Kitao O, Nakai H, Klene M, Li X, Knox JE, Hratchian HP, Cross JB, Bakken V, Adamo C, Jaramillo J, Gomperts R, Stratmann RE, Yazyev O, Austin AJ, Cammi R, Pomelli C, Ochterski JW, Ayala PY, Morokuma K, Voth GA, Salvador P, Dannenberg JJ, Zakrzewski VG, Dapprich S, Daniels AD, Strain MC, Farkas O, Malick DK, Rabuck AD, Raghavachari K, Foresman JB, Ortiz JV, Cui Q, Baboul AG, Clifford S, Cioslowski J, Stefanov BB, Liu G, Liashenko A, Piskorz P, Komaromi I, Martin RL, Fox DJ, Keith T, Al-Laham MA, Peng CY, Nanayakkara A, Challacombe M, Gill PMW, Johnson B, Chen W, Wong MW, Gonzalez C, Pople JA. Gaussian03. 2004.