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

Sebastião Gomes Silva; Renato Araújo da Costa; Mozaniel Santana de Oliveira; Jorddy Neves da Cruz; Pablo Luis B. Figueiredo; Davi do Socorro Barros Brasil; Lidiane Diniz Nascimento; Antônio Maia de Jesus Chaves Neto; Raul Nunes de Carvalho Junior; Eloisa Helena de Aguiar Andrade

doi:10.1371/journal.pone.0213393

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 (ΔG_bind) 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.

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

Editor: Marco Bonizzoni, The University of Alabama, UNITED STATES

Received: October 22, 2018; Accepted: February 19, 2019; Published: March 8, 2019

Copyright: © 2019 Silva et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Data Availability: All relevant data are within the manuscript and its Supporting Information files.

Funding: This work was supported by scholarships awarded to SGS. SGS thanks Secretaria Estadual de Educação do Pará.; Oliveira M. S (Process Number: 1662230) and Costa W. A. (Process Number: 1427204), thanks CAPES for the doctorate scholarship. The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Competing interests: The authors have declared that no competing interests exist.

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 Na₂SO₄ 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 Å⁻² 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⁻¹was 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 (ΔG_bind) were calculated according to the following equations: (1) Each free energy state is calculated by the following equations: (2) (3) (4)

The ΔG_bind values correspond to the sum of interaction energies in the gas phase between the protein and ligand (ΔE_MM), the desolvation free energy (ΔG_solv), and the system entropy (−TΔS). Here, ΔE_MM is the sum of internal energy (ΔE_internal), sum of bond length, angle, and dihedral energies), electrostatic contributions (ΔE_{electrostatic}), and van der Waals terms (ΔE_vdW). ΔG_solv is the sum of the polar contributions (ΔG_GB) and non-polar contributions (ΔG_SA). The external dielectric constant of the solute was defined as 80, while the internal dielectric constant was defined as 1. ΔG_SA 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 (ΔE_vdW) and electrostatic (ΔE_{electrostatic}) contributions in the gas phase, as well as the polar solvation (ΔG_pol) and nonpolar solvation (ΔG_nonpol) 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).

Download:

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.

Download:

Fig 1. Synthesis reaction of thymol acetate.

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

Download:

Fig 2. Thymol acetate chromatogram.

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

Download:

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 (IC₅₀) of 0.3 ± 0.02 μg/mL when physostigmine was used as positive control and exhibited an IC₅₀equal to 0.1 ± 0.003 μg/mL [43]. The essential oil of Thymus vulgaris L. (thymol: 12%) showed an IC₅₀ 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 IC₅₀: 0.74 mg/mL (thymol) and IC₅₀ equivalent to 0.063 mg/mL (carvacrol) [45]; 0.212 ± 0.011 mg/mL (thymol) and IC₅₀ 0.091 ± 0.011 mg/mL (carvacrol) [44]; 47.5 ± 1.08 mg/mL (thymol) and IC₅₀ 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.

Download:

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).

Download:

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.

Download:

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 ΔG_bind values of -18.49 and -26.88 kcal/mol, respectively. The binding free energy decomposition reveals that the van der Waals (ΔE_vdW) contribution represents a large proportion of the interactions of the ligands in AChE binding pocket. The electrostatic (ΔE_{electrostatic}) and non-polar (ΔG_SA) contributions were also favorable for the formation of the complexes. The ΔG values and their components are listed in Table 3.

Download:

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

ΔE_vdW represents the Van de Waals energy contribution, ΔE_ele represents the electrostatic energy, ΔG_GB 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.

Download:

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 (ΔG_bind = -26.88 kcal/mol) has a better binding energy than that ofthymol (ΔG_bind = -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 (ΔG_bind = -18.49 kcal/mol) and thymol acetate (ΔG_bind -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

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)

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.
- View Article
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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.
- View Article
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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.
- View Article
- Google Scholar
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.
- View Article
- Google Scholar
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.
- View Article
- Google Scholar
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.
- View Article
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
17. Thomsen R, Christensen MH. MolDock: A New Technique for High-Accuracy Molecular Docking. J Med Chem. 2006;49: 3315–3321. pmid:16722650
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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.
- View Article
- Google Scholar
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.
- View Article
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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.
- View Article
- Google Scholar
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.
- View Article
- Google Scholar
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.
- View Article
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
33. Izaguirre JA, Catarello DP, Wozniak JM, Skeel RD. Langevin stabilization of molecular dynamics. J Chem Phys. 2001;114: 2090–2098.
- View Article
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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.
- View Article
- Google Scholar
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.
- View Article
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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.
- View Article
- Google Scholar
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.
- View Article
- Google Scholar
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.
- View Article
- Google Scholar
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.
- View Article
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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.
- View Article
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
50. Quinn DM. Acetylcholinesterase: enzyme structure, reaction dynamics, and virtual transition states. Chem Rev. 1987;87: 955–979.
- View Article
- Google Scholar
51. Rosenberry TL. Acetylcholinesterase. Adv Enzym Relat Areas Mol Biol. 1975;43: 103–218.
- View Article
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar
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
- View Article
- PubMed/NCBI
- Google Scholar

[ref1] 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.
View Article
Google Scholar

[2] View Article

[3] Google Scholar

[ref2] 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
View Article
PubMed/NCBI
Google Scholar

[5] View Article

[6] PubMed/NCBI

[7] Google Scholar

[ref3] 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
View Article
PubMed/NCBI
Google Scholar

[9] View Article

[10] PubMed/NCBI

[11] Google Scholar

[ref4] 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
View Article
PubMed/NCBI
Google Scholar

[13] View Article

[14] PubMed/NCBI

[15] Google Scholar

[ref5] 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
View Article
PubMed/NCBI
Google Scholar

[17] View Article

[18] PubMed/NCBI

[19] Google Scholar

[ref6] 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.
View Article
Google Scholar

[21] View Article

[22] Google Scholar

[ref7] 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
View Article
PubMed/NCBI
Google Scholar

[24] View Article

[25] PubMed/NCBI

[26] Google Scholar

[ref8] 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.
View Article
Google Scholar

[28] View Article

[29] Google Scholar

[ref9] 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.
View Article
Google Scholar

[31] View Article

[32] Google Scholar

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

[ref11] 11. Funch LS, Harley RM, Funch R, Giulietti AM, Melo E. Plantas Úteis da Chapada Diamantina. São Carlos, Brazil: Rima; 2004.
View Article
Google Scholar

[35] View Article

[36] Google Scholar

[ref12] 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.
View Article
Google Scholar

[38] View Article

[39] Google Scholar

[ref13] 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
View Article
PubMed/NCBI
Google Scholar

[41] View Article

[42] PubMed/NCBI

[43] Google Scholar

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

[ref15] 15. NIST/EPANI Mass Spectral Database. 2011.

[ref16] 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
View Article
PubMed/NCBI
Google Scholar

[47] View Article

[48] PubMed/NCBI

[49] Google Scholar

[ref17] 17. Thomsen R, Christensen MH. MolDock: A New Technique for High-Accuracy Molecular Docking. J Med Chem. 2006;49: 3315–3321. pmid:16722650
View Article
PubMed/NCBI
Google Scholar

[51] View Article

[52] PubMed/NCBI

[53] Google Scholar

[ref18] 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
View Article
PubMed/NCBI
Google Scholar

[55] View Article

[56] PubMed/NCBI

[57] Google Scholar

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

[ref20] 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.
View Article
Google Scholar

[60] View Article

[61] Google Scholar

[ref21] 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.
View Article
Google Scholar

[63] View Article

[64] Google Scholar

[ref22] 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
View Article
PubMed/NCBI
Google Scholar

[66] View Article

[67] PubMed/NCBI

[68] Google Scholar

[ref23] 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
View Article
PubMed/NCBI
Google Scholar

[70] View Article

[71] PubMed/NCBI

[72] Google Scholar

[ref24] 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
View Article
PubMed/NCBI
Google Scholar

[74] View Article

[75] PubMed/NCBI

[76] Google Scholar

[ref25] 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
View Article
PubMed/NCBI
Google Scholar

[78] View Article

[79] PubMed/NCBI

[80] Google Scholar

[ref26] 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
View Article
PubMed/NCBI
Google Scholar

[82] View Article

[83] PubMed/NCBI

[84] Google Scholar

[ref27] 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
View Article
PubMed/NCBI
Google Scholar

[86] View Article

[87] PubMed/NCBI

[88] Google Scholar

[ref28] 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
View Article
PubMed/NCBI
Google Scholar

[90] View Article

[91] PubMed/NCBI

[92] Google Scholar

[ref29] 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.
View Article
Google Scholar

[94] View Article

[95] Google Scholar

[ref30] 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.
View Article
Google Scholar

[97] View Article

[98] Google Scholar

[ref31] 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.
View Article
Google Scholar

[100] View Article

[101] Google Scholar

[ref32] 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
View Article
PubMed/NCBI
Google Scholar

[103] View Article

[104] PubMed/NCBI

[105] Google Scholar

[ref33] 33. Izaguirre JA, Catarello DP, Wozniak JM, Skeel RD. Langevin stabilization of molecular dynamics. J Chem Phys. 2001;114: 2090–2098.
View Article
Google Scholar

[107] View Article

[108] Google Scholar

[ref34] 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
View Article
PubMed/NCBI
Google Scholar

[110] View Article

[111] PubMed/NCBI

[112] Google Scholar

[ref35] 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
View Article
PubMed/NCBI
Google Scholar

[114] View Article

[115] PubMed/NCBI

[116] Google Scholar

[ref36] 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
View Article
PubMed/NCBI
Google Scholar

[118] View Article

[119] PubMed/NCBI

[120] Google Scholar

[ref37] 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.
View Article
Google Scholar

[122] View Article

[123] Google Scholar

[ref38] 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.
View Article
Google Scholar

[125] View Article

[126] Google Scholar

[ref39] 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
View Article
PubMed/NCBI
Google Scholar

[128] View Article

[129] PubMed/NCBI

[130] Google Scholar

[ref40] 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.
View Article
Google Scholar

[132] View Article

[133] Google Scholar

[ref41] 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.
View Article
Google Scholar

[135] View Article

[136] Google Scholar

[ref42] 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.
View Article
Google Scholar

[138] View Article

[139] Google Scholar

[ref43] 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.
View Article
Google Scholar

[141] View Article

[142] Google Scholar

[ref44] 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
View Article
PubMed/NCBI
Google Scholar

[144] View Article

[145] PubMed/NCBI

[146] Google Scholar

[ref45] 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
View Article
PubMed/NCBI
Google Scholar

[148] View Article

[149] PubMed/NCBI

[150] Google Scholar

[ref46] 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.
View Article
Google Scholar

[152] View Article

[153] Google Scholar

[ref47] 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
View Article
PubMed/NCBI
Google Scholar

[155] View Article

[156] PubMed/NCBI

[157] Google Scholar

[ref48] 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
View Article
PubMed/NCBI
Google Scholar

[159] View Article

[160] PubMed/NCBI

[161] Google Scholar

[ref49] 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
View Article
PubMed/NCBI
Google Scholar

[163] View Article

[164] PubMed/NCBI

[165] Google Scholar

[ref50] 50. Quinn DM. Acetylcholinesterase: enzyme structure, reaction dynamics, and virtual transition states. Chem Rev. 1987;87: 955–979.
View Article
Google Scholar

[167] View Article

[168] Google Scholar

[ref51] 51. Rosenberry TL. Acetylcholinesterase. Adv Enzym Relat Areas Mol Biol. 1975;43: 103–218.
View Article
Google Scholar

[170] View Article

[171] Google Scholar

[ref52] 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
View Article
PubMed/NCBI
Google Scholar

[173] View Article

[174] PubMed/NCBI

[175] Google Scholar

[ref53] 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
View Article
PubMed/NCBI
Google Scholar

[177] View Article

[178] PubMed/NCBI

[179] Google Scholar

Figures

Abstract

Introduction

Materials and methods

Plant material and essential oil extraction

Hydrodistillation

Chemical constituent identification

Acetylcholinesterase assay

Semisynthesis of thymol acetate

Molecular docking

Molecular dynamics (MD) simulations

Binding free energy calculations using the MM-GBSA approach

Per-residue binding free energy decomposition

Results and discussion

Yield and composition of essential oil

Obtaining thymol acetate

Anticholinesterase activity

Interactions observed by molecular docking

Structural dynamics of the AChE systems

Binding free energy decomposition

Intermolecular interactions analyses

Conclusions

Supporting information

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

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

References