Abstract
Members of the Viola genus play important roles in traditional Asian herbal medicine. This study investigates the ability of Viola odorata L. extracts to inhibit Na+,K+-ATPase, an essential animal enzyme responsible for membrane potential maintenance. The root extract of V. odorata strongly inhibited Na+,K+-ATPase, while leaf and seeds extracts were basically inactive. A UHPLC-QTOF-MS/MS metabolomic approach was used to identify the chemical principle of the root extract’s activity, resulting in the detection of 35,292 features. Candidate active compounds were selected by correlating feature area with inhibitory activity in 14 isolated fractions. This yielded a set of 15 candidate compounds, of which 14 were preliminarily identified as procyanidins. Commercially available procyanidins (B1, B2, B3 and C1) were therefore purchased and their ability to inhibit Na+,K+-ATPase was investigated. Dimeric procyanidins B1, B2 and B3 were found to be inactive, but the trimeric procyanidin C1 strongly inhibited Na+,K+-ATPase with an IC50 of 4.5 µM. This newly discovered inhibitor was docked into crystal structures mimicking the Na3E1∼P·ADP and K2E2·Pi states to identify potential interaction sites within Na+,K+-ATPase. Possible binding mechanisms and the principle responsible for the observed root extract activity are discussed.
Similar content being viewed by others
Introduction
Na+,K+-ATPase (NKA; EC 7.2.2.13) is a transmembrane enzyme found in every animal cell. It is responsible for maintaining the plasma membrane potential and sodium gradient necessary for the function of many secondary active transporters1,2. Its importance in clinical practice is demonstrated by the frequent prescription of NKA inhibitors such as digoxin, which is used against cardiac insufficiency3. Moreover, rare mutations of NKA cause severe diseases including familial hemiplegic migraine type 2 and rapid-onset dystonia-parkinsonism4,5. Two hundred and thirty-six years after the first publication of William Withering’s book An Account of the Foxglove, and Some of Its Medical Uses, which describes medical uses of digitalis, the pharmacology of NKA remains interesting but incompletely understood. In addition to its role in heart failure management, recent high-throughput screening efforts suggest that NKA inhibition may have positive effects on hypercholesterolemia, several types of cancer, and viral infections6,7,8,9,10,11,12,13,14,15.
NKA consists of a ~ 112 kDa catalytic α subunit with ten transmembrane segments, and a ~ 55 kDa glycosylated β subunit with one transmembrane segment and three isoforms (β1-3) that acts as a chaperone and modulatory protein1,16,17. The catalytically active αβ heterodimer can also associate with one of the regulatory and stabilizing FXYD proteins18,19.
The catalytic cycle of NKA involves alternating-access ion transport and a series of large conformational changes20,21,22. As a translocase, NKA hydrolyzes ATP to transfer three Na+ ions to the extracellular space and two K+ ions to the cytoplasm23. The resulting ion gradient across the plasma membrane is central to the physiology of excitable cells.
Other important functions of NKA include thermogenesis24 and promotion of cell adhesion and interaction via its β1 subunit25,26,27. It was also shown that the β1 subunit can act as a general K+-dependent lectin28, and that its expression suppresses the motility and invasiveness of carcinoma cells29. In addition to these functions, NKA is part of a complex signaling network that is currently under investigation.
In signal transduction pathways, NKA acts as a receptor for cardiotonic steroids (CTSs), which are highly selective inhibitors that bind in a cavity accessible from the extracellular side30. Many studies have shown that upon CTS binding, NKA interacts with several downstream proteins via its intracellular side and triggers a signaling cascade31,32. The first report describing an endogenous CTS was published in 199133, and endogenous CTS are now regarded as a distinct class of hormones34,35. However, uncertainties remain concerning their biosynthesis and distribution36,37.
Viola is the largest genus of the Violaceae family, containing up to 600 species38. Members of this genus are used in traditional herbal medicine in many regions: V. odorata is used in Pakistan39, V. betonicifolia and V. canescens in the Himalayas40,41, V. hondoensis in Korea42, and V. yedoensis, V. kunawarensis, and V. tianschanica in Central and East Asia43,44.
Sweet violet (Viola odorata) is traditionally used against several diseases, including neurological disorders and hypertension45. Volatile components of V. odorata leaf extracts include nona-2,6-dienal and (Z)-hex-3-enal, which are the compounds primarily responsible for its odor. Most of the other volatile components are saturated or unsaturated aliphatic hydrocarbons and related oxidized species46. Another study identified phthalate, tetrahydrobenzofuranone derivative, and monoterpenoids as the main components of leaf extracts of an Iranian chemotype of V. odorata47. In addition, cyclic peptides called cyclotides (e.g. cycloviolacins) with diverse biological activities are found in various sweet violet organs48,49. The phytochemistry of the related species V. betonicifolia has been reviewed in detail50.
There is relatively little published information on the secondary metabolite content of root extracts of V. odorata, which are the subject of this work. Phenolic acids and the flavonoid compounds vitexin and rutin were quantified in the roots of V. tricolor during a study on mycorrhizal colonization that demonstrated the influence of mycorrhizal fungi on secondary metabolite levels51. Similar endophyte-induced changes in phytochemical composition are well documented52,53. In V. odorata, the diversity of endophytes in the roots was higher than in other organs54.
Procyanidins are polyphenolic compounds consisting of condensed flavan-3-ol subunits. Monomers of (–)-epicatechin or ( +)-catechin form oligomeric or polymeric structures that are classified based on the nature of their interflavan bonds. Specifically, they are categorized based on the stereochemistry of the interflavan bond (α or β) and the atoms linked by the bond: bonds between atoms C-4 and C-8 or C-4 and C-6 are classified as B-type linkages, while those between C-4 and C-8 along with C-2 → O → C-7 ether bond are A-type linkages55.
Here, we investigate the NKA-inhibiting activity of Viola odorata L. extracts and their components. Active compounds are identified by LC–MS correlation metabolomics56.
Results
Over 50 plant species were extracted and screened for NKA inhibiting activity (Supplementary Table S1 online). Individual parts of the most promising candidate, Viola odorata L., were then analyzed separately to evaluate their activity (Fig. 1A). Interestingly, although seed and leaf extracts showed no significant activity, root extracts inhibited NKA by over 90% at a concentration of 600 µg DW/mL. Therefore, the root extract was further investigated by separating it into fractions using C18, mixed-mode anion, and cation exchange SPE columns. This separatory process yielded 14 fractions, of which two (RP2 and MC2) were significantly active (Fig. 1B). A recently published correlation-metabolomic approach was then applied to these fractions to identify candidate features possibly responsible for the observed NKA-inhibiting activity56.
All isolated fractions were analyzed by non-targeted UHPLC-QTOF-MS/MS, which revealed 15,331 features in negative and 19,961 features in positive ion modes. Combined features from both modes were filtered based on their abundance and retention time, resulting in 2,748 features in total. These features were correlated with the inhibitory activity of each fraction and sorted according to the calculated correlation coefficients. The dimensionality of the data was further reduced by manual identification of the molecular ions of adducts, fragments, and isotopic peaks. Data from positive and negative modes were processed separately, which typically resulted in the presence of two rows independently referring to the same metabolite (see for example rows 1 and 2 in Table 1). The highly correlated candidate metabolites are shown in Table 1. Briefly, the higher the r value, the greater the likelihood that the row represents the active component. The majority of the candidate metabolites were annotated as B-type procyanidins, with the exception of one that was putatively identified as feruloyl putrescine. The commercially available procyanidins B1, B2, B3 and C1 were purchased and analyzed, allowing two candidates to be unambiguously identified as procyanidins C1 and B2. The NKA inhibitory activities of all available procyanidins (B1, B2, B3 and C1) were also determined (Fig. 2A; 2B), revealing that procyanidin C1 is a strong inhibitor with an IC50 of 4.5 ± 0.8 µM (Fig. 2C).
Due to the complexity of the V. odorata extracts, it was not possible to structurally identify the other B-type procyanidins from the candidate list. The diversity of B-type procyanidins in Viola odorata root is illustrated by the number of peaks with m/z 865 showing fragmentation patterns characteristic of B-type procyanidin trimers (Fig. 3). The procyanidin C1 (7.25 min) is a minor peak among detected procyanidin trimers in the V. odorata extract (Fig. 3B).
Molecular docking of procyanidin C1 to NKA was performed to identify possible interaction sites. Two crystal structures of NKA representing different catalytic intermediates were selected for the docking study: 3WGU (Na3E1[AlF4]–·ADP), which has a resolution of 2.80 Å and is an analogue of the Na3E1∼P·ADP state, and 2ZXE (K2E2·[MgF4]2–), which has a resolution of 2.40 Å and is analogous to the K2E2·Pi state.
Docking of procyanidin C1 into the structural analogue of the NKA of Na3E1∼P·ADP state (3WGU) yielded lower binding energies than those for docking into the K2E2·Pi state analogue. The binding site with the highest affinity for procyanidin C1 (binding energy: –11.5 kcal/mol) was located in close proximity to the nucleotide binding site (Fig. 4, left). For the K2E2·Pi state analogue (2ZXE), the lowest energy binding site of procyanidin C1 (binding energy: –9.9 kcal/mol) was in the extracellular part of the protein (Fig. 4, right). In this case, key protein–ligand interactions included two cation–π bonds between arginine guanidinium groups and the ortho-dihydroxyphenyl rings of procyanidin C1.
Discussion
We used a modern approach based on correlation metabolomics to identify active compounds in a plant extract. This procedure led us to a few candidate compounds and ultimately to the root extract constituent procyanidin C1, which strongly inhibited ATPase activity of NKA (IC50 4.5 µM). In contrast to the strong biological activity of the trimer procyanidin C1, commercially available procyanidin dimers (B1, B2 and B3) were inactive (IC50 > 100 µM). The NKA-inhibiting activity of procyanidin C1 was previously unknown, although a procyanidin-rich hawthorn extract WS® 1442 with a procyanidin content of 17.3–20.1% was previously shown to inhibit NKA58,59,60. Potential antiarrhythmic and positive ionotropic effects of WS® 1442 were also reported62. However, it remains unclear whether procyanidins or other components of hawthorn extract are primarily responsible for the activity of WS® 1442. For example, ursolic and oleanolic acids, which constitute roughly 0.6% of WS® 1442, are also moderate inhibitors of NKA, with IC50 values of 25 to 100 µM61,62. Our data on procyanidin C1 suggest that its strong NKA-inhibitory activity together with its high abundance in hawthorn extract may explain the cardiac glycoside-like properties of WS® 144263,64,65,66,67,68,69. Procyanidin C1 was also shown to inhibit gastric H+,K+-ATPase, which is another P-type ATPase with an E1–E2 catalytic mechanism like NKA. The gastric H+,K+-ATPase was inhibited by an extract from the aerial part of Cecropia glaziovii Snethl. as well as by its chemical constituents procyanidins C1 and B2, whose IC50’s were determined to be 46.5 and 40.6 µM, respectively70. Due to the presence of other phenolics and the relatively low activity of the identified constituents, the chemical principle of H+,K+-ATPase inhibition by Cecropia glaziovii extract remains unclear.
In this work we investigated the only commercially available trimer, procyanidin C1, which strongly inhibited NKA. Despite its strong activity, procyanidin C1 was present at relatively low concentrations in V. odorata when compared to other detected trimers (Fig. 3B),suggesting that they may also contribute to the extract’s activity. It should be noted that the top scoring metabolite given by correlation approach was an unknown trimeric B–type procyanidin, the isomer of procyanidin C1. Due to their similar chemical structures, the procyanidins inevitably coeluted from the SPE cartridges used in this work (C18, MAX and MCX). Therefore, their individual contributions to the overall extract activity are not readily distinguished by correlation alone. However, based on our metabolomic and biochemical data, we suggest that procyanidin C1 is at least partially responsible for the observed activity of V. odorata root extracts, acting most likely in combination with its isomer(s). Contributions of extract components other than procyanidin metabolites is unlikely because the only non-procyanidin compound correlated with NKA activity was annotated as feruloyl putrescine, a mono-conjugated phenolamide present in many species that have never been reported to inhibit NKA, including Salsola subaphylla, Zea mays, Oryza sativa, and Nicotiana tabacum71,72,73. It should also be noted that the presence of procyanidins in the root, the only active V. odorata part, might be associated with their role in the nitrogen cycle in the soil74,75.
In contrast to the strongly active trimeric procyanidin C1, the dimeric procyanidins B1, B2 and B3 were completely inactive towards NKA. Possible binding modes of procyanidin C1 were identified by molecular docking using the 3WGU and 2ZXE crystal structures of NKA proteins originating from Sus scrofa and Squalus acanthias, respectively. Both of these structures are available at considerably higher resolutions than other published NKA structures and were therefore the only Na3E1∼P·ADP and K2E2·Pi state analogues investigated in this work. Comparing the docking results for procyanidin B1, B2, B3, and C1, there is no specific binding site for procyanidin C1 exclusively. Therefore, the three catechin units in procyanidin C1 might be necessary to fulfil spatial requirements for the inhibitory effect. Two interesting bifacial cation–π interactions between arginine guanidinium groups and the ortho-dihydroxyphenyl rings of procyanidin C1 were found in the structure mimicking K2E2·Pi state. This interaction causes the ion pathway (which is formed by residues including one of the interacting arginines, R979) to become inaccessible from the extracellular side76. This binding mode may thus sterically impede ion exchange. It should be noted that residues R979 in loop L9-10 and D128 in loop L1-2 form a salt bridge in the E2P state, but separation of these residues is believed to be required for the motion of the TM2 transmembrane helix during the catalytic cycle77,78. Other residues from L1-2 are also important for the protein’s translocase function, so their interactions with the ligand could have additional effects on its conformational transitions77. In the best docking pose of the Na3E1∼P·ADP-mimicking structure (3WGU), the ligand plugged the whole hydrophilic cavity in front of the nucleotide-binding site in the cleft between the nucleotide-binding (N) and actuator (A) domains. In this case, the binding energy (–11.5 kcal/mol) was lower than that for the previously discussed 2ZXE structure (–9.9 kcal/mol). Due to the interdomain location of the binding site, residues from all three cytoplasmic domains are involved in procyanidin C1 binding. The interdomain space in the intracellular part of the protein was previously suggested to be a possible binding site for flavonolignans based on a docking study79. A subsequent fluorescence spectroscopy experiment proved that flavonolignans interact with the cytoplasmic segment connecting transmembrane helices TM4 and TM5, providing experimental evidence for a binding mode that may contribute to the inhibition of NKA by flavonolignans79. Therefore, our identification of a similarly located binding pose may help explain the inhibitory activity of procyanidin C1.
Procyanidins are highly abundant in human diet, with an estimated daily intake of 57.7 mg per person80. On the basis of cell monolayer penetration studies, they are believed to undergo paracellular absorption from the intestine81,82,83. In vivo animal and human studies have also shown that procyanidin dimers and trimers are stable under gastric conditions, are not degraded into monomers, and can be detected in rat plasma and urine after ingestion84,85,86,86,87,88,88,89,90,91. While quantitative data on the plasma concentrations of procyanidins in humans are limited, procyanidin B2 was detected in plasma at a concentration of 4.0 ± 0.6 nM after the consumption of 1.8 mg dimeric procyanidins per kg of body weight. Interestingly, treatment of the plasma samples with sulfatases and β-glucuronidases did not increase the measured procyanidin B2 concentration, suggesting that it was not conjugated90. Procyanidins were also detected in human plasma from individuals who had recently consumed procyanidin-rich plant foods/extracts. For example, administration of 0.375 g of cocoa per kg of body weight and 2 g of proanthocyanidin-rich grape seed extract per person resulted in the detection of procyanidins B2 at 41 ± 4 nM and B1 at 10.6 ± 2.5 nM in human plasma92,93. Ingested procyanidins thus seem to be readily absorbed in humans, resulting in detectable levels in human plasma. As such, they are both dietary compounds capable of inhibiting NKA and also potential sources of novel molecular structures that could be used as alternatives to cardiac glycosides in the treatment of congestive heart failure and cardiac arrhythmia.
Conclusions
Procyanidin C1 is a newly discovered NKA inhibitor whose molecular architecture could potentially be optimized to develop analogues with greater druglikeness. Its low micromolar IC50 makes all trimeric B-type procyanidins interesting targets for further mechanistic investigation and analysis of structure–activity relationships. Such studies could enable the development of novel NKA inhibitors with previously unexplored modes of action and/or binding sites. The successful application of correlation metabolomics in this work further demonstrates the high effectivity of this approach for identifying biologically active compounds in complex extracts.
Materials and methods
Chemicals
Procyanidin B1, B2, B3 and C1 standards (catalog numbers 89764, 89,552, 84,047 and 89,537, respectively) were purchased from PhytoLab (Germany). Stock solutions were prepared by dissolving procyanidins in 50% methanol to a final concentration of 10 mM and stored at –80 °C. The NKA inhibitor ouabain was purchased from Merck (Germany). All solutions were prepared using Milli-Q water with a resistivity above 18.2 MΩ·cm (25 °C).
Isolation of Na+,K+-ATPase
NKA was isolated using the protocol of Fedosova94 (Supplementary Fig. S1 online), with the following modifications: excess buffer I was poured out from the dissected outer medulla, and the soaked tissue was transferred into a blender. The outer medulla was then homogenized in buffer II with a tissue : buffer ratio of 1 : 3 while ¾ of buffer II was in the form of ice cubes. Aliquots of SDS-treated NKA isolate in buffer II were stored at –80 °C, and were determined to have a protein concentration of 0.44 mg/mL by the Lowry method.
Na+,K+-ATPase (NKA) activity assay
Authentic standards (procyanidins B1, B2, B3, and C1; ( +)-catechin and (–)-epicatechin) dissolved in 50% methanol with a stock concentration of 10 mM were diluted to a final concentration of 100 µM for use in assays. Crude extracts of Viola odorata L. seed, root, and leaf with concentration of 240 mg DW/mL in dimethyl sulfoxide were diluted to a final concentration of 600 µg DW/mL. Fractions of crude V. odorata L. root extract with stock concentrations of 240 mg DW/mL in 20% methanol with 0.1% HCOOH were diluted to a final concentration of 480 µg DW/mL.
Inorganic phosphate (Pi) quantification was conducted in microplates using a previously reported method79. Each ATPase reaction was performed in four wells without ouabain to determine total liberated Pi and in four wells with ouabain to quantify Pi not generated by ouabain-sensitive NKA activity.
The ouabain-sensitive ATPase activity of NKA (henceforth referred to as NKA activity) was measured by adding 30 µL of a 1.67 × concentrated buffered NKA solution to 10 µL of pre-concentrated inhibitor solution diluted in water. NKA was pre-incubated with the inhibitor for 5 min, then the ATPase reaction was started by adding 10 µL of 15 mM Na2ATP dissolved in 25 mM Tris base. The final buffered (pH 7.2) NKA solution contained 4 mM MgCl2, 20 mM KCl, 130 mM NaCl, 30 mM l-histidine, and 26 µg/mL of NKA protein with or without 1 mM ouabain. The ATPase reaction lasted for 6 min after which 75 µL of solution II was added to develop color (2.9% (w/v) sodium dodecyl sulfate, 0.5% (w/v) ammonium heptamolybdate tetrahydrate, 163 mM ascorbic acid, and 3.7% (w/v) HCl). After 8 min, color development was stopped by adding 125 µL of solution III (3.5% (w/v) bismuth citrate, 3.5% (w/v) trisodium citrate dihydrate, 3.7% (w/v) HCl). Finally, the absorbance was measured at 710 nm using an Infinite 200 microplate reader (Tecan, Switzerland).
The NKA activity observed in the presence of each inhibitor was normalized against a control activity determined by mixing the NKA solution with water instead of an inhibitor solution in order to obtain the relative NKA activity (%) for each inhibitor. NKA activity was calculated by subtracting the average absorbance of ouabain containing samples from that of samples without ouabain. Dose–response curves were generated by using a four-parameter logistic function (fixed 0–100% range) to fit the experimental data.
Plant collection and extraction
Viola odorata L., collected in Bělkovice Valley (Olomouc Region, Czech Republic), was identified by Michal Hroneš, Ph.D. (Department of Botany, Palacký University Olomouc). The collection complied with all applicable laws/guidelines, both institutional and national. According to IUCN guidelines, this taxon is not threatened being qualified for Least Concern category. As a common plant species, it’s collection for scientific purpose requires no permission/license. Plants were freeze-dried and homogenized using blade grinder. 30 mg of dried plant material (DW) was extracted in a microtube using 1 mL of 0.1% HCOOH in methanol. Glass balls were added to the microtube and the plant material was homogenized for 5 min using an MM 400 oscillatory ball mill (RETSCH, Germany) operating at 27 Hz. An ultrasound-assisted extraction (10 min, 45 kHz) in a USC100T ultrasonic bath (VWR, USA) and subsequent centrifugation (21,300 g, 25 °C, 10 min) followed. The supernatant was collected and evaporated at 37 °C in a TurboVap Classic LV nitrogen evaporator (Biotage, Sweden).
Plant extract fractionation
Methanol, aqueous ammonia and HCOOH used for fractionation were of LC–MS grade.
Reverse phase C18 SPE fractionation
A solid-phase extraction (SPE) C18 sorbent cartridge (Spe-ed 6 mL octadecyl/18; Applied Separations, USA) was conditioned with 12 mL of methanol and equilibrated with 2 mL of 0.1% HCOOH. Evaporated crude Viola odorata root extract (24 mg DW) dissolved in 2 mL of 0.1% HCOOH was loaded onto the column. Fractions RP0, RP1, RP2, RP3, RP4, RP5, RP6, and RP7 were obtained by eluting the crude extract from the sorbent with 2 mL of 0.1% HCOOH in aqueous methanol of gradually increasing concentration (0, 10, 20, 30, 50, 60, 70, and 100%).
Mixed-mode anion exchange SPE fractionation
A mixed-mode anion exchange (MAX) SPE sorbent cartridge (Oasis MAX, 150 mg, 6 mL; Waters, USA) was conditioned with 12 mL of methanol and equilibrated first with 6 mL of water followed by 6 mL of 5% aqueous ammonia. Evaporated crude V. odorata root extract (24 mg DW) dissolved in 5% aqueous ammonia was loaded onto the column and washed with 5 mL of the loading solvent (eluting fraction MA0). Fractions obtained by elution with 4 mL of methanol (MA1) and 4 mL of 2% methanolic HCOOH (MA2) were then collected.
Mixed-mode cation exchange SPE fractionation
A mixed-mode cation exchange (MCX) SPE sorbent cartridge (Oasis MCX, 150 mg, 6 mL; Waters, USA) was conditioned with 12 mL of methanol and equilibrated firstly with 6 mL of water, then with 6 mL of 2% aqueous HCOOH. Evaporated crude V. odorata root extract (24 mg DW) dissolved in 2% aqueous HCOOH was loaded onto the column and washed with 5 mL of loading solvent (eluting fraction MC0). Fractions obtained by elution with 4 mL of methanol (MC1) and 4 mL of 5% methanolic ammonia (MC2) were collected.
Each collected fraction was evaporated at 37 °C in a nitrogen evaporator and then redissolved in 100 µL of 20% methanol with 0.1% HCOOH. Aliquots were taken for use in LC–MS analyses and NKA inhibition assays.
UHPLC-QTOF-MS/MS
UHPLC-QTOF-MS/MS analyses were performed as described previously56. Briefly, samples dissolved in 20% methanol with 0.1% HCOOH were filtered through a 0.2 μm regenerated cellulose membrane microfilter (Grace, USA) prior to UHPLC analysis. 5 µL of each sample was injected onto a reversed-phase column (BEH C18, 1.7 μm, 2.1 × 150 mm, Waters, USA) maintained at 30 °C. Elution was performed over 22 min using the following binary gradient of acetonitrile (A) and 5 mM aqueous HCOOH (B): 0 min 5% A, 1 min 10% A, 12 min 35% A, 17 min 70% A, 17.5 min 100% A. The nebulizer gas pressure was set to 6 bar, desolvation temperature to 500 °C, and desolvation gas flow to 600 L/h. Analytes were ionized using an ESI source (120 °C, capillary 2 kV, cone voltage 25 V, cone gas flow 30 L/h) operating in both negative and positive ion modes. MS data were recorded in the m/z range of 70 − 1,500. MS/MS experiments were performed using a collision energy of 20 eV.
In-house developed MATLAB algorithms were applied to raw LC–MS data generated by MassLynx software to produce a list of features characterized by retention time, m/z, and peak area. Only features with retention times between 2 and 18 min and a peak area above 5,000 AU in at least one sample were considered. Pearson correlation coefficients were computed between feature area and inhibitory activity (expressed in µmol of ouabain equivalents per 1 g of DW). Features with the highest correlation coefficients were manually processed using Masslynx 4.1 software (Waters, USA) to identify candidate pseudo-molecular ions, their elemental compositions, and the corresponding mass accuracy. Procyanidins B1, B2, B3 and C1 were identified based on their retention times and comparisons to the fragmentation patterns of authentic standards.
Molecular docking
The structures of procyanidins B1, B2,B3 and C1 obtained from the PubChem database was optimized using Avogadro (Avogadro: an open-source molecular builder and visualization tool. Version 1.2. http://avogadro.cc/; https://doi.org/10.1186/1758-2946-4-17) using the energy minimizing tool with the General Amber Force Field. NKA crystal structures (PDB ID 3WGU, 2.80 Å109 and 2ZXE, 2.40 Å110) were downloaded from the Protein Data Bank (www.rcsb.org). Ligands and protein structures without waters, ions, and other ligands were prepared for docking using Autodock Tools (https://doi.org/10.1002/jcc.21256, version 1.5.6). Ligand bonds whose rotation generated different conformers were allowed to rotate freely in accordance with the software’s default settings. Docking was performed using Autodock Vina95, with a grid box covering the whole protein. Exhaustiveness was set to 100 and available 20 binding modes were obtained. Redocking was performed with grid box limited to frequently occupied binding sites.
Protein–ligand interactions of redocked binding modes were analyzed using the Protein–Ligand Interaction Profiler (PLIP)96. Residue numbering remained the same as in the original crystal structures. Figures for publication were prepared using PyMOL (The PyMOL Molecular Graphics System, Version 2.4.1, Schrödinger, LLC.).
Statistical analysis
All statistical analyses were performed using GraphPad Prism 8.0 software (GraphPad Software, Inc; USA). Data on NKA activity were analyzed by one-way ANOVA followed by Tukey’s post hoc test. Values are given as mean ± SD. Statistical significance is reported as **** p < 0.0001, *** p < 0.001, ** p < 0.01 and * p < 0.05.
Data availability
The datasets generated during and/or analysed during the current study are available from the corresponding author on reasonable request.
References
Kaplan, J. H. Biochemistry of Na, K-ATPase. Annu. Rev. Biochem. 71, 511–535 (2002).
Zdravkovic, I., Zhao, C., Lev, B., Cuervo, J. E. & Noskov, S. Y. Atomistic models of ion and solute transport by the sodium-dependent secondary active transporters. Biochim. Biophys. Acta Biomembr. 1818, 337–347 (2012).
Adams, K. F. et al. A perspective on re-evaluating digoxin’s role in the current management of patients with chronic systolic heart failure: Targeting serum concentration to reduce hospitalization and improve safety profile. Eur. J. Heart Fail. 16, 483–493 (2014).
De Carvalho Aguiar, P. et al. Mutations in the Na+/K+-ATPase α3 gene ATP1A3 are associated with rapid-onset dystonia parkinsonism. Neuron 43, 169–175 (2004).
Segall, L. et al. Alterations in the α2 isoform of Na, K-ATPase associated with familial hemiplegic migraine type 2. Proc. Natl. Acad. Sci. U. S. A. 102, 11106–11111 (2005).
Li, Z. et al. Dihydroouabain, a novel radiosensitizer for cervical cancer identified by automated high-throughput screening. Radiother. Oncol. 148, 21–29 (2020).
Wha Jun, D. et al. Ouabain, a cardiac glycoside, inhibits the fanconi anemia/BRCA pathway activated by DNA interstrand cross-linking agents. PLoS ONE 8, e75905 (2013).
Nilubol, N. et al. Four clinically utilized drugs were identified and validated for treatment of adrenocortical cancer using quantitative high-throughput screening. J. Transl. Med. 10, 198 (2012).
Guo, J. et al. Screening of natural extracts for inhibitors against Japanese encephalitis virus infection. Antimicrob. Agents Chemother. https://doi.org/10.1128/AAC.02373-19 (2020).
Prassas, I., Paliouras, M., Datti, A. & Diamandis, E. P. High-throughput screening identifies cardiac glycosides as potent inhibitors of human tissue kallikrein expression: Implications for cancer therapies. Clin. Cancer Res. 14, 5778–5784 (2008).
Rupaimoole, R., Yoon, B., Zhang, W. C., Adams, B. D. & Slack, F. J. A high-throughput small molecule screen identifies ouabain as synergistic with mir-34a in killing lung cancer cells. Iscience 23(2), 100878 (2020).
Zhang, L. et al. Quantitative high-throughput drug screening identifies novel classes of drugs with anticancer activity in thyroid cancer cells: Opportunities for repurposing. J. Clin. Endocrinol. Metab. 97, E319–E328 (2012).
Cayo, M. A. et al. A drug screen using human iPSC-derived hepatocyte-like cells reveals cardiac glycosides as a potential treatment for hypercholesterolemia. Cell Stem Cell 20, 478-489.e5 (2017).
Song, Y. et al. Inhibitors of Na+/K+ ATPase exhibit antitumor effects on multicellular tumor spheroids of hepatocellular carcinoma. Sci. Rep. 10, 1–16 (2020).
Simpson, C. D. et al. Inhibition of the sodium potassium adenosine triphosphatase pump sensitizes cancer cells to anoikis and prevents distant tumor formation. Cancer Res. 69, 2739–2747 (2009).
Blanco, G. & Mercer, R. W. Isozymes of the Na-K-ATPase: Heterogeneity in structure, diversity in function. Am. J. Physiol. Renal Physiol. 275, 633–650 (1998).
Toyoshima, C., Kanai, R. & Cornelius, F. First crystal structures of Na+, K+-ATPase: New light on the oldest ion pump. Structure 19, 1732–1738 (2011).
Clausen, M. V., Hilbers, F. & Poulsen, H. The structure and function of the Na, K-ATPase isoforms in health and disease. Front. Physiol. 8, 371 (2017).
Yap, J. Q., Seflova, J., Sweazey, R., Artigas, P. & Robia, S. L. FXYD proteins and sodium pump regulatory mechanisms. J. Gen. Physiol. https://doi.org/10.1085/jgp.202012633 (2021).
Monti, J. L. E., Montes, M. R. & Rossi, R. C. Steady-state analysis of enzymes with non-Michaelis-Menten kinetics: The transport mechanism of Na+/K+-ATPase. J. Biol. Chem. 293, 1373–1385 (2018).
Post, R. L., Kume, S., Tobin, T., Orcutt, B. & Sen, A. K. Flexibility of an active center in sodium-plus-potassium adenosine triphosphatase. J. Gen. Physiol. 54, 306–326 (1969).
Albers, R. W. Biochemical aspects of active transport. Annu. Rev. Biochem. 36, 727–756 (1967).
Kubala, M. ATP-binding to P-type ATPases as revealed by biochemical, spectroscopic, and crystallographic experiments. Proteins Struct. Funct. Bioinform. 64, 1–12 (2006).
Clarke, R. J., Catauro, M., Rasmussen, H. H. & Apell, H. J. Quantitative calculation of the role of the Na+,K+-ATPase in thermogenesis. Biochim. Biophys. Acta Bioenerg. 1827, 1205–1212 (2013).
Páez, O. et al. A Model for the homotypic interaction between Na+, K+-ATPase β1 subunits reveals the role of extracellular residues 221–229 in Its Ig-like domain. Int. J. Mol. Sci. 20(18), 4538 (2019).
Tokhtaeva, E. et al. Epithelial junctions depend on intercellular trans-interactions between the Na, K-ATPase β1 subunits. J. Biol. Chem. 286, 25801–25812 (2011).
Vagin, O., Tokhtaeva, E. & Sachs, G. The role of the β1 subunit of the Na, K-ATPase and its glycosylation in cell-cell adhesion. J. Biol. Chem. 281, 39573–39587 (2006).
Kitamura, N. et al. Mouse Na+/K+-ATPase β1-subunit has a K+-dependent cell adhesion activity for β-GlcNac-terminating glycans. Proc. Natl. Acad. Sci. U.S.A 102, 2796–2801 (2005).
Rajasekaran, S. A. et al. Na, K-ATPase β-subunit is required for epithelial polarization, suppression of invasion, and cell motility. Mol. Biol. Cell 12, 279–295 (2001).
Laursen, M., Yatime, L., Nissen, P. & Fedosova, N. U. Crystal structure of the high-affinity Na+, K+-ATPase–ouabain complex with Mg2+ bound in the cation binding site. Proc. Natl. Acad. Sci. U.S.A 110, 10958–10963 (2013).
Aperia, A., Akkuratov, E. E., Fontana, J. M. & Brismar, H. Na+-K+-ATPase, a new class of plasma membrane receptors. Am. J. Physiol.-Cell Physiol. 310(7), C491–C495 (2016).
Xie, Z. & Askari, A. Na+/K+-ATPase as a signal transducer. Eur. J. Biochem. 269, 2434–2439 (2002).
Hamlyn, J. M. et al. Identification and characterization of a ouabain-like compound from human plasma. Proc. Natl. Acad. Sci. U.S.A 88, 6259–6263 (1991).
Schoner, W. Endogenous cardiac glycosides, a new class of steroid hormones. Eur. J. Biochem. 269, 2440–2448 (2002).
Bagrov, A. Y., Shapiro, J. I. & Fedorova, O. V. Endogenous cardiotonic steroids: Physiology, pharmacology, and novel therapeutic targets. Pharmacol. Rev. 61, 9–38 (2009).
Lewis, L. K. et al. Ouabain is not detectable in human plasma. Hypertension 24, 549–555 (1994).
Baecher, S., Kroiss, M., Fassnacht, M. & Vogeser, M. No endogenous ouabain is detectable in human plasma by ultra-sensitive UPLC-MS/MS. Clin. Chim. Acta 431, 87–92 (2014).
Ballard, H. E., Sytsma, K. J. & Kowal, R. R. Shrinking the violets: Phylogenetic relationships of infrageneric groups in Viola (Violaceae) based on internal transcribed spacer DNA sequences. Syst. Bot. 23, 439–458 (1999).
Rahman, I. U. et al. Contributions to the phytotherapies of digestive disorders: Traditional knowledge and cultural drivers of Manoor Valley, Northern Pakistan. J. Ethnopharmacol. 192, 30–52 (2016).
Bhatt, V. P. & Negi, G. C. S. Ethnomedicinal plant resources of Jaunsari tribe of Garhwal Himalaya, Uttaranchal. Indian J. Tradit. Knowl. 5, 331–335 (2006).
Verma, G., Dua, V. K., Agarwal, D. D. & Atul, P. K. Anti-malarial activity of Holarrhena antidysenterica and Viola canescens, plants traditionally used against malaria in the Garhwal region of north-west Himalaya. Malar. J. 10, 1–5 (2011).
Chung, I. M., Kim, M. Y., Park, W. H. & Moon, H. I. Aldose reductase inhibitors from Viola hondoensis W. Becker et H Boss. Am. J. Chin. Med. 36, 799–803 (2008).
Zhu, Y., Zhao, L., Wang, X. & Li, P. Pharmacognostical and phytochemical studies of Viola tianschanica Maxim.—An Uyghur ethnomedicinal plant. J. Pharm. Pharmacogn. Res. 4, 95–106 (2016).
Wang, H., Cong, W. L., Fu, Z. L., Chen, D. F. & Wang, Q. Anti-complementary constituents of Viola kunawarensis. Nat. Prod. Res. 31, 2312–2315 (2017).
Feyzabadi, Z., Ghorbani, F., Vazani, Y. & Zarshenas, M. M. A critical review on phytochemistry, pharmacology of Viola odorata L. and related multipotential products in traditional Persian medicine. Phytother. Res. 31, 1669–1675 (2017).
Cu, J. Q., Perineau, F. & Gaset, A. Volatile components of violet leaves. Phytochemistry 31, 571–573 (1992).
Akhbari, M., Batooli, H. & Kashi, F. J. Composition of essential oil and biological activity of extracts of Viola odorata L. from central Iran. Nat. Prod. Res. 26, 802–809 (2012).
Parsley, N. C. et al. PepSAVI-MS reveals anticancer and antifungal cycloviolacins in Viola odorata. Phytochemistry 152, 61–70 (2018).
Narayani, M., Chadha, A. & Srivastava, S. Cyclotides from the indian medicinal plant Viola odorata (Banafsha): Identification and characterization. J. Nat. Prod. 80, 1972–1980 (2017).
Rizwan, K. et al. A comprehensive review on chemical and pharmacological potential of Viola betonicifolia: A plant with multiple benefits. Molecules 24, 3138 (2019).
Zubek, S., Rola, K., Szewczyk, A., Majewska, M. L. & Turnau, K. Enhanced concentrations of elements and secondary metabolites in Viola tricolor L. induced by arbuscular mycorrhizal fungi. Plant Soil 390, 129–142 (2015).
Welling, M. T., Liu, L., Rose, T. J., Waters, D. L. E. & Benkendorff, K. Arbuscular mycorrhizal fungi: effects on plant terpenoid accumulation. Plant Biol. 18, 552–562 (2016).
Harrison, M. J. Molecular and cellular aspects of the arbuscular mycorrhizal symbiosis. Annu. Rev. Plant Biol. 50(1), 361–389 (1999).
Katoch, M., Singh, A., Singh, G., Wazir, P. & Kumar, R. Phylogeny, antimicrobial, antioxidant and enzyme-producing potential of fungal endophytes found in Viola odorata. Ann. Microbiol. 67, 529–540 (2017).
Rauf, A. et al. Proanthocyanidins: A comprehensive review. Biomed. Pharmacother. 116, 108999 (2019).
Rárová, L. et al. Identification of narciclasine as an in vitro anti-inflammatory component of Cyrtanthus contractus by correlation-based metabolomics. J. Nat. Prod. 82, 1372–1376 (2019).
Enomoto, H., Takahashi, S., Takeda, S. & Hatta, H. Distribution of flavan-3-ol species in ripe strawberry fruit revealed by matrix-assisted laser desorption/ionization-mass spectrometry imaging. Molecules 25, 103 (2019).
Willer, E. A. et al. The vascular barrier-protecting hawthorn extract WS® 1442 raises endothelial calcium levels by inhibition of SERCA and activation of the IP3 pathway. J. Mol. Cell. Cardiol. 53, 567–577 (2012).
Schwinger, R. H. G., Pietsch, M., Frank, K. & Brixius, K. Crataegus special extract WS 1442 increases force of contraction in human myocardium cAMP-independently. J. Cardiovasc. Pharmacol. 35, 700–707 (2000).
Koch, E. & Malek, F. A. Standardized extracts from hawthorn leaves and flowers in the treatment of cardiovascular disorders—preclinical and clinical studies. Planta Med. 77, 1123–1128 (2011).
Yokomichi, T. et al. Ursolic acid inhibits Na+/K+-ATPase activity and prevents TNF-α-induced gene expression by blocking amino acid transport and cellular protein synthesis. Biomolecules 1, 32–47 (2011).
Chen, R. J. Y. et al. Steroid-like compounds in Chinese medicines promote blood circulation via inhibition of Na+/K+-ATPase. Acta Pharmacol. Sin. 31, 696–702 (2010).
Svedström, U. et al. High-performance liquid chromatographic determination of oligomeric procyanidins from dimers up to the hexamer in hawthorn. J. Chromatogr. A 968, 53–60 (2002).
Hellenbrand, N., Sendker, J., Lechtenberg, M., Petereit, F. & Hensel, A. Isolation and quantification of oligomeric and polymeric procyanidins in leaves and flowers of Hawthorn (Crataegus spp.). Fitoterapia 104, 14–22 (2015).
Svedström, U., Vuorela, H., Kostiainen, R., Laakso, I. & Hiltunen, R. Fractionation of polyphenols in hawthorn into polymeric procyanidins, phenolic acids and flavonoids prior to high-performance liquid chromatographic analysis. J. Chromatogr. A 1112, 103–111 (2006).
Cui, T. et al. Quantification of the polyphenols and triterpene acids in Chinese hawthorn fruit by high-performance liquid chromatography. J. Agric. Food Chem. 54, 4574–4581 (2006).
Cui, T., Nakamura, K., Tian, S., Kayahara, H. & Tian, Y.-L. Polyphenolic content and physiological activities of Chinese hawthorn extracts. Biosci. Biotechnol. Biochem. 70, 2948–2956 (2006).
Yang, B. & Liu, P. Composition and health effects of phenolic compounds in hawthorn (Crataegus spp.) of different origins. J. Sci. Food Agric. 92, 1578–1590 (2012).
Zumdick, S., Petereit, F., Luftmann, H. & Hensel, A. Preparative isolation of oligomeric procyanidins from Hawthorn (Crataegus spp.). Pharmazie 64, 286–288 (2009).
Souccar, C. et al. Inhibition of gastric acid secretion by a standardized aqueous extract of Cecropia glaziovii Sneth and underlying mechanism. Phytomedicine 15, 462–469 (2008).
Mizukasi, S., Tanabe, Y., Noguchi, M. & Tamaki, E. p-coumaroylputrescine, caffeoylputrescine and feruloylputrescine from callus tissue culture of Nicotiana tabacum. Phytochemistry 10, 1347–1350 (1971).
Li, Z. et al. Deep annotation of hydroxycinnamic acid amides in plants based on ultra-high-performance liquid chromatography-high-resolution mass spectrometry and its in silico database. Anal. Chem. 90, 14321–14330 (2018).
Ryabinin, A. A. & Il’Ina, E. M. The alkaloid of Salsola subaphylla. Dok. Akad. Nauk SSSR 67, 513 (1949).
Bardon, C. et al. Biological denitrification inhibition (BDI) with procyanidins induces modification of root traits, growth and N status in Fallopia x bohemica. Soil Biol. Biochem. 107, 41–49 (2017).
Kraus, T. E. C., Zasoski, R. J., Dahlgren, R. A., Horwath, W. R. & Preston, C. M. Carbon and nitrogen dynamics in a forest soil amended with purified tannins from different plant species. Soil Biol. Biochem. 36, 309–321 (2004).
Čechová, P., Berka, K. & Kubala, M. Ion Pathways in the Na+/K+-ATPase. J. Chem. Inf. Model. 56, 2434–2444 (2016).
Chebrolu, S., Ma, H. & Artigas, P. State-Dependent Movement between the First and Last External Loops of the Na/K Pump α Subunit. Biophys. J. 106, 582a (2014).
Young, V. C. & Artigas, P. Displacement of the Na+/K+-pump’s transmembrane domains demonstrate conserved conformational changes in P-type 2 ATPases. Proc. Natl. Acad. Sci. U.S.A 118, e2019317118 (2021).
Kubala, M. et al. Flavonolignans as a novel class of sodium pump inhibitors. Front. Physiol. 7, 115 (2016).
Gu, L. et al. Concentrations of proanthocyanidins in common foods and estimations of normal consumption. J. Nutr. 134, 613–617 (2004).
Kosińska, A. & Andlauer, W. Cocoa polyphenols are absorbed in Caco-2 cell model of intestinal epithelium. Food Chem. 135, 999–1005 (2012).
Mendoza-Wilson, A. M. et al. Absorption of dimers, trimers and tetramers of procyanidins present in apple skin by IEC-18 cell monolayers. J. Funct. Foods 27, 386–391 (2016).
Deprez, S., Mila, I., Huneau, J. F., Tome, D. & Scalbert, A. Transport of proanthocyanidin dimer, trimer, and polymer across monolayers of human intestinal epithelial Caco-2 cells. Antioxid. Redox Signal. 3, 957–967 (2001).
Serra, A. et al. Distribution of procyanidins and their metabolites in rat plasma and tissues in relation to ingestion of procyanidin-enriched or procyanidin-rich cocoa creams. Eur. J. Nutr. 52, 1029–1038 (2013).
Prasain, J. K. et al. Liquid chromatography tandem mass spectrometry identification of proanthocyanidins in rat plasma after oral administration of grape seed extract. Phytomedicine 16, 233–243 (2009).
Tsang, C. et al. The absorption, metabolism and excretion of flavan-3-ols and procyanidins following the ingestion of a grape seed extract by rats. Br. J. Nutr. 94, 170–181 (2005).
Serra, A. et al. Determination of procyanidins and their metabolites in plasma samples by improved liquid chromatography–tandem mass spectrometry. J. Chromatogr. B 877, 1169–1176 (2009).
Serra, A. et al. Bioavailability of procyanidin dimers and trimers and matrix food effects in in vitro and in vivo models. Br. J. Nutr. 103, 944–952 (2010).
Shoji, T. et al. Apple procyanidin oligomers absorption in rats after oral administration: Analysis of procyanidins in plasma using the porter method and high-performance liquid chromatography/tandem mass spectrometry. J. Agric. Food Chem. 54, 884–892 (2006).
Ottaviani, J. I., Kwik-Uribe, C., Keen, C. L. & Schroeter, H. Intake of dietary procyanidins does not contribute to the pool of circulating flavanols in humans. Am. J. Clin. Nutr. 95, 851–858 (2012).
Rios, L. Y. et al. Cocoa procyanidins are stable during gastric transit in humans. Am. J. Clin. Nutr. 76, 1106–1110 (2002).
Holt, R. R. et al. Procyanidin dimer B2 [epicatechin-(4β-8)-epicatechin] in human plasma after the consumption of a flavanol-rich cocoa. Am. J. Clin. Nutr. 76, 798–804 (2002).
Sano, A. et al. Procyanidin B1 is detected in human serum after intake of proanthocyanidin-rich grape seed extract. Biosci. Biotechnol. Biochem. 67, 1140–1143 (2003).
Fedosova, N. U. Purification of Na, K-ATPase from pig kidney in Methods in Molecular Biology 5–10 (Humana Press Inc., New York, 2016).
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. 31, 455–461 (2009).
Salentin, S., Schreiber, S., Haupt, V. J., Adasme, M. F. & Schroeder, M. PLIP: Fully automated protein-ligand interaction profiler. Nucleic Acids Res. 43, W443–W447 (2015).
Funding
This study was supported by ERDF project "Development of pre-applied research in nanotechnology and biotechnology" (No. CZ.02.1.01/0.0/0.0/17_048/0007323) and Palacky University (IGA_PrF_2022_007). Computational resources were supplied by the project "e-Infrastruktura CZ" (e-INFRA LM2018140) provided within the program Projects of Large Research, Development and Innovations Infrastructures.
Author information
Authors and Affiliations
Contributions
Conceptualization, J.G. and T.H.; methodology, T.H., J.G., M.Z. and M.K., writing-original draft preparation T.H., writing-review and editing, J.G., M.S. and M.K. All authors have read and agreed to the published version of the manuscript.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher's note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.
About this article
Cite this article
Heger, T., Zatloukal, M., Kubala, M. et al. Procyanidin C1 from Viola odorata L. inhibits Na+,K+-ATPase. Sci Rep 12, 7011 (2022). https://doi.org/10.1038/s41598-022-11086-y
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41598-022-11086-y
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.