Preprint
Article

Molecular and Pharmacokinetic Aspects of Acetylcholinesterase Inhibitory Potential of the Oleanane-type Triterpenes and Their Glycosides

Altmetrics

Downloads

105

Views

37

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

17 August 2023

Posted:

21 August 2023

You are already at the latest version

Alerts
Abstract
The acetylcholinesterase-inhibitory potential of the oleanane-type triterpenes and their glycosides from the Terminalia arjuna (Combreatceae) bark, i.e. arjunic acid, arjunolic acid, arjungenin, arjunglucoside I, sericic acid, and arjunetin, is presented. The studies are based on the in silico pharmacokinetic and biomimetic studies, the acetylcholinesterase (AChE) inhibitory activity tests, and the molecular docking research. Based on the calculated pharmacokinetic parameters, arjunetin and arjunglucoside I are indicated as able to cross the blood-brain barrier. The compounds of interest exhibit marked acetylcholinesterase inhibitory potential, which was tested in the TLC bioautography test. The longest time to reach brain equilibrium is observed for both the arjunic and arjunolic acids and the shortest one for arjunetin. All compounds exhibit high and relatively similar magnitude of binding energies, varying from ca. -15 to -13 kcal/mol. The superposition of the most favorable positions of all ligands interacting with AChE is analyzed. The correlation between the experimentally determined IC50 values and the steric parameters of the molecules is investigated. The inhibition of the enzyme by the analyzed compounds shows their potential to be used as cognition enhancing agents. For the most potent compound (arjunglucoside I; ARG), the kinetics of AChE inhibition is tested. The Michaelis–Menten constant (Km) for the hydrolysis of the acetylthiocholine iodide substrate was calculated to be 0.011 mM.
Keywords: 
Subject: Chemistry and Materials Science  -   Medicinal Chemistry

1. Introduction

Neurodegenerative diseases with memory impairment constitute a growing health issue in aging populations worldwide. Life extension in the last century resulting from the progress in medicine is related to the increasing neurodegenerative diseases morbidity, especially in developed countries. The World Health Organization recognized dementia as a public health priority. In 2017, the World Health Assembly endorsed “Global action plan on the public health response to dementia 2017-2025” [1]. According to the WHO report, in 2015, dementia affected 47 million people worldwide (or roughly 5% of the world’s elderly population). This number is predicted to increase to 75 million in 2030 and 132 million by 2050 [1]. Therefore, one of the key research requirements is to increase the effectiveness of existing and development of new therapeutic strategies as regards the treatment of neurodegenerative diseases with memory impairment [2].
Acetylcholinesterase (AChE) is an enzyme excreted to the synaptic cleft during the formation of action potential of a neuron in the brain. To memorize or resume memories the chemical activity of neurons is sustained by the excretion of acetylcholine (ACh) to the synaptic cleft. ACh is then targeting the post-synaptic membrane where it binds with cholinergic receptors to sustain the action potential. Later, an excessive amount of ACh is decomposed by the enzyme, namely AChE to prepare the neuron for the following activities. However, in elderly, a decreased volume of ACh is excreted to the synaptic cleft together with an elevated volume of AChE enzyme, which, in fact may lead to a significant reduction of the action potential number or strength [3]. The described changes on the chemical level in brain trigger memory impairment effects. That is why, the inhibition of AChE enzyme was found to stimulate neuronal functions and by this action led to an improved cognition. The inhibitors of AChE enzyme are now the first line drugs in the treatment of dementia and the Alzheimer’s disease (AD) [4,5,6]. Even if there are several registered drugs from the group of AChE inhibitors in the market, still there is a need for new, better tolerable, longer acting and stronger compounds.
Plant-derived secondary metabolites including oleanane-type triterpenes and their glycosides have demonstrated to be interesting sources of compounds with neuroprotective as well as memory enhancing effects e.g. medicagosides A-F from Medicago sativa L. [7,8], onjisaponins isolated from the roots of Polygala tenuifolia [9,10] or platycodins from Platycodi radix [11]. Also the Terminalia arjuna (Combreatceae), being a valued tree widely used in Indian traditional medicine, seems to be worthy of interest as a natural source of compounds with a cognition enhancing potential, e.g. arjunolic acid with neurons-protective potential from oxidative stress associated damage [12].
The herein described study will focus on the determination of AChE inhibitory potential by naturally occurring oleanane-type triterpenes and their glycosides as well as the elucidation the structural pattern of the interactions between the studied compounds and the AChE enzyme employing the ligand-protein docking methodology. The recognition of the AChE-inhibitory potential will be performed using the thin layer chromatography (TLC) bioautography method. In addition, both values of pharmacokinetic and physicochemical parameters connected with the blood-brain barrier (BBB) permeation will be determined using computational and biomimetic methods. They are used as screening tests in the drug discovery process [13] to study the broad-spectrum of biological activity, including the ability to cross specific biological barriers. The BBB is a selective barrier with the endothelium forming a much tighter interface than peripheral endothelia because the gaps between the capillary endothelial cells in most part of the brain are sealed by tight junctions and thus have a severely limited permeability [14]. The experimental determination of BBB permeability based on in vivo studies requires complex techniques which are usually time consuming and expensive [13,15]. Mainly due to ethical and economical reason in vivo experiment should be preceded by alternative tests including computational and biomimetic ones [16]. It is in line with Green Chemistry principles [17] and with the European Union Directive (Directive 2010/63/EU) which is based on the Three Rs principle of replace, reduce and refine the exploitation of animals for scientific purposes.

2. Materials and Methods

2.1. The analytes

The chemical structures of the investigated oleanane-type triterpenes and their glycosides of plant origin are presented in Table 1.

2.2. Chemicals

The pharmacopoeial standards of arjunic acid, arjunolic acid, arjungenin, arjunglucoside I, sericic acid and arjunetin were purchased from Sigma Aldrich (St. Louis, MO, USA; p.a.). The organic modifiers for micellar mobile phase i.e. acetonitrile and isopropanol as well as surfactants: polyoxyethylene (23) lauryl ether (Brij35) and dodecyl sodium sulfate (SDS), were purchased from Merck (Darmstadt, Germany; p.a.). The buffer components i.e. citric acid and disodium hydrogen phosphate (Na2HPO4) were purchased from Sigma Aldrich (Sigma Aldrich, St. Louis, MO, USA; p.a.). Distilled water was obtained from the Direct-Q3 UV apparatus (Millipore, Warsaw, Poland).

2.3. Chromatographic equipment

The Shimadzu Vp liquid chromatographic system (Shimadzu, Kyoto, Japan) equipped with LC 10AT pump, SPD 10A UV-Vis detector, SCL 10A system controller, CTO-10 AS chromatographic oven and Rheodyne injector valve with a 20 µL loop was applied in the HPLC measurements.

2.4. Chromatographic conditions

The solutions of pharmacopoeial standards of the studied compounds were prepared in methanol (Merck, Darmstadt, Germany; p.a.) at a concentration of 1 mg/mL. All the oleanane-type triterpenes and their glycosides proved to be in the neutral form in solution under experimental conditions. The optimization process of the chromatographic separation was made before the experiment. The flow rate of the mobile phases was established to 1 mL/min and the temperature was set at 25oC. The tested compounds were detected with the UV light at 210 nm.
The C18 encapped column (Purosphere; 125 × 4 mm i.d., 5 µm; Merck, Darmstadt, Germany) was used as the stationary phase while the buffered solutions of both Brij35 as SDS were used as mobile phases. The mobile phases compositions were as follows: Brij35 at the concentrations: 0.04; 0.06; 0.08; 0.10 mol/dm3 (pH 7.4) with the addition of acetonitrile (10% v/v); SDS at the concentrations: 0.06; 0.08; 0.10; 0.12 mol/dm3 (pH 7.4) with an addition of isopropanol as an organic modifier (7% v/v). The buffer was prepared from the solutions of both Na2HPO4 (0.02 mol/dm3) and citric acid (0.01 mol/dm3).
The dead time values were measured from the citric acid peaks. All the reported logarithms of the retention factor were measured three times. The values of peak asymmetry factor were in the acceptable range.

2.5. Pharmacokinetic in silico studies

All the BBB-pharmacokinetic descriptors were calculated using the ACD/Percepta software (version 2012, Advanced Chemistry Development, Inc., Toronto, ON, Canada).

2.6. TLC-based bioatographic assay towards the AChE inhibitory activity

Six standards of oleanane-type triterpenes: arjunic acid, arjunolic acid, arjungenin, arjunglucoside I, sericic acid, and arjunetin - purchased in Sigma Aldrich (St. Louis, MO, USA) were prepared at the concentration of 1 mg/mL in double-distilled water : methanol (50:50 v/v) and they were applied separately at the surface of the aluminium normal phase 10 cm x 10 cm TLC plate (Silica gel 60 F254, Merck, Darmstadt, Germany) with an autosampler (Camag, Muttenz, Switzerland) as the 6 mm zones, distant from one another of 1.5 cm horizontally and 2 cm vertically. Every reference solution was applied as 4, 6, 8 and 10 µL volume bands on three TLC plates.
The TLC plates were later subjected to the enzymatic assay according to the previously published protocol [18], with some modifications.
As the TLC plate was not developed in a TLC solvent system, but used directly in the TLC bioautographic assay, the authors modified the previously published protocol and sprayed the TLC with the substrate (2-naphtyl acetate) dissolved in distilled water at the quantity of 30 mg/20 mL. The TLC was dried in cold air and later the solution of AChE enzyme (AChE from electric eel type VI-S, Sigma Aldrich, St. Louis, CA, USA) dissolved in the aqueous solution of trizma buffer (pH 7.8) with bovine serum (500 mg/100 mL, Sigma Aldrich) at the quantity of 3 U/mL was sprayed on the TLC plate and incubated at the temperature of 37◦C for the following 20 min in a humid incubator. In the next step, the Fast Blue B solution (0.615 mg/mL) was sprayed on the plate and visualized active zones as white spots against violet background. The area of the discoloured zones was corresponding to the inhibitory strength of respective zones.
In the end the TLC plate was dried in the air, and analysed by Camag TLC visualizer at visible light. The peak areas of the discoloured zones were automatically calculated by the WinCats program (v. 1.4, Camag) and their size was compared to calculate the IC50 values that corresponded to the concentration of the standard giving half maximum inhibition of AChE enzyme.

2.7. Molecular docking procedure

The ligand molecules were obtained by the online SMILES translator (cactus.nci.nih.gov/translate) and subsequently optimized by using Avogadro 1.1.1 [19] and the UFF force field [20] (5000 steps, steepest descent algorithm). Flexible and optimized ligand molecules were docked into the binding pocket of the protein structure found in the PDB database (PDB:1EVE). Docking simulations were carried out in the AutoDockVina software [21]. The procedure was performed within the cuboid region of dimensions of 22 × 30 × 34 Å3 which covers the co-crystallized ligand present in the considered PDB record as well as the closest amino-acid residues that exhibit contact with this ligand. All the default procedures and algorithms implemented in AutoDockVina were applied during docking procedure. The rotatable torsional angles in both ligand molecules and the selected amino-acid sidechains within the binding cavity (Tyr334, Phe330, Phe75, Trp84, Glu199, Ser200, Tyr70, Tyr121, Trp279, Phe290, Phe331, Phe288, His440, Gln74, Leu282, Trp432, Asn85 and Asp285) were allowed to rotate. The visual inspections of each pose of the docked ligands were carried out in order to assure that the binding energies correspond to the structurally-analogous orientations. The procedure was validated in our previous work [22].

2.8. Kinetics of AChE inhibition

The samples of the most potent AChE inhibitor i.e., ARG, were prepared in 12 dilutions in the concentration range of 0.00045 - 0.92 mM in dimethyl sulfoxide (DMSO≥99.7 %; Sigma Aldrich). The Ellman’s colorimetric method [23] with some modifications [24] was applied. Each of the tested ARG samples (15 μL) was mixed with 20 μL of the AChE solution (from electric eel, Type VI-S; Sigma Aldrich; 0.28 U/mL) and after 5 minutes completed with 35 μL of acetylthiocholine iodide (ATChI; Sigma Aldrich; 1.5 mmol/L), 175 μL of 0.3 mmol/L 5,5′-dithiobis-2-nitrobenzoic acid (DTNB containing 10 mmol/L NaCl and 2 mmol/L MgCl2; Sigma Aldrich) and 100 μL with the Tris-HCl buffer (50 mmol/L, pH 8.0). The AChE, ATChI, and DTNB solutions were prepared in the Tris-HCl buffer. In order to eliminate the absorbance increase due to the spontaneous hydrolysis of the substrate, there were used “blank” samples composed of 15 μL of Tris-HCl buffer instead of ARG as well as of the above-mentioned compounds. The absorbance of the test samples was measured every minute for 32 minutes and it was subtracted from the absorbance of the “blank” sample. The background samples were prepared with 15 μL of each ARG solution and 330 μL of Tris-HCl buffer. The samples were incubated at room temperature for 30 minutes. The absorbance was measured at 412 nm (96-well microplate reader, Tecan Sunrise, Grödig, Austria). Each sample was analyzed in three repetitions. The linear regression analysis was conducted using the Minitab 18 Statistical Software (Minitab Inc., State College, PA, USA) and the values of the correlation coefficients, slopes, intercepts, and the standard errors were obtained.

2.9. Toxicity assay

To assess ARG toxic effect the ECOSAR (v. 1.11) free software was employed. Based on the ARG chemical structure both acute and chronic toxicity endpoints for fish, aquatic invertebrates (Daphnia), and green algae were measured.

3. Results

3.1. The BBB-pharmacokinetic in silico studies

The BBB pharmacokinetic descriptors were determined in silico using the ACD/Percepta software. The following parameters were calculated: logBB – the distribution of a substance in the blood-brain area (the BBB penetration descriptor), logPS – the rate of passive diffusion/permeability (the permeability-surface area product), logPS,Fu,brain – the brain/plasma equilibration rate, Fu – the fraction unbound in plasma and Fb – the fraction unbound in brain (Table 2).

3.2. The BBB-biomimetic studies

To determine the BBB permeability of the tested oleanane-type triterpenes and their glycosides, micellar chromatographic systems recognized as biomimetic ones were applied. For this purpose, both the Biopartitioning Micellar Chromatography (BMC) using non-ionic surfactant polyoxyethylene (23) lauryl ether (Brij35), and SDS-Micellar Chromatography using the anionic dodecyl sodium sulfate (SDS) were applied.
The relationship between the surfactant concentration in the effluent and the retention of analytes is described by Foley’s equation [25]:
1/k = (KMA/km)CM + 1/km
where k is the retention factor, CM is the total surfactant concentration in the mobile phase minus the critical micellization concentration (CMC), KMA is the analyte-micelle association constant, and km is the micellar retention factor at zero micelles concentration in the mobile phase which corresponds to the monomer surfactant concentration equal to CMC. These parameters describe in a simply way possible interactions in the micellar system that mimic biological environment. In this case, both Brij35 and SDS micelles can be treated like a simple BBB model.
To evaluate the KMA and km values, the relationship between the experimental 1/k and CM values should be plotted and then calculated from the slope and intercept of the plot.Very good linear relationships were obtained with the average values of R2 equal to 0.96 for BMC and 0.99 for SDS system. In Figure 1 1/k vs. CM relationships were presented.
Based on the Foley’s model, both log KMA and log km were calculated. These parameters are considered as lipophilicity descriptors due to the affinity to the surfactant-modified stationary phase (km) as well as binding to the micelles (KMA) [26]. The logarithm of the micellar retention factor, log km, is analogue to the logarithm of retention factor extrapolated to pure water (log kw) obtained in RP-LC system with the water-organic mobile phase. In this study, the log (km/KMA) values calculated from the slopes of Eq. 1 were taken as the micellar lipophilicity descriptors whereas the logarithm of the analyte–micelle association constant (logKMA) values obtained from micellar systems can be taken as an estimate of logBB values. The calculated logBB values based on the logKMA obtained from both BMC and SDS systems (logBB-BMC and logBB-SDS, respectively) as well as logBB in silico values are presented in Figure 2.

3.3. AChE inhibitory activity of the selected saponins in the TLC-bioautography assay

The selected assay is used to search for AChE inhibitory properties of single components or ingredients of mixtures that were introduced on a TLC plate. The TLC-bioautography assay was performed on a series of four dilutions of six reference solutions of arjunic acid, arjunolic acid, arjungenin, arjunglucoside I, sericic acid, and arjunetin provided evidence for the AChE – inhibitory properties of all selected triterpenes (Figure 3) that was dependant on the introduced concentration. As presented in Figure 3, the compounds were characterized by a similar inhibitory potential, in alike application volumes.
The imaging program (WinCats, Camag) enabled a relative quantitative analysis of the inhibition zones. As a result, from the zones of inhibition, peak areas were obtained. The transformation was necessary to calculate the IC50 values of every tested compound to compare their inhibitory potential towards the AChE enzyme (Table 3).

3.4. Molecular docking studies

All compounds exhibit high and relatively similar magnitude of binding energies, varying from ca. -15 to -13 kcal/mol. All energy values are negative which clearly speaks for strongly favorable binding in all considered cases. Contrary to our previous results, we did not observe any statistically-significant correlation between either experimentally-determined IC50 values or theoretically-predicted binding energies and molecular dimensions of the studied molecules (i.e. molecular volume and molecular surface area determined by using the 3vee.molmovdb.org online server with probe of 0.1 nm and high grid resolution; Figure 4 (C)). Such observation suggests that the intensity of binding to AChE is governed by those fragments of molecules which are common for all considered compounds. This is in line with experimentally-observed small scatter of the IC50 values.

3.5. The kinetics of AChE inhibition

The absorbance (A) vs. time [min] relationships were plotted (Figure 5) for each ARG concentration. The average value of the correlation coefficient was found to be 0.991. Due to the great linearity, further kinetic studies were carried out. There were calculated the following basic kinetic enzyme parameters: the Michaelis-Menten constant (Km) by means of the Lineweaver–Burk plot (Figure 6) and the maximum reaction velocity (Vmax).

3.6. In silico prediction of acute and chronic toxicity

Acute toxicity (short-term exposure) was assessed using lethal or effect concentration 50 (LC50 and EC50, respectively) whereas chronic toxicity (long-term exposure) was assessed using chronic values (ChV) obtained for fish, Daphnia, and green algae (Table 4).

4. Discussion

Natural products of plant origin can be interesting sources of compounds with neuroprotective properties. Among them the Traditional Chinese Medicine (TCM) herbs play an important role e.g. Ginkgo biloba L. [27], Panax ginseng [28,29,30,31,32], or Scutellaria baicalensis [33,34,35,36,37]. Other plants are also crucial source of a variety of compounds acting on the central nervous system (CNS) e.g. Olea europaea L. [37,38,39], Vitis vinifera L. [37,38,39], Salvia officinalis L. [40,41,42,43,44], Melissa parviflora [45], Berberis integerrima [46], or Carissa edulis [47]. The most important chemical groups of such compounds are saponins [48,49], tannins [50,51], flavonoids [52,53], alkaloids [54], etc. Among the above-mentioned groups of the CNS-active compounds are triterpenes and their glycosides [55,56,57,58]. These compounds can affect the CNS including nerve cells of the brain and spinal cord which control many direct body functions and behaviour. In the context of neuroprotective properties, firstly it is important to confirm the ability of a compound to cross the BBB.
The Terminalia arjuna accumulates bioactive triterpene glycosides (saponins) and aglycones (sapogenins), in a tissue-preferential manner [59]. Many triterpenes demonstrate therapeutic efficacy. In most cases, they can cross the BBB and may affect the CNS including nerve cells of the brain and spinal cord which control many direct body functions and behaviours. They may also affect the autonomic nervous system which includes the regulation of internal organs, heartbeat, circulation and breathing.
Oleanane triterpenoids/saponins (derived from β-amyrin) have also been reported to have mainly cardioprotective potential [59,60,61,62]. Moreover, numerous studies have confirmed their antioxidant [63], antimicrobial [64,65], anti-inflammatory [66], anticancer [67], precognitive [12], hepatoprotective [68], and others, activities.
It should be strongly emphasized that only the drug fraction unbound in media such as plasma can be transferred into body tissues. In vitro methods including ultrafiltration or equilibrium dialysis are most often used to measure the fraction unbound value of a drug. These in vitro obtained values are used not only for measurement of transfer rate into body tissues but also of the BBB permeability [69]. It should be remembered that research on the penetration of compounds through the biological barriers, including the BBB one, is carried out using in vivo methods in particular. However, for ethical and economic reasons, the need to use alternative methods to in vivo one, including the non-cell based-in vitro (biomimetic) and/or in silico (computational) has been emphasized in recent years.
Both biomimetic and computational BBB-pharmacokinetic studies are commonly used in the laboratory practice at the first stages of an experiment on biologically active compounds (potential drugs) and constitute an important stage of research in the drug design process. At the stage of the in silico studies the most important BBB-pharmacokinetic descriptors were calculated, i.e. the distribution of a substance in the blood-brain area, the rate of passive diffusion/permeability, the brain/plasma equilibration rate, the fraction unbound in plasma and the fraction unbound in brain. The blood-brain distribution (BB), frequently expressed as logBB, is defined as a ratio between the concentration in the brain and the concentration in the blood [70,71]. This experiment first identified 2 out of 6 tested compounds, i.e. arjunetin and arjunglucoside I, capable of crossing the BBB. However, it is commonly recognized that the most important parameter of the permeability through the BBB is the permeability – surface area product (PS) often expressed as logPS. These index is closely related to the cerebral blood flow (CBF) which is measured using various invasive as well as non-invasive techniques i.e. direct intravascular measurements, nuclear medicine, X-ray imaging, magnetic resonance imaging, ultrasound techniques, thermal diffusion, and optical methods. The most invasive methods require surgical access, arterial puncture, or catheterization while less invasive methods demand the intravenous injection of a contrast agent [72]. The CBF is a very important parameter for brain viability and its functions because it ensures proper delivery of oxygen which is necessary for the neuronal oxidative metabolism of energy substrates. It is defined as the blood volume that flows per unit mass per unit time in brain tissue and is typically expressed in units of mL blood∕(100 gtissue*min), or mL blood (100 mLtissue*min) [72] or in mL blood/(h*kg) [73]. Taking into account PS values calculated in this experiment, arjunetin exhibited the highest BBB-permeability potential, followed by arjungenin, arjunglucoside I and sericic acid (ex aequo) whereas both acids: arjunic and arjunolic one exhibited the lowest BBB permeability.
The scientific reports indicate that the time to reach brain equilibrium can be prolonged when the BBB permeability–surface area product (PS) or the fraction unbound in the brain decreases [74], therefore it can be noticed that the lower values of the PS or Fb, the longer the time to reach brain equilibrium is required [73]. In our experiment, no significant differences between Fb values were observed whereas the differences between the PS values are much greater (from 0.63 mL*h-1*kg-1 for both arjunic and arjunolic acids to 5 mL*h-1*kg-1 for arjunetin). Then, the longest time to reach brain equilibrium can be observed for the above-mentioned acids and the shortest for arjunetin. A high rate of penetration results from high BBB permeability as well as low brain tissue binding [73].
In addition, analyzing values from Table 2, it can be also seen that arjunetin and arjunglucoside I bind the least to blood plasma proteins (the highest value of free drug concentration, Fu) and shows the highest log BB value (0.73 and 0.12, respectively). The rest of the compounds have logBB values less than zero. Therefore, it can be presumed that among the tested compounds, arjunetin and arjunglucoside I are the substances that can penetrate the BBB to the greatest extent. However, the frequently used parameter for assessing the extent of the CNS distribution is also the ratio of brain/plasma partition coefficient, Kp,brain. This parameter – calculated for compounds that distribute solely by passive diffusion – is a function of the relative plasma and brain tissue unbound fractions at distribution equilibrium [74]. In our case, most substances i.e. arjunic acid, arjunolic acid, arjungenin, and sericic acid, have Kp,brain values less than 1 which can result from more extensive binding to proteins in plasma than those in brain tissues. Other explanation can be significant impairment the in CNS distribution such as the efflux transport at the BBB [74]. However, taking into account logBB values (Table 2) it can be assumed that these compounds have simply lower CNS-distribution potential contrary to arjunetin and arjunglucoside I with Kp,brain values 5.1 and 1.25, respectively.
There exists the free drug theory that postulates that all the distribution processes of the active substance within biological barriers depend on the unbound drug concentration [69,75]. It must be emphasized that drug in the blood is present both in unbound and bound form to plasma proteins and erythrocytes. In our experiment, two substances i.e. arjunetin and arjunglucoside I have the highest value of the fraction unbound in plasma (0.051 and 0.050, respectively) in contrast to other compounds with Fb values in the range of 0.012 to 0.016. This could confirm earlier suppositions that arjunetin has the greatest ability among the tested compounds to cross the blood-brain barrier. Nevertheless, it is also hypothesized that drugs binding to protein can rapidly dissociate and permeate in vivo through the BBB into the brain tissues [69]. Therefore there may exist some differences between drug concentration obtained in vivo in brain and that estimated in vitro based on the free drug concentration. Nevertheless, the ability of most drugs to cross the BBB is nowadays estimated using the free drug fraction theory with reasonably acceptable results [69].
The biomimetic studies were carried out to confirm (or not) the previously assumptions made based on the BBB-pharmacokinetic computational research. For this purpose, micellar liquid chromatography, using non-ionic Brij35(this type of chromatography is called the BMC) as well as anionic SDS surfactants, was applied. These methods are commonly used to assess the permeation of a substance through biological barriers [76,77,78,79,80]. The concentration of a surfactant in a micellar mobile phase must be above the critical micellization concentration (cmc) whereas the commonly used stationary phase is the octadecyl-modified silica gel [81,82,83]. Due to the wide application of micellar chromatography in the study of the penetration of compounds through biological barriers, it is a recognized technique in biomimetic studies on biologically active compounds.
Since the Brij35 micelle is assumed to be a kind of simple, chemical model of the biomembrane, the BMC technique can be useful in describing biological behaviours of different kinds of organic compounds. It can also mimic many biological processes such as BBB penetration, skin permeability, intestinal absorption and drug partitioning process in biological systems [81,82,83], and others. In our research, the logarithms of retention factor extrapolated to pure water (log km), for both the BMC and SDS systems, have been determined. This parameter is recognized to be alternative to the logarithm of n-octanol/water partition coefficient (logPo/w) lipophilicity descriptor.
In the research, each system was previously optimized by selecting the appropriate concentrations of surfactants, selecting the organic modifier and its concentration in the mobile phase. The surfactant solutions were buffered (pH 7.4). Moreover, according to the Foley’s equation [25], the interactions performed in the micellar systems have been characterized. Knowledge of the type of interactions between the analyte and the micelle, which in this case is a BBB model, can provide valuable information on the mechanism of interaction between a substance and a barrier. For this purpose, important physicochemical parameters such as KMA – the analyte-micelle association constant and PSW – the partition coefficient of an analyte between the stationary phase and water were calculated. Based on the above–mentioned parameters, one can conclude about the strength of analyte interaction with the biological membrane. Such studies can be very essential in the context of research on the biological activity of the tested compounds.
In the previous study [84], it was proved that the logarithm of the analyte–micelle association constant (logKMA) can characterize directly the passage of substances through the BBB comparable to the logBB pharmacokinetic parameter. Since the Foley’s model describes the retention behaviors of the analyte in the micellar system, which can be treated as a simple BBB model, therefore the parameters contained in it can characterize the biodistribution of the analyte in the BBB area. Very good linear relationships (R2 > 0.9) between 1/k and CM were obtained for all tested compounds (Figure 1), confirming that the Foley’s equation correctly describes the retention of solutes in the tested BMC and SDS chromatographic systems. LogKMA can be a useful tool for rapid assessment of the ability of a substance to cross the BBB, especially in the early stage of research. The obtained logKMA-BMC values confirmed that both compounds: arjunetin and arjunglucoside I interact the most with Brij micelles which is recognized as a simply biological membrane model. In the SDS system, the matter is more complicated. Due to probably electrostatic interactions between the analytes and anionic micelles as well as strong retention of compounds, the intercepts for 3 out of 6 equations are negative. Unfortunately, the intercepts less than zero have no physico-chemical sense because they are equal to reciprocal of km parameter being the retention factor in the system in which the concentration of free surfactant (CM) in the effluent is equal to zero. However, to eliminate the impact of possible electrostatic interactions, the log (km/KMA) values calculated from the slopes of Eq. 1, were taken into account. These values have been treated as logBB values (logBB-SDS; see Figure 2).
As shown in Figure 2, there are no significant differences between logBB values obtained using computational and biomimetic methods. The BBB-pharmacokinetic biomimetic studies confirmed that arjunetin and arjunglucoside I can cross the BBB and have therefore the greatest BBB-penetration potential among the tested compounds while arjunic and arjunolic acids have the smallest one.
Analyzing the IC50 values obtained in the TLC-bioautography assay towards the AChE inhibition, it can be stated that among the tested compounds arjunolic acid was found to be the strongest inhibitor, whereas arjunetin was the weakest one among the tested compounds. However, the differences in the obtained IC50 values are insignificant. Thus, it can be concluded that all compounds have a very similar affinity to AChE which was later confirmed by molecular docking studies.
Taking into account the biological potential of other metabolites of plant origin, the compounds tested in the study exhibit relatively strong inhibitory potential. Previous results on triterpenoids confirmed their AChE inhibitory potential. In the study on the metabolites from Centella asiatica, asiatic acid was found to be the strongest AChE inhibitor with the IC50 value of 15.05 ± 0.05 µM [85]. Also, the metabolites of Garcinia hombroniana delivered information on a high inhibitory potential of 2β-Hydroxy-3α-O-caffeoyltaraxar-14-en28-oic acid present in the plant [86]. In comparison to Amaryllidaceae alkaloids, like galanthamine that is registered as first line drug in the treatment of AD that was characterized by the IC50 value of 3.520 µM [87], the tested compounds seem to be promising.
As the ARG exhibits the lowest IC50 value [mM] among the tested compounds, it was applied in the AChE inhibition kinetic studies using the colorimentric Ellman’s test [23]. The Michaelis-Menten constant (Km) calculated based on the course of the curve (Figure 6) was found to be 0.000011 mol/L. The obtained relationships show that the rate of substrate-enzyme binding is concentration-dependent and reaches the maximum velocity equal 2.2 x 10-5.
The results of the docking study have been analyzed with respect to the mechanistic interaction pattern that may be significant in the context of interpretation of the obtained binding energies and recognizing the role of pharmacophore fragment. The summary given below relies on analyzing the ligand-protein contacts that take place if the distance between any corresponding atom pair is smaller than the arbitrarily accepted value of 0.4 nm. Figure 3 (A) shows the superposition of all most favorable ligand poses whereas Figure 3 (B) shows the most essential ligand-protein interactions. The extremely close match between superposed structures of all compounds agrees with the previous statement claiming that binding strength is determined by the pharmacophore-like, common molecular fragment of all compounds.
The detailed pattern of ligand-enzyme interactions is illustrated in Figure 3 (B), on the example of arjunglucoside I (i.e. the compound displaying the lowest IC50 value [mM]). However, due to similar orientations of all ligands in the binding cavity, the majority of conclusions can be transferred to remaining compounds. All ligands prefer roughly the same binding position in the enzyme cavity which enables them to block the catalytic site (the proximity to the catalytic histidine, His440, can be observed). The central fragment of ligand molecule (composed of aliphatic, cyclic moieties) interacts with aromatic cluster of sidechains, created by His440, Phe290, Trp84, Trp279, Phe288, Phe331, Tyr121, Tyr70, and Tyr334. Such contacts have a character of the CH-π interactions, supported (in some cases, e.g. His44) by hydrogen bonding with the neighboring fragments of ligand. The hydroxyl groups located at the edge of aliphatic, condensed fragment of the ligand, interact with Arg289 and Ser286. Both these contacts occur via hydrogen bonding and, surprisingly, involve backbone fragments of the protein (ligand can only be a hydrogen bonding donor). One can speculate about an analogous interaction in the case of Ile286 (also a backbone fragment) but, due to the lack of rotation around peptide bonds in the docking procedure, this was not explicitly observed. Interestingly, the ligand contacts with non-aromatic, hydrophobic sidechains are marginal and include only Ile287 and Leu127. Even in these cases, such proximitiesare rather an opportunistic consequence of much stronger interactions occurring with other, adjacent amino-acid residues.
The moiety of type and character varying between molecules (topologically-equivalent to the glucopyranose residue in the case of arjunglucoside I, illustrated in Figure 3) is located close to a set of polar amino-acid residues, including Asn85, Ser122, Gln69, and Asp72. The dominating character of involved interactions is hydrogen bonding, where the considered fragment of ligand molecule can play a role of both donor and acceptor. In spite of the presence of tryptophan and tyrosine sidechains in the close proximity of glucopyranosidic moiety, no CH-π stacking, characteristic for carbohydrate-protein binding, was observed. This may explain why this fragment of ligand molecule (or its lack) is not particularly crucial for binding strength; the hydrogen bond donors and acceptors present in this region of cavity can equally well be saturated by water molecules, providing roughly the same balance of energy.
Taking into account toxicity predicted values, it can be stated that ARG has a low potential for chronic and acute toxicity on fish, daphnia and green algae.

5. Conclusions

The results of the presented studies showed that naturally occurring oleanane-type triterpenes and their glycosides may by active on the AChE enzyme. Therefore these compounds, especially arjunetin and arjunglucoside I, can be novel drug candidates in the treatment of neurodegenerative diseases with memory impairment including the AD and can be an interesting source of further, deeper research in this aspect.

Author Contributions

Conceptualisation, K.S.; methodology, K.S., W.K.-K. W.P.; validation, K.S., W.K.-K., W.P.; formal analysis, K.S., W.K.-K., W.P.; investigation, K.S., W.K.-K., W.P.; resources, K.S., W.K.-K.; data curation, K.S.; writing—original draft preparation, K.S.; supervision, K.S.; project administration, K.S.; funding acquisition, K.S. and W.K.-K.. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Union of Lublin Universities in the framework of the program “Staż za miedzą” No. SzM/024/2022/I-N.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data sharing is not applicable to this article.

Conflicts of Interest

The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

References

  1. The WHO brochure “Global action plan on the public health response to dementia 2017 – 2025”, 2017; ISBN: 978-92-4-151348-7. Available online: https://www.who.int/publications/i/item/global-action-plan-on-the-public-health-response-to-dementia-2017---2025 (accessed on 3 October 2023).
  2. Wiendl, H.; Meuth, S.G. Pharmacological Approaches to Delaying Disability Progression in Patients with Multiple Sclerosis. Drugs 2015, 75, 947–977. [Google Scholar] [CrossRef] [PubMed]
  3. Zuin, M.; Cherubini, A.; Volpato, S.; et al. Acetyl-cholinesterase-inhibitors slow cognitive decline and decrease overall mortality in older patients with dementia. Sci Rep. 2022, 16, 12214. [Google Scholar] [CrossRef] [PubMed]
  4. Smyrska-Wieleba, N.; Mroczek, T. Natural Inhibitors of Cholinesterases: Chemistry, Structure–Activity and Methods of Their Analysis. Int. J. Mol. Sci. 2023, 24, 2722. [Google Scholar] [CrossRef]
  5. Coelho dos Santos, T.; Mota Gomes, T.; Araújo Serra Pinto, B.; et al. Naturally Occurring Acetylcholinesterase Inhibitors and Their Potential Use for Alzheimer's Disease Therapy. Front. Pharmacol. 2018, 9, 1192. [Google Scholar] [CrossRef] [PubMed]
  6. Berkov, S.; Atanasova, M.; Georgiev, B.; et al. The Amaryllidaceae alkaloids: An untapped source of acetylcholinesterase inhibitors. Phytochem. Rev. 2022, 21, 1415–1443. [Google Scholar] [CrossRef]
  7. Bora, K.S.; Sharma, A. Phytochemical and pharmacological potential of Medicago sativa: A review. Pharm. Biol. 2011, 49, 211–220. [Google Scholar] [CrossRef]
  8. Liu, X.-G.; Sun, Y.-Q.; Bian, J.; et al. Neuroprotective effects of triterpenoid saponins from Medicago sativa L. against H2O2-induced oxidative stress in SH-SY5Y cells. Bioorg. Chem. 2019, 83, 468–476. [Google Scholar] [CrossRef]
  9. Li, Z.; Liu, Y.; Wang, L.; Liu, X.; Chang, Q.; Guo, Z.; Liao, Y.; Pan, R.; Fan, T.P. Memory-Enhancing Effects of the Crude Extract of Polygala tenuifolia on Aged Mice. Evid Based Complement Alternat Med. 2014, 2014, 392324 htpps:. [Google Scholar] [CrossRef]
  10. Li, C.; Yang, J.; Yu, S.; et al. Triterpenoid Saponins with Neuroprotective Effects from the Roots of Polygala tenuifolia. Planta Med. 2008, 74, 133–141. [Google Scholar] [CrossRef]
  11. Son, I.H.; Park, Y.H.; Lee, S.I.; et al. Neuroprotective Activity of Triterpenoid Saponins from Platycodi radix Against Glutamate-induced Toxicity in Primary Cultured Rat Cortical Cells. Molecules 2007, 12, 1147–1152. [Google Scholar] [CrossRef]
  12. Yaidikar, L.; Thakur, S. Arjunolic acid, a pentacyclic triterpenoidal saponin of Terminalia arjuna bark protects neurons from oxidative stress associated damage in focal cerebral ischemia and reperfusion. Pharmacol. Rep. 2015, 67, 890–895. [Google Scholar] [CrossRef] [PubMed]
  13. Toropov, A.A.; Toropova, A.P.; Beeg, M.; et al. QSAR model for blood-brain barrier permeation. J. Pharmacol. Toxicol. Meth. 2017, 88, 7–18. [Google Scholar] [CrossRef] [PubMed]
  14. Van Bree, J.B.M.M.; De Boer, A.G.; Danhof, M.; et al. Drug transport across the blood-brain barrier. Pharm. Weekbl. 1992, 14, 305–310. [Google Scholar] [CrossRef]
  15. Iyer, M.; Mishra, R.; Han, Y.; et al. Predicting Blood–Brain Barrier Partitioning of Organic Molecules Using Membrane-Interaction QSAR Analysis. Pharm. Res. 2002, 19, 1611–1621. [Google Scholar] [CrossRef]
  16. Ekins, S.; Waller, C.L.; Swaan, P.W.; et al. Progress in predicting human ADME parameters in silico. J. Pharmacol. Toxicol. Meth 2000, 44, 251–272. [Google Scholar] [CrossRef]
  17. Anastas, P.; Eghbali, N. Green chemistry: Principles and practice. Chem. Soc. Rev. 2010, 39, 301–312. [Google Scholar] [CrossRef]
  18. Kukula-Koch, W.; Mroczek, T. Application of hydrostatic CCC-TLC-HPLC-ESI-TOF-MS for the bioguided fractionation of anticholinesterase alkaloids from Argemone mexicana L. roots. Anal. Bioanal. Chem. 2015, 407, 2581–2589. [Google Scholar] [CrossRef]
  19. Hanwell, M.D.; Curtis, D.E.; Lonie, D.C.; et al. Avogadro: An advanced semantic chemical editor, visualization, and analysis platform. J Cheminform. 2012, 13, 17 htpps. [Google Scholar] [CrossRef]
  20. Rappe, A.K.; Casewit, C.J.; Colwell, K.S.; et al. UFF, a full periodic table force field for molecular mechanics and molecular dynamics simulations. J. Am. Chem. Soc. 1992, 114, 10024–10035. [Google Scholar] [CrossRef]
  21. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2010, 30, 455–461. [Google Scholar] [CrossRef]
  22. Tarabasz, D.; Szczeblewski, P.; Laskowski, T.; et al. The Distribution of Glucosinolates in Different Phenotypes of Lepidium peruvianum and Their Role as Acetyl- and Butyrylcholinesterase Inhibitors-In Silico and In Vitro Studies. Int. J. Mol. Sci. 2022, 27, 4858. [Google Scholar] [CrossRef]
  23. Ellman, G.L.; Lourtney, D.K.; Andres, V.; Gmelin, G. A new and rapid colorimetric determination of acetylcholinesterase activity. Biochem. Pharmacol. 1961, 7, 88–95. [Google Scholar] [CrossRef] [PubMed]
  24. Szwajgier, D.; Baranowska-Wójcik, E. Terpenes and Phenylpropanoids as Acetyl- and Butyrylcholinesterase Inhibitors: A Comparative Study. Curr Alzheimer Res. 2019, 16, 963–973. [Google Scholar] [CrossRef] [PubMed]
  25. Foley, J.P. Critical compilation of solute-micelle binding constants and related parameters from micellar liquid chromatographic measurements. Anal. Chim. Acta 1990, 231, 237–247. [Google Scholar] [CrossRef]
  26. Janicka, M.; Śliwińska, A.; Sztanke, M.; et al. Combined Micellar Liquid Chromatography Technique and QSARs Modeling in Predicting the Blood–Brain Barrier Permeation of Heterocyclic Drug-like Compounds. Int. J. Mol. Sci. 2022, 23, 15887. [Google Scholar] [CrossRef] [PubMed]
  27. Birks, J.; Grimley, E.J. Ginkgo biloba for cognitive impairment and dementia. Cochrane Database Syst. Rev. 2009, 21, CD003120. [Google Scholar] [CrossRef] [PubMed]
  28. Joo, S.S.; Lee, D.I. Potential effects of microglial activation induced by ginsenoside Rg3 in rat primary culture: Enhancement of type a macrophage scavenger receptor expression. Arch. Pharm. Res. 2005, 28, 1164–1169. [Google Scholar] [CrossRef]
  29. Chen, F.; Eckman, E.A.; Eckman, C.B. Reductions in levels of the Alzheimer's amyloid β peptide after oral administration of ginsenosides. FASEB J. 2006, 20, 1269–1271. [Google Scholar] [CrossRef]
  30. Li, N.; Liu, B.; Dluzen, D.E.; Jin, Y. Protective effects of ginsenoside Rg2 against glutamateinduced neurotoxicity in PC12 cells. J. Ethnopharmacol. 2007, 111, 458–463. [Google Scholar] [CrossRef]
  31. Joo, S.S.; Yoo, Y.M.; Ahn, B.W.; et al. Prevention of inflammation-mediated neurotoxicity by Rg3 and its role in microglial activation. Biol. Pharm. Bull. 2008, 31, 1392–1396. [Google Scholar] [CrossRef]
  32. Shieh, P.C.; Tsao, C.W.; Li, J.S.; et al. Role of pituitary adenylate cyclase-activating polypeptide (PACAP) in the action of ginsenoside Rh2 against betaamyloid-induced inhibition of rat brain astrocytes. Neurosci. Lett. 2008, 434, 1–5. [Google Scholar] [CrossRef] [PubMed]
  33. Wakabayashi, I.; Yasui, K. Wogonin inhibits inducible prostaglandin E2 production in macrophages. Eur. J. Pharmacol. 2000, 406, 477–481. [Google Scholar] [CrossRef]
  34. Park, B.K.; Heo, M.Y.; Park, H.; et al. Inhibition of TPA-induced cyclooxygenase-2 expression and skin inflammation in mice by wogonin, a plant flavone from Scutellaria radix. Eur. J. Pharmacol. 2001, 425, 153–157. [Google Scholar] [CrossRef] [PubMed]
  35. Nakamura, N.; Hayasaka, S.; Zhang, X.Y.; et al. Effects of baicalein, baicalin, and wogonin on interleukin-6 and interleukin-8 expression, and nuclear factor-kb binding activities induced by interleukin-1beta in human retinal pigment epithelial cell line. Exp. Eye Res. 2003, 77, 195–202. [Google Scholar] [CrossRef] [PubMed]
  36. Suk, K.; Lee, H.; Kang, S.S.; et al. Flavonoid baicalein attenuates activationinduced cell death of brain microglia. J. Pharmacol. Exp. Ther. 2003, 305, 638–645. [Google Scholar] [CrossRef] [PubMed]
  37. La Vecchia, C.; Bosetti, C. Diet and cancer risk in Mediterranean countries. Hungar. Med. J. 2007, 1, 13–23. [Google Scholar] [CrossRef]
  38. Trichopoulou, A.; Lagiou, P.; Kuper, H.; et al. Cancer and Mediterranean dietary traditions. Cancer Epidemiol. Biomarkers Prev. 2000, 9, 869–873. [Google Scholar]
  39. Visioli, F.; Grande, S.; Bogani, P.; et al. The role of antioxidant in the Mediterranean diets: Focus on cancer. Eur. J. Cancer Prev. 2004, 13, 337–343. [Google Scholar] [CrossRef]
  40. Ninfali, P.; Mea, G.; Giorgini, S.; et al. Antioxidant capacity of vegetables, spices, and dressings relevant to nutrition. Brit. J. Nutr. 2005, 93, 257–266. [Google Scholar] [CrossRef]
  41. dos Santos-Neto, L.L.; de Vilhena Toledo, M.A.; Medeiros-Souza, P.; et al. The use of herbal medicine in Alzheimer’s disease--A systematic review. Evid.-Based Complement. Altern. Med. 2006, 3, 441–445. [Google Scholar] [CrossRef]
  42. Bozin, B.; Mika-Dukic, N.; Samojlik, I.; et al. Antimicrobial and antioxidant properties of rosemary and sage (Rosmarinus officinalis L. and Salvia officinalis L., Lamiaceae) essential oils. J. Agric. Food. Chem. 2007, 55, 7879–7885. [Google Scholar] [CrossRef] [PubMed]
  43. Cheung, S.; Tai, J. Anti-proliferative and antioxidant properties of rosemary Rosmarinus officinalis. Oncol. Rep. 2007, 17, 1525–1531. [Google Scholar] [CrossRef] [PubMed]
  44. Imanshahidi, M.; Hosseinzadeh, H. The pharmacological effects of Salvia species on the central nervous system. Phytother. Res. 2006, 20, 427–437. [Google Scholar] [CrossRef] [PubMed]
  45. Bhat, J.U.; Nizami, Q.; Asiaf, A.; et al. Anticonvulsant activity of methanolic and aqueous extracts of Melissa parviflora in experimentally induced Swiss albino mice. EXCLI J. 2012, 11, 1–6. [Google Scholar] [PubMed]
  46. Hosseinzadeh, H.; Ramezani, M.; Shafaei, H.; et al. Anticonvulsant effect of Berberis integerrima L. root extracts in mice. J. Acupunct. Meridian Stud. 2013, 6, 12–17. [Google Scholar] [CrossRef] [PubMed]
  47. Ya'u, J.; Yaro, A.H.; Malami, S.; et al. Anticonvulsant activity of aqueous fraction of Carissa edulis root bark. Pharm. Biol. 2015, 53, 1329–1338. [Google Scholar] [CrossRef] [PubMed]
  48. Sankari, M.; Chitra, V.; Silambujanaki, P.; et al. Anticonvulsant activity of ethanolic extract of Aegle marmelos (leaves) in mice. Int. J. Pharmtech Res. 2010, 2, 640–643. [Google Scholar]
  49. Anaka, O.N.; Ozolua, R.I.; Ikpefan, E.O.; et al. Anticonvulsant activity of the aqueous extract of Allium cepa L.(Amaryllidaceae) in rats and mice. J. Pharm. Biores. 2014, 11, 1–7. [Google Scholar] [CrossRef]
  50. Showraki, A.; Emamghoreishi, M.; Oftadegan, S. Anticonvulsant Effect of the Aqueous Extract and Essential Oil of Carum Carvi L. Seeds in a Pentylenetetrazol Model of Seizure in Mice. Iran. J. Med. Sci. 2016, 41, 200–208. [Google Scholar]
  51. Amabeoku, G.J.; Oluchi, N.M.; Davids, T.; et al. Evaluation of the anticonvulsant activity of the leaf methanol extract of Crassula arborescens (Mill.) Willd. (Crassulaceae) in mice. J. Pharm. Pharmacol. 2014, 2, 393–403. [Google Scholar]
  52. Chaulya, N.C.; Haldar, P.K.; Mukherjee, A. Antidiabetic activity of methanol extract of rhizomes of Cyperus tegetum Roxb (Cyperaceae). Acta Pol. Pharm. 2011, 68, 989–992. [Google Scholar] [PubMed]
  53. Chinchawade, A.B.; Deshmukh, D.B.; Gaikwad, D.D.; et al. Anticonvulsant Activity of Chloroform Extract of Bark and Root of Erythrina variegata L. Int. J. Pharm. Clin. Res. 2013, 5, 23–25. [Google Scholar]
  54. Nishanthi, N.; Mohanambal, E.; Devdass, G.; et al. Anticonvulsant Activity of Peperomia tetraphylla (G. Forst., Hook. & Arn.). Int. J. Nov. Trends Pharm. Sci. 2012, 2, 35–38. [Google Scholar]
  55. Barua, C.C.; Begum, S.A.; Barua, A.G.; et al. Anxiolytic and anticonvulsant activity of methanol extract of leaves of Alternanthera brasiliana (L.) Kuntze (Amaranthaceae) in laboratory animals. Indian J. Exp. Biol. 2013, 51, 450–457. [Google Scholar]
  56. Khan, I.; Karim, N.; Ahmad, W.; et al. GABA-A Receptor Modulation and Anticonvulsant, Anxiolytic, and Antidepressant Activities of Constituents from Artemisia indica Linn. Evid.-based Complement. Altern. Med. 2016, 2016, 1215393. [Google Scholar] [CrossRef]
  57. Li, D.-Q.; Zhou, L.; Wang, D.; et al. Neuroprotective oleanane triterpenes from the roots of Bupleurum chinense. Bioorg. Med. Chem. Lett. 2016, 26, 1594–1598. [Google Scholar] [CrossRef]
  58. Gao, Y.; Yan, H.; Jin, R.; et al. Antiepileptic activity of total triterpenes isolated from Poriacocos is mediated by suppression of aspartic and glutamic acids in the brain. Pharm. Biol. 2016, 54, 2528–2535. [Google Scholar] [CrossRef]
  59. Srivastava, G.; Sandeep, G.A.; et al. Transcriptome analysis and functional characterization of oxidosqualene cyclases of the arjuna triterpene saponin pathway. Plant Sci. 2020, 292, 110382. [Google Scholar] [CrossRef]
  60. Mohanty, I.R.; Borde, M.; Kumar, S.; et al. Dipeptidyl peptidase IV Inhibitory activity of Terminalia arjuna attributes to its cardioprotective effects in experimental diabetes: In silico, in vitro and in vivo analyses. Phytomedicine 2019, 57, 158–165. [Google Scholar] [CrossRef]
  61. Pawar, R.S.; Bhutani, K.K. Effect of oleanane triterpenoids from Terminalia arjuna – a cardioprotective drug on the process of respiratory oxyburst. Phytomedicine 2005, 12, 391–393. [Google Scholar] [CrossRef]
  62. Kapoor, D.; Vijayvergiya, R.; Dhawan, V. Terminalia arjuna in coronary artery disease: Ethnopharmacology, pre-clinical, clinical & safety evaluation. J. Ethnopharmacol. 2014, 155, 1029–1045. [Google Scholar] [PubMed]
  63. Pugazhendhi, A.; Shafreen, R.B.; Devi, K.P.; et al. Assessment of antioxidant, anticholinesterase and antiamyloidogenic effect of Terminalia chebula, Terminalia arjuna and its bioactive constituent 7-Methyl gallic acid – An in vitro and in silico studies. J. Mol. Liq. 2018, 257, 69–81. [Google Scholar] [CrossRef]
  64. Gupta, D.; Kumar, M. Evaluation of in vitro antimicrobial potential and GC-MS analysis of Camellia sinensis and Terminalia arjuna. Biotechnol. Rep. (Amst). 2016, 30, 19–25 hptts://. [Google Scholar] [CrossRef]
  65. Mandal, S.; Patra, A.; Samanta, A.; et al. Analysis of phytochemical profile of Terminalia arjuna bark extract with antioxidative and antimicrobial properties. Asian Pac. J. Trop. Biomed. 2013, 12, 960–966. [Google Scholar] [CrossRef] [PubMed]
  66. Dube, N.; Nimgulkar, C.; Bharatraj, D.K. Validation of therapeutic anti-inflammatory potential of Arjuna Ksheera Paka – A traditional Ayurvedic formulation of Terminalia arjuna. J. Tradit. Complement. Med. 2017, 7, 414–420. [Google Scholar] [CrossRef]
  67. Ahmad, M.S.; Ahmad, S.; Gautam, B.; et al. Terminalia arjuna, a herbal remedy against environmental carcinogenicity: An in vitro and in vivo study. Egypt. J. Medical Hum. Genet. 2014, 15, 61–67. [Google Scholar] [CrossRef]
  68. Bhattacharjee, B.; Pal, P.K.; Ghosh, A.K.; et al. Aqueous bark extract of Terminalia arjuna protects against cadmium-induced hepatic and cardiac injuries in male Wistar rats through antioxidative mechanisms. Food Chem. Toxicol. 2019, 124, 249–264. [Google Scholar] [CrossRef]
  69. Tanaka, H.; Mizojiri, K. Drug-protein binding and blood-brain barrier permeability. J Pharmacol. Exp. Ther. 1999, 288, 912–918. [Google Scholar]
  70. Platts, J.A.; Abraham, M.H.; Zhao, Y.H.; et al. Correlation and prediction of a large blood-brain distribution data set--an LFER study. Eur. J. Med. Chem. 2001, 36, 719–730. [Google Scholar] [CrossRef]
  71. Young, R.C.; Mitchell, R.C.; Brown, T.H.; et al. Development of a new physicochemical model for brain penetration and its application to the design of centrallyacting H2 receptor histamine antagonists. J. Med. Chem. 1988, 31, 656–671. [Google Scholar] [CrossRef]
  72. Fantini, S.; Sassaroli, A.; Tgavalekos, K.T.; et al. Cerebral blood flow and autoregulation: Current measurement techniques and prospects for noninvasive optical methods. Neurophotonics 2016, 3, 031411. [Google Scholar] [CrossRef]
  73. Xingrong, L.; Smith, B.J.; Chen, C.; et al. Use of a Physiologically Based Pharmacokinetic Model to Study the Time to Reach Brain Equilibrium: An Experimental Analysis of the Role of Blood-Brain Barrier Permeability, Plasma Protein Binding, and Brain Tissue Binding. J. Pharmacol. Exp. Ther. 2005, 313, 1254–1262. [Google Scholar] [CrossRef]
  74. Kalvass, J.C.; Maurer, T.S.; Pollack, G.M. Use of plasma and brain unbound fractions to assess the extent of brain distribution of 34 drugs: Comparison of unbound concentration ratios to in vivo p-glycoprotein efflux ratios. Drug Metab. Dispos. 2007, 35, 660–666. [Google Scholar] [CrossRef] [PubMed]
  75. Hammarlund-Udenaes, M.; Fridén, M.; Syvänen, S.; et al. On the rate and extent of drug delivery to the brain. Pharm. Res. 2008, 25, 1737–1750. [Google Scholar] [CrossRef] [PubMed]
  76. Escuder-Gilabert, L.; Molero-Monfort, M.; Villanueva-Camañas, R.M.; et al. Potential of biopartitioning micellar chromatography as an in vitro technique for predicting drug penetration across the blood–brain barrier. J. Chromatogr. B. 2004, 807, 193–201. [Google Scholar] [CrossRef]
  77. Martínez-Pla, J.J.; Martín-Biosca, Y.; Sagrado, S.; et al. Evaluation of the pH effect of formulations on the skin permeability of drugs by biopartitioning micellar chromatography. J. Chromatogr. A. 2004, 1047, 255–262. [Google Scholar] [CrossRef]
  78. Hadjmohammadi, M.; Salary, M. Biopartitioning micellar chromatography with sodium dodecyl sulfate as a pseudo α1-acid glycoprotein to the prediction of protein–drug binding. J. Chromatogr. B. 2013, 912, 50–55. [Google Scholar] [CrossRef]
  79. Tsopelas, F.; Danias, P.; Pappa, A.; et al. Biopartitioning micellar chromatography under different conditions: Insight into the retention mechanism and the potential to model biological processes. J. Chromatogr. A. 2020, 1621, 461027. [Google Scholar] [CrossRef]
  80. Dobričić, V.; Nikolic, K.; Vladimirov, S.; et al. Biopartitioning micellar chromatography as a predictive tool for skin and corneal permeability of newly synthesized 17β-carboxamide steroids. Eur. J. Pharm. Sci. 2014, 56, 105–112. [Google Scholar] [CrossRef]
  81. Molero-Monfort, M.; Escuder-Gilabert, L.; Villanueva-Camanas, R.M.; et al. Biopartitioning micellar chromatography: An in vitro technique for predicting human drug absorption. J. Chromatogr. B. 2001, 753, 225–236. [Google Scholar] [CrossRef]
  82. Escuder-Gilabert, L.; Martinez-Pla, J.J.; Sagrado, S.; et al. Biopartitioning micellar separation methods: Modelling drug absorption. J. Chromatogr. B. 2003, 797, 21–35. [Google Scholar] [CrossRef] [PubMed]
  83. Escuder-Gilabert, L.; Sanchis-Mallols, J.M.; Sagrado, S.; et al. Chromatographic quantitation of the hydrophobicity of ionic compounds by the use of micellar mobile phases. J. Chromatogr. A. 1998, 823, 549–559. [Google Scholar] [CrossRef]
  84. Stępnik, K. Biomimetic Chromatographic Studies Combined with the Computational Approach to Investigate the Ability of Triterpenoid Saponins of Plant Origin to Cross the Blood-Brain Barrier. Int. J. Mol. Sci. 2021, 30, 3573. [Google Scholar] [CrossRef]
  85. Jusril, N.; Juhari, A.M.; Abu Bakar, S.; et al. Combining In Silico and In Vitro Studies to Evaluate the Acetylcholinesterase Inhibitory Profile of Different Accessions and the Biomarker Triterpenes of Centellaasiatica. Molecules 2020, 25, 3353. [Google Scholar] [CrossRef] [PubMed]
  86. Jamila, N.; Khairuddean, M.; Yeong, K.K.; et al. Cholinesterase inhibitory triterpenoids from the bark of Garcinia hombroniana. J. Enzym. Inhib. Med. Chem. 2014, 30, 133–139. [Google Scholar] [CrossRef] [PubMed]
  87. Stavrakov, G.; Philipova, I.; Lukarski, A.; et al. Galantamine-curcumin hybrids as dual-site binding acetylcholinesterase inhibitors. Molecules 2020, 25, 3341. [Google Scholar] [CrossRef]
Figure 1. The 1/k vs. CM relationships obtained for BMC (above) and SDS (below) system (the numbers of analytes according to Table 1).
Figure 1. The 1/k vs. CM relationships obtained for BMC (above) and SDS (below) system (the numbers of analytes according to Table 1).
Preprints 82729 g001
Figure 2. LogBB values obtained using biomimetic and computational methods (the numbers of analytes according to Table 1).
Figure 2. LogBB values obtained using biomimetic and computational methods (the numbers of analytes according to Table 1).
Preprints 82729 g002
Figure 3. Result of the TLC-bioautography assay for the inhibition of acetylcholinesterase enzyme on the silica gel showing different concentrations of the tested triterpenes: arjunic acid (I), arjunolic acid (II), arjungenin (III), arjunglucoside I (IV), sericic acid (V), and arjunetin (VI).
Figure 3. Result of the TLC-bioautography assay for the inhibition of acetylcholinesterase enzyme on the silica gel showing different concentrations of the tested triterpenes: arjunic acid (I), arjunolic acid (II), arjungenin (III), arjunglucoside I (IV), sericic acid (V), and arjunetin (VI).
Preprints 82729 g003
Figure 4. (A) The superposition of the most favorable poses of all ligands interacting with AChE. (B) The most favorable location of the arjunglucoside I molecule bound to AChE. The ligand molecule is shown as thick sticks whereas all the closest amino-acid residues (of distance no larger than 0.4 nm) are represented by thin sticks. The description of the interaction types is given in the text. The residue numbering is compatible with the PDB:3EVE record. (C) The correlation between the experimentally-determined IC50 values (recalculated as ln(IC50)) and the molecular volume (blue points) or molecular area (red points) of studied compounds.
Figure 4. (A) The superposition of the most favorable poses of all ligands interacting with AChE. (B) The most favorable location of the arjunglucoside I molecule bound to AChE. The ligand molecule is shown as thick sticks whereas all the closest amino-acid residues (of distance no larger than 0.4 nm) are represented by thin sticks. The description of the interaction types is given in the text. The residue numbering is compatible with the PDB:3EVE record. (C) The correlation between the experimentally-determined IC50 values (recalculated as ln(IC50)) and the molecular volume (blue points) or molecular area (red points) of studied compounds.
Preprints 82729 g004
Figure 5. Absorbance (A) vs. time [min] for 0.00045 mM of ARG.
Figure 5. Absorbance (A) vs. time [min] for 0.00045 mM of ARG.
Preprints 82729 g005
Figure 6. Lineweaver–Burk plot for initial velocity (Vo) vs. substrate concentration [mmol/L].
Figure 6. Lineweaver–Burk plot for initial velocity (Vo) vs. substrate concentration [mmol/L].
Preprints 82729 g006
Table 1. The chemical structures of the tested compounds.
Table 1. The chemical structures of the tested compounds.
No. Name Chemical structure
1 Arjunic acid Preprints 82729 i001
2 Arjunolic acid Preprints 82729 i002
3 Arjungenin Preprints 82729 i003
4 Arjunglucoside I Preprints 82729 i004
5 Sericic acid Preprints 82729 i005
6 Arjunetin Preprints 82729 i006
Table 2. The BBB-pharmacokinetic and distribution parameters of the analyzed compounds calculated in silico (ACD/Percepta software).
Table 2. The BBB-pharmacokinetic and distribution parameters of the analyzed compounds calculated in silico (ACD/Percepta software).
Name logBB logPS logPS,Fu,brain Fu Fb
Arjunic acid -0.14 -3.2 -4.9 0.012 0.02
Arjunolic acid -0.13 -3.2 -5.0 0.012 0.02
Arjungenin -0.15 -2.9 -4.6 0.016 0.02
Arjunglucoside I 0.12 -2.9 -4.4 0.050 0.04
Sericic acid -0.15 -2.9 -4.6 0.016 0.02
Arjunetin 0.73 -2.3 -4.3 0.051 0.01
Table 3. The IC50 values calculated in the TLC-bioautography assay towards the AChE inhibition for the tested standards.
Table 3. The IC50 values calculated in the TLC-bioautography assay towards the AChE inhibition for the tested standards.
Name IC50 [mg/mL] IC50 [mM]
1. arjunic acid 10.12 0.0207
2. arjunolic acid 7.92 0.0162
3. arjungenin 9.83 0.0195
4. arjunglucoside I 8.78 0.0132
5. sericic acid 8.3 0.0164
6. arjunetin 12.23 0.0188
Table 4. Acute and chronic toxicity predicted in silico.
Table 4. Acute and chronic toxicity predicted in silico.
Organism Duration Endpoint Predicted value [mg/L]
Acute effects
Fish 96 h LC50 566.412
Daphnia 48 h 1364.175
Green algae 96 h EC50 721.889
Chronic effects
Fish ChV 56.863
Daphnia 1385.105
Green algae 118.412
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated