Subscriber access provided by UNIV OF NEBRASKA - LINCOLN
Letter
Active Components of Essential Oils as Anti-Obesity Potential Drugs investigated by in silico techniques Giosuè Costa, Maria Concetta Gidaro, Daniela Vullo, Claudiu T Supuran, and Stefano Alcaro J. Agric. Food Chem., Just Accepted Manuscript • DOI: 10.1021/acs.jafc.6b02004 • Publication Date (Web): 08 Jun 2016 Downloaded from http://pubs.acs.org on June 9, 2016
Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.
Journal of Agricultural and Food Chemistry is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.
Page 1 of 24
Journal of Agricultural and Food Chemistry
Active Components of Essential Oils as Anti-Obesity Potential Drugs investigated by in silico techniques Giosuè Costa†,, Maria Concetta Gidaro†*, Daniela Vullo‡, Claudiu T. Supuran‡, Stefano Alcaro† †
Dipartimento di Scienze della Salute, Università “Magna Græcia” di Catanzaro, Campus
Universitario “S. Venuta”, Viale Europa, Loc. Germaneto, 88100 Catanzaro, Italy. ‡
Laboratorio di Chimica Bioinorganica, Polo Scientifico, Università degli Studi di Firenze, Via
della Lastruccia 3, 50019 Sesto Fiorentino, Florence, Italy.
Corresponding Author *
Dr Maria Concetta Gidaro, Dipartimento di Scienze della Salute, Università “Magna Græcia” di
Catanzaro, Campus “S. Venuta”, Viale Europa, 88100 Catanzaro, Italy. E-mail:
[email protected] Phone: +39 0961 3694197
ACS Paragon Plus Environment
1
Journal of Agricultural and Food Chemistry
Page 2 of 24
1
ABSTRACT
2
In this study, for the first time, we have considered Essential Oils (EOs) as possible resources of
3
Carbonic Anhydrase inhibitors (CAIs), in particular against the mitochondrial isoform VA that,
4
actually, represents an innovative target for the obesity treatment. In silico structure-based virtual
5
screening was performed in order to speed up the identification of promising anti-obesity agents.
6
The potential hit compounds were submitted to in vitro assays and experimental results,
7
corroborated by molecular modeling studies, showed EOs components as a new class of CAIs with
8
a competitive mechanism of action due to the Zinc ion coordination within the active sites of these
9
metallo-enzymes.
10
Keywords
11
Essential Oils, Anti-obesity, Human Carbonic Anhydrase, Inhibition, Virtual Screening
ACS Paragon Plus Environment
2
Page 3 of 24
Journal of Agricultural and Food Chemistry
12
Introduction
13
Obesity is ranked by the World Health Organization (WHO) as one of the top 10 global health
14
problems, usually correlated with other pathologies such as cardiovascular and musculo-skeletal
15
disorders, type 2 diabetes, depression and several cancers including colon, breast and
16
endometrium.1 The pharmacotherapies for the obesity treatment play an important role for
17
increasing the weight loss established by diet and exercise alone, but actually, only a few anti-
18
obesity drugs are approved by the European Medicines Agency (EMA) and the US Food and Drug
19
Administration (FDA).2 In particular, five medications have been approved in the USA for chronic
20
weight
21
naltrexone/bupropion and liraglutide. Lorcaserin and phentermine/extended-release topiramate have
22
not been approved in the European Union, in fact, several controversies are still being debated and
23
one of the major uncertainties, for patients with obesity, is whether current pharmacological
24
treatments increase, reduce, or are neutral with respect to cardiovascular events.3 Moreover, drugs
25
acting on the metabolism have considerable side effects and, for these reasons, pharmacological
26
interventions for weight loss remain limited.4,5
27
Several studies have provided evidence that Carbonic Anhydrase inhibitors (CAIs) are emerging as
28
anti-obesity drugs.6-8 To date 16 α-CA isoforms were reported in mammals with different tissue
29
distribution and subcellular localization.9 These metallo-enzymes catalyze the reversible hydration
30
of carbon dioxide to afford the ionic species bicarbonate and protons as follows: CO2 + H2O ⇆
31
HCO3− + H+. The CA isoforms VA and VB, involved in lipogenesis, are expressed only in the
32
mitochondria and the clinical observation of body weight reduction, for patients treated with the
33
anticonvulsant topiramate and zonisamide (in phase II), is rationalized considering the potent
34
inhibitory effect of these mitochondrial enzymes.10,11
management:
orlistat,
lorcaserin,
phentermine/extended-release
ACS Paragon Plus Environment
topiramate,
3
Journal of Agricultural and Food Chemistry
Page 4 of 24
35
In this study, for the first time, Essential Oils (EOs) components were considered as a new class of
36
potential CAIs. Traditionally, EOs have been largely employed as antibacterial, antifungal and
37
insecticidal agents12,13 and several studies have confirmed these biological activities.14-19 In
38
addition, EOs have been also considered for their analgesic,20 sedative,21 anti-inflammatory,22-24
39
spasmolytic25 and locally anesthetic properties.26 Finally, EOs seem to have a great potential as
40
anticancer therapeutic agents but mechanisms of action are still unknown, considering that EO
41
extract is a very complex mixture of low weight molecules (20 to 70 components) at different
42
concentrations in the phytocomplex.27-31
43
Conversely, our study was focused on a specific interaction between the EOs content and the CAs
44
family. Considering EOs as new resources of CAIs, in particular of the mitochondrial isoforms VA,
45
a Structure-Based Virtual Screening (SBVS) was performed in order to early identify in silico new
46
potential anti-obesity drugs. Therefore, the Essential Oil University (EOU) chemoinformatics
47
database32 was used to build a chemical library of natural ligands that were virtually screened
48
against several crystallographic X-ray structures of CA, deposited into the bioinformatics Protein
49
Data Bank (PDB).33 The most relevant compounds selected on the basis of the predicted binding
50
affinity results, due to the ligand-target interactions, were submitted to in vitro assays and new
51
chemical scaffolds were discovered for CAIs.
52
Materials and Methods
53
In silico Structure-Based Virtual Screening
54
The 3D chemical structures of 2690 compounds, content in the EOU website,32 were downloaded
55
and imported into the graphic interface of Maestro 9.734 in order to create the ligands database used
56
for the SBVS. Duplicates were detected by the Canvas35 utility and Lipinski’s rule of five36,37 was
57
used to filter 1771 compounds. For each molecule protonated and tautomeric forms at pH 7.4 were
ACS Paragon Plus Environment
4
Page 5 of 24
Journal of Agricultural and Food Chemistry
58
calculated, by the “LigPrep”38 module, ending up with a library of 2304 compounds. The lowest
59
energy conformations, obtained by using the OPLS-2005 force field, represented the starting points
60
for docking simulations. The X-ray crystallographic structures of CA I, CA II and CA VA isoforms
61
were selected from the PDB33 with the respective entry codes 1AZM,39 4CQ040 and 1DMY,41 while
62
for the CA VB enzyme no crystallographic models are deposited in the PDB. After the preparation
63
procedure of each X-ray structure, the docking simulations were performed as reported in the
64
experimental section of our previous SBVS study.42 Regarding the CA VA isoform, the best docked
65
poses of compounds 1, 2, 6 and 11 were further submitted to a MD simulation of 100 ns, carried out
66
by Desmond 3.7.43 Therefore, in order to analyze ligand-target interactions involved in the
67
molecular recognition, we used the MD protocol already applied for investigating flavonoids
68
inhibition of the monoamine oxidases isoforms.44 Pymol version 1.745 was used for create Figure 1.
69
In vitro Carbonic Anhidrase Inibition assay
70
All CA isoforms (I, II and VA) used for in vitro assays were purified and obtained as previously
71
reported in other published works.46-48 An Applied Photophysics stopped-flow instrument has been
72
used for assaying the CA-catalyzed CO2 hydration activity.49 For compounds 1-11 four inhibitor
73
concentrations were used, starting from 10-5 M to 10-8 M, and testing them in triplicate. IC50 values
74
were obtained from dose-response curves working at seven different concentrations of test
75
compound, by fitting the curves using PRISM (www.graphpad.com) and nonlinear least squares
76
methods; values represent the mean of at least three different determinations as described by us
77
previously.50 The inhibition constants (Ki) were then derived by using the Cheng-Prusoff equation51
78
as follows: Ki = IC50/(1 + [S]/Km), where [S] represents the CO2 concentration at which the
79
measurement was carried out, and Km is the concentration of substrate at which the enzyme activity
80
is at half-maximal. The concentrations of enzymes used in the assay were as follows: hCA I: 11.7
ACS Paragon Plus Environment
5
Journal of Agricultural and Food Chemistry
Page 6 of 24
81
nM, hCA II: 6.9 nM, hCA VA 10.5 nM. The kcat/KM values for these enzymes were as follows: 5.0
82
x 107 M-1 x s-1 (hCA I); 1.5 x 108 M-1 x s-1 (hCA II); 2.9 x 107 M-1 x s-1 (hCA VA), respectively.
83
Results and discussion
84
Theoretical studies and in silico methods, by means of high performing computers, recently
85
represent new innovative approaches in the drug discovery process. In fact, excluding the clinical
86
trials, chemoinformatics and bioinformatics are playing important roles in every stage of the drug
87
discovery pipeline.52 The EOU website is an example of a chemoinformatics database, providing
88
accurate information about the world of EOs.
89
In this study, we have considered for the first time EOs as resources of anti-obesity agents by
90
inhibition of the mitochondrial CA VA/VB isoforms. With the aim of performing a SBVS, the 3D
91
compounds used to build the ligands database were downloaded from the EOU website and
92
prepared as reported in the experimental section. Regarding the X-ray crystal structures of the
93
enzymes, we selected the PDB models of ubiquitous CA I, CA II and mitochondrial CA VA as
94
previously reported for another SBVS study, performed in our laboratory for identifying new
95
natural CAIs.42 Unfortunately, in the PDB there are no crystal structures for the CA VB isoform
96
and, therefore, it was not considered for the SBVS. Molecular docking simulations and the virtual
97
screening procedure, already published, allowed us to select the new hit compounds on the basis of
98
the theoretical affinity values,42 reported in the Table 1S.
99
Ten molecules, identified in silico as promising CA VA inhibitors, were submitted to in vitro
100
assays using the human enzymes in order to confirm the ability to exert the predicted biological
101
activity. In the Table 1 for compounds 1-10 and the positive control acetazolamide, also called
102
AZA (11), are reported inhibition constant (Ki) values, calculated as displayed in the experimental
103
section.
ACS Paragon Plus Environment
6
Page 7 of 24
104
Journal of Agricultural and Food Chemistry
Table 1 should be positioned here
105 106
Biological results showed all compounds as week inhibitors of the slow cytosolic isozyme CA I,
107
demonstrated by Ki values in the range of 3.44-8.81 µM, and ineffective inhibitors of the rapid
108
cytosolic isozyme CA II with Ki >100 µM, except for compounds 1, 3, 7 and 8. As regard the
109
mitochondrial CA VA, it was poorly inhibited by compounds 8 and 10, with Ki values of 43.10 µM
110
and 44.00 µM, while the other compounds showed a Ki in the range of 3.80-9.69 µM. The 2-
111
hydroxyisobutyric acid (2) proved to be the best CA VA inhibitor, with a Ki value equal to 3.80
112
µM, followed by 3Z-nonenoic acid (6) and β-thujaplicin (1) with Ki values 4.52 µM and 7.50 µM,
113
respectively. Compound 9 has been already tested against several CA isoforms.53
114
In order to analyze the ligand-target interactions, involved in the molecular recognition, and with
115
the aim of taking into account the mechanism of action of new scaffolds, the best docked poses of
116
compounds 1, 2 and 6, in complex with the CA VA X-ray structure 1DMY,41 were further
117
investigated by a Molecular Dynamics (MD) simulation of 100 ns. Moreover, the reference 11 co-
118
crystallized with the CA VA PDB model, was also used as positive control for the in silico
119
simulations. In Figure 1, it is displayed a competitive mechanism of action for volatile compounds
120
of EOs, coordinating the Zinc ion and establishing other non-covalent ligand-target interactions
121
with amino-acid residues of the CA VA catalytic site.
122 123
Figure 1 should be positioned here
124 125
In particular, considering the different types of CAIs binding modes,54 all selected compounds
126
showed metal-coordination binding to Zinc ion (Figures 1 and 1S). Docking poses were submitted
ACS Paragon Plus Environment
7
Journal of Agricultural and Food Chemistry
Page 8 of 24
127
to molecular dynamics simulations (MDs) in order to rationalize the physical basis of the
128
complexes, as reported in the plots of Figure 2S. For compound 2, the carboxyl group in the
129
dissociated form and the hydroxyl substituent at C2 are key moieties in the metal coordination that
130
was maintained for all MD simulation time (Figure 2).
131 132
Figure 2 should be positioned here
133 134
Moreover, stabilization of the best pose is attributed to H-bond contacts with THR199 and
135
ARG209, hydrophobic interactions with TRP209 and HIS94/96/119, also involved in ionic and π-π
136
contacts (Figures 1, 1S and 2S). The theoretical results showed a stabilization of compound 2,
137
better than 1 and 6, but not than AZA (11), due to the great number of H-bond contacts. These
138
findings are summarized in Figure 3, as timeline representation of H-bond contacts established
139
between enzyme and ligand over the course of the MD trajectory.
140 141
Figure 3 should be positioned here
142 143
Enzyme activity analysis is only a preliminary investigation but it is the first step, in the drug
144
discovery process, for identifying new potential anti-obesity agents targeting the CA V isoforms.
145
Actually, there are several patents with mitochondrial CAIs, also in combination therapy,
146
confirming that inhibition of enzymes involved in de novo lipogenesis represents a new valid
147
approach for the obesity treatment.1In addition, from a polypharmacological prospective, EOs
148
contents and their derivatives may be considered Multi-Target Ligands (MuTaLig), in fact, other
149
experimental results reported in literature have shown EOs compounds as potential Histone
ACS Paragon Plus Environment
8
Page 9 of 24
Journal of Agricultural and Food Chemistry
150
Deacetylase inhibitors (HDACIs), promising targets in the cancer therapy, by binding the Zinc ion
151
also present in this metallo-enzyme.55,56
152
In conclusion, in this work, docking and MDs enabled us to investigate the binding modes of
153
ligands against the CA VA target, confirming the important role of molecular modeling techniques
154
to speed up the identification of bioactive compounds, reducing time and costs of research
155
procedures and taking into account multi target interactions. Moreover, EOs and other natural
156
products could be investigated in silico in order to early identify new potent anti-obesity drugs
157
and/or new scaffolds to focus the lead optimization studies, considering that it has been
158
demonstrated that CA VA inhibition with sulfonamides changes the metabolism of pyruvate,
159
acetate, and succinate in mitochondria and reduces dramatically lipogenesis.57
160
ACS Paragon Plus Environment
9
Journal of Agricultural and Food Chemistry
Page 10 of 24
161
ABBREVIATIONS
162
EOs, Essential Oils; CAIs, Carbonic Anhydrase inhibitors; WHO, World Health Organization;
163
EMA, European Medicines Agency; FDA, Food and Drug Administration; SBVS, Structure-Based
164
Virtual Screening; EOU, Essential Oil University; PDB, Protein Data Bank; AZA, acetazolamide;
165
Ki, inhibition constant; MD, Molecular Dynamics; MDs, Molecular Dynamics simulations; hCA,
166
human Carbonic Anhydrase; MuTaLig, Multi-Target Ligands; HDACIs, Histone Deacetylase
167
inhibitors.
168
ACKNOWLEDGMENTS
169
This work was supported by Interregional Research Center for Food Safety and Health at the
170
Magna Græcia University of Catanzaro (MIUR PON a3_00359). The authors thank Teresa Piperno,
171
biotech student at Magna Græcia University of Catanzaro, for her availability and support to the
172
experimental section and Dr. Simon Cross, Molecular Discovery Ltd, for the corrections in the
173
language style of manuscript. The authors also thank the COST action CA15135 (Multi-target
174
paradigm for innovative ligand identification in the drug discovery process MuTaLig) for the
175
support.
176
ASSOCIATED CONTENT
177
Supporting Information. In the Table 1S are reported the theoretical affinity results of molecular
178
docking for compounds 1-11. 2D structures of 1, 2, 6 and 11 with ligand-target interactions of the
179
best docked poses are shown in Figure 1S. These contacts were also monitored during MDs and
180
reported in Figure 2S. This material is available free of charge via the Internet at
181
http://pubs.acs.org.”
182
Notes
ACS Paragon Plus Environment
10
Page 11 of 24
Journal of Agricultural and Food Chemistry
183
This project won the first national prize of the Italian Society for the Research of Essential Oils
184
(SIROE), presented from the Dr. Maria Concetta Gidaro during the 3rd Congress to Rome on
185
November 2015.
186
ACS Paragon Plus Environment
11
Journal of Agricultural and Food Chemistry
187
REFERENCES
188
1.
189 190
2.
Haslam, D. Weight management in obesity - past and present. Int. J. Clin. Pract. 2016, 70, 206-217.
3.
193 194
Scozzafava, A.; Supuran, C. T.; Carta, F. Antiobesity carbonic anhydrase inhibitors: a literature and patent review. Expert Opin. Ther. Pat. 2013, 23, 725-735.
191 192
Page 12 of 24
Bray, G.A.; Frühbeck, G.; Ryan, D.H.; Wilding, JP. Management of obesity. Lancet. 2016, Epub ahead of print. doi: 10.1016/S0140-6736(16)00271-3.
4.
195
Rodríguez, J.E.; Campbell, K.M. Past, Present, and Future of Pharmacologic Therapy in Obesity. Prim Care. 2016, 43, 61-67.
196
5.
Rankin, W.; Wittert, G. Anti-obesity drugs. Curr. Opin. Lipidol. 2015, 26, 536-543.
197
6.
Supuran, C. T. Carbonic anhydrase inhibitors as emerging drugs for the treatment of obesity.
198 199
Expert Opin. Emerg. Drugs. 2012, 17, 5-11. 7.
200 201
for the treatment of obesity. Expert Opin. Emerg. Drugs 2008, 13, 383-392. 8.
202 203
9.
208
Supuran, C. T. Drug interaction considerations in the therapeutic use of carbonic anhydrase inhibitors. Expert Opin. Drug. Metab. Toxicol. 2016, 3, 1-9.
10.
206 207
Supuran, C. T. Carbonic anhydrases: novel therapeutic applications for inhibitors and activators. Nat. Rev. Drug Discov. 2008, 7, 168-181.
204 205
Supuran, C. T.; Di Fiore, A; De Simone, G. Carbonic anhydrase inhibitors as emerging drugs
Antel, J.; Hebebrand, J. Weight-reducing side effects of the antiepileptic agents topiramate and zonisamide. Handb Exp Pharmacol. 2012, 209, 433-66.
11.
Vitale, R.M.; Pedone, C.; Amodeo, P.; Antel, J.; Wurl, M.; Scozzafava, A.; Supuran, C. T.; De Simone, G. Molecular modeling study for the binding of zonisamide and topiramate to the
ACS Paragon Plus Environment
12
Page 13 of 24
Journal of Agricultural and Food Chemistry
209
human mitochondrial carbonic anhydrase isoform VA. Bioorg. Med. Chem. 2007, 15, 4152-
210
4158.
211
12.
Kiran, S.; Prakash, B. Assessment of Toxicity, Antifeedant Activity, and Biochemical
212
Responses in Stored-Grain Insects Exposed to Lethal and Sublethal Doses of Gaultheria
213
procumbens L. Essential Oil. J. Agric. Food Chem. 2015, 63, 10518-10524.
214
13.
Yeom, H. J.; Jung, C. S.; Kang, J.; Kim, J.; Lee, J. H.; Kim, D. S.; Kim, H. S.; Park, P. S.;
215
Kang, K. S.; Park, I. K. Insecticidal and acetylcholine esterase inhibition activity of
216
Asteraceae plant essential oils and their constituents against adults of the German cockroach
217
(Blattella germanica). J. Agric. Food Chem. 2015, 63, 2241-2248
218
14.
219 220
Bakkali, F.; Averbeck, S.; Averbeck, D.; Idaomar, M. Biological effects of essential oils--a review. Food Chem. Toxicol. 2008, 46, 446-475.
15.
Di Vito, M.; Mattarelli, P.; Modesto, M.; Girolamo, A.; Ballardini, M.; Tamburro, A.;
221
Meledandri, M.; Mondello, F. In Vitro Activity of Tea Tree Oil Vaginal Suppositories against
222
Candida spp. and Probiotic Vaginal Microbiota. Phytother. Res. 2015, 29, 1628-1633.
223
16.
Rosato, A.; Piarulli, M.; Corbo, F.; Muraglia, M.; Carone, A.; Vitali, M. E.; Vitali, C. In vitro
224
synergistic antibacterial action of certain combinations of gentamicin and essential oils. Curr.
225
Med. Chem. 2010, 17, 3289-3295.
226
17.
Diao, W. R.; Hu, Q. P.; Feng, S. S.; Li, W. Q.; Xu, J. G. Chemical composition and
227
antibacterial activity of the essential oil from green huajiao (Zanthoxylum schinifolium)
228
against selected foodborne pathogens. J. Agric. Food Chem. 2013, 61, 6044-6049.
229
18.
Mimica-Dukic, N.; Bozin, B.; Sokovic, M.; Simin, N. Antimicrobial and antioxidant activities
230
of Melissa officinalis L. (Lamiaceae) essential oil. J. Agric. Food Chem. 2004, 52, 2485-
231
2489.
ACS Paragon Plus Environment
13
Journal of Agricultural and Food Chemistry
232
19.
Page 14 of 24
Seo, S. M.; Jung, C. S.; Kang, J.; Lee, H. R.; Kim, S. W.; Hyun, J.; Park, I. K. Larvicidal and
233
acetylcholinesterase inhibitory activities of apiaceae plant essential oils and their constituents
234
against aedes albopictus and formulation development. J. Agric. Food Chem. 2015, 63, 9977-
235
9986.
236
20.
Friedman, M. Chemistry and multibeneficial bioactivities of carvacrol (4-isopropyl-2-
237
methylphenol), a component of essential oils produced by aromatic plants and spices. J.
238
Agric. Food Chem., 2014, 62, 7652-7670.
239
21.
de Sousa, D. P.; de Almeida Soares Hocayen, P.; Andrade, L.N.; Andreatini, R. A Systematic
240
Review of the Anxiolytic-Like Effects of Essential Oils in Animal Models. Molecules 2015,
241
20, 18620-18660.
242
22.
Rufino, A. T.; Ribeiro, M.; Judas, F.; Salgueiro, L.; Lopes M. C.; Cavaleiro, C.; Mendes, A.
243
F. Anti-inflammatory and chondroprotective activity of (+)-α-pinene: structural and
244
enantiomeric selectivity. J. Nat. Prod. 2014, 77, 264-269.
245
23.
Altinier, G.; Sosa, S.; Aquino, R. P.; Mencherini, T.; Della Loggia, R.; Tubaro, A.
246
Characterization of topical antiinflammatory compounds in Rosmarinus officinalis L. J.
247
Agric. Food Chem., 2007, 55, 1718-1723.
248
24.
Chao, L. K.; Hua, K. F.; Hsu, H. Y.; Cheng, S. S.; Liu, J. Y.; Chang, S. T. Study on the
249
antiinflammatory activity of essential oil from leaves of Cinnamomum osmophloeum. J.
250
Agric. Food Chem. 2005, 53, 7274-7278.
251
25.
Rivero-Cruz, B.; Rivero-Cruz, I.; Rodríguez, J. M.; Cerda-García-Rojas, C. M.; Mata, R.
252
Qualitative and quantitative analysis of the active components of the essential oil from
253
Brickellia veronicaefolia by nuclear magnetic resonance spectroscopy. J. Nat. Prod. 2006, 69,
254
1172-6.
ACS Paragon Plus Environment
14
Page 15 of 24
255
Journal of Agricultural and Food Chemistry
26.
Toni, C.; Martos-Sitcha, J. A.; Baldisserotto, B.; Heinzmann, B. M.; de Lima Silva, L.;
256
Martínez-Rodríguez, G.; Mancera, J. M. Sedative effect of 2-phenoxyethanol and essential oil
257
of Lippia alba on stress response in gilthead sea bream (Sparus aurata). Res. Vet. Sci. 2015,
258
103, 20-27.
259
27.
Russo, R.; Corasaniti, M. T.; Bagetta, G.; Morrone, L. A. Exploitation of cytotoxicity of some
260
essential oils for translation in cancer therapy. Evid. Based Complement. Alternat. Med. 2015,
261
2015, 1-9.
262
28.
Guerrini, A.; Lampronti, I.; Bianchi, N.; Zuccato, C.; Breveglieri, G.; Salvatori, F.; Mancini,
263
I.; Rossi, D.; Potenza, R.; Chiavilli, F.; Sacchetti, G.; Gambari, R.; Borgatti, M. Bergamot
264
(Citrus bergamia Risso) fruit extracts as γ-globin gene expression inducers: phytochemical
265
and functional perspectives. J. Agric. Food Chem. 2009, 57, 4103-4111.
266
29.
Sun, H.; Ni, H.; Yang, Y.; Wu, L.; Cai, H. N.; Xiao, A. F.; Chen, F. Investigation of sunlight-
267
induced deterioration of aroma of pummelo (Citrus maxima) essential oil. J. Agric. Food
268
Chem. 2014, 62, 11818-11830.
269
30.
Siani, A. C.; Tappin, M. R.; Ramos, M. F.; Mazzei, J. L.; Ramos, M. C.; De Aquino Neto, F.
270
R.; Frighetto, N. Linalool from Lippia alba: study of the reproducibility of the essential oil
271
profile and the enantiomeric purity. J. Agric. Food Chem. 2002, 50, 3518-3521.
272
31.
273
Schaneberg, B. T.; Khan, I. A. Comparison of extraction methods for marker compounds in the essential oil of lemon grass by GC. J. Agric. Food Chem. 2002, 50, 1345-1349.
274
32.
http://essentialoil.university/
275
33.
The Research Collaboratory for Structural Bioinformatics (RCSB) Protein Data Bank (PDB),
276 277
available online at http://www.rcsb.org. 34.
Maestro, version 9.7, Schrödinger, LLC, New York, NY, 2014.
ACS Paragon Plus Environment
15
Journal of Agricultural and Food Chemistry
Page 16 of 24
278
35.
Canvas, version 2.2, Schrödinger, LLC, New York, NY, 2014.
279
36.
Veber, D.; Johnson, S.; Cheng, H.; Smith, B.; Ward, K.; Kopple, K. Molecular properties that
280 281
influence the oral bioavailability of drug candidates. J. Med. Chem. 2002, 45, 2615-2623. 37.
282
Walters, W.; Murcko, M. Prediction of 'drug-likeness'. Adv. Drug Deliv. Rev. 2002, 54, 255271.
283
38.
LigPrep, version 2.9, Schrödinger, LLC, New York, NY, 2014.
284
39.
Chakravarty, S.; Kannan, K. K. Drug-protein interactions. Refined structures of three
285
sulfonamide drug complexes of human carbonic anhydrase I enzyme. J. Mol. Biol. 1994, 243,
286
298-309.
287
40.
Moeker, J.; Peat, T. S.; Bornaghi, L. F.; Vullo, D.; Supuran, C. T.; Poulsen, S. A. Cyclic
288
secondary sulfonamides: unusually good inhibitors of cancer-related carbonic anhydrase
289
enzymes. J. Med. Chem. 2014, 57, 3522-3531.
290
41.
Boriack-Sjodin, P. A.; Heck, R. W.; Laipis, P. J.; Silverman, D. N.; Christianson, D. W.
291
Structure determination of murine mitochondrial carbonic anhydrase V at 2.45-A resolution:
292
implications for catalytic proton transfer and inhibitor design. Proc. Natl. Acad. Sci. 1995, 92,
293
10949-10953.
294
42.
Gidaro, M. C.; Alcaro, F.; Carradori, S.; Costa, G.; Vullo, D.; Supuran, C. T.; Alcaro, S.
295
Eriocitrin and Apigenin as New Carbonic Anhydrase VA Inhibitors from a Virtual Screening
296
of Calabrian Natural Products. Planta Medica. 2015, 81, 533-540.
297
43.
Desmond Molecular Dynamics System, version 3.7, D. E. Shaw Research, New York, NY,
298
2014. Maestro-Desmond Interoperability Tools, version 3.7, Schrödinger, New York, NY,
299
2014.
ACS Paragon Plus Environment
16
Page 17 of 24
300
Journal of Agricultural and Food Chemistry
44.
Gidaro, M. C.; Astorino, C.; Petzer, A.; Carradori, S.; Alcaro, F.; Costa, G.; Artese, A.;
301
Rafele, G.; Russo, F. M.; Petzer, J. P.; Alcaro, S. Kaempferol as Selective Human MAO-A
302
Inhibitor: Analytical Detection in Calabrian Red Wines, Biological and Molecular Modeling
303
Studies. J. Agric. Food. Chem. 2016, 64, 1394-400.
304
45.
The PyMOL Molecular Graphics System, Version 1.7.0.0 Schrödinger, LLC
305
46.
De Simone, G.; Supuran, C. T. Antiobesity carbonic anhydrase inhibitors. Curr. Top. Med.
306 307
Chem. 2007, 7, 879-884 47.
Poulsen, S. A., Wilkinson, B. L.; Innocenti, A.; Vullo, D.; Supuran, C. T. Inhibition of human
308
mitochondrial carbonic anhydrases VA and VB with para-(4-phenyltriazole-1-yl)-
309
benzenesulfonamide derivatives. Bioorg. Med. Chem. Lett. 2008, 18, 4624-4627.
310
48.
D’Ascenzio, M.; Carradori, S.; De Monte, C.; Secci, D.; Ceruso, M.; Supuran, C. T. Design,
311
synthesis and evaluation of N-substituted saccharin derivatives as selective inhibitors of
312
tumor-associated carbonic anhydrase XII. Bioorg. Med. Chem. 2014, 22, 1821-1831.
313
49.
Khalifah, R. G. The carbon dioxide hydration activity of carbonic anhydrase. I. Stop-flow
314
kinetic studies on the native human isoenzymes B and C. J. Biol. Chem. 1971, 246, 2561-
315
2573.
316
50.
317 318
Lopez, M.; Drillaud, N.; Bornaghi, L. F.; Poulsen, S. A. Synthesis of S-glycosyl primary sulfonamides. J. Org. Chem. 2009, 74, 2811-2816.
51.
Cheng, H. C. The power issue: determination of KB or Ki from IC50. A closer look at the
319
Cheng-Prusoff equation, the Schild plot and related power equations. J. Pharmacol. Toxicol.
320
Methods 2001, 46, 61-71.
ACS Paragon Plus Environment
17
Journal of Agricultural and Food Chemistry
321
52.
Page 18 of 24
Baig, M. H.; Ahmad, K.; Adil, M.; Khan, Z. A.; Khan, M. I.; Lohani, M.; Khan, M. S.;
322
Kamal, M. A. Drug Discovery and In Silico Techniques: A Mini-Review. Enz. Eng. 2014, 4,
323
123-133.
324
53.
Innocenti, A.; Beyza Oztürk Sarikaya, S.; Gülçin, I.; Supuran, C. T. Carbonic anhydrase
325
inhibitors. Inhibition of mammalian isoforms I-XIV with a series of natural product
326
polyphenols and phenolic acids. Bioorg. Med. Chem. 2010, 18, 2159-2164.
327
54.
Ekinci, D.; Karagoz, L.; Ekinci, D.; Senturk, M.; Supuran, C. T. Carbonic anhydrase
328
inhibitors: in vitro inhibition of α isoforms (hCA I, hCA II, bCA III, hCA IV) by flavonoids.
329
J. Enzyme Inhib. Med. Chem. 2013, 28, 283-288.
330
55.
Wang, F.; Lu, W.; Zhang, T.; Dong, J.; Gao, H.; Li, P.; Wang, S.; Zhang, J. Development of
331
novel ferulic acid derivatives as potent histone deacetylase inhibitors. Bioorg. Med. Chem.
332
2013, 21, 6973-6980.
333
56.
Ononye, S. N.; VanHeyst, M. D.; Oblak, E. Z.; Zhou, W.; Ammar, M.; Anderson, A. C.;
334
Wright, D. L. Tropolones as lead-like natural products: the development of potent and
335
selective histone deacetylase inhibitors. ACS Med. Chem. Lett. 2013, 4, 757-761.
336
57.
Arechederra, R.L.; Waheed, A.; Sly, W.S.; Supuran, C.T.; Minteer, S.D. Effect of
337
sulfonamides as carbonic anhydrase VA and VB inhibitors on mitochondrial metabolic
338
energy conversion. Bioorg. Med. Chem. 2013, 21, 1544-1548.
ACS Paragon Plus Environment
18
Page 19 of 24
Journal of Agricultural and Food Chemistry
Figure captions Figure 1. Best docked poses into the catalytic site of CA VA Xray model (PDB code 1DMY41) with metal coordination of the Zinc ion (Zn) and H-bond contacts for: acetazolamide 11 (a), cocrystallized inhibitor of the PDB model, and compounds 1 (d), 2 (b) and 6 (c). All ligands are rendered as cyan carbon sticks, THR199 amino-acid residue involved in the H-bond is shown as yellow carbon sticks, while histidines coordinating Zn sphere as yellow carbon lines. The enzyme is represented as a grey cartoon with transparent surface and all non-carbon atoms are colored according to atom types. Figure 2. A schematic representation of the metal coordination between compound 2 and Zn ion within the CA VA catalytic site. The ligand-target interaction was maintained for all MD simulation time, indicated by the percentage value of 100. Figure 3. Timeline representation of H-bond contacts that CA VA target makes with ligands 1 (green), 2 (red) and 6 (yellow) compared to AZA 11 (blue) over the course of the MD simulations.
ACS Paragon Plus Environment
19
Journal of Agricultural and Food Chemistry
Page 20 of 24
Table 1. 2D chemical structures of ten hit compounds (1-10) selected in silico, by a SBVS study, for the inhibition assays against CA I, II and VA isoforms. The traditional names of the new hit compounds are: 1. β-thujaplicin; 2. 2-hydroxyisobutyric acid; 3. 4-isopropylbenzoic acid; 4. methyl geranate; 5. 3-phenylpropyl benzoate; 6. 3Z-nonenoic acid; 7. Z-geranic acid; 8. 2-methylhexanoic acid; 9. ferulic acid; 10. 5-methylfuran-2-carboxylic acid. 11. Acetazolamide is the CA inhibitor used as positive control. Biological data are reported as inhibition constant (Ki) values in µM. TARGET COMPOUND CAI CAII CA VA 1
HO
4.98
90.60
7.50
5.18
>100
3.80
3.44
85.40
7.38
8.27
>100
8.36
3.75
>100
7.83
6.09
>100
4.52
8.81
18.90
9.07
5.06
3.71
43.10
6.82
>100
9.69
4.89
>100
44.00
0.25
0.01
0.06
O O OH
2 HO
HO
3 O
O
4
O
O
5 O OH
6 O O
7 OH
8
HO O O
9
O
HO
OH
HO
10
O O N N
O
11 N H
S
O S O NH2
ACS Paragon Plus Environment
20
Page 21 of 24
Journal of Agricultural and Food Chemistry
Figure 1.
ACS Paragon Plus Environment
21
Journal of Agricultural and Food Chemistry
Page 22 of 24
Figure 2.
ACS Paragon Plus Environment
22
Page 23 of 24
Journal of Agricultural and Food Chemistry
Figure 3.
ACS Paragon Plus Environment
23
Journal of Agricultural and Food Chemistry
Page 24 of 24
For Table of Contents Use Only
Active Components of Essential Oils as Anti-Obesity Potential Drugs investigated by in silico techniques Giosuè Costa† , Maria Concetta Gidaro†,*, Daniela Vullo‡, Claudiu T. Supuran‡, Stefano Alcaro†
ACS Paragon Plus Environment
24