Search tips
Search criteria 


Logo of nihpaAbout Author manuscriptsSubmit a manuscriptNIH Public Access; Author Manuscript; Accepted for publication in peer reviewed journal;
Chem Biol Interact. Author manuscript; available in PMC Oct 23, 2011.
Published in final edited form as:
PMCID: PMC3199116
Efficacy ranking of triterpenoids as inducers of a cytoprotective enzyme and as inhibitors of a cellular inflammatory response via their electron affinity and their electrophilicity index
René V. Bensasson,a Vincent Zoete,b Gaston Berthier,c Paul Talalay,d and Albena T. Dinkova-Kostovade*
aMuséum National d'Histoire Naturelle, Molécules de Communication et Adaptation des Microorganismes, Paris, France
bSwiss Institute of Bioinformatics, Molecular Modeling Group, Lausanne, Switzerland
cLaboratoire de Chimie Théorique, Université Pierre et Marie Curie, Paris Cedex, France
dDepartment of Pharmacology and Molecular Sciences, Johns Hopkins University School of Medicine, Baltimore, MD, USA
eBiomedical Research Institute, University of Dundee, Dundee, United Kingdom
*Corresponding author : Albena T. Dinkova-Kostova, Biomedical Research Institute, University of Dundee, Dundee DD1 9SY, United Kingdom; Tel. 44 1382 740045; Fax: 44 1382 669993; a.dinkovakostova/at/
Electron affinity (EA) and electrophilicity index (ω) of 16 synthetic triterpenoids (TP), previously identified as inducers of cytoprotective enzymes and as inhibitors of cellular inflammatory responses, have been calculated by the molecular orbital method. Linear correlations were obtained by plotting the values of EA, as well as those of ω versus (i) the potencies of induction of NAD(P)H quinone reductase (NQO1, EC, a cytoprotective enzyme, expressed via the concentration of TP required to double the specific activity of NQO1 (CD value) and (ii) the values of their anti-inflammatory activity expressed via the IC-50 of TP for suppression of upregulation of inducible nitric oxide synthase (iNOS, EC, both previously experimentally determined. The observed correlations demonstrate quantitatively for a series of triterpenoids that their electrophilicity is a major factor determining their potency as inducers of the cytoprotective phase 2 response and as inhibitors of inflammatory processes.
Keywords: electron affinity, inflammation, NAD(P)H-quinone oxidoreductase 1, quantum chemical calculations, structure-activity relationships
Recent observations have contributed to a new understanding of the 19th century Rudolf Virchow's hypothesis of a link between inflammatory processes induced by chronic injury and a predisposition to certain cancers [1]. Indeed, a chronic injury can promote a cascade of events leading to mutations and carcinogenesis via (i) release of leukocytes, (ii) production by leukocytes of reactive oxygen and nitrogen species (ROS, RNS), which can lead to oxidative damage and mutations, (iii) induction of cell responses involving wound healing, angiogenesis and tissue remodeling [2,3]. A linear correlation has been observed between the anti-inflammatory potencies of a series of synthetic triterpenoid (TP) analogues of oleanolic acid and their ability to induce the phase 2 enzyme NAD(P)H-quinone reductase (NQO1, EC [4]. It has also been shown that this correlation extends to many other inducers that belong to structurally distinct chemical classes [5]. NQO1 is an obligatory two electron-reducing enzyme. Its gene expression is regulated coordinately with a battery of >100 cytoprotective enzymes via the Keap1/Nrf2/ARE pathway that protects against damaging electrophilic species, including endogenous ROS and RNS arising from aerobic metabolism, toxic electrophilic exogenous compounds, and electrophilic metabolites. The induction of NQO1 by test agents is a reliable biomarker of the ability of these agents to protect against tumor development in many animal models [6,7]. In our previous study [4] two bioassays were used: the first expressed the induction of the phase 2 enzyme NQO1 via the concentration of TP (called CD value) required to double the specific activity of NQO1 in the murine hepatoma cell line Hepa1c1c7, the second quantified the inhibition of a cellular inflammatory process via the IC-50 value of TP for suppression of induction of nitric oxide synthase (iNOS, EC by interferon-gamma (IFN-γ) in primary mouse macrophages. This study demonstrated that (i) a Michael reaction acceptor functionality, i.e., olefinic function conjugated to electron-withdrawing group(s) in the TP, is necessary for inducer activity of NQO1 and for blocking inflammation, (ii) induction of phase 2 enzymes by TP is selective and independent of phase 1 enzymes, (iii) induction requires the presence of both Keap1, the protein sensor for inducers, and transcription factor Nrf2, (iv) the most potent TP (TP-225) reacts with thiol groups of the Keap1 sensor; such reaction renders Keap1 unable to target transcription factor Nrf2 for ubiquitination and proteasomal degradation, which allows Nrf2 to bind to the antioxidant response element (ARE), and activate the transcription of NQO1 and other ARE-regulated genes. These studies also suggested that both abilities of TP to induce NQO1 and to inhibit inflammatory processes, which stimulate the inducible form of nitric oxide (NO) synthase (iNOS), could be contributing to the potent chemoprotective effects of TP, recently observed in several animal models [8-10]. Increased levels of NO can react with superoxide anions and produce peroxynitrite anion (ONOO-), an oxidant and a nitrating species, which can damage DNA and increase cancer risk [11].
The aim of the present investigation is to search for quantitative structure-activity relationships between the physico-chemical properties of TP and their potencies as inducers of NQO1 and as inhibitors of inflammation. Since these TP are electrophiles and interact, via their Michael acceptor functionalities, with sulfhydryl groups of Keap1, the property involved in this interaction should be the reduction potential E (TP/TP•-). A quantum mechanical calculation of the energy of the lowest unoccupied molecular orbital E (LUMO) represents a satisfactory measure of the electron-acceptor properties of the molecule in question [12]. This E (LUMO) is linearly correlated with the reduction potential E (TP/TP•-) in solution and with the electron affinity (EA) of TP in the gas phase as already demonstrated for a series of aromatic hydrocarbons [13].
In the present work, three possible theoretical approaches for estimating the electron affinity (EA) of our triterpenoid compounds are tested. One uses directly the E (LUMO) value of the neutral species TP, the second and the third consider a system where an electron is removed from the open-shell structure of the TP•- anion, using either the energy of the highest occupied molecular orbital, E (HOMO), or the energy of the singly occupied molecular orbital, E (SOMO), located between E (HOMO) and E (LUMO) on the axis of the orbital energies of TP•-. For ranking the efficacy of our 16 TP, another descriptor has been used, the electrophilicity index ω. Further comments concerning the determinations and usefulness of these two descriptors, electron affinity (EA) and electrophilicity index (ω), will be given in the section below dealing with computational methods.
In search of correlations, the three different values of EA, as well as those of ω of the triterpenoid compounds were plotted versus: (i) the induction potencies of NQO1 expressed via the concentration of TP (CD) required to double the specific activity of NQO1, and (ii) their anti-inflammatory activity expressed via the half maximal inhibitory concentration, IC-50, of TP for suppression of induction of nitric oxide synthase (iNOS), both previously experimentally determined [4,14,15].
In addition to the purely electronic factor, another physico-chemical property will be considered, the hydrophobic character of the TP expressed by log P, P representing the relative distribution in an octanol/water mixture obtained from the ratio of the concentration of a compound in octanol over its concentration in aqueous solution.
2.1 Electron affinity (EA)
It is customary to justify the theoretical procedures of electron affinity (EA) determination by referring to the so-called Koopmans theorem, specific of self-consistent field (SCF) treatments for closed-shell systems [16]. Formal relationships between E (HOMO), or the E (LUMO) energy levels of a closed-shell system and their energy change when an electron is respectively removed from the highest occupied molecular orbital (HOMO) or attached to the lowest unoccupied molecular orbital (LUMO) are given in classical textbooks of Quantum Chemistry [17]. Linear correlations between (i) E (HOMO) of aromatic hydrocarbons and their ionization potentials in the gas phase as well as their electrochemical oxidation potentials in solution, and (ii) E (LUMO) of aromatic hydrocarbons and their electron affinities (EA) in the gas phase and their electrochemical reduction potentials in solution are reported by Pullman and Pullman [18]. Moreover, the Koopmans theorem, valid for closed-shell structures, without unpaired electrons, is extended to radicalic open-shell systems, using E (SOMO) [19]. Some overlooked details of this Koopmans theorem have been clearly stated by Angeli [20].
Given the structural complexity of the TP molecules, self-consistent field (SCF) calculations by the AM1 quantum mechanical procedure have been carried out with the Hyperchem 7.51 program. The restricted Hartree-Fock (RHF) formalism, using standard AM1 basis sets of atomic orbitals was retained for computing the orbital energy levels of TP and TP•- systems. A set of filled and empty molecular orbitals is generated for each of them, from which were selected the relevant HOMO, LUMO and SOMO levels linearly correlated with electron affinities EA of (TP) and reduction potentials E (TP/TP•-). Let us present now some details on the validity of our approach from the theoretical point of view.
In the survey of Mulliken on molecular orbital theory [21], it was pointed out that the success of predictions based on Koopmans theorem for ionization potentials comes from a balance between two sources of error in these calculations, namely lack of self-consistency of the electronic and geometrical structure assigned to the final cation and neglect of correlation energies in the picture of both initial and final states. It can be shown without any calculation that these two terms have opposite signs, hence the Mulliken conclusion. Extending the argument of Mulliken to non-closed shell configurations is not immediate, but Koopmans theorem is probably valid, although to a minor extent, for the SOMO orbital of radicals as discussed by Dodds and McWeeny [22]. It is not possible to argue in the same way about the method using E (LUMO) of neutral TP for estimating electron affinities EA, because there is no reason why the two sources of error mentioned above should have opposite signs in this case, in contrast with the use of E (HOMO) or E (SOMO) for the estimation of ionization potentials. Thus, EA evaluations from orbital energies only lie finally in a comparison a posteriori between the calculation results and the experimental data. A good linear relationship observed between them in a series of chemically related molecules means that there is probably a quantity transferable from a compound to another, compensating all the deficiencies of the treatment. Let us add finally that the evaluation of electron affinities of molecules M via a calculation of the ionization energies of the charged species M•- obtained by putting an extra electron in a SOMO or a LUMO level, would be a valuable alternative [23]. However, the application of Koopmans theorem to the resulting anion TP•- implies that a SCF treatment of the charged species has been made effectively, in addition to the treatment carried out from the neutral TP, as it has been done in the present work by the AM1 method. The electron affinities (EA) of our TP reflect the capability of these agents to accept one electron from a donor.
2.2 Electrophilicity index (ω)
In the frame of the molecular orbital method, it is possible to describe the efficacy of TP compounds by means of another reactivity descriptor expressed in terms of E (HOMO) and E (LUMO) orbital energies of the neutral TP, namely the electrophilicity index (ω)
equation M1
where (χ) is the electronegativity and (η) the chemical hardness of the molecule. Use of ω was suggested empirically by Maynard et al. in their study of the reaction between the human immunodeficiency virus type 1 (HIV-1) nucleocapsid p7 (NCp7) zinc finger thiolates and different electrophilic ligands [24]. This was justified later by Parr et al. in the frame of the density functional theory (DFT) [25]. The alternative of using standard DFT calculations would be possible if theoretical values in agreement with experimental ionization potentials could be obtained in a simple way. Now the E (HOMO) DFT values are generated by the so-called non-interacting electron model and therefore are lower than the exact ones by a deviation ΔE, the evaluation of which would necessitate an intricated linear response procedure [26]. Moreover, the determination of realistic E (LUMO) values via DFT with respect to experimental EA is still debated [27]. Let us add that, in connection with the anti-inflammatory properties of some non-steroidal drugs, we have recently pointed out the fact that correlation coefficients between E (HOMO) and biological efficacy computed by our AM1 calculations are better than those computed by corresponding DFT or even ab initio methods [28]. Equations for the calculations of χ and η are :
equation M2
where I and EA are respectively the ionization potential and electron affinity. In terms of orbital energies, these equations can be written as:
equation M3
Thus, the electrophility index becomes
equation M4
This concept of electrophilicity index, its usefulness as well as the electrophilicity scales based on different physico-chemical properties have been reviewed by Chattaraj et al. [29]. In the section “Results”, we will briefly compare the electrophilicity scale, ω observed for our TP with the ranking based on their electron affinities EA.
2.3 Hydrophobic character
The calculation of logP was carried out using the XlogP program, which is based on the summation of atomic contributions and includes correction factors for some intramolecular interactions [30].
Statistical calculations and graph plotting were performed using the Kaleidagraph software (version 3.6). This software calculates the linear curve fits, using the least-squared error method.
3.1 Potencies of synthetic TP analogues of oleanolic acid as inducers of NQO1 versus their electron affinity (EA)
The molecular structures of the TP under study are shown in Fig.1
Figure 1
Figure 1
Structures of the triterpenoids under study
Induction potencies are assessed by using the CD value, the Concentration of each TP required to Double the NQO1 specific activity in Hepa 1c1c7 murine hepatoma cells. These CD values have been previously experimentally determined by us [4]. For convenience, the induction potencies are now reported as pCD values (Table 1), where pCD is the logarithm of the reciprocal of CD (expressed in M units). Notably, these vary by four orders of magnitude from 9.55 for the most active TP-225 to 5.41 for the least active TP-46.
Table 1
Table 1
Energy values in eV of E (LUMO) of TP, E (HOMO) of TP·-, E (SOMO) of TP·-, E (HOMO) of TP, CD values in μM unit, the corresponding pCD values and log P values of the TP
In Table 1 are also reported the energy values in eV of the E (LUMO) of neutral TP, the E (HOMO) of TP•- anions, the E (SOMO) of TP•- anions, the E (HOMO) of neutral TP, the CD values, and the pCD values of the TP. The E (LUMO) of TP, and E (HOMO) and E (SOMO) of the anion TP•- reported in Table 1, have a sign which is the inverse of those of the corresponding reduction potentials E (TP/TP•-). Thus, the most electrophilic TP is 225 and the least electrophilic is 46. Of note, TP-46, as well as TP-69 and TP-156 are the least potent inducers and have one Michael acceptor function (α,β-unsaturated carbonyl group), in contrast with the other 13 TP, which have two Michael acceptor functions. In addition, inactive oleanolic acid has no Michael acceptor function (Fig.1).
An absence of correlation of the biological activities of these TP with their electron-donating property is observed when plotting the pCD values versus the E (HOMO) of the TP, which represents the ease of a TP to donate an electron and is linearly correlated with their oxidation potential E (TP•+/TP). This absence of correlation is consistent with the fact that the TP act as electron-acceptors via their electrophilic Michael acceptor groups in their direct interaction with highly reactive thiols of the sensor protein Keap1 [4]. Plots of the pCD versus E (LUMO), versus E (HOMO) of the anion TP•- and versus the energy of the singly occupied molecular orbital E (SOMO) for the 16 active TP (Fig.1) give respectively linear correlations expressed by equations (1), (2) and (3) with correlation coefficients r:
equation M5
equation M6
equation M7
The plots of correlations (1), (2), and (3) are rather similar to one another with the same main outliers TP-233 and 191.
The linear plot pCD versus E (SOMO) corresponding to equation (3) has the best correlation coefficient r and is shown in Fig.2.
Figure 2
Figure 2
Plot of pCD (logarithm of the reciprocal of CD, expressed in M units) of the 16 triterpenoid inducers of NQO1 versus their E (SOMO), corresponding to equation (3)
Outlier TP-233 is dissimilar to the other compounds in the series by its large hydrophobic substituent -CO2(CH2)7CH3 in position 17 (Fig.1), which is the cause of its high logP value (8.65) relatively to the logP values of the other TP (Table 1). In contrast with TP-233, outlier TP-191 is dissimilar by the presence of two hydrophilic substituents –CO2H in positions 2 and 17, cause of the lowest logP value (1) in the TP series. These special structures can be considered as the cause of the presence of TP-233 and TP-191 relatively far from the linear fit of the correlations (1), (2), and (3) (Fig.2), and (5), (6), and (7) (Fig.3). If the TP-233 and TP-191 are taken out, the correlation coefficients r are significantly improved:
equation M8
equation M9
equation M10
Figure 3
Figure 3
Plot of pIC-50 (logarithm of the reciprocal of IC-50, expressed in M units) of the 16 triterpenoid inhibitors of induction of iNOS expressing the anti-inflammatory activities versus their E (SOMO), corresponding to equation (6)
The value of the correlation coefficients might still be improved by taking into account other physico-chemical properties such as the steric factors of the TP approach to their binding site, which are more difficult to quantify. The present correlations demonstrate an unquestionable general quantitative tendency: the higher is the electron affinity of a TP [i.e: the higher is its reduction potential E (TP/TP•-)], the greater its potency as NQO1 inducer. The importance of the electron affinity (EA) of the TP for their biological efficacy is consistent with the observations that the introduction of a second Michael acceptor functionality increases their potency and that a reaction of the TP occurs with the thiol groups of the Keap 1 sensor [4]. Reactions of Michael acceptor-containing TP and their tricyclic derivatives with nucleophile thiols have been reported previously [31-33].
3.2 Potencies of synthetic TP for suppression of iNOS induction versus their electron affinity (EA)
Table 2 reports the IC-50 and the pIC-50 values, cologarithm of IC-50, determined for the suppression of induction of nitric oxide synthase (iNOS) (M) in mouse macrophages, expressing the anti-inflammatory activities of our series of TP.
Table 2
Table 2
IC-50 values taken from references (14,15) and pIC-50 values for inhibition of iNOS activation by IFN-gamma by TP in primary mouse macrophages
The plot of pIC-50 versus the E (LUMO), versus E (HOMO) of the anion TP•- and versus the energy of the singly occupied molecular orbital E (SOMO) for the 16 active TP of Fig.1, give respectively linear correlations expressed by equations (4), (5) and (6) with correlation coefficients r :
equation M11
equation M12
equation M13
The plot pIC-50 versus E (SOMO) corresponding to equation (6) has the best correlation coefficient r, and is shown in Fig. 3.
Again, TP-233 and TP-191 are outliers for reasons related to their structure. If they are taken out of the plot, the correlation coefficients r are improved:
equation M14
equation M15
equation M16
New equations expressing the correlations for all compounds as a function of both parameters EA and logP are reported below. These new expressions lead to two remarks: (i) the EA parameters are by far the major factor controlling the biological properties, (ii) the new correlation coefficients are identical to those expressing pCD and pIC-50 as an exclusive function of the EA parameters (see equations 1, 2, 3, 4, 5, 6)
equation M17
equation M18
The two remarks above are also valid for the equations 1', 2', 3', 4', 5', 6', where the two outliers TP-233 and TP-191 have been taken out.
3.3 Potencies of synthetic TP analogues of oleanolic acid as inducers of NQO1 and as inhibitors of iNOS induction versus their electrophilicity index (ω)
As mentioned above, the electrophilicity index, ω, of every TP, reported in Table 3 is calculated in terms of E (HOMO) and E (LUMO) orbital energies of the neutral TP by means of the formula: equation M19
Table 3
Table 3
Electronegativity, χ, chemical hardness, η, and electrophilicity index, ω, of our triterpenoids
Plotting this elecrophilicity index, ω, versus E(SOMO) gives a linear correlation expressed by:
ω = 1.90 – 0.30 E (SOMO) with a correlation coefficient r = 0.92 This indicates that both scales, based on distinct physico-chemical parameters, are rather similar. Compared with the plots of pCD and pIC50 versus E(SOMO) expressed by equations (3), (3'), (6), and (6'), plots of pCD and pIC50 as a function of the electrophilicity index, ω, showed rather similar correlation coefficients expressed by the following equations: for all our 16 TP
equation M20
equation M21
and without the outliers TP 233 and 191
equation M22
equation M23
The efficacy ranking observed by using the electrophilicity index, ω, is in good agreement with that obtained by means of EA expressed by the E (SOMO) for our 16 TP.
Our approach could be used to analyze the efficacy ranking of other electrophilic ligands reacting with the highly reactive thiolates of Keap1 as well as those of other proteins.
More than 20 years ago, Talalay et al. provided the first insight into the structural and electronic requirements for phase 2 inducer activity by demonstrating that many inducers are Michael reaction acceptors [34]. Our calculations on triterpenoid electrophiles, acting as Michael acceptors via their α,β-unsaturated ketone functionalities, rationalize these observations. The present findings complement earlier studies on other classes of NQO1 inducers, considered as “anti-oxidants”, including diphenols [35], phenylpropenoids [36], and flavonoids [37] for which linear correlations between their NQO1 induction potencies and their E (HOMO) are observed. These linear correlations are due to the fact that these classes of electron-donating molecules (diphenols, phenylpropenoids, and flavonoids) act via a two-step mechanism involving (i) oxidation of these inducers to their quinone or semiquinone derivatives, which are Michael acceptors and (ii) oxidation or alkylation of highly reactive thiols of Keap1 by these electrophilic quinones or semiquinones. The transformation of these diphenols, phenylpropenoids, and flavonoids into oxidant species can occur by various routes, including their oxidation by reactive oxygen and nitrogen species (ROS, RNS), by autooxidation or by non-enzymatic (e.g., metal-catalyzed) as well as enzymatic mechanisms [35]. This necessary transformation of the diphenols, phenylpropenoids, and flavonoids into oxidant species is consistent with an important observation of Prochaska et al. [38, 39] that among diphenol inducers, only those with oxidizable phenolic hydroxyl groups (i.e., ortho- or para-, but not meta- diphenols) are inducers. Very recently, we have obtained experimental evidence that under aerobic conditions, Cu++ catalyzes the oxidation of ortho- and para-hydroquinones to their corresponding quinones, which subsequently react with cysteine residues of Keap1, resulting in activation of Nrf2, and that such oxidation to quinones is both necessary and sufficient for activation of the Keap1/Nrf2 pathway by oxidizable diphenols [40].
The correlations between the electron affinity (EA), and their electrophilicity index, ω, of the TP and their potencies as inducers of the cytoprotective phase 2 response and as inhibitors of inflammatory processes provide insights into the mechanism(s) by which the TP is recognized by the intracellular sensor molecules. In the case of the Keap1/Nrf2/ARE-dependent upregulation of the cytoprotective phase 2 response, it is now widely accepted that highly reactive cysteine residues of Keap1 are the sensors for inducers. Similarly, cysteine 179 of IκB kinase β (IKKβ) has been implicated as the intracellular target of TP which triggers events that ultimately lead to inhibition of the IFNγ-mediated upregulation of iNOS [41, 42]. Up to date, adducts of TP and cysteine residues of either Keap1 or IKKβ have not been identified despite of the strong spectroscopic (UV and NMR) evidence that the Michael addition reaction does occur when TP are incubated with thiols [4, 43], suggesting that such adducts do form, but are unstable, and reversible. Indeed, NMR spectroscopic evaluation has demonstrated the nucleophilic addition of dithiothreitol to TP-151 and the decomposition of the resulting adduct to regenerate the TP [43].
One may then speculate that the ultimate result of such a reversible addition could be a change in the oxidation state of the sensor cysteine(s), e.g, formation of disulfide bond(s) within Keap1 and IKKβ. Indeed, we have previously detected intermolecular disulfide-linked dimers of Keap1 in extracts of cells that were treated with inducers of three different chemical types, i.e., sulforaphane, 1,2-dithiole-3-thione, and bis(2-hydroxybenzylidene)acetone, and their facile conversion to reduced monomers by treatment with dithiothreitol [44]. Likewise, treatment with either TP-155 ot TP-151 resulted in formation of a species, corresponding to a dimer of IKKβ based on its molecular mass, that reacts with anti-IKKβ antibody in cells ectopically expressing wild-type-, but not C179A mutant IKKβ [42]. These findings imply that following chemical modification of specific cysteine thiols, the modifying agents are displaced by intermolecular sulfhydryl disulfide interchange to lead to a covalent disulfide-linked dimer of Keap1 or IKKβ as shown in Scheme 1.
Scheme 1
Scheme 1
Scheme 1
A nucleophilic attack on the electrophilic carbon (C1) on ring A of a Michael-acceptor bearing triterpenoid by a cysteine thiolate ion of a protein sensor leads initially to adduct formation that in turn is rapidly attacked by a second cysteine thiolate (more ...)
Notably, there is now experimental evidence for the formation of both intra- and intermolecular disulfide bridges in Keap1. Thus, very recently, Fourquet et al. [45] reported that exposure to reactive oxygen (ROS) and nitrogen (RNS) species leads to formation of an intramolecular disulfide linking Cys226 and Cys613, as well as an intermolecular disulfide bridging two Keap1 molecules through Cys151 of each monomer, and that even in untreated cells a fraction of Keap1 carries the intramolecular disulfide. In the case of a triterpenoid reacting with Keap1, such speculated outcome as shown in Scheme 1 will allow for the possibility of regeneration of the sensor proteins upon reduction of the disulfide bridge(s), unlike an irreversible inactivation by a stable adduct formation that may result in toxicity and will require de novo protein synthesis for regain of function.
In our present calculations, the TP outer energy levels computed by the molecular orbital method are correlated to (i) their electron affinity (EA) in the gas phase, and (ii) to their reduction potential E (TP/TP•-) in solution. Our quantum chemical study establishes quantitative structure-activity relationships between EA and biological efficacies. The electrophilicity index ω of the TP calculated using electronegativity and chemical hardness parameters, gives correlations and efficacy ranking analagous to those obtained by the efficacy scaling via EA, expressed by E (SOMO) calculations. These correlations demonstrate quantitatively, for this series of triterpenoids, that their electron-acceptor property is by far the major factor that determines their potency as inducers of the cytoprotective phase 2 response and as inhibitors of inflammatory processes. Our results illustrate the Szent-Györgyi concept of “submolecular electronic processes”, in particular those involving electron-donor or electron-acceptor properties of exogenous substances in biological systems [46]. More specifically, in the present study, the efficacy of the Michael addition of a thiol nucleophile to an electron deficient carbon of a triterpenoid has been interpreted using the electron affinity of the acceptor as a global reactivity descriptor. Moreover, the energy change predicted from orbital energies gives information on the most favored attack position of the electrophilic compound considered, as pointed out by A. Pullman [47].
The authors are grateful for financial support from the American Cancer Society (Grant RSG-07-157-01-CNE), the National Institutes of Health (Grants CA06793 and CA93780), Research Councils UK, Cancer Research UK (C20953/A10270), the Royal Society, the Anonymous Trust, Tenovus Scotland, the American Institute for Cancer Research, the Lewis B. and Dorothy Cullman Foundation, and W. Patrick McMullan and the McMullan Family Fund. The triterpenoids were kindly supplied by M.B. Sporn, T. Honda and G. Gribble, Dartmouth School of Medicine (Hanover, NH).
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.
All authors contributed equally to this work
Conflict of interest
The authors declare that there are no conflicts of interest.
1. Virchow R. Cellular Pathology as Based upon Physiological and Pathological Histology. J. B. Lippincott; Philadelphia: 1863. [PubMed]
2. Christen S, Hagen TM, Shigenaga MK, Ames BN. Chronic inflammation, mutation and cancer. In: Parsonnet J, editor. Microbes and Malignancy: Infection as a Cause of Human Cancer. Oxford University Press; New York: 1999. pp. 35–88.
3. Coussens LM, Werb Z. Inflammation and cancer. Nature. 2002;420:860–867. [PMC free article] [PubMed]
4. Dinkova-Kostova AT, Liby KT, Stephenson KK, Holtzclaw WD, Gao X, Suh N, Williams C, Risingsong R, Honda, T. T, Gribble GW, Sporn MB, Talalay P. Extremely potent triterpenoid inducers of the phase 2 response: correlations of protection against oxidant and inflammatory stress. Proc. Natl. Acad. Sci. USA. 2005;102:4584–4589. [PubMed]
5. Liu H, Dinkova-Kostova AT, Talalay P. Coordinate regulation of enzyme markers for inflammation and for protection against oxidants and electrophiles. Proc. Natl. Acad. Sci. USA. 2008;105:15926–15931. [PubMed]
6. Fahey JW, Dinkova-Kostova AT, Stephenson KK, Talalay P. The “Prochaska” microtiter plate bioassay for inducers of NQO1. In: Sies H, Packer L, editors. Quinones and quinone enzymes, part B. Methods in Enzymology. Vol. 382. Academic Press; San Diego: 2004. pp. 243–258. [PubMed]
7. Kang Y-H, Pezzuto JM. Induction of Quinone Reductase as a Primary Screen for Natural Product Anticarcinogens. In: Sies H, Packer L, editors. Quinones and quinone enzymes, part B. Methods in Enzymology. Vol. 382. Academic Press; San Diego: 2004. pp. 380–414. [PubMed]
8. Liby KT, Yore MM, Sporn MB. Triterpenoids and rexinoids as multifunctional agents for the prevention and treatment of cancer. Nat. Rev. Cancer. 2007;7:357–369. [PubMed]
9. Liby K, Royce DB, Williams CR, Risingsong R, Yore MM, Honda T, Gribble GW, Dmitrovsky E, Sporn TA, Sporn MB. The synthetic triterpenoids CDDO-methyl ester and CDDO-ethyl amide prevent lung cancer induced by vinyl carbamate in A/J mice. Cancer Res. 2007;67:2414–2419. [PubMed]
10. Liby K, Risingsong R, Royce DB, Williams CR, Yore MM, Honda T, Gribble GW, Lamph WW, Vannini N, Sogno I, Albini A, Sporn MB. Prevention and treatment of experimental estrogen receptor-negative mammary carcinogenesis by the synthetic triterpenoid CDDO-methyl Ester and therexinoid LG100268. Clin Cancer Res. 2008;14:4556–4563. [PubMed]
11. Ohshima H, Bartsch H. Chronic infections and inflammatory processes as cancer risk factors: possible role of nitric oxide in carcinogenesis. Mutat Res. 1994;305:253–264. [PubMed]
12. Maccoll A. Reduction potential of conjugated systems. Nature. 1949;163:178–179.
13. Matsen FA. Electron affinities, methyl affinities and ionization energies of condensed ring aromatic hydrocarbons. J. Chem. Phys. 1956;24:602–606.
14. Honda T, Rounds BV, Bore L, Finlay HJ, Favaloro FG, Jr, Suh N, Wang Y, Sporn MB, Gribble GW. Synthetic oleanane and ursane triterpenoids with modified rings A and C: a series of highly active inhibitors of nitric oxide production in mouse macrophages. J. Med. Chem. 2000;43:4233–4246. [PubMed]
15. Honda T, Honda Y, Favaloro FG, Jr, Gribble GW, Suh N, Place AE, Rendi MH, Sporn MB. A novel dicyanotriterpenoid, 2-cyano-3,12-dioxooleana-1,9(11)-dien-28-onitrile, active at picomolar concentrations for inhibition of nitric oxide production. Bioorg. Med. Chem. Lett. 2002;12:1027–1030. [PubMed]
16. Koopmans T. Ordering of Wave Functions and Eigenvalues to the Individual Electrons of an Atom. Physica. 1933;1:104–113.
17. Daudel R, Lefebvre R, Moser C. Quantum Chemistry, Methods and applications. Interscience Publishers; New York & London: 1959.
18. Pullman B, Pullman A. Quantum Biochemistry. Interscience Publishers; New York, USA: 1963. pp. 132–135.
19. Moffit W. Term values in hybrid states (see Eqs A1, A2) Proc.Roy.Soc A. 1950;202:534–548.
20. Angeli C. Physical interpretation of Koopmans’ theorem: a criticism of the current didactic presentation. J. Chem. Educ. 1998;75:1494–1497.
21. Mulliken RS. Molecular Orbital Method and Molecular Ionization Potentials. Phys. Rev. 1948;74:736–738.
22. Dodds JL, McWeeny R. Orbital energies and Koopmans’ theorem in opn-shell Hartree-Fock theory. Chem. Phys. Letters. 1972;13:9–12.
23. Heinrich N, Koch W, Frenking G. On the use of Koopmans’ theorem to estimate negative electron affinities. Chem. Phys. Letters. 1986;124:20–25.
24. Maynard AT, Huang M, Rice WG, Covell DG. Reactivity of the HIV-1 nucleocapsid protein p7 zinc finger domains from the perspective of density-functional theory. Proc. Natl. Acad. Sci. U.S.A. 1998;95:11578–11583. [PubMed]
25. Parr RG, Szentpàly L. v., Liu S. Electrophilicity index. J. Am. Chem. Soc. 1999;121:1922–1924.
26. Chong DP, Gritsenko OV, Baerends EJ. Interpretation of the Kohn-Sham orbital energies as approximate vertical ionization potentials. J. Chem. Physics. 2002;116:1760–1772.
27. Filipetti A. Electron affinity in density functional theory in the local spin-density approximation. Physical Rev. A. 1998;57:914–919.
28. Bensasson RV, Berthier G, Zoete V. “Comment on “Assessing the Efficacy of Nonsteroidal Anti-Inflammatory Drugs Through the Quantum Computation of Molecular Ionization Energies” Text and Supporting Information available as an ASAP article on the Journal of Physical Chemistry A website. in press. Publication Date (Web): October 9, 2009, DOI: 10.1021/jp9074652. [PubMed]
29. Chattaraj PK, Sarkar U, Roy DR. Electrophilicity Index. Chem. Rev. 2006;106:2065–2091. [PubMed]
30. Wang R, Fu Y, Lai L. A New Atom-Additive Method for Calculating Partition Coefficients. J. Chem. Inf. Comput. Sci. 1997;37:615–621.
31. Brookes PS, Morse K, Ray D, Tompkins A, Young SM, Hilchey S, Salim S, Konopleva M, Andreeff M, Phipps R, Bernstein SH. The triterpenoid 2-cyano-3,12-dioxooleana-1,9-dien-28-oic acid and its derivatives elicit human lymphoid cell apoptosis through a novel pathway involving the unregulated mitochondrial permeability transition pore. Cancer Res. 2007;67:1793–1802. [PubMed]
32. Ahmad R, Raina D, Meyer C, Kufe D. Triterpenoid CDDO-methyl ester inhibits the Janus-activated kinase-1 (JAK1)-signal transducer and activator of transcription-3 (STAT3) pathway by direct inhibition of JAK1 and STAT3. Cancer Res. 2008;68:2920–2926. [PMC free article] [PubMed]
33. Liby K, Yore MM, Roebuck BD, Baumgartner KJ, Honda T, Sundararajan C, Yoshizawa H, Gribble GW, Williams CR, Risingsong R, Royce DB, Dinkova-Kostova AT, Stephenson KK, Egner PA, Yates MS, Groopman JD, Kensler TW, Sporn MB. A novel acetylenic tricyclic bis-(cyano enone) potently induces phase 2 cytoprotective pathways and blocks liver carcinogenesis induced by aflatoxin. Cancer Res. 2008;68:6727–6733. [PMC free article] [PubMed]
34. Talalay P, De Long MJ, Prochaska HJ. Identification of a common chemical signal regulating the induction of enzymes that protect against chemical carcinogenesis. Proc Natl Acad Sci USA. 1988;85:8261–8265. [PubMed]
35. Bensasson RV, Zoete V, Dinkova-Kostova AT, Talalay P. Two-Step Mechanism of Induction of the Gene Expression of a Prototypic Cancer-Protective Enzyme by Diphenols. Chem. Res. Toxicol. 2008;21:805–812. [PubMed]
36. Zoete V, Rougée M, Dinkova-Kostova AT, Talalay P, Bensasson RV. Redox ranking of inducers of a cancer protective enzyme via the energy of their highest occupied molecular orbital. Free Radic. Biol. Med. 2004;36:1418–1423. [PubMed]
37. Lee-Hilz YY, Boerboom AM, Westphal AH, Berkel WJ, Aarts JM, Rietjens IM. Pro-oxidant activity of flavonoids induces EpRE-mediated gene expression. Chem. Res. Toxicol. 2006;19:1499–1505. [PubMed]
38. Prochaska HJ, De Long MJ, Talalay P. On the mechanisms of induction of cancer-protective enzymes: a unifying proposal. Proc Natl Acad Sci USA. 1985;82:8232–8236. [PubMed]
39. Prochaska HJ, Talalay P. Huang MT, Ho CT, Lee CY, editors. Phenolic antioxidants as inducers of anticarcinogenic enzymes. Phenolic compounds in food in their effects on health II, ACS Symposium Series No.507. 1992. pp. 150–159.
40. Wang XJ, Hayes JD, Higgins LG, Wolf CR, Dinkova-Kostova AT. Activation of the Nrf2 signalling pathway by copper-mediated redox cycling of para- and ortho–hydroquinones. Chem Biol. 2010;17:75–85. [PubMed]
41. Yore MM, Liby KT, Honda T, Gribble GW, Sporn MB MB. The synthetic triterpenoid 1-[2-cyano-3,12-dioxooleana-1,9(11)-dien-28-oyl]imidazole blocks nuclear factor-κB activation through direct inhibition of IκB kinase β Mol Cancer Ther. 2006;5:3232–3239. [PubMed]
42. Ahmad R, Raina D, Meyer C, Kharbanda S, Kufe D. Triterpenoid CDDO-Me blocks the NF-κB pathway by direct inhibition of IKKκ on Cys-179. J Biol Chem. 2006;281:35764–35769. [PubMed]
43. Couch RD, Browning RG, Honda T, Gribble GW, Wright DL, Sporn MB, Anderson AC. Studies on the reactivity of CDDO, a promising new chemopreventive and chemotherapeutic agent: implications for a molecular mechanism of action. Bioorg Med Chem Lett. 2005;15:2215–2219. [PubMed]
44. Wakabayashi N, Dinkova-Kostova AT, Holtzclaw WD, Kang MI, Kobayashi A, Yamamoto M, Kensler TW, Talalay P. Protection against electrophile and oxidant stress by induction of the phase 2 response: fate of cysteines of the Keap1 sensor modified by inducers. Proc Natl Acad Sci USA. 2004;101:2040–2045. [PubMed]
45. Fourquet S, Guerois R, Biard D, Toledano MB. Activation of NRF2 by nitrosative agents and H2O2 involves KEAP1 disulfide formation. J Biol Chem. 2010;285:8463–8471. [PubMed]
46. Szent-Györgyi A. Submolecular Biology. Radiation Research. 1960;(Suppl.2):4–18.
47. Pullman A. On the possibilities of correlating the redox potential of reversible systems to the molecular orbitals involved in the elecron transfer. J. Theoret. Biol. 1962;2:259–262.