Flavonoids and terpenoids with PTP-1B inhibitory properties from the infusion of salvia amarissima ortega

Article abstract of DOI:10.3390/molecules25153530

An infusion prepared from the aerial parts of Salvia amarissima Ortega inhibited the enzyme protein tyrosine phosphatase 1B (PTP-1B) (IC₅₀~88 and 33 μg/mL, respectively). Phytochemical analysis of the infusion yielded amarisolide (1), 5,6,4⁰-trihydroxy-7,3⁰-dimethoxyflavone (2), 6-hydroxyluteolin (3), rutin (4), rosmarinic acid (5), isoquercitrin (6), pedalitin (7) and a new neo-clerodane type diterpenoid glucoside, named amarisolide G (8a,b). Compound 8a,b is a new natural product, and 2–6 are reported for the first time for the species. All compounds were tested for their inhibitory activity against PTP-1B; their IC₅₀ values ranged from 62.0 to 514.2 μM. The activity was compared to that of ursolic acid (IC₅₀ = 29.14 μM). The most active compound was pedalitin (7). Docking analysis predicted that compound 7 has higher affinity for the allosteric site of the enzyme. Gas chromatography coupled to mass spectrometry analyses of the essential oils prepared from dried and fresh materials revealed that germacrene D (15) and β-selinene (16), followed by β-caryophyllene (13) and spathulenol (17) were their major components. An ultra-high performance liquid chromatography coupled to mass spectrometry method was developed and validated to quantify amarisolide (1) in the ethyl acetate soluble fraction of the infusion of S. amarissima.

Full text of DOI:10.3390/molecules25153530

molecules
Article
Flavonoids and Terpenoids with PTP-1B Inhibitory
Properties from the Infusion of
Salvia amarissima Ortega
Eric Salinas-Arellano 1 , Araceli Pérez-Vásquez , Isabel Rivero-Cruz , Rafael Torres-Colin ,
1
1
2
Martín González-Andrade 3 , Manuel Rangel-Grimaldo 1 and Rachel Mata *
1,
1
Facultad de Química, Universidad Nacional Autónoma de México, Mexico City 04510, Mexico;
ersalinass@hotmail.com (E.S.-A.); perezva@unam.mx (A.P.-V.); riveroic@unam.mx (I.R.-C.);
manuel_erg_p9@hotmail.com (M.R.-G.)
Instituto de Biología, Universidad Nacional Autónoma de México, Mexico City 04510, Mexico;
2
rafael.torres@ib.unam.mx
Facultad de Medicina, Universidad Nacional Autónoma de México, Mexico City 04510, Mexico;
3
martin@bq.unam.mx
*
Correspondence: rachel@unam.mx; Tel.: +52-55-56225289
ꢀ
ꢁꢂꢀꢃꢄꢅꢆꢇ
Academic Editor: Raffaele Capasso
Received: 17 June 2020; Accepted: 30 July 2020; Published: 1 August 2020
ꢀ
ꢁꢂꢃꢄꢅꢆ
Abstract: An infusion prepared from the aerial parts of Salvia amarissima Ortega inhibited the enzyme
protein tyrosine phosphatase 1B (PTP-1B) (IC ~88 and 33
µ
g/mL, respectively). Phytochemical
50
0
0
analysis of the infusion yielded amarisolide (
-hydroxyluteolin ( ), rutin ( ), rosmarinic acid (
neo-clerodane type diterpenoid glucoside, named amarisolide G (8a
natural product, and are reported for the first time for the species. All compounds were tested for
their inhibitory activity against PTP-1B; their IC50 values ranged from 62.0 to 514.2 M. The activity
was compared to that of ursolic acid (IC50 = 29.14 M). The most active compound was pedalitin
). Docking analysis predicted that compound has higher affinity for the allosteric site of the
enzyme. Gas chromatography coupled to mass spectrometry analyses of the essential oils prepared
from dried and fresh materials revealed that germacrene D (15) and -selinene (16), followed by
-caryophyllene (13) and spathulenol (17) were their major components. An ultra-high performance
1
), 5,6,4 -trihydroxy-7,3 -dimethoxyflavone (
2
),
) and a new
). Compound 8a is a new
6
3
4
5), isoquercitrin (6), pedalitin (7
,b
,b
2–6
µ
µ
(
7
7
β
β
liquid chromatography coupled to mass spectrometry method was developed and validated to
quantify amarisolide (1) in the ethyl acetate soluble fraction of the infusion of S. amarissima.
Keywords: Salvia amarissima; PTP-1B activity; amarisolide G; diabetes
1
. Introduction
Type 2 diabetes mellitus is a metabolic disease characterized by chronic hyperglycemia due to
insulin resistance, or the relative absence of the hormone. The prevalence of the disease is continuously
increasing, with approximately 463 million people living with diabetes nowadays. Mexico is one of the
countries more affected by type 2 diabetes mellitus, with more than 12 million cases. The population
sector more affected is the indigenous people owing to variations in its traditional way of life and
the effects of industrial developments [
to treat the symptoms of diabetes; in some cases, the patients combine allopathic therapies with the
botanical remedies [ ]. These plants are an essential part of the country’s alternative medical care,
1]. Mexican population employs more than 300 plant species
2
and the best testimony of their efficacy is their persistence in Mexican markets and other places for
crude or fresh drug selling. Therefore, it is crucial to analyze these plants to establish their composition,
Molecules 2020, 25, 3530; doi:10.3390/molecules25153530
www.mdpi.com/journal/molecules

Molecules 2020, 25, 3530
2 of 18
security, efficacy, and to develop a suitable methodology for quality control of the crude drug following
good practice guidelines.
Quality control of herbal drugs is the base for their efficacy and safety. Quality control
of herbal drugs aims to define their identity, purity and content of active principles or marker
compounds. The chemical composition of plants, and hence of their therapeutic preparations, is
variable, so standardization is necessary to guarantee comparable therapeutic effects. To prove the
constant composition of herbal preparations, there are appropriate standard analytical methods to
establish relevant criteria for uniformity. Standard analytical techniques include, among others,
high-performance liquid chromatography. For many years, the World Health Organization (WHO) has
encouraged all its country members to elaborate pharmacopeic monographs providing comprehensive
scientific information on the quality of their most commonly used medicinal plants. Following WHO
guidelines, Mexico has developed the Mexican Herbal Pharmacopeia, which contained monographs
with definitions, analytical techniques for identity and composition, as well as storage regulations of
the most widely used Mexican Herbal drugs [2,3].
Like other pharmaceutical products, herbal drugs should fulfill the basic requirements of being
efficacious and safe. To establish herbal drugs’ efficacy and safety is necessary to perform preclinical
and clinical assays, including those of the healer or medical doctor in rural communities who apply
locally produced herbal preparations. When assessing the efficacy of the plants, it is essential to
study both the traditional preparations and their components; this is because sometimes the efficacy
is attained with the combinations of compounds in the preparations, which may be acting through
synergy, network pharmacology or by targeting several nonrelated proteins involved in the pathology
of a disease. Altogether, these studies can lead to the development of standardized phytomedicines of
good quality and discover good drug candidates or molecules useful for lead optimization or even
fragment-based drug discovery [3].
Among the species highly valued in Mexico for treating diabetes is Salvia amarissima Ortega (syn.
Salvia circinata Cav.) belonging to the mint family. It is a perennial aromatic shrub native to Mexico,
listed as medicinal in the catalog of plants from the Royal Botanical Expeditions to New Spain. Like
many other New World Salvia species, S. amarissima is melittophilous (bee-pollinated). A tea brewed
from dried aerial parts of the plant is useful in Mexican folk medicine for treating diabetes, ulcers and
helminthiases [
to improve the diabetic condition [
neo-clerodane diterpenoids, including amarisolide (1) [
some flavonoids [6,8]. The traditional preparation of the plant collected in Puebla, Mexico, as well as
amarisolide (1), and some of the flavonoids showed inhibitory activity against mammal -glucosidases
in vitro and in vivo. The preparation and compounds were hypoglycemic and reduced the postprandial
peak significantly during an oral sucrose tolerance test in healthy mice [ ]. Some of the seco-clerodane
diterpenoids were cytotoxic against a few human cancer cell lines, had modulatory activity in a
breast cancer cell line resistant to vinblastine and exhibited antiprotozoal action [ 11]. Furthermore,
the traditional preparation lack toxicity when tested according to the Lorke criteria. [ ]. More recently,
4,
5
]. The species is commonly regarded as “insulina” (insulin), referring to its efficacy
]. Previous phytochemical studies allowed the isolation of some
], a few seco-clerodane diterpenoids [ 10] and
5
6
6–
α
8
7,9–
8
the antinociceptive properties of an aqueous extract of the plant, 1 and 7 were demonstrated [12].
Based on the above considerations, this investigation aimed: (i) to determine the effect of the
traditional preparation (infusion) of the plant collected in Oaxaca and their components on the activity
of the protein tyrosine phosphatase (PTP-1B) in order to assess a new molecular target, and get a
better insight in the in vivo hypoglycemic effect previously demonstrated [8]. This target was chosen,
considering that this enzyme acts as a negative regulator of insulin and leptin dependent signal cascades
holding therapeutic utility in type 2 diabetes mellitus and obesity [13]. (ii) To analyze the chemical
composition of the essential oil. (iii) To set up an appropriate procedure using Ultra-High-Performance
Liquid Chromatography (UHPLC) to quantify one active or marker compound of the plant’s infusion.
The chemical composition of the plant’s essential oil and the UHPLC procedure will allow developing

Molecules 2020, 25, 3530
3 of 18
a pharmacopeic monograph of S. amarissima, as they represent valuable identity and composition
tests, respectively.
2
. Results and Discussion
Scheme 1 summarizes the whole work.
Scheme 1. Workflow diagram of this study. * Efficacy studies. ** Quality control analyses. *** Discovery
of a lead molecule.
2
.1. Chemical Constituents of the Aqueous Extract
An aqueous extract (AE) from S. amarissima inhibited the PTP-1B activity significantly, with an
IC50 value of 88.6
afforded eight compounds (Figure 1), namely amarisolide (
), 6-hydroxyluteolin ( ), rutin ( ), rosmarinic acid ( ), isoquercitrin (
type diterpenoid glycoside, named amarisolide G (8a ). Compound 8a,b
±
5.4 µg/mL. Workup of AE by solvent partitioning and repeated chromatography
0
0
1
), 5,6,4 -trihydroxy-7,3 -dimethoxyflavone
), pedalitin ( ) and a neo-clerodane
, is a new natural product
(2
3
4
5
6
7
,b
characterized by conventional spectroscopic and spectrometric techniques. The known compounds
were identified by comparing their spectroscopic data with those previously described (Figures S1–S11,
Supplementary Material) [8,14,15]. Compounds 2–6 are reported for the first time for this species.
In addition, thin layer chromatographic (TLC) analysis revealed that these compounds were present in
the AE prepared from the fresh material.
Product 8a,b was obtained as a white solid optically active. DART-HRMS (Direct Analysis
In Real Time-High Resolution Mass Spectrometry) established its molecular formula as C H O .
2
6
36 11
−
1
−1
The IR spectrum included bands for hydroxyl (3365 cm ) and
functionalities (Figure S12, Supplementary Material) [ ]. The NMR spectra of 8a
S13–S18, Supplementary Material) were closely similar to those of amarisolide D (8c), a neo-clerodane
type of diterpenoid with an -unsaturated- -lactone at C-4/C-5, a five-membered keto- -lactol
methyl ether and a -D-glucopyranosyloxy moiety at C-2 [ ]. Thus, in compound 8a the signals
for the five-membered keto- -lactol methyl ether were replaced by those of a keto- -lactol moiety.
α
,
β-unsaturated-
γ
-lactone (1749 cm
)
8
,b (Table 1; Figures
α
,β
γ
γ
β
8
,b
γ
γ

Molecules 2020, 25, 3530
4 of 18
Accordingly, the NMR spectra showed signals at
δ 5.90, 6.05 (brs, H-14)/δC 99.4, 99.9 (C-14), δ
H H
5
.96 (brs, H-16)/δ 116.1 (C-16), δ 171.4 (C-13) and δ 171.9 (C-15) (Table 1). Since the resonances
C C C
for H-14 appeared as two separate signals, the intensity of each corresponding to one-half proton,
13
compound 8a,b is a 1:1 mixture of C-14 epimers. The doubling of C-14 signal in the C-NMR spectra
13
provided further evidenced (Table 1). The C-NMR chemical shifts (Table 1) of C-11-C-16 and the
key HMBC (Heteronuclear Multiple Bond Correlation) correlations of H-12 and H-14 established
the point of attachment of the ethyl fragment (C-11-C-12) to the keto−γ-lactol ring. The most
relevant HMBC correlations were H-12a ( H
δ
1.72) with C-14 ( C
5.90 and 6.05) with C-12 ( C
H
δ 99.4 and 99.9); H-12b (δ 1.56) with
H
C-14 (
δ
C
99.4 and 99.9); and H-14 (
δ
δ
33.8). The NOESY (Nuclear
Overhauser Effect Spectroscopy) interactions revealed that the relative configuration at the stereogenic
centers of 8a was identical to that of amarisolides A–D [ ]. The electronic circular dichroism
spectrum of 8a showed negative Cotton effects at ~212 and ~250 nm due to the electronic transitions
π→π* and n→π*, respectively, of the -unsaturated- -lactone. The latter data indicated that
the absolute configuration at the stereogenic centers C-2, C-5, C-8, C-9 and C-10 of compound
was S, S, R, R and R, respectively. The D configuration of the -glucopyranosyloxy moiety
was established as previously described [ ]. On the basis of these evidences, compound 8a
was characterized as (2S,5S,8R,9R,10R,14R,S)-2-(O- -d-glucopyranosyl)-neo-clerodan-14-hydroxy-3,
3-diene-14,15;18,19-diolide (8a,b) and was designated with the trivial name of amarisolide G.
,
b
8
,
b
α
,
β
γ
8
a
,b
β
8
,b
β
1
Figure 1. Structures of compounds (1–8a,b) from S. amarissima.
2
.2. Chemical Constituents of the Essential Oil
Dried and fresh plant materials were analyzed to assess any change during the drying process.
Since it is an aromatic plant, the chemical profile of the essential oil is valuable as an identity test.
The essential oil content of distilled aerial parts was 0.02% and 0.03% for fresh and dried material,
respectively. In each case, eight major compounds were identified, representing 99.96 and 99.97%
of the composition of the essential oil, respectively. As observed in Table 2 and Figures S19 and
S20 (Supplementary Material), the major components in both samples were germacrene D (15) and

Molecules 2020, 25, 3530
5 of 18
β-selinene (16), followed by
β-caryophyllene (13) and spathulenol (17). The only mutually exclusive
constituents were 3-methoxy-p-cymene (9) present in the dried material, and δ-elemene (11) found
only in the fresh plant. These differences are not due to seasonal changes because the plant material
was from the same batch.
Table 1. H and 13C NMR spectroscopic data for compound 8a,b in DMSO-d6.
1
a
b
Position
δ
(J/Hz)
δ
C
H
1
α: 1.35 dd (13.6, 3.2)
β: 1.81 brd (13.7)
4.44–4.49 m
26.3
2
3
4
5
6
70.0
131.1
143.0
45.4
6.66 d (6.4)
α: 1.72 td (12.7,3.3)
β: 1.22–1.26 m
α: 1.53–1.58 m
β: 1.65−1.70 m
1.65–1.70 m
33.8
7
27.6
8
9
36.2
37.8
40.5
27.6
33.8
1
1
1
0
1
2
2.21 d (15.1)
1.59–1.64 m
α: 1.71–1.74 m
β: 1.53–1.64 m
1
1
3
4
171.4
99.4
99.9
α: 5.90 brs
β: 6.05 brs
1
1
1
1
5
6
7
8
171.9
116.1
15.9
169.2
70.9
5.96 brs
0.80 d (6.5)
1
9a
4.39 d (8.2)
4.02 d (8.2)
0.57 s
1
9b
2
1
2
3
4
5
0
17.7
102.8
74.1
77.2
70.7
77.4
0
0
0
0
0
4.31 d (7.7)
2.96–3.05 m
3.13–3.16 m
2.96–3.05 m
3.13–3.16 m
α: 3.43–3.47 m
β: 3.63–3.69 m
7.75 brs
0
6
61.6
1
2
3
4
6
4-OH
0
0
0
0
-OH
-OH
-OH
-OH
4.92 brs
4.92 brs
4.92 brs
3.33 s
Assignments based on the analysis of the HSQC (Heteronuclear single quantum coherence spectroscopy) and
HMBC experiments; J values (Hz) in parentheses. 700 MHz. b 175 MHz.
a
It is worth mentioning that the oils from other Salvia species analyzed also contains sesquiterpenes.
In this context, Salvia ceratophylla, S. aethiopis L., S. palaestina Bentham and S. xanthocheila Boiss. ex
Benth., are rich sources of germacrene D (15) [16]; β-caryophyllene (13) is the major component of
S. nemorosa L., S. verticillata L., S. virgata Ortega and S. hydrangea DC. ex Benth. Finally, germacrene B
was the major compound of S. syriaca L. [17,18].

Molecules 2020, 25, 3530
6 of 18
Table 2. GC-MS identified compounds from the essential oils of S. amarissima.
Peak Area (%)
Sa-Batch 1
Compound
CAS Number
R [a]
I
Sa-Batch 2
3
-Methoxy-p-cymene (9)
1076-56-8
1686-15-3
20307-84-0
5208-58-2
87-44-5
1139-30-6
37839-63-7
17066-67-0
6750-60-3
1219
1313
1329
1378
1413
1447
1476
1491
1576
4.40
2.10
(
E)-Pinocarvyl acetate (10)
δ-Elemene (11)
5.98
1.99
4.24
15.05
7.68
25.09
28.35
11.59
99.97
α-Bourbonene (12)
β-Caryophyllene (13)
α-Caryophyllene (14)
Germacrene D (15)
β-Selinene (16)
4.11
21.27
4.21
24.26
30.16
9.45
Spathulenol (17)
Total peak area (%)
99.96
[a] RI retention indices relative to series of n-alkanes (C8-C20) on a DB-5 column.
2
.3. Evaluation of Compounds 1–8a,b on the PTP-1B Inhibitory Activity
One of the major causes of type 2 diabetes mellitus is insulin resistance, which occurs when
the hormone cannot activate signaling pathways in major metabolic tissues (muscles, fat and liver).
Insulin resistance involves several inhibitory molecules that interfere with tyrosine phosphorylation
of the insulin receptor. Among those, PTP-1B is a promising therapeutic target since it inactivates
crucial signaling effectors in the insulin- and leptin-signaling cascades by dephosphorylating their
tyrosine residues [13]. Therefore, natural products
the flavonoids tested, the most active component was pedalitin (
19.0 g/mL), whereas of the diterpenoids was amarisolide ( ) (279.9
cases, the inhibitory effect was concentration-dependent. The inhibitory activity against PTP-1B
reported for compounds and correlated well with that found in this work [19 20]. On the other hand,
rosmarinic acid ( ) is an ingredient of an active extract of Rosmarinus officinalis against PTP-1B [21].
However, in this study we report for the first time the effect of the pure . The activity of the isolates
8a,b was compared to that of ursolic acid (UA, IC50 = 28.1 1.2 M), which in other studies displayed
lower IC50 values (~3.1 µM) [22].
1–
8a,b were tested against PTP-1B (Table 3). Among
) with an IC50 of 62.0 4.1
26.0 M; 137.9 g/mL). In both
7
±
±
µM
(
µ
1
µ
µ
4
6
,
5
5
1
–
±
µ
Table 3. Protein tyrosine phosphatase 1B (PTP-1B) inhibitory activity of compounds (1–8a,b) from
S. amarissima expressed as IC .
50
µg/mL ± SD
µM ± SD
1
2
3
4
5
137.9 ± 12.8
62.5 ± 2.3
24.2 ± 0.7
197.3 ± 10.1
49.4 ± 6.7
120.7 ± 14.7
19.6 ± 1.3
269.7 ± 3.6
-
279.9 ± 26.0
189.4 ± 7.0
80.1 ± 2.2
323.4 ± 16.6
137.1 ± 18.7
259.9 ± 31.7
62.0 ± 4.1
6
7
8
a,b
UA [a]
514.2 ± 6.8
28.1 ± 1.2
[
a] Positive control. Values are representative of three independent determinations.
The results of the PTP-1B are significant. They suggest that the traditional preparation
of S. amarissima, with hypoglycemic and antihyperglycemic properties demonstrated in vivo [8]
(
(
i.e., the overall action), contains compounds such as
i.e., PTP-1B and others) within the same signaling network thus shutting insulin signaling cascade
and 8a , with weaker activity
, altogether put forth a biochemical effect by synergism (i.e., a synergy between weakly
1–8a,b that might weakly target different proteins
process by network pharmacology. It is also possible that compounds
than compound
1
–6
,b
7

Molecules 2020, 25, 3530
7 of 18
active compounds against PTP-1B). Finally, molecules like compounds
1
,
4,
5
and
7
can exert their
action binding different targets such as PTP-1B and
The fact that rutin ( ) [23] and rosmarinic acid ( ) [21] are multitarget antidiabetic compounds, and
compounds 1 and 7 inhibited α-glycosidases in vivo [8] strengthen any of these possibilities.
α-glycosidases, among others (polypharmacology).
4
5
2
.4. Docking Study
To predict the preferred binding orientation of compounds 1, UA and
7 into PTP-1B, we
performed a docking analysis. Compound
included for comparative purposes. These substances
1
was not among the most active compounds but
and UA were docked with the co-crystallized
1, 7
structure of PTP-1B at the catalytic (PDB ID: 1G7F) and allosteric (PDB ID: 1T49) sites. The results in
Figure 2 and Table 4 indicate that all ligands bind at the same site as the co-crystallized ligands 892
(
3-(3,5-dibromo-4-hydroxy-benzoyl)-2-ethyl-benzofuran-6-sulfonic acid (4-sulfamoyl-phenyl)-amide)
and INZ (2-{4-[(2s)-2-[({[(1s)-1-carboxy-2-phenylethyl]amino}carbonyl)amino]-3-oxo-3-(pentylamino)
propyl]phenoxy}malonic acid) [24 25]. The estimated energy binding is different for each compound.
has a higher affinity for the catalytic site, which is more and conserved site, while
and UA targeted its more hydrophobic and less conserved allosteric site. The RMSD
and are in the range of 2 to 3 Å, which indicates that the ligands do not precisely
overlap but maintain the correct orientation [26 27]. The amino acids interacting with compounds
and , as well as AU are similar to those previously reported for other inhibitors (Table 4; Figures S30
and S31 of the Supplementary Material) [24 25]. At the allosteric site, compound has hydrophobic
interactions with Ala189, Leu192, Phe280 and Phe196 while AU with Leu192, Phe280 and Phe196.
Therefore, pedalitin ( ) behaves as UA regarding its higher affinity for the allosteric site. For UA, the
key structural feature is the pentacyclic core displaying a nonpolar characteristic, which interacts with
nonpolar residues in the allosteric site [28]. For compound , apparently beside the tricyclic structure,
the lack of oxygen at C-3 of the flavonoid core seems to be essential. It will be necessary to pursue a
kinetic analysis to determine if experimentally compound is an allosteric inhibitor of PTP-1B. It is
,
Compound
compound
values for ligands
1
7
1
7
,
1
7
,
7
7
7
7
essential to mention, that the less-conserved PTP-1B allosteric site is an ideal target for a compound to
inhibit PTP-1B activity because the problems associated with inhibition at catalytic site will disappear.
Thus, this study may provide an important lead among flavones for the development of novel PTP-1B
allosteric inhibitors.
Figure 2. Structural models of the binding sites of the PTP-1B ligand complexes. In cyan cartoons
(PTP-1B protein), red sticks (
1), green sticks (7) purple sticks (UA) at the catalytic (A) and allosteric
sites (B). Images created with Pymol.

Molecules 2020, 25, 3530
8 of 18
Table 4. Results of the dockings analyses.
Catalytic Site
Allosteric Site
RMSD (Å) 3
Ki
BE
(kcal/mol)
Ki
BE
(kcal/mol)
RMSD (Å) 3
Residues 2
Residues 2
µM) 1
(µM) 1
(
Tyr46, Asp48, Val49, Phe182, Gly183,
Cys215, Ser216, Ala217, Gly220,
Arg221, Gln262, Thr263 and Gln266
Phe196, Lys197, Arg199,
Glu200, Gly277, Phe280
and Ile281
1
7
0.55
9.94
2.56
−8.5
−6.8
−7.6
2.89
8.39
2.16
0.78
−6.8
−7.7
−8.3
2.15
2.27
1.73
Tyr46, Asp48, Val49, Trp179, Asp181,
Phe182, Gly183, Ala217, Ile219,
Arg221, Gln262 and Gln 266
Ala189, Leu192, Glu276,
Gly277, Phe280, Phe196,
Lys197 and Glu200
2.73
Phe182, Gly183, Cys215, Ala217,
Gly218, Ile219, Gly220, Arg221,
Gln262, Thr263 and Gln266
Asn193, Phe196,
Asp263, Phe280, Leu192
and Ile281
UA
1.35
1
3
i i
values were calculated using the equation ∆G = RT lnK [
26] 2 Residues interacting 4 Å of the compound;
K
RMSD were calculated from the co-crystallized ligands for UA; and for 1 and 7 from UA. BE, binding energy.
2
.5. Molecular Dynamics of PTP-1B-Compounds 1, 7 and UA Complexes
Molecular dynamics (MD) studies of the complexes PTP-1B-compounds 1, 7 and UA were carried
out to evaluate the stability of the docked complexes illustrate in Figure 3. Table 5 shows the theoretical
parameters calculated from 100 ns of MD. All PTP-1B-compound complexes show negative G (affinity
parameter) consistent with their stability. Compound has a G similar to UA at the catalytic site,
while compound has a lower G than AU at the same site. Figures 3 and 4 show the structural models,
RMSD and Root Mean Square Fluctuations (RMSF) of the molecular trajectories for compounds and
, respectively. The RMSD of the complexes is lower with respect to PTP-1B in both the catalytic and
∆
1
∆
7
∆
1
7
allosteric MDs, which indicates the conformational stability of the complexes. In the RMSF analysis,
it is observed an area between residues 27–50 (corresponding to a loop), which is stabilized with the
ligands. The data obtained with the MDs are complementary and conscious with the docking data and
experimental results.
Figure 3. Molecular dynamics at the catalytic site of PTP-1B for compound
1. (A) Structural models,
(
B) RMSD and (C) Root Mean Square Fluctuations (RMSF) from molecular trajectories. AU = UA.

Molecules 2020, 25, 3530
9 of 18
Table 5. Calculations of the theoretical energy components from the trajectories of molecular dynamics.
Complex.
∆Evdw
∆Eele
20.02 1.08
−1.21 ± 0.46
18.53 0.77
−4.27 ± 0.31
∆GGB
∆GNP
∆Gsolv
∆Gbind
1
2
PTP-1B-UA
PTP-1B-UA
PTP-1B-1
PTP-1B-7
1
−
−
−
−
21.88
±
±
±
±
0.44
−
±
27.60 ± 1.63
10.89 ± 0.42
41.62 ± 0.55
12.90 ± 0.33
−
−
−
−
41.91
±
±
±
±
1.26
25.39 ± 1.04
7.73 ± 4.21
36.84 ± 0.53
9.06 ± 0.33
−
−
−
−
16.51
20.05
17.92
31.16
±
±
±
±
0.36
0.21
0.61
0.46
28.67
36.22
35.95
0.23
0.47
0.45
27.78
54.76
40.22
0.48
1.00
0.55
−
±
Catalytic site; 2 Allosteric site;
∆
G
E
vdw, contributions by van der Waals interactions;
np, nonpolar solvation energy; solv, desolvation free energy (
Gnonpol); ∆Gbind, binding affinity.
∆Eele, electrostatic energy;
∆
∆
G
GB, polar solvation energy;
∆
∆
G
solv GB
∆G = ∆G +
Figure 4. Molecular dynamics at the allosteric site of PTP-1B for compound
B) RMSD and (C) RMSF from molecular trajectories.
7. (A) Structural models,
(
2
.6. Drug Likeness for Compounds 1, 7 and UA
According to the SwissTargetPrediction (http://www.swisstargetprediction.ch/index.php),
and Molinspiration (http://www.molinspiration.com/cgi-bin/properties) databases which predict
the most probable targets of small bioactive molecules, compounds and UA 1 could target any
1
,
7
, 7
protein with percentages of 26.7, 40 and 40 %, respectively. SwissTargetPrediction predicted PTP-1B
inhibition for these compounds with probabilities of 0, 0.1266 and 0.95, respectively. According to these
predictions, UA should have been 7.5 times more active than compound
UA was only 2.2 times more than compound 7 (Figures S31–S34).
7; however, experimentally,
Next, using the Osiris Property Explorer server (http://cheminformatics.ch/propertyExplorer),
relevant properties for compounds and UA were calculated and summarized in Table 6.
1, 7

Molecules 2020, 25, 3530
10 of 18
These properties indicate whether a molecule is a potential drug. The logP value is a measure
of a compound’s hydrophilicity. Low hydrophilicity and, therefore, high logP values cause poor
absorption. For a compound being well absorbed, its logP value must not be greater than 5.0; logP
values between 1.35 and 1.8 indicate perfect oral and intestinal absorption. Thus, compounds
1 and 7
could have proper absorption, but not UA.
Table 6. Physicochemical properties of compounds 1, 7 and UA.
1
7
UA
LogP
1.01
−3.48
492.56
−2.96
9
4
10
6
0.40
1.92
−2.58
316.26
1.8
6
−6.11
456.71
−3.66
3
2
10
1
0.26
Solubility (LogS)
Molecular weight
Druglikeness
H bond acceptor
H bond donor
Nb stereocenters
Nb rotable bonds
Drug score
7
4
0
2
0.52
The data was calculated using the OSIRIS Property Explorer server (http://www.cheminfo.org/Chemistry/
Cheminformatics/Property_explorer/index.html).
The aqueous solubility of a compound influences its absorption and distribution characteristics;
a low solubility goes along with inadequate absorption. More than 80% of the drugs on the market
have a logS value higher than
−
4. Compounds
1
and UA present
−
3.48 and
−
6.11, respectively.
The drug-likeness parameter is a complex balance of various molecular properties and structural
features that determine whether a molecule is similar to the known drugs. These properties, mainly
hydrophobicity, electronic distribution, hydrogen bonding characteristics, molecule size and flexibility,
and of course presence of various pharmacophoric features influence the behavior of molecule in a
living organism, including bioavailability, transport properties, affinity to proteins, reactivity, toxicity,
metabolic stability and many others. A positive value indicates that a molecule contains predominantly
fragments, which are frequently present in commercial drugs [29]; compound
value of 1.8. The H bond acceptor and H bond donor’s parameters indicate a molecule’s ability to
interact to a greater or lesser degree with a protein; compounds and present a higher number of
possible hydrogen bridges than UA. Finally, the drug-score is an indicator that qualifies the potential
7 has a drug-likeness
1
7
of a compound for being a drug based on all the calculated parameters; compound
drug-score (0.52), which is in harmony with the experimental and theoretical data.
7 has the best
2
.7. Development and Validation of an UHPLC-MS Method for Quantifying 1
Initial assessments about the complexity of samples of AE were based on visual comparisons
of their chromatographic profiles. The ethyl acetate soluble fraction of AE yielded the best
®
profile. Chromatographic separation was performed on Acquity UHPLC BEH Shield C column
18
(
2.1 × 100 mm, 1.7 µm) applying a binary gradient elution of water (0.1% formic acid) and MeCN.
The total run time was 10 min. As illustrated in Figure 5 compounds and are present in the
1
,
2
4–7
chromatogram. All compounds were identified by their retention times and m/z values corresponding
−
to [M
−
H] ions. Compounds
1
,
2
and
4
–
7
showed an effective baseline resolution. The pseudo
−
−
molecular ions of these compounds appeared at m/z 609.54 [M
; R 1.60 min), 359.23 [M
R 3.65 min) and 491.23 [M
was amarisolide (
isolated from this species, which makes it an excellent marker compound for quality control. It is
worth mentioning that the neo-clerodanes type compounds detected in the infusion were and 8a,b
but not the minor diterpenoids we previously isolated from the organic extract of the plant [8].
−
H] (
4
; R 1.36 min), 463.46 [M − H
]
T
−
−
−
(6
−
H] (5; R 2.46 min), 315.48 [M
−
−
H] (7; R 2.82 min), 329.70 [M
−
H] (
2
;
T
T
T
−
H] (
1; R 4.89 min). The main component of the ethyl acetate fraction
T
T
1), then selected as a marker for validation. So far amarisolide (1) has been only
1
,

Molecules 2020, 25, 3530
11 of 18
Figure 5. Representative LC chromatogram of the ethyl acetate soluble fraction prepared from the
infusion of S. amarissima (Detection: 270 nm). Peak identification: : R 1.36 min; : R 1.60 min; : RT
.46 min; 7: R 2.82 min; 2: R 3.65 min; 1: R 4.89 min.
4
6
5
T
T
2
T
T
T
The analytical method was validated in terms of precision, accuracy, linearity and recovery
according to the Q2 (R1) guideline published by the International Conference on Harmonisation
(
ICH) [30]. The linearity of the system was tested using a concentration range of
1
between 5
2
to 100
µg/mL and was found to be linear (R = 0.9994 (UV) and 0.9921 (Electrospray Ionization
Mass Spectrometry, ESI-MS) in the concentration range used. The CV was less than 0.13% at each
concentration level analyzed. The limit of identification (LOD) and quantification (LOQ) values
were 1.22 and 3.70
ESI-MS detection.
µg/mL, respectively for UV detection; and 0.60 and 1.82 µg/mL, respectively for
The linearity of the method was tested by recovery assay; the linear regression equation were
found to be y = 5636.35x 4337.55 (UV detection) and y = 5252.80x + 64,101.76 (ESI-MS detection).
−
The recovery ranges for the standard were expressed as the concentration detected as a percentage of the
expected concentration and were found in the ranges of 100.7–101.7 % for UV detection and 83.7–95.7%
for ESI-MS detection. The reproducibility and repeatability of the analytical method were evaluated in
terms of the intermediate precision by analyzing 6 replicates of the stock solution (50
different days. The relative SD (RSD; n = 6) was calculated for each sample evaluated. The results
indicated that their chromatographic patterns were similar showing the presence of amarisolide ( ).
was quantified, and the
µg/mL) in two
1
The CV values for accuracy were less than 0.11%. Subsequently, compound
1
mean concentration calculated was 116 mg/g in dry matter.
3
. Materials and Methods
3
.1. General Procedures
IR spectra were recorded using a Bruker Tensor 27 FT-IR spectrophotometer (Bruker Corp.,
Billerica, MA, USA). Optical rotations were recorded at the sodium d-line wavelength using a Perkin
◦
Elmer model 343 polarimeter at 20 C (Perkin Elmer, MA, USA). NMR spectra were registered
1
on a Bruker AVANCE III HD with TCI CryoProbe 700 H-C spectrometer at 700 MHz ( H) or 175
13
MHz ( C), using TMS as an internal standard (Bruker Corp., Billerica, MA, USA). DARTHRMS
were acquired with a JEOL AccuTOF-DART JMS-T100LC (JEOL Ltd., Tokyo, Japan) spectrometer in
positive mode. For GC-MS analyses, an Agilent 6890N series gas chromatograph coupled to a LECO
(Laboratory Equipment Corporation) time-of-flight mass spectrometer detector (MS-TOF; Agilent
Technology, Santa Clara, CA, USA) was used. UHPLC-MS analyses were performed on a Waters

Molecules 2020, 25, 3530
12 of 18
Acquity UHPLC-H
®Class system (Waters, Darmstadt, Germany) equipped with a quaternary pump,
sample manager, column oven and photodiode array detector (PDA) interfaced with an SQD2 single
mass spectrometer detector with an electrospray ion source. Column chromatography (CC) was carried
out on Sephadex LH-20 (GE Healthcare, IL, USA). Thin layer chromatographic (TLC) analyses were
performed on silica gel 60 F254 plates (Merck, Darmstadt, Germany), or C -silica gel matrix plates
18
Analtech plates (Merck, Darmstadt, Germany), visualization of the plates was carried out using an
(
NH ) Ce(SO ) (10%) solution in H SO . Reagent-grade EtOAc, CHCl , CH Cl and MeOH were
4 4 4 4 2 4 3 2 2
purchased from J.T. Baker (Avantor Performance Materials, PA, USA). MeCN, MeOH and water LC-MS
or HPLC grades were purchased from Honeywell Burdick & Jackson (Morristown, NJ, USA). All other
analytical grade solvents and reagents were obtained from various commercial sources. Amarisolide
(1
) was isolated from the species S. amarissima in the present study. The purity was determined to be
more than 98% by UHPLC-MS.
3
.2. Plant Material
Salvia amarissima was collected in Capulálpam de Méndez, Ixtlán de Juárez, Oaxaca, in January
2
019 (Sa-Batch 1 (fresh) and Sa-Batch 2 (air-dried)). A voucher specimen (Number 1502277) was
a, UNAM. R. Torres-Colin
achieved the botanical identification of the plant. The plant was air-dried and ground to a fine powder
2 mm) in a Thomas Wiley Model 4 Mill.
deposited at the National Herbarium of Mexico (MEXU), Instituto de Biolog
í
(
3
.3. Extracts and Essential Oils Preparation
AE from S. amarissima (dried aerial parts) was prepared with 250 mL of boiling water and 12.5 g
of the crude drug for 30 min. After filtration, the aqueous extract was concentrated in vacuo to obtain
.1 g of a green residue. This process was repeated as necessary to prepare 10 g of AE. The ethyl acetate
soluble fraction was prepared via partitioning with EtOAc (3 250 mL) from the aqueous extract.
0
×
The resulting organic phase was dried over anhydrous sodium sulfate and concentrated in vacuo to
yield 130 mg of a brown residue (yield 1.0%).
EOs were prepared from fresh (Sa-batch 1) and air-dried (Sa-batch 2) plant material (200 g in 1.5 L
of distilled water) by hydrodistillation in a modified Clevenger type apparatus for 3 h. In both cases,
the hydrodistilled was extracted with CH Cl (3
×
2 L). The resulting organic phases were dried over
2
2
Na SO and concentrated in vacuo to yield an oily yellow residue (0.040 g, yield 0.02% in the case of
2
4
◦
4 C until
the fresh material, and 0.063 g, yield 0.03% for the dried plant). All samples were stored at
chemical analysis.
−
3
.4. Separation of Active Compounds from the Ethyl Acetate Soluble Fraction
The ethyl acetate soluble fraction (100 mg) was subjected to CC on Sephadex LH-20 (400 g)
using MeOH as eluent; fractions were pooled into 20 secondary fractions (F –F ) according to their
1
20
TLC profiles. From fraction F16 (73 mg) crystallized 60 mg of amarisolide (
of fraction F20 (10 mg) yielded 1.2 mg of 5,6,4 -trihydroxy-7,3 -dimethoxyflavone (2) and 4.2 mg of
1
). Preparative RP-TLC
0
0
6
-hydroxyluteolin (3).
3
.5. Separation of Active Compounds from AE
AE (4.3 g) was fractionated via CC on Sephadex LH-20 using a gradient system of methanol–water
water 40–100%); this process gave 12 secondary fractions (AE –AE ). From fraction AE (35 mg)
(
1
12
9
◦
crystallized 30 mg of rutin (
4; m.p. 241–242 C). From fraction AE11 (15 mg) crystallized 5 mg of
isoquercitrin (
to yield 16 mg of rosmarinic acid (
fraction AE (11 mg) afforded 1 mg of pedalitin (
6). Fraction AE (220 mg) was further purified on a Sephadex CC, eluting with MeOH,
6
5). Preparative TLC on silica gel [ethyl acetate-methanol (85:15)] of
7
). Finally, preparative RP-TLC of AE (10 mg, MeOH)
5
6
afforded 4 mg of an epimeric mixture of amarisolide G (8a,b).

Molecules 2020, 25, 3530
13 of 18
◦
20
Amarisolide G (8a,b): White solid; m.p. 133–135 C. [α] = –149 (c = 1 mg/mL, MOH). UV
D
−1
(
(
MeOH): λmax (log ε) 206 (0.612) nm. IR (KBr): νmax 3365, 1749 cm . ECD (c 0.2 mM, MeOH): λmax
∆ε) 212 (–5.21), 250 (–4.65) nm. H and C-NMR: see Table 1. HRESIMS: m/z 525.2317 [M + H] (calc.
1
13
+
5
25.2330 for C H O ).
26 37 11
3
.6. Enzymatic Hydrolysis of 8a,b
Compound 8a (2 mg) was mixed with
,b
β
-glucosidase (2 mg, Sigma-Aldrich, MO, USA) in
◦
phosphate buffer solution (2 mL, 100 mM at pH 7); and kept at 40 C for 15 days. Subsequently,
the reaction mixture was extracted with CHCl and the aqueous phase was concentrated to dryness
3,
and subjected to TLC analysis. d-Glucose was identified by comparison of the retention factor and
optical rotation value with those of the authentic sample.
3
.7. UHPLC-MS Analysis and Method Validation
®
The analytical method (Figure 5) was developed using an Acquity UHPLC BEH Shield C18
◦
column (2.1
×
100 mm, 1.7
µ
m) at 40 C. The mobile phase consisted of (A) water (0.1% formic acid)
and (B) acetonitrile with a linear gradient elution program: 0–10 min, 20–100% (B); 10–10.5 min, 20%
B); 10.5–13 min, 20% (B). The flow rate was set to 0.3 L/min, and the sample injection volume was
.0 L; detection was achieved with a PDA detector at 270 nm. For the identification of compounds,
each sample was analyzed with the electrospray ion source operating in both positive (ESI ) and
negative (ESI ) ionization modes. The ESI-MS conditions consisted of capillary voltage at 3.0 or 2.5 kV
(
µ
3
µ
+
−
◦
in positive and negative ion modes, respectively; dry heater temperature 150 C; and nitrogen as the
sheath gas flow. MS spectra were obtained within a mass range of m/z 100–1000 using nitrogen as the
collision gas. The MassLynk software (version 4.1) was used to control of the UHPLC-MS system and
for data acquisition and processing.
The method was validated according to the ICH guidelines [30]. For linearity, amarisolide (
accurately weighed and dissolved in dioxane-methanol (v/v, 1:1) to prepare stock solution at a final
concentration of 1 mg/mL. Six working solutions in the range of 5–100 g/mL for the standard were
1) was
µ
prepared from serial dilutions from the stock solution. Each concentration was prepared in sextuplicate.
2
The linearity was assessed estimating the slope, y-intercept and coefficient of determination (R ) using
the least squares method. Limits of detection (LOD) and quantification (LOQ) for the standard were
determined at signal-to-noise (S/N) ratios of 3 and 10, respectively. Recovery experiments were carried
out to evaluate the accuracy, assaying independently three amounts equivalent to 50 (ca. 10
00 (ca. 50 g/mL) or 125% (ca. 75 g/mL). At each level, compound was added simultaneously to
the ethyl acetate soluble fraction (50 g/mL). Each sample was injected twice and analyzed according
µg/mL),
1
µ
µ
1
µ
to the method previously described. The mean percentage recovery for the standard was found to
be between 98 and 102% by means of Fisher’s F test [30]. Finally, the repeatability and inter-day
precision was evaluated by testing six identical samples according to the above described method on
two consecutive days and by two different analysts by triplicate. The relative standard deviation (RSD)
was calculated for each determination as a measure of precision and repeatability.
3
.8. GC-MS Analysis of the Essential Oils
For GC-MS analyses, compounds were separated on a DB-5 capillary column (Supelco, Bellefonte,
PA, USA) with the following temperature program: oven temperature was programmed from 40 to
◦
◦
◦
2
60 C at 4 C/min during 20 min, and finally up to 340 C for 20 min isothermally; injector and MS
◦
transfer line temperatures were set at 200 and 300 C, respectively; Helium was used as the carrier gas
at a constant flow rate of 1 mL/min; split ratio, 1:20. A mixture of the homologous series of n-alkanes
(
C –C ) in CH Cl was directly injected into the GC under the above temperature program, in order
8 20 2 2
to calculate the linear retention indices (R ). All mass spectra were acquired in EI mode (scan range m/z
I
4
0–400, ionization energy 70 eV). The components were identified using retention index (R ) of peaks
I
in the chromatogram [31,32] and by comparison of their mass spectra with those of standard library

Molecules 2020, 25, 3530
14 of 18
data (NIST) of the GC-MS system and literature data or with those of authentic samples available
commercially. All determinations were performed in triplicate.
3
.9. Protein Tyrosine Phosphatase 1B Inhibition Assay
The expression and purification of hPTP-1B was performed as previously described [33]. Aqueous
extract (AE), ethyl acetate soluble fraction, pure compounds and positive control were dissolved in
DMSO or MeOH or Tris buffer solution (Tris-HCl, 20 mM, pH 7). Aliquots of 0–10 L of testing
materials (triplicated) were incubated for 5 min with 20 L of enzyme stock solution in Tris-HCl
L of p-nitrophenylphosphate (pNPP; 5 mM) was added and further
µ
µ
(
22 nM). After incubation, 10
µ
◦
incubated for 15 min at 25 C; then, the absorbance was determined (λmax 415 nm). For all samples,
the inhibitory activity was determined as a percentage in comparison to the blank (Tris-HCl) according
to the following equation:
!
A
4
15t
%
PTP1B = 1 − A
× 100
(1)
4
15C
where % PTP-1B is the percentage of inhibition, A415t is the corrected absorbance of the extracts, fraction,
or compounds under testing (A415end A415initial) and A415C is the absorbance of the blank (A415end blank
A415initial blank). The IC50 was calculated by regression analysis, using the following equation:
−
−
A
1
00
%
Inhibition =
ꢀ
ꢁ
(2)
s
I
1
+
IC50
where A100 is the maximum inhibition, I is the inhibitor concentration, IC is the concentration required
50
to inhibit the activity of the enzyme by 50% and S is the cooperative degree.
3
.10. Docking Studies
To perform the docking at the catalytic site, the PTP1B-INZ complex corresponding to the PDB
1
G7F was used, which has a resolution of 1.8 A. For the allosteric site, the PDB 1T49 (resolution of 1.9)
was used. The two PDBs used were selected considering the resolution and that they had co-crystallized
ligands at the sites of interest. All compounds were built using the HyperChem 8.0 release program
and optimized geometrically using the Gaussian 09 program, revision A.02 (Gaussian Inc., Wallingford,
CT, USA) at DFT B3LYP/3-21G level of theory. The protein and ligands were further prepared using
the utilities implemented by AutoDockTools 1.5.4 (http://mgltools.scripps.edu/). The protein was
adding polar hydrogen atoms, Kollman united-atom partial charges, and to the ligands computing
Gasteiger–Marsilli formalism charges, rotatable groups which were assigned automatically as were the
active torsions. Blind docking was carried out using AutoDock Vina version 2.0 [34]. The root mean
square deviation (RMSD) values were obtained by comparing the best pose generated in AutoDock.
The initial parameters used for the active site were a grid box size was 42 Å
and z dimensions and grid center 9.73, 18.00, 971 to 1G7F.PDB. For the allosteric site docking were
a grid box size was 42 Å 40 Å 40 Å in the x, y and z dimensions and grid center 9.73, 18.00, 971
×
40 Å
×
40 Å in the x, y
×
×
used 1T49.pdb. For both sites, the exhaustiveness was 25, and the ten best poses were obtained. The
analysis of the docking was made with PyMol (Maestro, Schrödinger, LLC, New York, NY) [35].
3
.11. Molecular Dynamics Simulation
All the structural complexes were verified, cleaned and ordered with the pdb4amber scrip before
starting the preparation in order to generate suitable topologies from the LEaP module of AMBER
9 [36 37]. Each structure and complex was subjected to the following protocol: hydrogens and
other missing atoms were added using the LEaP module with the leaprc.protein.ff19SB parameter
1
,
−
+
set; Cl or K counterions were added to neutralize the system; the complexes were then solvated
in an octahedral box of explicit TIP3P model water molecules localizing the box limits at 12 Å from
the protein surface. Molecular dynamic simulations were performed at 1 atm and 315 K, maintained

Molecules 2020, 25, 3530
15 of 18
with the Berendsen barostat and thermostat, using periodic boundary conditions and particle mesh
Ewald sums (grid spacing of 1 Å) for treating long-range electrostatic interactions with a 10 Å cutoff for
computing direct interactions. The SHAKE algorithm was used to satisfy bond constraints, allowing
the employment of a two fs time step for the integration of Newton’s equations as recommended
in the Amber package [36
for all residues. All calculations were made using Graphics Processing Units (GPU) accelerated MD
engine in AMBER (pmemd.cuda), a program package that runs entirely on CUDA (Compute Unified
,38]. Amber leaprc.protein.ff19SB force field [39] parameters were used
®
Device Architecture)-enabled GPUs [40]. The protocol consisted of performing a minimization of the
initial structure, followed by 50ps heating and pressure equilibration at 315 K and 1.0 atm pressure,
respectively. Finally, the system is equilibrated with 500ps before starting the production of MD.
The production of the MD consisted of 100 ns for each complex. Frames were saved at ten ps intervals
for subsequent analysis. All analyses were done using CPPTRAJ [41] part of AMBER19 utilities and
OriginPro 9.1. The calculations of RMSD and Root Mean Square Fluctuations (RMSF) were made,
considering the C, CA and N. The charts were built with OriginPro 2018 SR1, and the trends were
adjusted with the function processing smooth (method lowess span). VMD and PyMOL [35] were
used to visualize and create the images from the MD.
3
.12. Chemoinformatic Properties of Compounds 1, 7 and AU
The biological and chemoinformatic properties of compounds 1, 7 and AU were explored using
the servers Swiss TargetPrediction (http://www.swisstargetprediction.ch/index.php), Molinspiration
http://www.molinspiration.com/cgi-bin/properties) and Osiris Property Explorer server (http://che
minformatics.ch/propertyExplorer/) [29,42].
(
4
. Conclusions
The AE of S. amarissima contains rutin (
4) and rosmarinic acid (5), which inhibit intestinal glucose
absorption, promote glucose uptake in muscle cells and suppress insulin-resistance, among other
effects. On the other hand, pedalitin (7), which behave in silico as an allosteric inhibitor of PTP-1B,
could contribute to its overall antidiabetic action via α-glucosidase and selective PTP-1B inhibition,
and other mechanisms yet to be determined. The overall action of AE could be attained via network
pharmacology, synergism and or polypharmacology. Altogether, our studies on S. amarissima tend to
support its medicinal use for the treatment of diabetes in Mexican folk medicine. The chromatographic
analyses developed and validated in this study will allow the development of a pharmacopeic
monograph, and generate standardized preparations of this very Mexican plant. The analytical UHPLC
method was suitable for its intended purpose, the quantification of amarisolide (1) according to the Q2
(R1) guideline. Overall, the scientific information generated for this plant will contribute to its rational
use in Mexican folk medicine. Like many other New World melittophilous Salvia species, S. amarissima
is a rich source of bioactive compounds.
Supplementary Materials: The following are available online. NMR, IR and MS spectrum of the isolated or
identified compounds, TIC chromatograms of the essential oils and Residues of interaction to 4 Å of compounds 1,
7
and UA at the catalytic and allosteric sites of PTP-1B are available online.
Author Contributions: This work is part of the PhD thesis of E.S.-A. who performed the experiments; analysis of
spectral data, A.P.-V.; UHPLC-MS method validation, I.R.-C.; docking analysis, M.R.-G.; expression and purification
of hPTP-1B, docking, molecular dynamics and parameters calculations, M.G.-A.; botanical characterization of the
species, R.T.-C.; conception, funding acquisition, structure elucidation, writing and editing, R.M. All authors have
read and agreed to the published version of the manuscript.
Funding: The research was supported by grants from CONACyT CB A1-S-11226 and DGAPA IN 217320 awarded
to R.M.

Molecules 2020, 25, 3530
16 of 18
Acknowledgments: The authors recognize the valuable support of M.I. Vel
á
squez-L
ó
pez from Facultad de
Medicina; and Georgina Duarte from Facultad de Química. We are indebted to Biol. Itzi Fragoso and Martha
Mart
We also recognize the support of Direcci
DGTIC, UNAM) for the resources to carry out computational calculations through the Miztli supercomputing
í
nez Gordillo (Facultad de Ciencias, UNAM) for their valuable help in the identification of the plant material.
ó
n General de C mputo y de Tecnolog as de Informaci n y Comunicaci
ó
í
ó
ón
(
system (LANCAD-UNAM-DGTIC-313). E.S.-A. acknowledges the fellowship from CONACyT (289212) to pursue
graduate studies.
Conflicts of Interest: The authors declare no conflict of interest.
References
1
.
International Diabetes Federation (IDF). Diabetes Atlas, 9th ed. Available online: http://www.diabetesatlas.org
accessed on 13 March 2020).
(
2
.
Mata, R.; Figueroa, M.; Navarrete, A.; Rivero-Cruz, I. Chemistry and Biology of Selected Mexican Medicinal
Plants. In Progress in the Chemistry of Organic Natural Products; Kinghorn, A.D., Falk, H., Gibbons, S.,
Kobayashi, J., Asakawa, Y., Liu, J.K., Eds.; Springer Nature Switzerland AG: Cham, Switzerland, 2019;
Volume 108, pp. 1–142.
3
4
.
.
Efferth, T.; Greten, H.J. Quality control for medicinal plants. Med. Aromat. Plants 2012, 1, 1000e131.
Calder
Comisi
pp. 1–1406.
pez, F.C.E.; S
ó
ó
n, R.G.; Rzedowski, J. Flora Fanerog
n Nacional para el Conocimiento y Uso de la Biodiversidad, P
á
mica del Valle de M
é
xico, 2nd ed.; Instituto de Ecolog
í
a, A.C. y
átzcuaro: Mexico City, Mexico, 2010;
5
6
7
.
.
.
L
ó
ánchez, D.G.; Arrieta, B.D.; Rom
án, G.J.H. Estudio preliminar fitoquímico y de la actividad
antimicrobiana de Salvia amarissima Ortega. Cienc. y Tecnol. 2010, 9, 67–76.
Maldonado, E.; C rdenas, J.; Boj rquez, H.; Escamilla, E.M.; Ortega, A. Amarisolide, a neo-clerodane
á
ó
diterpene glycoside from Salvia amarissima. Phytochemistry 1996, 42, 1105–1108. [CrossRef]
Bautista, E.; Fragoso-Serrano, M.; Ortiz-Pastrana, N.; Toscano, R.A.; Ortega, A. Structural elucidation and
evaluation of multidrug-resistance modulatory capability of amarissinins A-C, diterpenes derived from
Salvia amarissima. Fitoterapia 2016, 114, 1–6. [CrossRef] [PubMed]
8
.
.
Flores-Bocanegra, L.; Gonz
Salvia circinata. J. Nat. Prod. 2017, 80, 1584–1593. [CrossRef] [PubMed]
Fragoso-Serrano, M.; Ortiz-Pastrana, N.; Luna-Cruz, N.; Toscano, R.A.; Alpuche-Solis, A.G.; Ortega, A.;
Bautista, E. Amarisolide F, an acylated diterpenoid glucoside and related terpenoids from Salvia amarissima.
J. Nat. Prod. 2019, 82, 631–635. [CrossRef]
ález-Andrade, M.; Bye, R.; Linares, E.; Mata, R. α-Glucosidase inhibitors from
9
1
0. Bautista, E.; Fragoso-Serrano, M.; Toscano, R.A.; García-Peña, M.R.; Ortega, A. Teotihuacanin, a diterpene
with an unusual spiro-10/6 system from Salvia amarissima with potent modulatory activity of multidrug
resistance in cancer cells. Org. Lett. 2015, 17, 3280–3282. [CrossRef]
1
1. Calzada, F.; Bautista, E. Plants used for the treatment of diarrhoea from Mexican Flora with amoebicidal and
giadicidal activity, and their phytochemical constituents. J. Ethnopharmacol. 2020, 253, 1–49. [CrossRef]
1
2. Moreno-P
é
rez, G.F.; Gonz
ález-Trujano, M.E.; Martínez-Gordillo, M.J.; San Miguel-Chávez, R.;
Basurto-Peña, F.A.; Dorazco-Goz
á
lez, A.; AguirreHern ndez, E. Amarisolide A and pedalitin as bioactive
á
compounds in the antinociceptive effects of Salvia circinata. Bot. Sci. 2019, 97, 355–365. [CrossRef]
3. Klaman, L.D.; Boss, O.; Peroni, O.D.; Kim, J.K.; Martino, J.L.; Zabolotny, J.M.; Moghal, N.; Lubkin, M.;
Kim, Y.B.; Sharpe, A.H. Increased energy expenditure, decreased adiposity, and tissue-specific insulin
sensitivity in protein-tyrosine phosphatase 1B-deficient mice. Mol. Cell. Biol. 2000, 20, 5479–5489. [CrossRef]
1
1
1
4. Voirin, B.;
Viricel, M.R.;
Favre-Bonvin, J.;
Van den Broucke, C.O.;
Lemli, J.
0
0
5,6,4 -Trihydroxy-7,3 -dimethoxyflavone and other methoxylated flavonoids isolated from Thymus
satureiodes. Planta Med. 1985, 51, 523–525. [CrossRef] [PubMed]
5. Airoldi, C.; Sironi, E.; Dias, C.; Marcelo, F.; Martins, A.; Rauter, A.P.; Nicotra, F.; Jimenez-Barbero, J. Natural
compounds against Alzheimer’s disease: Molecular recognition of A
extract and its major component, rosmarinic acid, as investigated by NMR. Chem. Asian, J. 2013, 8, 596–602.
CrossRef] [PubMed]
β1-42 peptide by Salvia sclareoides
[
1
6. Bagci, E.; Kocak, A. Essential oil composition of the aerial parts of two Salvia, L. (S. multicaulis Vahl. Enum
and S. tricochlada Bentham) species from east Anatolian region (Turkey). Int. J. Sci. Tech. 2008, 3, 13–18.

Molecules 2020, 25, 3530
17 of 18
1
1
1
7. Salehi, P.; Tolami, L.B.; Sefidkon, F. Essential oil composition of Salvia xanthocheila Boiss. ex Benth. from Iran.
J. Essent. Oil Res. 2005, 17, 442–443. [CrossRef]
8. Sonboli, A.; Babakhani, B.; Mehrabian, A.R. Antimicrobial activity of six constituents of the essential oil from
Salvia. Z. Naturforsch. 2006, 61, 160–164. [CrossRef]
9. Zhao, B.T.; Le, D.D.; Nguyen, P.H.; Ali, M.Y.; Choi, J.S.; Min, B.S.; Shin, H.M.; Rhee, H.I.; Woo, M.H. PTP1B,
α
glucosidase, and DPP-IV inhibitory effects for chromene derivatives from the leaves of Smilax china L. Chem.
Biol. Interact. 2016, 253, 27–37. [CrossRef]
2
0. Xu, J.; Wang, X.; Yue, J.; Sun, Y.; Zhang, X.; Zhao, Y. Polyphenols from a corn leaves (Quercus liaotungensis)
protect pancreatic beta cells and their inhibitory activity against
α-glucosidase and protein tyrosine
phosphatase 1B. Molecules 2018, 23, 2167. [CrossRef]
2
2
1. Ngoa, Y.L.; Lau, C.H.; Chua, L.S. Review on rosmarinic acid extraction, fractionation and its anti-diabetic
potential. Food Chem. Toxicol. 2018, 121, 687–700. [CrossRef]
2. Oh, J.; Choe, S.; Kim, D.W.; Lee, S.H.; Na, M.K. Rhododendric acid A, a new ursane-type PTP1B inhibitor
from the endangered plant Rhododendron brachycarpum G. Don. Bioorg. Med. Chem. Lett. 2012, 22,
6
116–6119.
2
2
3. Al-Ishaq, R.K.; Abotaleb, M.; Kubatka, P.; Kajo, K.; Busselberg, D. Flavonoids and their anti-diabetic effects:
Cellular mechanisms and effects to improve blood sugar levels. Biomolecules 2019, 9, 430. [CrossRef]
4. Wiesmann, C.; Barr, J.J.; Kung, J.; Zhu, J.; Erlanson, D.A.; Shen, W.; Fahr, B.J.; Zhong, M.; Taylor, L.; Randal, M.;
et al. Allosteric Inhibition of Protein Tyrosine Phosphatase 1B. Nat. Struct. Mol. Biol. 2004, 11, 730–737.
[CrossRef] [PubMed]
2
2
2
2
5. Shrestha, S.; Seong, S.H.; Park, S.G.; Min, B.S.; Jung, A.H.; Chi, J.S. Insight into the PTP1B inhibitory activity
of arylbenzofurans: An in vitro and in silico study. Molecules 2019, 24, 2893. [CrossRef] [PubMed]
6. Gohlke, H.; Klebe, G. Approaches to the description and prediction of the binding affinity of small-molecule
ligands to macromolecular receptors. Angew. Chem. Int. Ed. Engl. 2002, 41, 2644–2676. [CrossRef]
7. Ramírez, D.; Caballero, J. Is it reliable to take the molecular docking top scoring position as the best solution
without considering available structural data? Molecules 2018, 23, 1038. [CrossRef] [PubMed]
8. Zhang, X.-Y.; LI, W.; Wang, J.; Li, N.; Ccheng, M.S.; Koike, K. Protein tyrosine phosphatase 1B inhibitory
activities of ursane-type triterpenes from Chinese raspberry, fruits of Rubus chingii. Chin. J. Nat. Med. 2019
7, 15–21.
,
1
2
3
9. Khan, T.; Lawrence, A.J.; Azad, I.; Raza, S.; Joshi, S.; Khan, A.R. Computational drug designing and prediction
of important parameters using in silico methods- A review. Curr. Comput. Aided Drug. Des. 2019, 15, 384–397.
[CrossRef]
0. International Conference on Harmonization (ICH) Q2(R1): Validation of Analytical Procedures: Text
and Methodology. In ICH: International Council for Harmonisation of Technical Requirements for
Pharmaceuticals for Human Use. Geneva, Switzerland, 1 November 2005. Available online: http:
//academy.gmp-compliance.org/guidemgr/files/Q2(R1).pdf (accessed on 17 June 2020).
3
3
3
1. Van den Dool, H.; Kratz, P.D. Generalization of the retention index system including linear temperature
programmed Gas-Liquid partition Chromatography. J. Chromatogr. A. 1963, 11, 463–471. [CrossRef]
2. Adams, R.P. Identification of Essential Oil Components by Gas Chromatography/Mass Spectrometry, 4th ed.; Allured
Publishing Corporation: Carol Stream, IL, USA, 2007; pp. 1–804.
3. Rangel-Grimaldo, M.; Mac
-Glucosidase and protein tyrosine phosphatase 1B inhibitors from Malbranchea circinata. J. Nat. Prod. 2020
3, 675–683. [CrossRef]
ías-Rubalcava, M.L.; González-Andrade, M.; Raja, H.; Figueroa, M.; Mata, R.
α
,
8
3
3
3
4. 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, 31, 455–461. [CrossRef]
5. DeLano, W.L. Use of PYMOL as a communications tool for molecular science. Abstr. Pap. A.C.S. 2004, 228,
U313–U314.
6. Case, D.A.; Betz, R.M.; Cerutti, D.S.; Cheatham, T.E., III; Darden, T.A.; Duke, R.E.; Giese, T.J.; Gohlke, H.;
Goetz, A.W.; Homeyer, N.; et al. AMBER. 2015. Available online: http://ambermd.org/doc12/Amber16.pdf
(accessed on 31 July 2020).
3
7. Case, D.A.; Cheatham, T.E., III; Darden, T.; Gohlke, H.; Luo, R.; Merz, K.M., Jr.; Onufriev, A.; Simmerling, C.;
Wang, B.; Woods, R.J. The Amber biomolecular simulation programs. J. Comp. Chem. 2005, 26, 1668.
[CrossRef] [PubMed]

Molecules 2020, 25, 3530
18 of 18
3
3
8. Walker, R.C.; Crowley, M.F.; Case, D.A. The implementation of a fast and accurate QM/MM potential method
in Amber. J. Comp. Chem. 2008, 29, 1019. [CrossRef] [PubMed]
9. Tian, C.; Kasavajhala, K.; Kellon, A.; Belfon, K.; Raguette, L.; Huang, H.; Migues, A.N.; Bickel, J.; Wang, Y.;
Pincay, J.; et al. ff19SB: Amino-acid specific protein backbone parameters trained against quantum mechanics
energy surfaces in solution. J. Chem. Theory Comput. 2020, 16, 528–552. [CrossRef] [PubMed]
0. Salomon-Ferrer, R.; Gotz, A.W.; Poole, D.; Le Grand, S.; Walker, R.C. Routine microsecond molecular
dynamics simulations with AMBER on GPUs. 2. Explicit solvent particle mesh Ewald. J. Chem. Theory
Comput. 2013, 9, 3878–3888. [CrossRef] [PubMed]
4
4
4
1. Roe, D.R.; Cheatham, T.E., III. PTRAJ and CPPTRAJ: Software for processing and analysis of molecular
dynamics trajectory data. Cheatham. J. Chem. Theory Comput. 2013, 9, 3084–3095. [CrossRef] [PubMed]
2. Bhal, S.K. LogP-making sense of the value; Advanced Chemistry Development, Inc. (ACD/LABS): Toronto, ON,
Canada, 2007; pp. 1–4.
Sample Availability: Samples of the compounds are not available from the authors.
©
2020 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access
article distributed under the terms and conditions of the Creative Commons Attribution
(CC BY) license (http://creativecommons.org/licenses/by/4.0/).