Chemical profile of Lippia thymoides, evaluation of the acetylcholinesterase inhibitory activity of its essential oil, and molecular docking and molecular dynamics simulations

Abstract

The essential oils of the fresh and dry flowers, leaves, branches, and roots of Lippia thymoides were obtained by hydrodistillation and analyzed using gas chromatography (GC) and GC–mass spectrometry (MS). The acetylcholinesterase inhibitory activity of the essential oil of fresh leaves was investigated on silica gel plates. The interactions of the key compounds with acetylcholinesterase were simulated by molecular docking and molecular dynamics studies. In total, 75 compounds were identified, and oxygenated monoterpenes were the dominant components of all the plant parts, ranging from 19.48% to 84.99%. In the roots, the main compounds were saturated and unsaturated fatty acids, having contents varying from 39.5% to 32.17%, respectively. In the evaluation of the anticholinesterase activity, the essential oils (detection limit (DL) = 0.1 ng/spot) were found to be about ten times less active than that of physostigmine (DL = 0.01ng/spot), whereas thymol and thymol acetate presented DL values each of 0.01 ng/spot, equivalent to that of the positive control. Based on the docking and molecular dynamics studies, thymol and thymol acetate interact with the catalytic residues Ser203 and His447 of the active site of acetylcholinesterase. The binding free energies (ΔGbind) for these ligands were -18.49 and -26.88 kcal/mol, demonstrating that the ligands are able to interact with the protein and inhibit their catalytic activity.

About the Authors:

Sebastião Gomes Silva

Roles Conceptualization

Affiliation: Program of Post-Graduation in Chemistry, Federal University of Pará, Belém, PA, Brazil

Renato Araújo da Costa

Roles Conceptualization, Formal analysis

Affiliation: Program of Post-Graduation in Chemistry, Federal University of Pará, Belém, PA, Brazil

Mozaniel Santana de Oliveira

Roles Investigation, Visualization

Affiliation: LABEX/FEA (Faculty of Food Engineering), Program of Post-Graduation in Food Science and Technology, Federal University of Para, Belém, PA, Brazil

ORCID logo http://orcid.org/0000-0002-4076-2443

Jorddy Neves da Cruz

Roles Visualization

Affiliation: Laboratory of Preparation and Computation of Nanomaterials, Federal University of Pará, Belém, PA, Brazil

ORCID logo http://orcid.org/0000-0003-0529-3714

Pablo Luis B. Figueiredo

Roles Investigation, Methodology, Validation, Visualization

Affiliations Program of Post-Graduation in Chemistry, Federal University of Pará, Belém, PA, Brazil, Department of Natural Sciences, State University of Pará, Belém, PA, Brazil

ORCID logo http://orcid.org/0000-0003-1365-3513

Davi do Socorro Barros Brasil

Roles Conceptualization

Affiliation: Program of Post-Graduation in Chemistry, Federal University of Pará, Belém, PA, Brazil

Lidiane Diniz Nascimento

Roles Investigation, Writing – original draft

Affiliations Program of Post-Graduation in Engineering of Natural Resources of Amazon, Federal University of Pará, Belém, PA, Brazil, Adolpho Ducke Laboratory, Botany Coordinating, Museu Paraense Emílio Goeldi, Belém, PA, Brazil

Antônio Maia de Jesus Chaves Neto

Roles Resources, Software

Affiliations Laboratory of Preparation and Computation of Nanomaterials, Federal University of Pará, Belém, PA, Brazil, Program of Post-Graduation in Engineering of Natural Resources of Amazon, Federal University of Pará, Belém, PA, Brazil

ORCID logo http://orcid.org/0000-0002-9730-3512

Raul Nunes de Carvalho Junior

Roles Conceptualization, Funding acquisition, Methodology, Project administration, Visualization, Writing – review & editing

* E-mail: [email protected]

Affiliations LABEX/FEA (Faculty of Food Engineering), Program of Post-Graduation in Food Science and Technology, Federal University of Para, Belém, PA, Brazil, Program of Post-Graduation in Engineering of Natural Resources of Amazon, Federal University of Pará, Belém, PA, Brazil

ORCID logo http://orcid.org/0000-0002-4433-6580

Eloisa Helena de Aguiar Andrade

Roles Conceptualization, Funding acquisition, Investigation, Methodology, Project administration, Writing – review & editing

Affiliations Program of Post-Graduation in Chemistry, Federal University of Pará, Belém, PA, Brazil, Adolpho Ducke Laboratory, Botany Coordinating, Museu Paraense Emílio Goeldi, Belém, PA, Brazil

Introduction

Alzheimer's disease is considered one of the major public health problems worldwide and one of the main complications of this pathology is the activity deficit of cholinergic neurons. This fact can be reversed and/or attenuated by elevating the levels of the neurotransmitter acetylcholine in the neuronal synaptic area. The use of cholinesterase inhibitors is an effective therapeutic approach [1]. The inhibitors increase the availability of neurotransmitters by inhibiting their main catalytic enzymes, acetyl- and butyrylcholinesterase, thus diminishing the cholinergic deficit and relieving the symptoms of Alzheimer's patients [2]. The oldest inhibitor of these enzymes is physostigmine, an alkaloid of the shrub Physostigma venenosum Balf [3]. Synthetic and semisynthetic inhibitors, such as galantamine, donepezil, tacrine, and rivastigmine, can also be used, but these drugs have disadvantages such as short half-lives and adverse side effects including hepatotoxicity and gastrointestinal irritation [1,4,5]. This has encouraged a search for new inhibitors from natural sources, and some examples are the secondary metabolites present in essential oils [1,4,6–8].

Lippia thymoides is a native and endemic Brazilian species with distribution in the states of Bahia and Minas Gerais in the Caatinga and Cerrado types of vegetation [9]. Folk medicine makes use of this plant for the treatment of wounds, and the leaves are used as an antipyretic and digestive, as well as in the treatment of bronchitis and rheumatism [10,11].

The objective of this study was to obtain and analyze the chemical composition of the essential oils from different parts of L. thymoides and to evaluate the anticholinesterase potential of the oil of the fresh leaves and its main constituents(thymol and thymol acetate), as well as to evaluate the interactions of acetylcholinesterase (AChE) with thymol and thymol acetate by molecular docking and molecular dynamics simulations and free energy calculations using the molecular mechanics generalized Born surface area (MM/GBSA) method.

Materials and methods

Plant material and essential oil extraction

The authors declare that no specific permissions were required for these locations/activities; and we confirm that the field studies did not involve endangered or protected species.

L. thymoides was collected in the Municipality of Abaetetuba, Eastern Amazon, State of Pará, Brazil. The botanical identification was made by comparison with authentic samples and incorporated into the "João Murça Pires" herbarium of the Museu Paraense Emilio Goeldi (Belém, Pará, Brazil) under catalog number MG 213373.

The essential oils were obtained from fresh and dried parts of L. thymoides.The drying of the plant material was conducted in a forced air convection oven for two days at 34°C. After that the material was ground and submitted to hydrodistillation.

Hydrodistillation

The extraction of the essential oils of fresh flowers (FFL, 10g), leaves (FLE, 40g), branches (FB, 40g), and roots (FR, 40g), and dried flowers (DFL, 10g), leaves (DLE, 40g), branches (DB, 40g), and roots (DR, 40g) was carried out by hydrodistillation in a Clevenger apparatus for 3 h. Subsequently, the oils were submitted to centrifugation, dehydrated with anhydrous Na2SO4 and stored at 8°C in a freezer. The moisture content was determined using a moisture analyzer (ID50, Marte) at the moment of extraction. The yield of essential oil (%) was calculated and is expressed as milliliters / 100-g of dried material [12].

Chemical constituent identification

The chemical composition was determined according to the approach described by Da Silva et al. [13], where the qualitative analyses of essential oils was carried using a gas chromatography–mass spectrometry (GC-MS) Thermo Focus DSQ-II system under the following the operating conditions: DB-5MS silica capillary column (30 m × 0.25 mm; 0.25 μm film thickness), temperature program: 60–250°C with a gradient of 3°C/min); injector temperature: 240°C, helium carrier gas (linear velocity of 32 cm/s, measured at 100°C), injection type: split less (0.1 μL of a 2:1000 n-hexane solution),and ion-source temperature and other parts: 200°C. The ionization was achieved by electron impact at 70 eV. Quantitative sample data were obtained using a GC with a flame ionization detector in a Focus GC-FID, which was operated under the same conditions as the GC-MS, except for the carrier gas, which was nitrogen. The identification of volatile components was based on the linear retention index (IR), which was calculated in relation to the retention times of a homologous series of n-alkanes and the fragmentation pattern observed in the mass spectra by comparison with authentic samples from the libraries of the data system and the literature [14,15].

Acetylcholinesterase assay

The AChE assay was conducted according to the method reported by Marston et al. [16] The enzyme AChE (500 U), from Electrophorus electric us (electric eel, Sigma Aldrich, Missouri, EUA, E.C. 3.1.1.7), was dissolved in tris-hydrochloric acid buffer (pH 7.8) and stabilized by the addition of bovine serum albumin fraction V (0.1%). Thymol, thymol acetate, and the essential oil of L. thymoides fresh leaves were applied to thin-layer chromatography (TLC) plates to obtain spots with concentrations from 0.01 to 1000 ng/spot. Physostigmine was used as a positive control. The plates were sprayed with the AChE solution (3.33 U/mL), dried, and incubated at 37°C for 20 min. The enzyme activity was detected by spraying with a solution of 0.25% of 1-naphtyl acetate in ethanol and a 0.25% aqueous solution of Fast Blue B salt (20 mL). Potential acetylcholinesterase inhibitors appeared as clear zones on a purple background.

Semisynthesis of thymol acetate

To obtain thymol acetate, thymol (Sigma–Aldrich), acetic anhydride, and pyridine, which acts as a catalyst, were used. Thymol (4g) was acetylated with acetic anhydride in the presence of pyridine for 24 h at 25°C. The excess acetic anhydride was removed by storage of the sample in a desiccator for 12 h. The reaction mixture was partitioned with dichloromethane and water for the removal of acetic acid. The semi synthesis and identification of the pure compound was monitored by GC and GC-MS analysis.

Molecular docking

To analyze the interactions between the ligands (thymol and thymol acetate) and AChE, molecular docking simulations were carried out using the Molegro Virtual Docker (MVD) 5.5 program [17]. The crystallographic structure of AChE was obtained from the Protein Data Bank (https://www.rcsb.org/) (PDB code: 1C2B) [18]. Both the structures of the enzyme and the ligands were prepared using the MVD module. Before the docking simulation, for each complex, hydrogen atoms and partial atomic charges were added. The active site of AChE was positioned in a spherical grid of 10-Å diameter, and all residues of AChE binding site were included using the following spatial coordinates of the central cavity: x = 26.40, y = 79.03, and z = 20.20. A grid resolution equal to 0.3 Å was used. Molecular docking was performed using the standard MolDock algorithm in MVD.

Molecular dynamics (MD) simulations

Molecular dynamics (MD) simulation for all ligand–protein systems obtained after docking were performed to evaluate possible conformational changes in the protein structure. First, the partial atomic charges were calculated using Gaussian09 [19], using the restrained electrostatic potential (RESP) protocol at the HF/6-31G* level of theory [20,21]. The antechamber [22] module of Amber 16 package was used to parameterize the ligand [23,24], which was described by the General Amber Force Field (GAFF) [25].

The protonation state of all ionizable residues was determined by the PDB2PQR server [26,27], and the protein structure was treated with the ff14SB force-field [28] in all MD simulations. The protein system was solvated in octahedral periodic box with the TIP3P explicit solvation model [29]. Counterions were added to neutralize the system charges. The system contained approximately 16,200 water molecules and 55,500 atoms. A cut-off distance of 12 Å and the particle mesh Ewald method [30] was used for electrostatic calculations, and the SHAKE algorithm was used to keep all the hydrogen bonds at their pre-defined equilibrium distances during minimization [31]. Before the MD simulation, production systems underwent simulations for energy minimization, heating, and equilibration that were performed with the Sander module and pmemd. CUDA [32].

Energy minimization was performed in two steps. In each step, 1,000 steepest descent and conjugate gradient algorithm cycles were applied. In these steps, bad contacts and possible steric conflicts were removed, and the protein system acquired the most energetically favorable conformational state. First, the protein structure was restrained with a harmonic force constant of 100 kcal/mol Å−2, and the water molecules and counterions were not treated with harmonic restraints. Then, the harmonic constraint was removed to perform the MD run of the protein system (protein, water, and counterions). This system was then heated to 300 K in five steps for 500 ps. In the first four steps, we applied a harmonic force constant of 25 kcal/mol Å−2 to restrain all heavy atoms. In the last heating step, the harmonic force constraint was removed. The Langevin thermostat [33] within a collision frequency of 2 ps−1was used to maintain the temperature. Thus, 5 ns of MD simulation without harmonic restraints at 300 K was carried out to obtain system equilibrium. Finally, a 100-ns MD simulation was performed for all protein systems.

Binding free energy calculations using the MM-GBSA approach

The molecular mechanics–generalized Born surface area (MM-GBSA) method [34,35] was used to measure the binding free energies of the ligands, thymol and thymol acetate. The binding free energies (ΔGbind) were calculated according to the following equations:(1)Each free energy state is calculated by the following equations:(2)(3)(4)

The ΔGbind values correspond to the sum of interaction energies in the gas phase between the protein and ligand (ΔEMM), the desolvation free energy (ΔGsolv), and the system entropy (−TΔS). Here, ΔEMM is the sum of internal energy (ΔEinternal), sum of bond length, angle, and dihedral energies), electrostatic contributions (ΔEelectrostatic), and van der Waals terms (ΔEvdW). ΔGsolv is the sum of the polar contributions (ΔGGB) and non-polar contributions (ΔGSA). The external dielectric constant of the solute was defined as 80, while the internal dielectric constant was defined as 1. ΔGSA was determined from the solvent accessible surface area (SASA) estimated using the LCPO algorithm. For the free energy calculations, 1,000 MD frames were used, corresponding to the last 5 ns of the simulations.

Per-residue binding free energy decomposition

To analyze the free energy contributions of each amino acid residue of protein pocket for ligand interaction, the binding free energy was decomposed into the van der Waals (ΔEvdW) and electrostatic (ΔEelectrostatic) contributions in the gas phase, as well as the polar solvation (ΔGpol) and nonpolar solvation (ΔGnonpol) contributions [36], using the following equation:(5)

Results and discussion

Yield and composition of essential oil

The yield of L. thymoides essential oil varied in quantity in the different plant parts and with treatment (fresh and dried). The flowers produced the highest essential oil yields, both fresh and dried (FFL 5.8%; DFL 7.3%). DB presented an oil yield equal to 0.14%, whereas FB, FR, and DR contained only traces of essential oil (Table 1).

[Figure omitted. See PDF.]

Table 1. Essential oil composition of different organs of Lippia thymoides.

https://doi.org/10.1371/journal.pone.0213393.t001

Generally, the essential oil yield of a plant varies depending on the part, seasonality, and geographical distribution, among other factors. For example, samples of L. citriodora produced different oil yields depending on the part extracted, and it was verified that the highest yields were obtained from the flowers [37]. Low oil yields from branches and roots (0.1% for both) were also registered in a species of domesticated L. muliflora [38]. The yield of dried leaves of L. thymoides from state of Bahia (Brazil) showed seasonal variation, and its values varied between 2.14% and 2.93% [39].

Seventy-five components were identified in the oil, comprising (95.53–99.54%) of the total composition (Table 1). Oxygenated monoterpenes were identified in all plant organs, and the percentage contents of this group varied from 19.48% (FR) to 84.99% (FFL)

The monoterpene and sesquiterpene hydrocarbons varied respectively from 0.26% (FFL) to 28.78% (DFL) and 2.41% (FB) to 12.95% (FFL), but these were absent in the essential oils of the roots. Saturated fatty acids were also identified, whose content varied from 0.06% (DFL) to 41.62% (FR), but these were absent in the leaves; unsaturated fatty acids were detected only in the branches and roots, ranging from 6.26% (FB) to 36.37% (FR).

Thymol was the only constituent present in all the organs of L. thymoides and was the main component of the oils of the flowers, leaves, and branches. This oxygenated monoterpenoid varied among all the plant organs from 19.34% (FR) to 66.33% (FLE). In the three parts where thymol was the main component. In addition, for the fourth plant sample, quantitative variations in p-cymene (0.07% FFL and 8.36% DLE), β-caryophyllene (1.29% FB and 9.55% FFL), γ-terpinene (0.15% FFL and 15.06% DFL), and thymol acetate (5.07% FB and 33.81% FFL) were observed.

The essential oils from the roots were characterized, and hexadecanoic (palmitic) acid was found to be the main component (38.02% DR; 40.92% FR), followed by (9Z)-octadecenoic (oleic) acid (27.4% DR; 28.21% FR), thymol (19.34% FR; 22.18% DR), (9Z,12Z)-9,12-octadecadienoic (linoleic) acid (4.49% FR), and (11Z)-11-hexadecenoic acid (2.73% DR; 3.67% FR). It is important to note that hexadecanoic (palmitic) acid was also present in the oils obtained from the branches (2.77% DB; 6.11% FB) and flowers (0.06% DFL; 0.16% FFL), and(9Z)-octadecenoic (oleic) acid was also identified in the oils of the branches (1.58% DB; 6.26% FB). Please refer to S1 and S2 Figs with the chromatogram ions for the chemical composition of the fractions of essential oils.

Silva et al. [39]reported that the essential oil from the leaves of a specimen of L. thymoidescontains sesquiterpeneβ-caryophyllene (17.22–26.27%) as the main component, followed by borneol (4.45–7.36%), camphor (3.22–8.61%), camphene (2.64–5.66%), and germacrene-D (4.72–6.18%). This chemical profile is different from described in the present study. This is possible related to biotic and abiotic factors of the collection site of these two specimens; according to Ribeiro et al. [40], these factors qualitatively and quantitatively affect the yield and composition of secondary metabolites.

Studies with species of the same genus, such as L. multiflora Moldenke, showed that the essential oils obtained from the leaves, flowers, branches, and roots do notvary in chemical profile, and the composition was essentially dominated by monoterpenes [38]. Similarly, the volatiles of the vegetative parts (leaves and branches, flowers, and fruit) of Lippia citriodora were characterized as having the same constituents in all parts of the plant: geranial (30.67–36.87%), neral (21.71–28.33%), and limonene (6.07–7.27%) [37].

Species of other genera, such as Hertia cheirifolia, showed significant variation in the composition of the essential oil between the different parts of the plant, where the flower buds and flowers were characterized by a higher content of drimane-type sesquiterpenes (34.08% and 32.55%, respectively), while the leaves and fruits predominantly contained α-pinene (35.63% and 33.17%, respectively) [41]. In addition, the essential oils of the roots, leaves, branches, inflorescences, and fruits of Kelussia odoratissima Mozaff. allcontained the same major compounds: (Z)-ligustilide (54.0–86.0%) and (2e)-decen-1-ol (2.0–12.3%). However, among the samples of these parts, a complex mixture of up to thirty-two different chemical compounds was identified [42].

Obtaining thymol acetate

Thymol acetate was obtained by the acetylation of thymol (Fig 1) in 99.38% purity. The chromatogram and mass spectrum of thymol acetate are shown in Figs 2 and 3.

[Figure omitted. See PDF.]

Fig 1. Synthesis reaction of thymol acetate.

https://doi.org/10.1371/journal.pone.0213393.g001

[Figure omitted. See PDF.]

Fig 2. Thymol acetate chromatogram.

https://doi.org/10.1371/journal.pone.0213393.g002

[Figure omitted. See PDF.]

Fig 3. Thymol acetate mass spectrum.

https://doi.org/10.1371/journal.pone.0213393.g003

Anticholinesterase activity

The aim of this study is to contribute to the search for new inhibitors of AChE from natural sources. Thus, the anticholinesterase activity of the essential oil from the leaves of L. thymoides was investigated. In addition, we determined the activity of two of the main constituents of the oil, thymol and thymol acetate (99.38% purity).

The essential oil of the leaves of L. thymoides (thymol: 66.33%) showed a detection limit (DL) at a concentration of 0.1 ng/spot, about ten times less active than physostigmine (DL = 0.01 ng/spot). Because thymol and thymol acetate presented a DL values of 0.01 ng/spot, equivalent to that of physostigmine, the alkaloid was used as a positive control. The higher detection concentration of thymol and thymol acetate in the oil may be because these components have a lower total content in the sample mixture compared to the isolate. Thus, the anticholinesterase action of the essential oil of L. thymoides may be tied to these two components. These results are in agreement with data from the literature, which indicates that the essential oil of Origanum ehrenbergii (thymol: 19.6%) shows strong inhibitory activity against AChE, having a 50% inhibitory concentration (IC50) of 0.3 ± 0.02 μg/mL when physostigmine was used as positive control and exhibited an IC50equal to 0.1 ± 0.003 μg/mL [43]. The essential oil of Thymus vulgaris L. (thymol: 12%) showed an IC50 of 0.216 ± 0.011 mg/mL [44]. O. vulgare subsp. vulgare essential oil (thymol: 58.31%) was found to have a concentration of1.64 ± 0.002 mg galantamine equivalents per gram (GALAEs/g) [1].

Thymol and carvacrol have anticholinesterase activity reported by several authors, with different values of IC50: 0.74 mg/mL (thymol) and IC50 equivalent to 0.063 mg/mL (carvacrol) [45]; 0.212 ± 0.011 mg/mL (thymol) and IC50 0.091 ± 0.011 mg/mL (carvacrol) [44]; 47.5 ± 1.08 mg/mL (thymol) and IC50 182 ± 1.32 mg/mL (carvacrol) [46]. These studies have shown that the anticholinesterase activity of thymol is promising. Thus, the specimen of L. thymoides described in this work can be a natural source with potential for the development of new phytotherapics for the treatment of Alzheimer's disease or even be included in the diet of people with Alzheimer's. It is important to emphasize that in our previous work [47] it was shown that this specimen presents thymol as the majority constituent of essential oil throughout its seasonal cycle, which reinforces its viability as a source of obtaining this component.

Interactions observed by molecular docking

Previously, molecular docking has been used to elucidate the interactions that occur between the ligands with the active site residues of AChE [4,48,49]. In the present study, molecular docking simulations were used to evaluate the complementarity and interactions between thymol and thymol acetate with the AChE binding site.

The binding energies, expressed by the MolDock scores obtained from molecular docking simulations, are shown in Table 2. The differences in energies can be explained by the difference in the molecular volume between thymol and thymol acetate.

[Figure omitted. See PDF.]

Table 2. Results of docking energies obtained by the MolDock score.

https://doi.org/10.1371/journal.pone.0213393.t002

Thymol and thymol acetate form hydrogen bonding (H-bond) interactions with the Ser203 catalytic residue in the AChE binding pocket (Fig 4).

[Figure omitted. See PDF.]

Fig 4. The result of molecular docking conformation obtained.

(A) Molecular interactions for Thymol and (b) Thymol acetate in AChE binding pocket.

https://doi.org/10.1371/journal.pone.0213393.g004

The O1 atom of the hydroxyl group, which belongs to thymol, forms a H-bond interaction with Ser203 at a distance of 3.14 Å. In contrast, thymol acetate forms two H-bond interactions with the Ser203 residue. The interactions with oxygen atom O2 of the Ser203 residue have distances of 2.66 and 2.98 Å.

Considering the catalytic mechanism of AChE, residues Ser203 and His447 are directly involved in covalent bond formation and breakage. Residue Ser203 is a nucleophilic site, whereas the His447 imidazole group acts as a catalytic base, accepting a proton transferred from Ser203 [50,51]. Thus, because molecular docking calculations have shown that thymol and thymol acetate interact with the same catalytic site residues of the enzyme, these inhibitors may show similar anti-AChE activity.

Structural dynamics of the AChE systems

The structural dynamics of AChE complexed with the ligands thymol and thymol acetate during MD simulations were analyzed using the root mean square deviation (RMSD). The RMSD plots of the AChE structure were analyzed based on the Cα backbone, whereas the ligand structure was analyzed using the heavy atoms alone.

The RMSD plot of AChE complexed to thymol and thymol acetate reaches a plateau (Fig 5), indicating that AChE structure is stable, showing few conformational changes over the 100-ns MD simulation. The ligands showed high stability at the AChE binding site. Over the 100-ns MD run, the ligands remained bound to AChE and did not undergo drastic conformational changes that altered their interactions with the catalytic residues.

[Figure omitted. See PDF.]

Fig 5. RMSD values for 100 ns of MD for AChE systems.

The protein backbone is colored in black, thymol in blue and thymol acetate in green. (a) RMSD sistem plot Thymol- AChE e (b) RMSD sistem plot Thymol acetate-AChE.

https://doi.org/10.1371/journal.pone.0213393.g005

Binding free energy decomposition

Binding energy calculations show that thymol and thymol interact with high affinity in the AChE binding pocket, having ΔGbind values of -18.49 and -26.88 kcal/mol, respectively. The binding free energy decomposition reveals that the van der Waals (ΔEvdW) contribution represents a large proportion of the interactions of the ligands in AChE binding pocket. The electrostatic (ΔEelectrostatic) and non-polar (ΔGSA) contributions were also favorable for the formation of the complexes. The ΔG values and their components are listed in Table 3.

[Figure omitted. See PDF.]

Table 3. Free energy variation (ΔGbind) and its components.

ΔEvdW represents the Van de Waals energy contribution, ΔEele represents the electrostatic energy, ΔGGB polar contribution and ΔGSA non-polar contribution. All values are expressed in kcal/mol.

https://doi.org/10.1371/journal.pone.0213393.t003

Intermolecular interactions analyses

To explore the relative energy contribution to the overall ligand binding energy for each residue in the AChE binding pocket, the per-residue binding free energy decomposition was analyzed in Amber16. The results obtained are shown in Fig 6.

[Figure omitted. See PDF.]

Fig 6.

Binding free energy decomposition for each residue that interact with (A) thymol and (B) thymol acetate. The bars represent the energy contribution values: van der Waals contributions (blank bars), electrostatic contributions (striped bars) and total energy contribution for each residue (black bars). Residues highlighted in red are from catalytic site, blue are from anionic subsite and green are from oxyanion hole.

https://doi.org/10.1371/journal.pone.0213393.g006

The AChE structure contains different binding sites [52]. The catalytic residues Ser203, Glu337, and His447 are located in the same cavity of the enzyme, which also contains other important residues for enzymatic activity, for example, Gly121, Gly122, and Ala204 (oxyanion hole), Trp86, Tyr133, Tyr337, and Phe338 (anionic subsite), and Phe295 and Phe297 (acyl pocket) [53].

Thymol acetate (ΔGbind = -26.88 kcal/mol) has a better binding energy than that ofthymol (ΔGbind = -18.49 kcal/mol); that is, the binding energy of thymol acetate is greater than that of thymol, and both the van der Waals and electrostatic energies are important for ligand binding in AChE cavity.

Thymol formed intermolecular interactions with different residues of the AChE binding pocket, in particular, with residue His447, which belongs to the catalytic triad (His 447, Tyr337, and Trp86) located on the anionic subsite. Residues His447 and Gly122 belong to the oxyanion hole. In contrast, thymol acetate formed different intermolecular interactions with AChE binding pocket and the electrostatic and van der Waals energy contributions were found to be significant in the interactions with residues Ser203 and His447 (catalytic triad), thus stabilizing ligand binding. Thymol acetate also formed important interactions with residues from two other sites, Trp86 and Tyr133 (anionic subsite) and Gly121 (oxyanion hole).

Conclusions

Thymol is the main constituent in the essential oil from the flowers (37.86% DFL; 48.04% FFL), leaves (58.9% DLE; 66.33% FLE), and branches (63.59% FB; 66.2% DB) of L. thymoides, and the essential oil of the fresh leaves had the highest content of this oxygenated monoterpene. In the roots, hexadecanoic acid (palmitic acid) was the main constituent (38.02% DR; 40.92% DF) of the oil. The essential oil of the fresh leaves (DL = 0.1 ng/spot) and two of its mains components, thymol (DL = 0.01 ng/spot) and thymol acetate (DL = 0.01 ng/spot), showed inhibitory activity against acetylcholinesterase on TLC layers.

The MolDock scores obtained for thymol and thymol acetate were favorable. The poses obtained from molecular docking simulations showed that the ligands forms H-bond interactions with the Ser203 residue, which belongs to the catalytic triad. The RMSD values obtained over 100 ns of MD simulation showed that the ligands are stable in the AChE binding pocket. The binding free energies for thymol (ΔGbind = -18.49 kcal/mol) and thymol acetate (ΔGbind -26.88 kcal/mol) indicate the ligands are stable and bind with affinity to AChE. The per-residue binding free energy decomposition revealed that the ligands form interactions with residues that are important for catalytic activity. Some residues belonging to the catalytic triad (Ser203 and His447), anionic subsite (Trp86, Tyr333, and Tyr337), and oxyanion hole (Gly121 and Gly122) form intermolecular interactions that stabilize ligand binding.

Supporting information

[Figure omitted. See PDF.]

S1 Fig. ion chromatogram relative to the essential oil composition of leaves and flowers of Lippia thymoides.

https://doi.org/10.1371/journal.pone.0213393.s001

(DOCX)

S2 Fig. ion chromatogram relative to the essential oil composition of branches and roots of Lippia thymoides.

https://doi.org/10.1371/journal.pone.0213393.s002

(DOCX)

Citation: Silva SG, da Costa RA, de Oliveira MS, da Cruz JN, Figueiredo PLB, Brasil DdSB, et al. (2019) Chemical profile of Lippia thymoides, evaluation of the acetylcholinesterase inhibitory activity of its essential oil, and molecular docking and molecular dynamics simulations. PLoS ONE 14(3): e0213393. https://doi.org/10.1371/journal.pone.0213393

References

1. Sarikurkcu C, Zengin G, Oskay M, Uysal S, Ceylan R, Aktumsek A. Composition, antioxidant, antimicrobial and enzyme inhibition activities of two Origanum vulgare subspecies (subsp. vulgare and subsp. hirtum) essential oils. Ind Crops Prod. 2015;70: 178–184.

2. Kurt BZ, Gazioglu I, Dag A, Salmas RE, Kayık G, Durdagi S, et al. Synthesis, anticholinesterase activity and molecular modeling study of novel carbamate-substituted thymol/carvacrol derivatives. Bioorg Med Chem. 2017;25: 1352–1363. pmid:28089589

3. Docherty JR. The pharmacology of α1- and α2-adrenoceptors: Evidence for and against a further subdivision. Pharmacol Ther. 1989;44: 241–284. pmid:2577511

4. Silva NNS, Silva JRA, Alves CN, Andrade EHA, da Silva JKR, Maia JGS. Acetylcholinesterase Inhibitory Activity and Molecular Docking Study of 1-Nitro-2-Phenylethane, the Main Constituent of Aniba canelilla Essential Oil. Chem Biol Drug Des. 2014;84: 192–198. pmid:24661632

5. Wszelaki N, Kuciun A, Kiss A. Screening of traditional European herbal medicines for acetylcholinesterase and butyrylcholinesterase inhibitory activity. Acta Pharm. 2010;60: 119–128. pmid:20228046

6. da Silva JKR, Pinto LC, Burbano RMR, Montenegro RC, Guimarães EF, Andrade EHA, et al. Essential oils of Amazon Piper species and their cytotoxic, antifungal, antioxidant and anti-cholinesterase activities. Ind Crops Prod. 2014;58: 55–60.

7. Jukic M, Politeo O, Maksimovic M, Milos M, Milos M. In Vitroacetylcholinesterase inhibitory properties of thymol, carvacrol and their derivatives thymoquinone and thymohydroquinone. Phyther Res. 2007;21: 259–261. pmid:17186491

8. Orhan I, Şener B, Kartal M, Kan Y. Activity of Essential Oils and Individual Components against Acetyl-and Butyrylcholinesterase. Zeitschrift fur Naturforsch—Sect C J Biosci. 2008;63: 547–553.

9. Salimena FRG, Mulgura M 2015. Lippia in Lista de Espécies da Flora do Brasil. Jardim Botânico do Rio de Janeiro. [Internet]. 2015.

10. Almeida MZ De. Plantas Medicinais. 3a. Edufba. Salvador-BA; 2011.

11. Funch LS, Harley RM, Funch R, Giulietti AM, Melo E. Plantas Úteis da Chapada Diamantina. São Carlos, Brazil: Rima; 2004.

12. Figueiredo PLB, Silva RC, Kelly J, Silva R, Suemitsu C, Mourão RH V, et al. Chemical variability in the essential oil of leaves of Araçá (Psidium guineense Sw.), with occurrence in the Amazon. Chem Cent J. Springer International Publishing; 2018; 1–11.

13. da Silva JKR, Silva NNS, Santana JFS, Andrade EHA, Maia JGS, Setzer WN. Phenylpropanoid-rich essential oils of Piper species from the Amazon and their antifungal and anti-cholinesterase activities. Nat Prod Commun. 2016;11: 1907–1911. pmid:30508363

14. ADAMS RP. Identification of Essential Oil Components by Gas Chromatography/Mass Spectrometry. 4th ed. Carol Stream, IL: Allured Publishing Corporation; 2007.

15. NIST/EPANI Mass Spectral Database. 2011.

16. Marston A, Kissling J, Hostettmann K. A rapid TLC bioautographic method for the detection of acetylcholinesterase and butyrylcholinesterase inhibitors in plants. Phytochem Anal. 2002;13: 51–54. pmid:11899607

17. Thomsen R, Christensen MH. MolDock: A New Technique for High-Accuracy Molecular Docking. J Med Chem. 2006;49: 3315–3321. pmid:16722650

18. Bourne Y, Grassi J, Bougis PE, Marchot P. Conformational Flexibility of the Acetylcholinesterase Tetramer Suggested by X-ray Crystallography. J Biol Chem. 1999;274: 30370–30376. pmid:10521413

19. Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, et al. Gaussian 09. Pittsburgh, Wallingford CT; 2009.

20. Cornell WD, Cieplak P, Bayly CI, Kollmann PA. Application of RESP charges to calculate conformational energies, hydrogen bond energies, and free energies of solvation. J Am Chem Soc. 1993;115: 9620–9631.

21. Wang J, Cieplak P, Kollman PA. How well does a restrained electrostatic potential (RESP) model perform in calculating conformational energies of organic and biological molecules? J Comput Chem. 2000;21: 1049–1074.

22. Wang J, Wang W, Kollman PA, Case DA. Automatic atom type and bond type perception in molecular mechanical calculations. J Mol Graph Model. 2006;25: 247–260. pmid:16458552

23. Case DA, Cheatham TE, Darden T, Gohlke H, Luo R, Merz KM, et al. The Amber biomolecular simulation programs. J Comput Chem. 2005;26: 1668–1688. pmid:16200636

24. Cruz JN, Costa JFS, Khayat AS, Kuca K, Barros CAL, Neto AMJC. Molecular dynamics simulation and binding free energy studies of novel leads belonging to the benzofuran class inhibitors of Mycobacterium tuberculosis Polyketide Synthase 13. J Biomol Struct Dyn. Taylor & Francis; 2018;1102: 1–12. pmid:29633908

25. Wang J, Wolf RM, Caldwell JW, Kollman PA, Case DA. Development and testing of a general amber force field. J Comput Chem. 2004;25: 1157–1174. pmid:15116359

26. Dolinsky TJ, Nielsen JE, McCammon JA, Baker NA. PDB2PQR: an automated pipeline for the setup of Poisson-Boltzmann electrostatics calculations. Nucleic Acids Res. 2004;32: W665–W667. pmid:15215472

27. Li H, Robertson AD, Jensen JH. Very fast empirical prediction and rationalization of protein pKa values. Proteins Struct Funct Bioinforma. 2005;61: 704–721. pmid:16231289

28. Maier JA, Martinez C, Kasavajhala K, Wickstrom L, Hauser KE, Simmerling C. ff14SB: Improving the Accuracy of Protein Side Chain and Backbone Parameters from ff99SB. J Chem Theory Comput. 2015;11: 3696–3713. pmid:26574453

29. Jorgensen WL, Chandrasekhar J, Madura JD, Impey RW, Klein ML. Comparison of simple potential functions for simulating liquid water. J Chem Phys. 1983;79: 926–935.

30. Darden T, York D, Pedersen L. Particle mesh Ewald: An N ⋅log (N) method for Ewald sums in large systems. J Chem Phys. 1993;98: 10089–10092.

31. Ryckaert J-P, Ciccotti G, Berendsen HJC. Numerical integration of the cartesian equations of motion of a system with constraints: molecular dynamics of n-alkanes. J Comput Phys. 1977;23: 327–341.

32. Salomon-Ferrer R, Götz AW, Poole D, Le Grand S, Walker RC. Routine Microsecond Molecular Dynamics Simulations with AMBER on GPUs. 2. Explicit Solvent Particle Mesh Ewald. J Chem Theory Comput. 2013;9: 3878–3888. pmid:26592383

33. Izaguirre JA, Catarello DP, Wozniak JM, Skeel RD. Langevin stabilization of molecular dynamics. J Chem Phys. 2001;114: 2090–2098.

34. Genheden S, Ryde U. The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin Drug Discov. 2015;10: 449–461. pmid:25835573

35. Kollman PA, Massova I, Reyes C, Kuhn B, Huo S, Chong L, et al. Calculating Structures and Free Energies of Complex Molecules: Combining Molecular Mechanics and Continuum Models. Acc Chem Res. 2000;33: 889–897. pmid:11123888

36. Gohlke H, Kiel C, Case DA. Insights into Protein–Protein Binding by Binding Free Energy Calculation and Free Energy Decomposition for the Ras–Raf and Ras–RalGDS Complexes. J Mol Biol. 2003;330: 891–913. pmid:12850155

37. Shahhoseini R, Hosseini N, Ghorbanpour M. Study of Essential Oil Content and Composition of Different Parts of Lemon verbena (Lippia citriodora) Grown in Iran. J Essent Oil Bear Plants. 2014;17: 120–125.

38. Tsiba G, Célestine NL, Yaya M, Jean O, Antoine AA, Jean-claude C, et al. Variation in the chemical composition of the essential oils of different organs of domesticated Lippia multiflora Moldenke. 2010;9: 7009–7013.

39. Silva FS, Menezes PMN, Sá PGS De, Oliveira ALDS, Souza EAA, Almeida JRG da S, et al. Chemical composition and pharmacological properties of the essential oils obtained seasonally from Lippia thymoides. Pharm Biol. 2016;54: 25–34. pmid:25856708

40. Ribeiro AF, Andrade EHA, Salimena FRG, Maia JGS. Circadian and seasonal study of the cinnamate chemotype from Lippia origanoides Kunth. Biochem Syst Ecol. Elsevier Ltd; 2014;55: 249–259.

41. Rahali N, Mehdi S, Younsi F, Boussaid M, Messaoud C. Antioxidant, α-amylase, and acetylcholinesterase inhibitory activities of Hertia cheirifolia essential oils: Influence of plant organs and seasonal variation. Int J Food Prop. 2017;2912: 1–15.

42. Raeisi S, Mirjalili MH, Nadjafi F, Hadian J. Variability in the essential oil content and composition in different plant organs of Kelussia odoratissima Mozaff. (Apiaceae) growing wild in Iran. J Essent Oil Res. 2015;27: 283–288.

43. Loizzo MR, Menichini F, Conforti F, Tundis R, Bonesi M, Saab AM, et al. Chemical analysis, antioxidant, antiinflammatory and anticholinesterase activities of Origanum ehrenbergii Boiss and Origanum syriacum L. essential oils. Food Chem. 2009;117: 174–180.

44. Aazza S, Lyoussi B, Miguel MG. Antioxidant and Antiacetylcholinesterase Activities of Some Commercial Essential Oils and Their Major Compounds. Molecules. 2011;16: 7672–7690. pmid:21900869

45. Jukic M, Politeo O, Maksimovic M, Milos M, Milos M. In Vitroacetylcholinesterase inhibitory properties of thymol, carvacrol and their derivatives thymoquinone and thymohydroquinone. Phyther Res. 2007;21: 259–261. pmid:17186491

46. Öztürk M. Anticholinesterase and antioxidant activities of Savoury (Satureja thymbra L.) with identified major terpenes of the essential oil. Food Chem. 2012;134: 48–54.

47. Silva SG, Figueiredo PLB, Nascimento LD, da Costa WA, Maia JGS, Andrade EHA. Planting and seasonal and circadian evaluation of a thymol-type oil from Lippia thymoides Mart. & Schauer. Chem Cent J. Springer International Publishing; 2018;12: 113. pmid:30421173

48. Nemukhin A V., Lushchekina S V., Bochenkova A V., Golubeva AA, Varfolomeev SD. Characterization of a complete cycle of acetylcholinesterase catalysis by ab initio QM/MM modeling. J Mol Model. 2008;14: 409–416. pmid:18343962

49. Pilger C, Bartolucci C, Lamba D, Tropsha A, Fels G. Accurate prediction of the bound conformation of galanthamine in the active site of torpedo californica acetylcholinesterase using molecular docking. J Mol Graph Model. 2001;19: 288–296. pmid:11449566

50. Quinn DM. Acetylcholinesterase: enzyme structure, reaction dynamics, and virtual transition states. Chem Rev. 1987;87: 955–979.

51. Rosenberry TL. Acetylcholinesterase. Adv Enzym Relat Areas Mol Biol. 1975;43: 103–218.

52. Xu Y, Colletier J-P, Weik M, Jiang H, Moult J, Silman I, et al. Flexibility of Aromatic Residues in the Active-Site Gorge of Acetylcholinesterase: X-ray versus Molecular Dynamics. Biophys J. Elsevier; 2008;95: 2500–2511. pmid:18502801

53. Ordentlich A, Barak D, Kronman C, Flashner Y, Leitner M, Segall Y, et al. Dissection of the human acetylcholinesterase active center determinants of substrate specificity. Identification of residues constituting the anionic site, the hydrophobic site, and the acyl pocket. J Biol Chem. 1993;268: 17083–95. pmid:8349597