Next Article in Journal
Determination of Volatile Organic Compounds and Endogenous Extracts and Study of Expression Patterns of TPS and BSMT in the Flowers of Seven Lilium Cultivars
Previous Article in Journal
Mechanism and Progress of Natural Products in the Treatment of NAFLD-Related Fibrosis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discovering Potential Compounds for Venous Disease Treatment through Virtual Screening and Network Pharmacology Approach

by
Oscar Salvador Barrera-Vázquez
,
Juan Luis Escobar-Ramírez
,
Jacinto Santiago-Mejía
,
Omar Francisco Carrasco-Ortega
and
Gil Alfonso Magos-Guerrero
*
Department of Pharmacology, Faculty of Medicine, University National Autonomous of Mexico (UNAM), Mexico City 04510, Mexico
*
Author to whom correspondence should be addressed.
Molecules 2023, 28(24), 7937; https://doi.org/10.3390/molecules28247937
Submission received: 31 October 2023 / Revised: 28 November 2023 / Accepted: 30 November 2023 / Published: 5 December 2023
(This article belongs to the Section Natural Products Chemistry)

Abstract

:
Peripheral venous hypertension has emerged as a prominent characteristic of venous disease (VD). This disease causes lower limb edema due to impaired blood transport in the veins. The phlebotonic drugs in use showed moderate evidence for reducing edema slightly in the lower legs and little or no difference in the quality of life. To enhance the probability of favorable experimental results, a virtual screening procedure was employed to identify molecules with potential therapeutic activity in VD. Compounds obtained from multiple databases, namely AC Discovery, NuBBE, BIOFACQUIM, and InflamNat, were compared with reference compounds. The examination of structural similarity, targets, and signaling pathways in venous diseases allows for the identification of compounds with potential usefulness in VD. The computational tools employed were rcdk and chemminer from R-Studio and Cytoscape. An extended fingerprint analysis allowed us to obtain 1846 from 41,655 compounds compiled. Only 229 compounds showed pharmacological targets in the PubChem server, of which 84 molecules interacted with the VD network. Because of their descriptors and multi-target capacity, only 18 molecules of 84 were identified as potential candidates for experimental evaluation. We opted to evaluate the berberine compound because of its affordability, and extensive literature support. The experiment showed the proposed activity in an acute venous hypertension model.

Graphical Abstract

1. Introduction

Venous disease (VD) is defined as a condition in which the veins cannot transport blood unidirectionally to the heart with a flow adapted to the needs of tissue drainage, temperature regulation, and hemodynamic reserve, regardless of their position and activity. The first manifestation of VD is increased venous pressure (venous hypertension or high blood pressure in the veins) with or without reflux [1]. In the United States, it affects over 25 million adults and over 6 million have advanced disease. The total cost of medical care is estimated at more than USD 3 billion per year [2]. In Latin America, it is estimated that the prevalence of VD is 68.11% [3], and in Mexico, the incidence rate in 2022 was 142.27/100,000 inhabitants [4].
The treatments of choice for VD are surgery, sclerotherapy, and mechanical compression. However, pharmacological treatment with drugs referred to as phlebotonics is frequently used due to their easy administration. In the Anatomical Therapeutic Chemical (ATC) system, phlebotonics are categorized as vasoprotective agents, specifically within the group of capillary stabilizing agents (ATC 2022) [5]. Most of these agents are derived from natural flavonoids extracted from plants [6].
Although phlebotonics are known as vasoactive drugs, their mechanism of action is not scientifically well established despite the lack of studies examining their pharmacological and clinical properties. These drugs have been found to impact macrocirculation, such as enhancing venous tone [7], as well as microcirculatory parameters such as capillary hyper-permeability [8]. Controversy exists regarding the clinical relevance of the efficacy and benefit–risk balance of phlebotonics. A meta-analysis study suggests a slight reduction in leg swelling, but it may not significantly affect quality of life or ulcer healing. Phlebotonics probably increase adverse events such as gastrointestinal disorders [6]. The therapeutic benefits of phlebotonics may not be fully realized in the VD because of their limitations and mechanisms of action.
Computer-aided drug design (CADD) is a discipline that uses a range of chemical-molecular and quantum methodologies. CADD allows us to investigate, develop, and synthesize medicinal chemical compounds [9]. One computational approach used in computer-aided drug design involves the application of chemoinformatics. This discipline specifically concentrates on extracting, processing, and extrapolating significant data from chemical structures [9]. Chemoinformatics is used to compare drugs, natural products (NPs), and semi-synthetic molecules (SMSs) to find potential candidates for pre-clinical research.
CADD has been promoted by various disciplines, such as chemoinformatics and network pharmacology, which were employed in the present investigation to find new drugs for VD treatment. CADD has proven to be effective in finding new candidates from NPs and SMSs that succeed better and have fewer side effects [10].
We tested a candidate found through virtual screening to address increased venous pressure, which is the first manifestation of VD. This evaluation was performed in a rat model of acute venous hypertension.

2. Results

2.1. Creation of a Reference Compounds Dataset

During the first phase of our research, MEDLINE, EMBASE, and Scopus servers enabled us to generate a dataset of reference compounds with phlebotonic, anti-inflammatory, and antioxidant properties. We found 145 reference compounds. Osiris DataWarrior software V5.2.1 and SwissADME server were used to determine their one-dimensional (1D) and two-dimensional (2D) molecular descriptors. Through the use of hierarchical clustering (K-means) and principal component analysis (PCA) via RStudio, an analysis was conducted on the molecular descriptors to determine the most significant reference compounds obtained from the literature (Figure 1). Supplementary Figure S1 contains the PCA results. The molecular descriptors include Fragments, Rotatable Bonds, Veber #violations, Molecular Flexibility, Fraction Csp3 Stereo Centers, sp3-Atoms, Saturated Rings, Non-Aromatic Rings, Molecular Complexity, MR, Total Surface Area, #Heavy atoms, Non-H Atoms, Molweight, Monoisotopic Mass, Total Molweight, Synthetic Accessibility, Hetero-Rings, Small Rings, Rings Closures, H-Donors, Polar Surface Area, TPSA, H-Acceptors, Non-C/H Atoms, Electronegative Atoms, and Muegge #violations. Cluster one showed the greatest similarity between the following 23 molecules: escin, crocin, echinomycin, cyclosporine A, amphotericin B, everolimus, rapamycin, chetomin, astragaloside IV, 20(R)-ginsenoside Rh2, epigallocatechin 3-gallate, troxerutin, O-(beta-hydroxyethyl)-rutoside, salvianolic acid B, keracyanin chloride, rutin, hidrosmin, diosmin, hesperidin, linarin, isorhoifolin, naringin, and betanin. Because of their great similarity in their molecular descriptors, these compounds are marked in green (Figure 1). The outstanding compounds were classified as reference compounds for the rest of the study and comprise various chemical families, such as saponins, carotenoids, cyclic peptides, macrolides, and alkaloids. Many flavonoids, including flavonoid glycosides, flavones, flavanols (catechins), flavonols, and anthocyanins, are present. Escin, diosmin, and hesperidin are used as first-line treatments for VD [11]. Table 1 presents a summary of the characteristics of the compounds in cluster one, including their activity and K-means coefficient. Escin, crocin, echinomycin, astragaloside IV, epigallocatechin 3-gallate, cyclosporin A, chetomin, and troxerutin were the most similar because of their high K-means coefficient.
Table 1. Main characteristics of the molecules of cluster one.
Table 1. Main characteristics of the molecules of cluster one.
Reference CompoundsIUPAC NamePubChem CIDPharmacological ActivityK-Means CoefficientReferences
Escin(2S,3S,4S,5R,6R)-6-[[(4S,6bS,8R,9R)-9-acetyloxy-8-hydroxy-4,8a-bis(hydroxymethyl)-4,6a,6b,11,11,14b-hexamethyl-10-[(E)-2-methylbut-2-enoyl]oxy-1,2,3,4a,5,6,7,8,9,10,12,12a,14,14a-tetradecahydropicen-3-yl]oxy]-4-hydroxy-3-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-(hydroxymethyl)oxan-2-yl]oxy-5-[(2R,3R,4S,5S,6R)-3,4,5-trihydroxy-6-(hydroxymethyl)oxan-2-yl]oxyoxane-2-carboxylic acid6476031Anti-edematous, anti-inflammatory, and venotonic agent2.03[12]
Crocinbis[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-[[(2R,3R,4S,5S,6R)-3,4,5-trihydroxy-6-(hydroxymethyl)oxan-2-yl]oxymethyl]oxan-2-yl] (2E,4E,6E,8E,10E,12E,14E)-2,6,11,15-tetramethylhexadeca-2,4,6,8,10,12,14-heptaenedioate5281233Inhibitor inhibits STAT3 activation induced by IL-61.50[13]
EchinomycinN-[2,4,12,15,17,25-hexamethyl-27-methylsulfanyl-3,6,10,13,16,19,23,26-octaoxo-11,24-di(propan-2-yl)-20-(quinoxaline-2-carbonylamino)-9,22-dioxa-28-thia-2,5,12,15,18,25-hexazabicyclo[12.12.3]nonacosan-7-yl]quinoxaline-2-carboxamide3197Anti-cancer agent and inhibitor of HIF1.39[14]
Cyclosporin A(3S,6S,9S,12R,15S,18S,21S,24S,30S,33S)-30-ethyl-33-[(E,1R,2R)-1-hydroxy-2-methylhex-4-enyl]-1,4,7,10,12,15,19,25,28-nonamethyl-6,9,18,24-tetrakis(2-methylpropyl)-3,21-di(propan-2-yl)-1,4,7,10,13,16,19,22,25,28,31-undecazacyclotritriacontane-2,5,8,11,14,17,20,23,26,29,32-undecone5284373Immunosuppressant agent and inhibitor of TCR signaling via NFAT-independent1.20[15]
Amphotericin B(1R,3S,5R,6R,9R,11R,15S,16R,17R,18S,19E,21E,23E,25E,27E,29E,31E,33R,35S,36R,37S)-33-[(2R,3S,4S,5S,6R)-4-amino-3,5-dihydroxy-6-methyloxan-2-yl]oxy-1,3,5,6,9,11,17,37-octahydroxy-15,16,18-trimethyl-13-oxo-14,39-dioxabicyclo[33.3.1]nonatriaconta-19,21,23,25,27,29,31-heptaene-36-carboxylic acid5280965Antibiotic agent for the treatment of life-threatening fungal infections and modulator of the immune system0.97[16]
Everolimus(1R,9S,12S,15R,16E,18R,19R,21R,23S,24E,26E,28E,30S,32S,35R)-1,18-dihydroxy-12-[(2R)-1-[(1S,3R,4R)-4-(2-hydroxyethoxy)-3-methoxycyclohexyl]propan-2-yl]-19,30-dimethoxy-15,17,21,23,29,35-hexamethyl-11,36-dioxa-4-azatricyclo[30.3.1.04,9]hexatriaconta-16,24,26,28-tetraene-2,3,10,14,20-pentone6442177Anticancer agent and mTOR inhibitor0.99[17]
Rapamycin(1R,9S,12S,15R,16E,18R,19R,21R,23S,24E,26E,28E,30S,32S,35R)-1,18-dihydroxy-12-[(2R)-1-[(1S,3R,4R)-4-hydroxy-3-methoxycyclohexyl]propan-2-yl]-19,30-dimethoxy-15,17,21,23,29,35-hexamethyl-11,36-dioxa-4-azatricyclo[30.3.1.04,9]hexatriaconta-16,24,26,28-tetraene-2,3,10,14,20-pentone5284616Inhibitor of mTOR complex 1 (mTORC1), which phosphorylates substrates including S6 kinase 1 (S6K1), eIF4E-binding protein 1 (4E-BP1), transcription factor EB (TFEB), unc-51-like autophagy-activating kinase 1 (Ulk1), and growth factor receptor-bound protein 10 (GRB-10)0.91[18]
Chetomin14-(hydroxymethyl)-3-[3-[[4-(hydroxymethyl)-5,7-dimethyl-6,8-dioxo-2,3-dithia-5,7-diazabicyclo[2.2.2]octan-1-yl]methyl]indol-1-yl]-18-methyl-15,16-dithia-10,12,18-triazapentacyclo[12.2.2.01,12.03,11.04,9]octadeca-4,6,8-triene-13,17-dione10417379Inhibitor of HIF-1α/p300 interaction, and antitumor agent1.16[19]
Astragaloside IV(2R,3R,4S,5S,6R)-2-[[(1S,3R,6S,8R,9S,11S,12S,14S,15R,16R)-14-hydroxy-15-[(2R,5S)-5-(2-hydroxypropan-2-yl)-2-methyloxolan-2-yl]-7,7,12,16-tetramethyl-6-[(2S,3R,4S,5R)-3,4,5-trihydroxyoxan-2-yl]oxy-9-pentacyclo[9.7.0.01,3.03,8.012,16]octadecanyl]oxy]-6-(hydroxymethyl)oxane-3,4,5-triol13943297Anti-inflammatory agent and inhibitor of NF-kappaB activation and adhesion molecule expression1.29[20]
20(R)-ginsenoside Rh2(2R,3R,4S,5S,6R)-2-[[(3S,5R,8R,9R,10R,12R,13R,14R,17S)-12-hydroxy-17-[(2R)-2-hydroxy-6-methylhept-5-en-2-yl]-4,4,8,10,14-pentamethyl-2,3,5,6,7,9,11,12,13,15,16,17-dodecahydro-1H-cyclopenta[a]phenanthren-3-yl]oxy]-6-(hydroxymethyl)oxane-3,4,5-triol54580480Antitumoral agent0.63[21]
Epigallocatechin 3-gallate(2S,3S,4S,5R,6S)-6-[2,6-dihydroxy-4-[[(2R,3R)-5-hydroxy-7-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-(hydroxymethyl)oxan-2-yl]oxy-2-(3,4,5-trihydroxyphenyl)-3,4-dihydro-2H-chromen-3-yl]oxycarbonyl]phenoxy]-3,4,5-trihydroxyoxane-2-carboxylic acid102025303Anti-carcinogen, anti-tumorigenesis, anti-proliferation, anti-angiogenesis, and antioxidant agent, cell death inductor1.27[22]
Troxerutin2-[3,4-bis(2-hydroxyethoxy)phenyl]-5-hydroxy-7-(2-hydroxyethoxy)-3-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-[[(2R,3R,4R,5R,6S)-3,4,5-trihydroxy-6-methyloxan-2-yl]oxymethyl]oxan-2-yl]oxychromen-4-one5486699Radioprotective and antioxidant agent1.12[23]
O-(beta-hydroxyethyl)-rutoside2-(3,4-dihydroxyphenyl)-5-hydroxy-7-(2-hydroxyethoxy)-3-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-[[(2R,3R,4R,5R,6S)-3,4,5-trihydroxy-6-methyloxan-2-yl]oxymethyl]oxan-2-yl]oxychromen-4-one9852585An agent used as a treatment for disorders of the venous and microcirculatory systems0.95[24]
Salvianolic acid B(2R)-2-[(E)-3-[(2S,3S)-3-[(1R)-1-carboxy-2-(3,4-dihydroxyphenyl)ethoxy]carbonyl-2-(3,4-dihydroxyphenyl)-7-hydroxy-2,3-dihydro-1-benzofuran-4-yl]prop-2-enoyl]oxy-3-(3,4-dihydroxyphenyl)propanoic acid6451084Anti-inflammatory and antioxidant agent0.98[25]
Keracyanin chloride(2R,3R,4R,5R,6S)-2-[[(2R,3S,4S,5R,6S)-6-[2-(3,4-dihydroxyphenyl)-5,7-dihydroxychromenylium-3-yl]oxy-3,4,5-trihydroxyoxan-2-yl]methoxy]-6-methyloxane-3,4,5-triol chloride29231Anti-inflammatory and antioxidant agent0.85[26]
Rutin2-(3,4-dihydroxyphenyl)-5,7-dihydroxy-3-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-[[(2R,3R,4R,5R,6S)-3,4,5-trihydroxy-6-methyloxan-2-yl]oxymethyl]oxan-2-yl]oxychromen-4-one5280805Anticancer, antidiabetic, antimicrobial, anticoagulant, antioxidant, cytoprotective, vasoprotective, anticarcinogenic, neuroprotective and cardioprotective agent0.85[27]
Hidrosmin5-(2-hydroxyethoxy)-2-(3-hydroxy-4-methoxyphenyl)-7-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-[[(2R,3R,4R,5R,6S)-3,4,5-trihydroxy-6-methyloxan-2-yl]oxymethyl]oxan-2-yl]oxychromen-4-one3087722Venoactive agent and post-thrombotic syndrome protector0.88[28]
Diosmin5-hydroxy-2-(3-hydroxy-4-methoxyphenyl)-7-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-[[(2R,3R,4R,5R,6S)-3,4,5-trihydroxy-6-methyloxan-2-yl]oxymethyl]oxan-2-yl]oxychromen-4-one5281613Agent for treatment of chronic venous insufficiency and varicose veins, with antioxidant, anticancer activities0.82[29,30]
Hesperidin(2S)-5-hydroxy-2-(3-hydroxy-4-methoxyphenyl)-7-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-[[(2R,3R,4R,5R,6S)-3,4,5-trihydroxy-6-methyloxan-2-yl]oxymethyl]oxan-2-yl]oxy-2,3-dihydrochromen-4-one53477767Antioxidant, neuroprotective, and anti-inflammatory agent0.82[31]
Linarin5-hydroxy-2-(4-methoxyphenyl)-7-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-[[(2R,3R,4R,5R,6S)-3,4,5-trihydroxy-6-methyloxan-2-yl]oxymethyl]oxan-2-yl]oxychromen-4-one5317025Antioxidant and anti-inflammatory agents0.76[32]
Isorhoifolin5-hydroxy-2-(4-hydroxyphenyl)-7-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-[[(2R,3R,4R,5R,6S)-3,4,5-trihydroxy-6-methyloxan-2-yl]oxymethyl]oxan-2-yl]oxychromen-4-one9851181Antioxidant agent0.75[33]
Naringin(2S)-7-[(2S,3R,4S,5S,6R)-4,5-dihydroxy-6-(hydroxymethyl)-3-[(2S,3R,4R,5R,6S)-3,4,5-trihydroxy-6-methyloxan-2-yl]oxyoxan-2-yl]oxy-5-hydroxy-2-(4-hydroxyphenyl)-2,3-dihydrochromen-4-one442428Antioxidant, antitumor, antiviral, antibacterial, anti-inflammatory, anti-adipogenic, and cardioprotective agent0.72[34]
Betanin1-[(2E)-2-(2,6-dicarboxy-2,3-dihydro-1H-pyridin-4-ylidene)ethylidene]-6-hydroxy-5-[(2S,3R,4S,5S,6R)-3,4,5-trihydroxy-6-(hydroxymethyl)oxan-2-yl]oxy-2,3-dihydroindol-1-ium-2-carboxylate12300103Antioxidant and anti-inflammatory agents0.46[35]

2.2. Selection of Drug-like NPs and SMSs Based on QED

We searched candidates compounds with potential usefulness in VD from NP and SMSs molecules collected from several databases. At the end of the search, 41,655 molecules derived from InflamNat (657 molecules), BIOFACQUIM (422 molecules), NuBBE (155 molecules), and AC Discovery (7320 NP molecules and 33,063 SMS molecules) were obtained. The molecular descriptors of 41,655 molecules and 145 reference molecules were analyzed to examine their distribution regarding the QED (Quantitative Estimate of Druglikeness) index. These results are shown in Figure 2 and Figure 3, where it is observed that the values of molecular descriptors for NPs, SMSs, and reference compounds are distributed within a similar range. This tendency in the distribution is attributed to the fact that most reference compounds are derived from natural sources.
At the end of QED analysis, molecules with a score greater than 0.5 were considered drug-like by their potential use of being administered orally. From 41,655 molecules, 27,868 molecules passed the rule (QED > 0.5), of which 2964 were NP molecules (10.1%) and 24,904 were SMS molecules (96.51%). The total number of molecules that violated the rule (QED 0.5) was 12,670 molecules (4511 were NP molecules and 8159 were SMS molecules).
To analyze and remove compounds that have non-desirable effects, we used the Osiris Data Warrior software V5.2.1. This analysis identified 26,530 molecules (1626 NP molecules and 24,904 SMS molecules) free from potential tumorigenic and deleterious functions, including in the reproductive system. The chemical space of our molecules obtained using three-dimensional principal component analysis (3D-PCA) allowed us to identify drug-like NPs and SMSs that were similar to our reference compounds. At the end of the analysis, 1953 similar molecules (1293 NPs and 660 SMSs) were identified by their chemical space (Figure 4).

2.3. Selection of NPs and SMSs Structurally Similar to Reference Compounds

To further analyze the obtained result via 3D-PCA, we performed a fingerprint analysis. Of the 1953 drug-like molecules obtained via 3D-PCA, 1846 were found in PubChem because this server contains various types of structure searches in 2D, which are important for fingerprint analysis. To achieve this aim, a comparison was made between the chemical structure of 1846 drug-like molecules and the 23 reference compounds obtained from cluster 1. Figure 5 presents an illustrative example of fingerprint analysis using the reference compound diosmin. The elbow method shows that the optimal number of clusters for clusterization is three (Figure 5A). Figure 5B shows the comparison of fingerprints between the 1846 drug-like molecules and the reference compounds through the Tanimoto coefficient and Ward’s method clustering. Figure 5C,D illustrate the Dunn index and silhouette coefficient (Sc), which confirm the clusters obtained using Ward’s method for the 1846 drug-like molecules. The number of drug-like molecules similar to each one from the 23 reference compounds is shown in Table 2. The highest Sc value is found in the cluster that contains the phlebotonic reference compound. This result suggests a strong similarity between 1846 drug-like molecules and 23 reference compounds. Finally, the mean of the Sc of the clusterization performed with the 23 reference compounds was 0.35. Supplementary Figure S2 shows the complete analysis for each of the 23 reference compounds. Supplementary Materials show the clusters of those compounds that were structurally similar to the reference compounds. This information is shown in these sections since it is difficult to present it clearly in Figure 5 because of the large amount of data.

2.4. Network Analysis to Obtain the Best Molecules with Potential Usefulness in VD through Their Multi-Target Capacity

An important requirement for the construction of structural networks is that analyzed molecules interact with pharmacological targets. From the 1846 molecules obtained from fingerprint analysis, only 229 showed pharmacological targets on the PubChem server. To enhance compound interaction with the VD networks, a macro instruction was used to match all the involved genes. At the end of the macro, only 84 out of 229 molecules interacted directly with the genes involved in the VD networks (Table 3). A compound–target network studied the interactions between the genes and molecules with potential usefulness in VD. Our compound–target network (Figure 6) for VD comprised 1430 nodes and 2232 edges. At the end of the network analysis, 23 molecules were identified with the best multi-target capacity. However, only 18 molecules have the largest number of targets involved in more than one pathway of the VD network. So, we identified these 18 molecules with a multi-target and multi-pathway capacity (Figure 6).
The 18 potential phlebotonic molecules are represented in Figure 6 as yellow nodes and their pharmacological targets as orange nodes. These potential candidates are (E)-3-phenyl-N-(2-phenylethyl)prop-2-enamide; 3-hydroxy-1-(4-hydroxyphenyl)propan-1-one; 2-(4-methoxyphenyl)-2,3-dihydrochromen-4-one; 5-hydroxy-7-methoxy-3-(4-methoxyphenyl)chromen-4-one; 7-hydroxy-2-phenylchromen-4-one; 5,7-dihydroxy-2-(4-methoxyphenyl)chromen-4-one; 2-hydroxy-6-(4-hydroxy-2-methoxy-6-methoxycarbonylphenoxy)-4-methylbenzoic acid; 16,17-dimethoxy-5,7-dioxa-13-azoniapentacyclo[11.8.0.02,10.04,8.015,20]henicosa-1(13),2,4(8),9,14,16,18,20-octaene chloride; 5,7-dihydroxy-3-(4-methoxyphenyl)chromen-4-one; (E)-N-[(4-hydroxy-3-methoxyphenyl)methyl]-8-methylnon-6-enamide; (E)-1-(2,4-dihydroxy-6-methoxyphenyl)-3-phenylprop-2-en-1-one; 5,7-dihydroxy-2-phenylchromen-4-one; 7-hydroxy-3-(4-methoxyphenyl)chromen-4-one; 9-methoxy-2,2-dimethyl-6H-pyrano[3,2-c]quinolin-5-one; (2S)-5,7-dihydroxy-2-(4-hydroxyphenyl)-2,3-dihydrochromen-4-one; 1-[(3,4-dimethoxyphenyl)methyl]-6,7-dimethoxyisoquinoline;hydrochloride; 3-(4-hydroxyphenyl)-1-(2,4,6-trihydroxyphenyl)propan-1-one; and (2Z)-2-[(3,4-dihydroxyphenyl)methylidene]-6-hydroxy-1-benzofuran-3-one. Table 4 summarizes the characteristics of compounds with potential usefulness in VD derived from the compound–target network.
Table 4. Main characteristics of the best candidate compounds with usefulness in VD obtained via the multi-target capacity.
Table 4. Main characteristics of the best candidate compounds with usefulness in VD obtained via the multi-target capacity.
CompoundsIUPAC NameStructureSourcePharmacological ActivityNumber of TargetsReferences
N-Phenethylcinnamamide
CID 795855
(E)-3-phenyl-N-(2-phenylethyl)prop-2-enamide Molecules 28 07937 i001Acmella radicansInhibitor agent of prostaglandin and leukotriene biosyntheses4[36,37]
3,4′-Dihydroxypropiophenone
CID 638759
3-hydroxy-1-(4-hydroxyphenyl)propan-1-one Molecules 28 07937 i002The roots of the Carissa edulisAChE inhibitor, cytotoxic and antioxidant agent.7[38]
4′-Methoxyflavanone
CID 102928
2-(4-methoxyphenyl)-2,3-dihydrochromen-4-one Molecules 28 07937 i003Isolated from
Isaria fumosorosea
Inhibitor agent of glycation and aldose reductase activity4[39]
7,4′-dimethoxy-5-hydroxy isoflavone
CID 5386259
5-hydroxy-7-methoxy-3-(4-methoxyphenyl)chromen-4-one Molecules 28 07937 i004Isolated from the roots of Lotus polyphyllosAnticancer agent6[40,41]
7-hydroxyflavone
CID 5281894
7-hydroxy-2-phenylchromen-4-one Molecules 28 07937 i005Isolated from M. indicaVasorelaxant and anti-inflammatory agent5[42,43,44]
Acacetin
CID 5280442
5,7-dihydroxy-2-(4-methoxyphenyl)chromen-4-one Molecules 28 07937 i006Robinia pseudoacaciaAntiviral agent19[45]
Asterric acid
CID 3080568
2-hydroxy-6-(4-hydroxy-2-methoxy-6-methoxycarbonylphenoxy)-4-methylbenzoic acid Molecules 28 07937 i007Aspergillus terreusAntibiotic agent4[46]
Berberine chloride
CID 12456
16,17-dimethoxy-5,7-dioxa-13-azoniapentacyclo[11.8.0.02,10.04,8.015,20]henicosa-1(13),2,4(8),9,14,16,18,20-octaene chlorideMolecules 28 07937 i008Rhizoma coptidisAntibacterial and anti-inflammatory agents4[47]
Biochanin A
CID 5280373
5,7-dihydroxy-3-(4-methoxyphenyl)chromen-4-one Molecules 28 07937 i009Red clover, soy, alfalfa sprouts, peanuts, chickpea (Cicer arietinum), and other legumesAntiviral agent9[48]
Capsaicin
CID 1548943
(E)-N-[(4-hydroxy-3-methoxyphenyl)methyl]-8-methylnon-6-enamideMolecules 28 07937 i010Chili peppersThe agent involved in the treatment of obesity, diabetes, cardiovascular conditions, cancer, airway diseases, itch, gastric, and urological disorders4[49]
Cardamonin
CID 641785
(E)-1-(2,4-dihydroxy-6-methoxyphenyl)-3-phenylprop-2-en-1-one Molecules 28 07937 i011Alpinia blepharocalyx, Alpinia gagnepainii, Alpinia conchigera, Alpinia hainanensis, Alpinia malaccensis, Alpinia mutica, Alpinia pricei, Alpinia rafflesiana, Alpinia speciosa, Amomum subulatum, Artemisia absinthium, Boesenbergia pandurata, Boesenbergia rotunda, Carya cathayensis, Cedrelopsis grevei, Combretum apiculatum, Comptonia peregrina, Desmos cochinchinensis, Elettaria cardamomum, Helichrysum forskahlii, Kaempferia parviflora, Morella pensylvanica, Piper dilatatum, Piper hispidum, Polygonum ferrugineum, Polygonum lapathifolium, Polygonum persicaria, Populus fremontii, Populus × euramericana, a hybrid between Populus deltoides and Populus nigra, Syzygium samarangense, Vitex leptobotrys, and Woodsia scopulina.Anti-inflammatory, antineoplastic, and antioxidant agent4[50]
Chrysin
CID 5281607
5,7-dihydroxy-2-phenylchromen-4-one Molecules 28 07937 i012Isolated from Pyrus pashia fruitVasorelaxant agent18[42,51]
Formononetin
CID 5280378
7-hydroxy-3-(4-methoxyphenyl)chromen-4-one Molecules 28 07937 i013Red cloverAnti-cancer agent4[52]
Haplamine
CID 648601
9-methoxy-2,2-dimethyl-6H-pyrano[3,2-c]quinolin-5-one Molecules 28 07937 i014Isolated from Haplophyllum perforatuAnti-cancer agent9[53]
Naringenin
CID 439246
(2S)-5,7-dihydroxy-2-(4-hydroxyphenyl)-2,3-dihydrochromen-4-one Molecules 28 07937 i015Citrus fruits, bergamot, tomatoes, and other fruitsAnti-Hepatitis C virus, antiaging, anti-Alzheimer, antiasthma, anticancer, anti-Chikungunya virus, anticonvulsant, anti-dengue virus, antidiabetic, anti-Edwardsiellosis, anti-hyperlipidemic, anti-inflammatory, antimicrobial, antioxidant, antiplatelet, anti-stroke damage, cardioprotective, chronic kidney disease, expectorant, eye-protective, fertility, immunomodulatory, laxative, hepatoprotective, pregnancy, radioprotective, and weight-loss agent4[34]
Papaverine hydrochloride
CID 6084
1-[(3,4-dimethoxyphenyl)methyl]-6,7-dimethoxyisoquinoline;hydrochloride Molecules 28 07937 i016Papaver somniferum (Opium poppy)Vasodilator agent4[54]
Phloretin
CID 4788
3-(4-hydroxyphenyl)-1-(2,4,6-trihydroxyphenyl)propan-1-one Molecules 28 07937 i017AppleAntioxidative, anti-inflammatory, anti-microbial, anti-allergic, anticarcinogenic, anti-thrombotic, and hepatoprotective agent4[55]
Sulfuretin
CID 5281295
(2Z)-2-[(3,4-dihydroxyphenyl)methylidene]-6-hydroxy-1-benzofuran-3-one Molecules 28 07937 i018The stem bark of Albizzia julibrissin and heartwood of Rhus vernicifluaAdipogenesis inhibitor7[56]

2.5. Selection of Berberine for Experimental Validation

Once we found the most promising candidates with potential usefulness in the VD (18 compounds), we selected the following seven compounds as the most suitable to be tested experimentally: 5-hydroxy-7-methoxy-3-(4-methoxyphenyl)chromen-4-one; 5,7-dihydroxy-2-(4-methoxyphenyl)chromen-4-one; 16,17-dimethoxy-5,7-dioxa-13-azoniapentacyclo[11.8.0.02,10.04,8.015,20]henicosa-1(13),2,4(8),9,14,16,18,20-octaene chloride; 5,7-dihydroxy-3-(4-methoxyphenyl)chromen-4-one; 5,7-dihydroxy-2-phenylchromen-4-one; 7-hydroxy-3-(4-methoxyphenyl)chromen-4-one; and (2Z)-2-[(3,4-dihydroxyphenyl)methylidene]-6-hydroxy-1-benzofuran-3-one. Because of the availability, accessibility, administration, and literature background, 16,17-dimethoxy-5,7-dioxa-13-azoniapentacyclo[11.8.0.02,10.04,8.015,20]henicosa-1(13),2,4(8),9,14,16,18,20-octaene chloride, commonly known as berberine, was selected to be evaluated in a model of acute venous hypertension [57,58,59].

2.6. Effects of Ligation of the Veins on the Venous Pressure of the Left Hind Limb

Peripheral acute venous hypertension was induced in anesthetized rats by ligation of the right hind limb (RHL). After 15 min of stabilization, the mean arterial pressure (MAP), heart rate (HR), and venous pressure (VP) readings were taken. Figure 7 illustrates the comparison of VP values of the right hind limb (RHL) and left hind limb (LHL) before and after the ligation of the LHL veins. There were no significant differences in the baseline MAP and HR values before and after ligation (Table 5). In contrast, the ligation of four veins (superficial epigastric, proximal caudal, medial proximal geniculate, and popliteal) in the LHL increased the VP significantly both for rats that will be administered with berberine from a baseline of 11.42 mm Hg to 28.28 mm Hg (Figure 7A) and β-escin baseline of 11.4 mm Hg to 31.48 mm Hg (Figure 7B).

2.7. Effects of Berberine, and β-Escin on MAP, HR, and VP of the Anesthetized Rat

Successive injections of berberine (0.001, 0.003, 0.01, 0.03, 0.1, 0.3, 1 mg/kg, i.v.) and β-escin (0.56, 1, 17, 31, 56, and 10 mg/kg) were administered independently in anesthetized rats with venous ligation in the LHL. The last doses of berberine reduce the MAP with an increase, probably with a reflecting increase in HR (Figure 8A). In contrast, the administration of β-escin significantly reduces the HR and MAP with a higher dose than that employed in the experiments performed with berberine (Figure 8B). The administration of increased berberine doses resulted in a significant reduction in acute venous hypertension induced by the ligation of different veins in the LHL (Figure 9A). A similar effect but of lesser intensity was produced by the administration of increasing doses of β-escin (Figure 9B). The VP of the RHL used as a control does not show significant differences with the administration of berberine or β-escin. Scheme 1 summarizes the steps followed in this work.

3. Discussion

Drug discovery from natural products has a long and successful history. Natural product research through computer-aided drug design techniques has been broadly successful, particularly when we explore natural products’ chemical space to identify bioactive compounds, with an emphasis on drug discovery [60]. In this study, we employed the application of chemoinformatics and network pharmacology techniques to identify compounds that potentially exhibit usefulness in VD. The selection of natural products as the basis of our investigation was motivated by the prevalence of natural-origin compounds called phlebotonics in VD treatment, which have strong associations with important signaling pathways [61,62,63,64,65,66,67,68].
In this study, we found 145 compounds reported with usefulness in the VD in three databases: MEDLINE, EMBASE, and Scopus. After clustering analysis, we discovered three groups with similar molecular descriptors. The main cluster of similar drugs involved in VD includes escin, crocin, echinomycin, cyclosporin A, amphotericin B, everolimus, rapamycin, chetomin, astragaloside IV, 20(R)-ginsenoside Rh2, epigallocatechin 3-gallate, troxerutin, O-(beta-hydroxyethyl)-rutoside, salvianolic acid B, keracyanin chloride, rutin, hidrosmine, diosmine, hesperidin, linarin, isorhoifolin, naringin, and betanin. The cluster that contained the most representative reference compounds was corroborated by the PCA analysis.
Our dataset shows a strong clustering tendency in flavonoid-related molecular descriptors. Flavonoids provide many health benefits, such as protecting against cardiovascular diseases, allergies, diabetes, inflammation, and oxidative damage [11].
Interestingly, all of them have different action mechanisms with potential usefulness in VD. For instance, escin is the active component of Aesculus hippocastanum, the horse chestnut, which has been used in traditional medicine for centuries and is currently used to treat hemorrhoids, varicose veins, hematoma, and venous congestion; this compound has shown anti-edematous, anti-inflammatory and venotonic activity. Its chemical structure contains a trisaccharide linked to the 3-OH residue, such as glucose, xylose, and galactose. Moreover, esterified domains in the C21 and C22 with an organic acid such as angelic, tiglinic, or acetic acid are also present [12]. Crocin, a water-soluble carotenoid, has been reported to demonstrate strong antioxidant activity against reactive oxygen species. Its mechanism is defined by the NF-κB pathway and NF-κBp65 translocation. This translocation inhibits the secretion of pro-inflammatory cytokines. Echinomycin is a cyclic peptide with antibacterial activity from the quinoxaline family produced by a strain of Streptomyces echinatus. This compound inhibits a factor involved in the signaling pathways in varicose veins, hypoxia-inducible factor-1 (HIF-1) [69]. Cyclosporine A is a cyclic non-ribosomal peptide of eleven amino acids. This peptide has been reported to have both immunosuppressive and mitochondrial inhibition activity. These characteristics have been given to Cyclosporine A to be used as a treatment for inflammatory pathologies, for instance, coronary artery disease (CAD), and acute myocardial infarction (AMI). Also, the pro-inflammatory genes IL-1β, IL-6, and tumor necrosis factor (TNF)-α have been targeted by this cyclic non-ribosomal peptide [70], which is also involved in VD [61]. Amphotericin B, a polyene antifungal with reported antibiotic activity, down-regulates the upregulation of inflammatory-related genes such as IL-6 and IL-8, reported in TNF-α-stimulated gingival epithelial cells with the involvement of p38 MAP kinase and ERK [71]. Both interleukins are targeted in the process of VD. Rapamycin and everolimus, which are mammalian targets of rapamycin (mTOR) inhibitors, have been reported to have significant anti-inflammatory activity. Their mechanism is comprised of avoiding the release of IL-8 and vascular endothelial growth factor (VEGF), allowing the preservation of the anti-inflammatory cytokine interleukin-1 receptor antagonist (IL-1RA). Pre-incubation with everolimus has shown a significative potent effect in TNF-α-treated neutrophils, which play a role in VD. This pre-incubation reduces the release of VEGF, IL-8, and IL-1RA [72]. Chetomin is an HIF-1α/p300 interaction inhibitor, a factor involved in VD [19]. The astragaloside IV (AS-IV), a 3-O-beta-D-xylopyranosyl-6-O-beta-D-glucopyranosylcycloastragenol extracted from the Chinese medical herb Astragalus membranaceous (Fisch), which has been reported to have in vivo anti-inflammatory activity, has been shown to inhibit the TNFalpha-induced specific mRNA levels for E-selectin and VCAM-1 [20]. 20(R)-ginsenoside Rh2, a minor stereoisomer of ginsenoside Rh2, has been reported to be an outstanding phlebotonic reference drug because of its several activities involved in VD such as matrix metalloproteinase inhibitory, anti-inflammatory, and anti-oxidative activity [73]. Epigallocatechin-3-gallate (EGCG) is a type of catechin found in green tea with reported anti-inflammatory and antioxidant properties [22], while troxerutin, is a flavonoid with radioprotective, antioxidant, and other diverse pharmacological activities [23]. O-(beta-hydroxyethyl)-rutosides (HR) is employed to treat chronic venous disease, as well as the indications of chronic venous insufficiency (VD), varicose veins, and deep venous disease [74]. Salvianolic acid B has been proposed to be an excellent drug candidate for the treatment and prevention of cardiovascular diseases [75]. Keracyanin chloride has been reported to reduce the expression of endothelial inflammatory antigens [26]. Rutin, hidrosmin, diosmin, hesperidin, and linarin have been used as pharmacological treatments for VD [6]. Linarin is a component of a micronized purified flavonoid fraction (MPFF) called Daflon. MPFF mainly contains diosmin (90%) and other active flavonoids like hesperidin, diosmetin, linarin, and isorhoifolin. It is used to treat inflammation and chronic venous insufficiency [76]. Naringenin is a flavonoid belonging to the flavanones subclass. It is widely distributed in several citrus fruits, bergamot, tomatoes, and other fruits. Several biological activities have been ascribed to this phytochemical, among them antioxidant and anti-inflammatory activity [34]. Betanin is widely found in red beet and is the most common beta cyanin pigment that acts as a stimulator of antioxidant defense mechanisms, has considerable free-radical scavenger activity, and is an anti-inflammatory agent [35].
Lipinski’s rule of five is commonly employed to assess the drug-likeness of a compound. Nevertheless, the rule does not extend to substrates of biological transporters or natural products [77]. According to research [77], 16% of approved oral drugs do not adhere to at least one of Lipinski’s criteria, while 6% fail to meet two or more. Because of this rationale, we opt to substitute Lipinski’s rule of five with the Quantitative Estimation of Druglikeness (QED) as the benchmark for assessing the drug-likeness of NP and SMSs [77]. This estimation enabled us to enhance our acquisition of NPs and SMSs that could potentially serve as drug candidates. Within this context, the QED displayed greater leniency in the selection criteria. QED allowed us to avoid bias and enrich our selection of NPs and SMSs, obtaining a more realistic and graded result. In reality, many approvals violate the criteria of this rule [77]. The molecular descriptors’ distribution showed resemblance in both datasets, even with the presence of molecules of different characteristics, encompassing both natural and synthetic products, in the phlebotonic dataset. Consequently, these data suggest that natural products could serve as a dependable resource for drug discovery.
Although there are many clustering methods, only some are used in practice; two were used in this study: hierarchical grouping and the K-means. They have been used successfully in fingerprint analysis [78]. Within the scope of this study, the utilization of fingerprint analysis facilitated the acquisition of 1846 molecules because they resemble reference molecules found in cluster 1. This finding suggests that the physicochemical characteristics of these 1846 compounds exhibit a higher level of structural similarity to reference compounds. Various clustering methods were utilized to validate these analyzes. The establishment of a compound–target network facilitated the identification of compounds with potential usefulness in VD based on their therapeutic targets. The collected data aided us in confirming that the similarity in chemical structure can lead to similar biological functions. This is clear as certain NPs exhibited targets associated with various pathways related to VD, including the inflammatory response pathway, cellular response to hypoxia, PI3K-AKT-NFKB pathway, ion channel transport, Il-18 signaling pathway, VEGF-VEGFR-2, senescence-associated secretory phenotype (SASP), oxidative stress, and intrinsic pathways for apoptosis. However, few SMSs from the AC discovery database had reported targets in PubChem and none were related to VD pathways. Hence, it is crucial to conduct additional bioassays to identify more drug targets and document them in these databases. It is crucial to consider compounds that may not have direct associations with VD targets to prevent them from being disregarded, as some of these compounds might display phlebotonic activity because of structural similarities with compounds identified in cluster one of phlebotonics.
Through this study, compounds with a multi-target capacity were identified, which had been previously selected based on fingerprints. Thus, the development of multi-target drugs is essential as a promising strategy to address complex, multifactorial disorders [79], including VD [61]. Through the analysis of the structure and functionality of multi-target molecules, this study has identified three compounds that exhibit an enhanced likelihood of selectively acting on VD. These new molecules open a pre-clinic study field of great interest considering the knowledge regarding the NPs that must be confirmed in preclinical studies.
Traditional Chinese medicine has been widely used for the treatment of various diseases, obtaining great relevance within clinical applications. VDs that affect the vasculature of the heart, cerebrovascular disease, atherosclerosis, and diabetic complications have harmed the quality of life of patients and present an increase in the burden of health care services. Berberine, an isoquinoline alkaloid from Rhizoma coptidis, has wide application in traditional Chinese medicine because of its antibacterial and anti-inflammatory properties. With this background, a greater number of studies have been generated that have allowed us to identify several cellular and molecular targets for berberine, indicating its potential as an alternative therapeutic strategy for vascular diseases, in addition to providing novel evidence supporting the potential therapeutic use of berberine to combat VD. Currently, berberine displays remarkable anti-inflammatory, antioxidant, antiapoptotic, and anti-autophagic activity through the regulation of multiple signaling pathways, including AMP-activated protein kinase (AMPK), nuclear factor κB (NF-κB), mitogen-regulatory Silent Information System Activated Protein Kinase 1 (SIRT-1), Hypoxia-Inducible Factor 1α (HIF-1α), Phosphoinositide Vascular Endothelial Growth Factor 3-Kinase (PI3K), Protein Kinase B (Akt), Janus Kinase 2 (JAK-2), Ca2+ channels and endoplasmic reticulum stress, in addition to modulating Na+, Ca2+ concentration and lipid metabolism in vascular smooth muscle cells [47]. All these actions increase the possibility of obtaining therapeutic activity in venous activity because of its similarity with the structure of tested phlebotonic drugs and its multi-target ability as a protector in various vascular according to recent in vitro and in vivo experimental reports.
In our study, berberine was utilized due to its potential venous antihypertensive effect and its ability to protect vascular endothelial cells, enhance vascular remodeling, and inhibit inflammation, oxidative stress, autophagy, and apoptosis. The results of our experimental model of peripheral venous hypertension confirm a significant venous antihypertensive effect. However, the evaluation conducted so far focused on acute venous hypertension, necessitating an assessment of a chronic hypertension model to determine the efficacy of berberine in mitigating microcirculatory changes, safeguarding endothelial cells, reducing inflammation, and preventing ischemia by enhancing endothelial function.
It is widely recognized that phlebotonic activity extends beyond the mere elevation of venous wall tone. In the modern context, veno-active drugs should possess the ability to stimulate lymphatic drainage and enhance microcirculation [80], employing various mechanisms [81,82,83,84]. Endothelial dysfunction holds significance within the venous system as it may lead to a decrease in venodilation, resulting in an increase in venous tone and a decrease in venous compliance. This, in turn, enhances the mean circulatory filling pressure, leading to the maintenance or alteration of cardiac workload and contributing to the development of cardiovascular diseases. Changes in vein function may manifest early on, even before the onset of these diseases. However, is not yet fully understood if the venous endothelium dysfunction is involved in these alterations, so further studies are required [85]. This hypothesis may support our idea that berberine could be a more efficient phlebotonic than β-escin, because berberine generates vasorelaxation, avoiding endothelial damage and thus preventing the progression of VD. A similar process has been reported during the treatment with the phlebotonic rutin [86]. Our discovery supports the notion that berberine stimulates vasorelaxation mediated by the endothelium and enhances vasodilation in VSMCs by partially reducing oxidative stress [87], so we expect that these positive effects may present in the venous environment.
The effects of vasorelaxation may be involved because the targets of berberine found in our compound–target network were the KCNH2 and NFE2L2 genes as the most relevant. KCNH2 is a voltage-activated potassium channel found in cardiac muscle, nerve cells, and microglia. This gene encodes the pore-forming subunit of a rapidly activating-delayed rectifier potassium channel that plays an essential role in the final repolarization [88]. NFE2L2 is a transcription factor that shapes the antioxidant response in driving cancer progression, metastasis, and resistance to therapy [89].
The use of virtual screening and network pharmacology allowed us, for the first time, to find that berberine produces venous antihypertension activity. Our finding, coupled with the background of berberine activities, makes it a better candidate with usefulness in VD because it is a natural product with therapeutic potential in vascular diseases [47]. In addition, it has been shown to improve other conditions that are associated with VD, such as diabetes mellitus [47].
These experimental results belong only to an evaluation of one of many other candidates found in this work, which will be tested in a model of chronic venous hypertension that is currently under development.

4. Materials and Methods

To achieve the proposed aims, several were created:

4.1. In Silico Phase

4.1.1. Creation of a Reference Compound Dataset

Collection of Useful Compounds in VD

We created a dataset by reviewing literature from MEDLINE, EMBASE, and Scopus up to March 2023, focusing on useful compounds in VD. The search descriptors in this work were “venous disease”, “phlebotonics”, “venotonics”, “venous insufficiency treatment”, “venous disease treatment”, “chronic venous disease”, and “pathways”.

Selection of Reference Compounds from the Dataset Based on Hierarchical Analysis and Chemical Space Analysis

The molecular structure of each compound found in the literature (reference compounds dataset) was translated into a Simplified Molecule Input Line Entry System (SMILES) obtained from PubChem (https://pubchem.ncbi.nlm.nih.gov/) [90] (accessed on 3 April 2023). Subsequently, the SwissADME server (http://www.swissadme.ch) [91] (accessed on 10 April 2023) and Osiris Data Warrior software V5.2.1 were used to analyze the physicochemical properties and molecular descriptors. Total Molweight, Molweight, Monoisotopic Mass, cLogP, cLogS, H-Acceptors, H-Donors, Total Surface Area, Relative PSA, Polar Surface Area, Druglikeness, Shape Index, Molecular Flexibility, Molecular Complexity, Fragments, Non-H Atoms, Non-C/H Atoms, Electronegative Atoms, Stereo Centers, Rotatable Bonds, Rings Closures, Aromatic Atoms, sp3-Atoms, Symmetric atoms, Small Rings, Carbo-Rings, Hetero-Rings, Saturated Rings, Non-Aromatic Rings, Aromatic Rings, #Heavy atoms, #Aromatic heavy atoms, Fraction Csp3, MR, TPSA, log Kp (cm/s), Lipinski #violations, Ghose #violations, Veber #violations, Egan #violations, Muegge #violations, Bioavailability Score, and Synthetic Accessibility were the molecular descriptors for this work [92].
By using the K-means algorithm and distance matrix, similarities in the molecular descriptors of the reference compounds were identified. This analysis was performed with the complexheatmap package [93] and R-studio version 3.4. PCA was used to analyze the chemical space and data distribution to validate these results. PCA is defined as an orthogonal linear transformation technique that can transform the data into a new coordinate system. We utilized a two-dimensional system in our analysis [94]. Briefly, the variance of the data maximized on the first coordinate was called the first principal component. The second variance was maximized on the second coordinate; we used the factoextra package with R-Studio version 3.4.

4.1.2. Drug-Likeness of NPs and SMSs from BIOFACQUIM, Inflamnat, NuBBE, and AC Discovery Databases Was Determined Using the QED Index

Data Collection of NPs and SMSs

To obtain compounds with potential usefulness in VD, in this work, we used different databases that contained both NPs and SMSs. NPs were obtained from databases like InflamNat, BIOFACQUIM, NuBBE, and AC Discovery. SMS were only obtained from the AC Discovery database.
We analyzed only the molecules that were available on that server.

Determination of the Drug-Likeness of NP and SMSs

To determine the drug-likeness of NP and SMS compounds, a drug-likeness estimator based on the Quantitative Estimation of Drug-likeness (QED) was used. This one is part of RDKit, a piece of open-source cheminformatics software (version 2023.03) [95] (https://www.rdkit.org/docs/source/rdkit.Chem.QED.html (accessed on 10 April 2023)). The QED is represented as an integrative score to evaluate the compounds’ preference to be considered a hit [77]. This method quantifies the drug-likeness, considering the following molecular descriptors: molecular weight (Da), Ghose–Crippen–Viswanadhan octanol–water partition coefficient (AlogP), the number of H-acceptors (HBA), the number of H-donors (HBD), the number of rotatable bonds, total polar surface area (TPSA), and aromatic bound account. We analyzed the molecular descriptors of the reference compounds, NPs, and SMSs to identify similar patterns in the datasets. A histogram and box plot were accomplished to represent graphically this distribution. The score of QED uses the molecular descriptors mentioned above. We selected the NPs and SMSs with high QED scores (over 0.5).

Filtering of the NP and SMS Dataset

To eliminate NPs and SMSs with potential non-desirable effects, we identified the functions related to tumorigenesis, reproduction, and other effects [96]. Osiris Datawarrior software V5.2.1 was employed to determine these potential non-desirable effects. Because of the large number of drug-like molecules available, we used PCA to select molecules that resembled phlebotonics in their chemical space. We utilized both two- and three-dimensional systems in our analysis comparing the NPs and SMSs versus reference compounds [94,97]. We used Osiris Datawarrior software V5.2.1 for the determination of PCA. Only the molecules with similar chemical space to reference compounds were considered for the next step of the work.

4.1.3. The Selection of Compounds with Potential Usefulness in VD Based on Comparison of the Fingerprints of the Drug-like NPs, SMSs, and Reference Compounds

Fingerprints of Drug-like NPs, SMSs, and Reference Compounds

We used fingerprint analysis to identify molecules similar to the reference compounds present in our drug-like NPs in the SMS dataset. Molecular fingerprints are one of the most common representations of chemical structures in chemoinformatics. They represent chemical information in any chemical entity through binary vectors. Commonly, the binary fingerprint is the representation of a chemical structure. The vector’s positions show the presence (1) or absence (0) of predetermined features in the fingerprint design [98].They are essential cheminformatic tools for virtual screening and mapping chemical space [99]. This method was supplemented with clustering of chemical compounds by the similarity of their molecular fingerprints to identify similar structures in an extensive set of similar data, as previously reported [78]. Reference compounds were compared to drug-like molecules using a fingerprint analysis. To achieve this structural comparison, we used ChemmineR and rcdk in RStudio version 3.4 [78]. Because the molecular structures from all datasets were aromatic compounds, the “extended” method is preferred. An extended number of fingerprints with a default length of 1024 (number of bits) were used. The preferred analysis gave us a higher level of detail in the structure analysis compared to the MACCS-type fingerprints. By using the Tanimoto coefficient, we compared the similarity of NPs, SMSs, and reference compounds through cluster analysis. Later, the molecules were classified into three categories based on their distances using Ward’s clustering method, a widely used algorithm in drug discovery. To confirm the clustering results, the Dunn index and the silhouette coefficient were used [78]. The cluster, which had a reference compound, was employed to search for compounds with potential usefulness in VD and generate a compound–target network.

4.1.4. Selection of the Best Compounds with Potential Usefulness in VD through Their Multi-Target Capacity

Target Search of the Compounds Obtained Using Fingerprint Analysis

The pharmacological networks were constructed by searching for the reported targets of each drug-like compound in the PubChem server. Only drug-like molecules with reported targets were considered for this phase of the study.

Compound–Target Network Generation

We used network pharmacology to select compounds with potential usefulness that interact with VD networks. By linking drug-like compounds and their reported targets with genes from VD-related pathways in Wikipathways, this network was formed. The pathways used in this analysis were the inflammatory response pathway, cellular response to hypoxia, PI3K-AKT-NFKB pathway, ion channel transport, IL-18 signaling pathway, VEGF-VEGFR-2, senescence-associated secretory phenotype (SASP), oxidative stress, intrinsic pathways for apoptosis, and angiogenesis [61,62,63,64,65,66,67,68]. For more detail on the targets, pharmacological functions, and genes data from the VD network, see the Supplementary materials.
Cytoscape software version 3.8.2 was used to construct the compound–target network [100]. Only multi-target compounds participating in over three pathways were considered as compounds with potential usefulness in VD.

Selection of the NPs and SMSs through Their Availability, Accessibility(cost), and Ease of Administration

The best candidates with potential usefulness in VD obtained using network pharmacology were evaluated in an in vivo model. We started with various criteria based on their availability, accessibility, and administration. The information was collected from the literature and commercial sites. The Scheme 1 shows a summary of the strategy followed in this work.

4.2. Experimental Phase

4.2.1. Materials

Berberine chloride (Sigma Aldrich Cat. B3251-10G) (St. Louis, MO, USA) was dissolved in normal saline-dimethyl sulfoxide (DMSO) (9:1). β-escin (Sigma Aldrich Cat. E-1378 10G) was dissolved in phosphates buffer to a pH of 7. Urethane (Sigma Aldrich Cat. U-2500 500G)/chloralose (Sigma Aldrich Cat. C-0128 500G) (1000/100 mg) was dissolved in water.

4.2.2. Animals

Male Wistar rats weighing 300–350 g were used for all experiments. The animals were maintained at room temperature (21–23 °C) on a 12 h light:12 h dark cycle and fed with pellets food (5001 Rodent Laboratory Chow) ad libitum. Their care was in line with Mexican standards [101,102] and international official guidelines [103]. The procedures were approved by the Research Ethics Committee of the Faculty of Medicine, UNAM (Project No. 008-CIC-2022).

4.2.3. Register of PAM, HR, and VP in Anesthetized Rats

Wistar rats were anesthetized with urethane/chloralose (1000/100 mg/kg i.p. 100 g/0.2 mL). Using extreme care to prevent skin irritation, the limbs and throat were shaved. Through vertical incisions, the veins of the thigh, such as the superficial epigastric, proximal caudal, medial proximal geniculate, and popliteal [104], were surgically exposed and subsequently dissected. These veins of the left hind limb were ligated with microvascular clips (11 mm) (George. Tiemann Co., Long Island City, New York, NY, USA). The veins of the right hind limb were not ligated (sham control group).
In all animals, PE 50 cannulas (BD IntramedicTM, Sparks, MD, USA) were inserted into the right carotid artery, left jugular vein, and right and left femoral veins to measure medium arterial pressure, heart rate, and venous pressure in both hind limbs, drug administration and recording of the venous pressure, respectively. In the trachea, a PE 240 cannula was inserted to facilitate spontaneous breathing. Model P231D pressure transducers (Gould-Statham Instruments, Oxnard, CA, USA) were connected to the cannulas inserted in the carotid artery and femoral vein to record MAP, HR, and PV [105,106]. The signal from the transducer was electronically damped and recorded using a Model 79 Grass polygraph (Grass Instrument, Quincy, MA, USA). HR was recorded through another channel from the polygraph with a Grass 7P4 tachograph triggered by the pulse waves from the unfiltered transducer signal, and the pressure change was transmitted to The software Lab-ViewTM 21 SP1 (National Instruments, Austin, TX, USA) via an interface NI USB-6009 multifunction DAQ (National Instruments, USA). The animals’ temperature was regulated at 38 °C through the utilization of a heating table connected to a regulation unit. All experiments used an electronic thermometer to regulate the temperature probe control of the heating table. The animals were stabilized for 15 min after the insertion of the cannula.

4.2.4. Dose–Response Curves of the Cardiovascular System to Berberine, and β-Escin

After 15 min of stabilization, MAP, HR, and VP readings were taken. Successive injections independent of berberine (0.001, 0.003, 0.01, 0.03, 0.1, 0.3, 1 mg/kg, i.v., and β-escin (0.56, 1, 17, 31, 56 and 10 mg/Kg) were administered. A time interval separated each bolus injection to allow recovery of the stability of MAP, HR, and PV in both hind limbs. The group comprised eight animals receiving berberine and six receiving β-escin treatment. The obtained results were used to construct dose–response curves.
For more details, see the section before the Introduction. All the acronyms/initialisms/abbreviations are deciphered within it.

5. Conclusions

Through our analysis of chemical structures and compound–target networks, we have identified drug-like compounds with potential usefulness in VD. We identified 18 compounds with potential usefulness in VD. Only the following compounds have a previously reported method for experimental use: 7-hydroxyflavone, acacetin, berberine chloride, biochanin A, chrysin, formononetin, and sulfuretin. Berberine, our chosen candidate, significantly reduced peripheral acute venous hypertension in anesthetized rats. This effect shows the start of confirmation to our approach for discovering new drugs using virtual screening and network pharmacology. Additional preclinical experiments are needed to determine whether berberine is effective for treating chronic peripheral venous hypertension.
Our results show that testing the most promising candidates in future experiments could lead to the discovery of new VD drugs. VD is a complex illness for which there is a lack of experimental models.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28247937/s1, Figure S1: Principal component analysis (PCA) scoreplot of the analyzed reference compounds useful in VD with respect to their molecular descriptors; Figure S2: Illustrative representation of the fingerprint analyzes of each reference compound.

Author Contributions

Conceptualization, G.A.M.-G., J.S.-M. and O.S.B.-V.; methodology, O.S.B.-V.; software, O.S.B.-V.; validation, J.L.E.-R. and O.S.B.-V.; formal analysis, O.S.B.-V.; investigation, O.S.B.-V.; resources, O.F.C.-O.; data curation, O.S.B.-V.; writing—original draft preparation, O.S.B.-V.; writing—review and editing, G.A.M.-G.; visualization, J.L.E.-R. and O.F.C.-O.; supervision, G.A.M.-G.; project administration, G.A.M.-G.; funding acquisition, O.F.C.-O. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by División de Investigación, Facultad de Medicina, UNAM. Project number: 008-CIC-2022.

Institutional Review Board Statement

This study complied with Mexican standards NOM-062-ZOO-1999; NOM-087-SEMARNAT-SSA1-2002 and Guide for the Care and Use of Laboratory Animals, revised in 2011. Animal use and management align with the approved project by the Internal Committee for the Care and Use of Laboratory Animals (CICUAL), School of Medicine, National Autonomous University of Mexico (UNAM) (Project No. 008-CIC-2022) for the Phytopharmacology Laboratory, Department of Pharmacology, School of Medicine, National Autonomous University of Mexico (UNAM).

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available in the Supplementary Materials.

Acknowledgments

UNAM’s School of Medicine assists with implementing project 008-CIC-2022. Oscar Salvador Barrera-Vázquez is grateful to the Dirección General de Asuntos del Personal Académico (DGAPA), Universidad Nacional Autónoma de México, for the Subprograma de Incorporación de Jóvenes Académicos de Carrera (SIJA).

Conflicts of Interest

The authors have no conflict of interest to declare.

Abbreviations

Venous disease (VD); computer-aided drug design (CADD); natural products (NPs); semi-synthetics (SMSs); Quantitative Estimate of Druglikeness (QED); micronized purified flavonoid fraction (MPFF); three-dimensional principal component analysis (3D-PCA); principal component analysis (PCA); one-dimensional (1D); two-dimensional (2D); IUPAC (International Union of Pure and Applied Chemistry); venous pressure (VP); right hind limb (RHL); left hind limb (LHL); mean arterial pressure (MAP); heart rate (HR); acute venous hypertension (AVH); molecular weight range (Da); Ghose–Crippen–Viswanadhan octanol–water partition co-efficient (AlogP); number of H-acceptors (HBA); number of H-donors (HBD); number of rotatable bonds, total polar surface area (TPSA; standard error of the mean (SEM); analysis of variances (ANOVA); AMP-activated protein kinase (AMPK); nuclear factor κB (NF-κB); Vascular Endothelial Growth Factor Receptor-2 (VEGF-VEGFR-2); Senescence-Associated Secretory Phenotype (SASP); mitogen-regulatory Silent Information System Activated Protein Kinase 1 (SIRT-1); Hypoxia-Inducible Factor 1α (HIF-1α); Phosphoinositide Vascular Endothelial Growth Factor 3-Kinase (PI3K); Protein Kinase B (Akt); Janus Kinase 2 (JAK -2); vascular smooth muscle cells (VSMCs); potassium voltage-gated channel subfamily H member 2 (KCNH2); nuclear factor, erythroid 2-like 2 (NFE2L2); mammalian target of rapamycin (mTOR); vascular endothelial growth factor (VEGF); interleukin-1 receptor antagonist (IL-1RA); silhouette coefficient (Sc).

References

  1. Kurz, X.; Kahn, S.R.; Abenhaim, L.; Clement, D.; Norgren, L.; Baccaglini, U.; Berard, A.; Cooke, J.P.; Cornu-Thenard, A.; Depairon, M.; et al. Chronic venous disorders of the leg: Epidemiology, outcomes, diagnosis and management. Summary of an evidence-based report of the VEINES task force. Venous Insufficiency Epidemiologic and Economic Studies. Int. Angiol. J. Int. Union. Angiol. 1999, 18, 83–102. [Google Scholar]
  2. Kim, Y.; Png, C.Y.M.; Sumpio, B.J.; DeCarlo, C.S.; Dua, A. Defining the human and health care costs of chronic venous insufficiency. Semin. Vasc. Surg. 2021, 34, 59–64. [Google Scholar] [CrossRef] [PubMed]
  3. Javier, J.J.; Ortiz, P. Treatment of chronic venous insufficiency in Latin America. J. Vasc. Surg. Venous Lymphat. Disord. 2020, 8, 667–675. [Google Scholar] [CrossRef] [PubMed]
  4. Gobierno de México. Anuario de Morbilidad 1984–2022. Incidencia de Enfermedad por Grupo de Edad y Entidad Federativa. Available online: https://epidemiologia.salud.gob.mx/anuario/html/incidencia_enfermedad.html (accessed on 11 June 2023).
  5. WHO. Anatomical Therapeutic Chemical (ATC) System Classification. Available online: http://www.whocc.no/atc_ddd_index/?code=C05 (accessed on 10 November 2023).
  6. Martinez-Zapata, M.J.; Vernooij, R.W.; Simancas-Racines, D.; Uriona Tuma, S.M.; Stein, A.T.; Moreno Carriles, R.M.M.; Vargas, E.; Bonfill Cosp, X. Phlebotonics for venous insufficiency. Cochrane Database Syst. Rev. 2020, 11, Cd003229. [Google Scholar] [CrossRef] [PubMed]
  7. Tsouderos, Y. Venous tone: Are the phlebotonic properties predictive of a therapeutic benefit? A comprehensive view of our experience with Daflon 500 mg. Z. Fur Kardiol. 1991, 80 (Suppl. S7), 95–101. [Google Scholar]
  8. Behar, A.; Nathan, P.; Lavieuville, M.; Allaert, F.A. Effect of veinotonyl 75 on the capillary permeability test using technetium albumin in cyclic orthostatic edemas. Phlebologie 1993, 46, 721–731. [Google Scholar]
  9. Prieto-Martínez, F.D.; López-López, E.; Eurídice Juárez-Mercado, K.; Medina-Franco, J.L. Chapter 2—Computational Drug Design Methods—Current and Future Perspectives. In In Silico Drug Design; Roy, K., Ed.; Academic Press: Cambridge, MA, USA, 2019; pp. 19–44. [Google Scholar]
  10. Barrera-Vázquez, O.S.; Magos-Guerrero, G.A.; Escobar-Ramírez, J.L.; Gomez-Verjan, J.C. Natural Products as a Major Source of Candidates for Potential Senolytic Compounds obtained by in silico Screening. Med. Chem. 2023, 19, 653–668. [Google Scholar] [CrossRef]
  11. Cazaubon, M.; Benigni, J.P.; Steinbruch, M.; Jabbour, V.; Gouhier-Kodas, C. Is There a Difference in the Clinical Efficacy of Diosmin and Micronized Purified Flavonoid Fraction for the Treatment of Chronic Venous Disorders? Review of Available Evidence. Vasc. Health Risk Manag. 2021, 17, 591–600. [Google Scholar] [CrossRef]
  12. Gallelli, L. Escin: A review of its anti-edematous, anti-inflammatory, and venotonic properties. Drug Des. Devel Ther. 2019, 13, 3425–3437. [Google Scholar] [CrossRef]
  13. Sumaiya, S.; Naved, T.; Sharma, A.; Sarwat, M. Chapter 1—Amelioration of Liver Ailments by Saffron (Crocus sativus) and Its Secondary Metabolites. In Saffron; Sarwat, M., Sumaiya, S., Eds.; Academic Press: Cambridge, MA, USA, 2020; pp. 1–20. [Google Scholar]
  14. Riboldi, E.; Sica, A. Chapter 4—Modulation of Innate Immunity by Hypoxia. In The Innate Immune Response to Noninfectious Stressors; Amadori, M., Ed.; Academic Press: Cambridge, MA, USA, 2016; pp. 81–106. [Google Scholar]
  15. Otsuka, S.; Melis, N.; Dutta, D.; Weigert, R.; Ashwell, J.D. Immunosuppressive activity of cyclosporin A in vivo via NFAT-independent inhibition of TCR signaling. J. Immunol. 2020, 204, 161.113. [Google Scholar] [CrossRef]
  16. Mesa-Arango, A.C.; Scorzoni, L.; Zaragoza, O. It only takes one to do many jobs: Amphotericin B as antifungal and immunomodulatory drug. Front. Microbiol. 2012, 3, 286. [Google Scholar] [CrossRef] [PubMed]
  17. Su, Y.-C.; Liao, H.F.; Yu, C.C.; Lin, Y.C. Anticancer activity of the mTOR inhibitor (everolimus) and dual mTORC1/mTORC2 Inhibitor (AZD2014) on mouse lymphocytic leukemia both in vitro and in vivo. Ann. Oncol. 2016, 27, vi8. [Google Scholar] [CrossRef]
  18. Lamming, D.W. Inhibition of the Mechanistic Target of Rapamycin (mTOR)-Rapamycin and Beyond. Cold Spring Harb. Perspect. Med. 2016, 6, a025924. [Google Scholar] [CrossRef] [PubMed]
  19. Viziteu, E.; Grandmougin, C.; Goldschmidt, H.; Seckinger, A.; Hose, D.; Klein, B.; Moreaux, J. Chetomin, targeting HIF-1α/p300 complex, exhibits antitumour activity in multiple myeloma. Br. J. Cancer 2016, 114, 519–523. [Google Scholar] [CrossRef] [PubMed]
  20. Zhang, W.J.; Hufnagl, P.; Binder, B.R.; Wojta, J. Antiinflammatory activity of astragaloside IV is mediated by inhibition of NF-kappaB activation and adhesion molecule expression. Thromb. Haemost. 2003, 90, 904–914. [Google Scholar] [CrossRef]
  21. Lv, Q.; Rong, N.; Liu, L.J.; Xu, X.L.; Liu, J.T.; Jin, F.X.; Wang, C.M. Antitumoral Activity of (20R)- and (20S)-Ginsenoside Rh2 on Transplanted Hepatocellular Carcinoma in Mice. Planta Medica 2016, 82, 705–711. [Google Scholar] [CrossRef]
  22. Min, K.J.; Kwon, T.K. Anticancer effects and molecular mechanisms of epigallocatechin-3-gallate. Integr. Med. Res. 2014, 3, 16–24. [Google Scholar] [CrossRef]
  23. Panat, N.A.; Singh, B.G.; Maurya, D.K.; Sandur, S.K.; Ghaskadbi, S.S. Troxerutin, a natural flavonoid binds to DNA minor groove and enhances cancer cell killing in response to radiation. Chem. Biol. Interact. 2016, 251, 34–44. [Google Scholar] [CrossRef]
  24. Rehn, D.; Hennings, G.; Nocker, W.; Diebschlag, W. Time course of the anti-oedematous effect of O-(beta-hydroxyethyl)-rutosides in healthy volunteers. Eur. J. Clin. Pharmacol. 1991, 40, 625–627. [Google Scholar] [CrossRef]
  25. Wang, S.X.; Hu, L.M.; Gao, X.M.; Guo, H.; Fan, G.W. Anti-inflammatory activity of salvianolic acid B in microglia contributes to its neuroprotective effect. Neurochem. Res. 2010, 35, 1029–1037. [Google Scholar] [CrossRef]
  26. Blando, F.; Calabriso, N.; Berland, H.; Maiorano, G.; Gerardi, C.; Carluccio, M.A.; Andersen, Ø.M. Radical Scavenging and Anti-Inflammatory Activities of Representative Anthocyanin Groupings from Pigment-Rich Fruits and Vegetables. Int. J. Mol. Sci. 2018, 19, 169. [Google Scholar] [CrossRef] [PubMed]
  27. Ganeshpurkar, A.; Saluja, A.K. The Pharmacological Potential of Rutin. Saudi Pharm. J. 2017, 25, 149–164. [Google Scholar] [CrossRef] [PubMed]
  28. Monreal, M.; Callejas, J.M.; Martorell, A.; Sahuquillo, J.C.; Contel, E. Prevention of Post-Thrombotic Syndrome with Hidrosmina: A Prospective Pilot Study. Phlebology 1997, 12, 21–24. [Google Scholar] [CrossRef]
  29. Naso, L.; Martínez, V.R.; Lezama, L.; Salado, C.; Valcarcel, M.; Ferrer, E.G.; Williams, P.A.M. Antioxidant, anticancer activities and mechanistic studies of the flavone glycoside diosmin and its oxidovanadium(IV) complex. Interactions with bovine serum albumin. Bioorganic Med. Chem. 2016, 24, 4108–4119. [Google Scholar] [CrossRef] [PubMed]
  30. Zheng, Y.; Zhang, R.; Shi, W.; Li, L.; Liu, H.; Chen, Z.; Wu, L. Metabolism and pharmacological activities of the natural health-benefiting compound diosmin. Food Funct. 2020, 11, 8472–8492. [Google Scholar] [CrossRef]
  31. Hajialyani, M.; Hosein Farzaei, M.; Echeverría, J.; Nabavi, S.M.; Uriarte, E.; Sobarzo-Sánchez, E. Hesperidin as a Neuroprotective Agent: A Review of Animal and Clinical Evidence. Molecules 2019, 24, 648. [Google Scholar] [CrossRef]
  32. Kim, B.; Lee, J.H.; Seo, M.J.; Eom, S.H.; Kim, W. Linarin down-regulates phagocytosis, pro-inflammatory cytokine production, and activation marker expression in RAW264.7 macrophages. Food Sci. Biotechnol. 2016, 25, 1437–1442. [Google Scholar] [CrossRef]
  33. Sroka, Z.; Fecka, I.; Cisowski, W. Antiradical and anti-H2O2 properties of polyphenolic compounds from an aqueous peppermint extract. Z. Fur Naturforschung. C J. Biosci. 2005, 60, 826–832. [Google Scholar] [CrossRef]
  34. Salehi, B.; Fokou, P.V.T.; Sharifi-Rad, M.; Zucca, P.; Pezzani, R.; Martins, N.; Sharifi-Rad, J. The Therapeutic Potential of Naringenin: A Review of Clinical Trials. Pharmaceuticals 2019, 12, 11. [Google Scholar] [CrossRef]
  35. Ahmadi, H.; Nayeri, Z.; Minuchehr, Z.; Sabouni, F.; Mohammadi, M. Betanin purification from red beetroots and evaluation of its anti-oxidant and anti-inflammatory activity on LPS-activated microglial cells. PLoS ONE 2020, 15, e0233088. [Google Scholar] [CrossRef]
  36. Cortez-Espinosa, N.; Aviña-Verduzco, J.A.; Ramírez-Chávez, E.; Molina-Torres, J.; Ríos-Chávez, P. Valine and phenylalanine as precursors in the biosynthesis of alkamides in Acmella radicans. Nat. Prod. Commun. 2011, 6, 857–861. [Google Scholar] [CrossRef]
  37. Tseng, C.F.; Iwakami, S.; Mikajiri, A.; Shibuya, M.; Hanaoka, F.; Ebizuka, Y.; Padmawinata, K.; Sankawa, U. Inhibition of in vitro prostaglandin and leukotriene biosyntheses by cinnamoyl-beta-phenethylamine and N-acyldopamine derivatives. Chem. Pharm. Bull. 1992, 40, 396–400. [Google Scholar] [CrossRef] [PubMed]
  38. Achenbach, H.; Waibel, R.; Addae-Mensah, I. Lignans and other constituents from Carissa edulis. Phytochemistry 1983, 22, 749–753. [Google Scholar] [CrossRef]
  39. Patil, K.K.; Gacche, R.N. Inhibition of glycation and aldose reductase activity using dietary flavonoids: A lens organ culture studies. Int. J. Biol. Macromol. 2017, 98, 730–738. [Google Scholar] [CrossRef] [PubMed]
  40. Tamaki, H.; Satoh, H.; Hori, S.; Ohtani, H.; Sawada, Y. Inhibitory effects of herbal extracts on breast cancer resistance protein (BCRP) and structure-inhibitory potency relationship of isoflavonoids. Drug Metab. Pharmacokinet. 2010, 25, 170–179. [Google Scholar] [CrossRef] [PubMed]
  41. Abdel-Kader, M.S.; Amer, M.E.; Tang, S.; Kingston, D.G. Two new isoflavone derivatives from the roots of an Egyptian collection of Lotus polyphyllos. Nat. Prod. Res. 2006, 20, 922–926. [Google Scholar] [CrossRef] [PubMed]
  42. Singh, M.; Kaur, M.; Silakari, O. Flavones: An important scaffold for medicinal chemistry. Eur. J. Med. Chem. 2014, 84, 206–239. [Google Scholar] [CrossRef]
  43. Sengupta, B.; Sahihi, M.; Dehkhodaei, M.; Kelly, D.; Arany, I. Differential roles of 3-Hydroxyflavone and 7-Hydroxyflavone against nicotine-induced oxidative stress in rat renal proximal tubule cells. PLoS ONE 2017, 12, e0179777. [Google Scholar] [CrossRef]
  44. Vinutha, K.; Pavan, G.; Pattar, S.; Kumari, N.S.; Vidya, S.M. Aqueous extract from Madhuca indica bark protects cells from oxidative stress caused by electron beam radiation: In vitro, in vivo and in silico approach. Heliyon 2019, 5, e01749. [Google Scholar] [CrossRef]
  45. Adhikari, B.; Marasini, B.P.; Rayamajhee, B.; Bhattarai, B.R.; Lamichhane, G.; Khadayat, K.; Adhikari, A.; Khanal, S.; Parajuli, N. Potential roles of medicinal plants for the treatment of viral diseases focusing on COVID-19: A review. Phytother. Res. PTR 2021, 35, 1298–1312. [Google Scholar] [CrossRef]
  46. Lee, H.J.; Lee, J.H.; Hwang, B.Y.; Kim, H.S.; Lee, J.J. Fungal metabolites, asterric acid derivatives inhibit vascular endothelial growth factor (VEGF)-induced tube formation of HUVECs. J. Antibiot. 2002, 55, 552–556. [Google Scholar] [CrossRef] [PubMed]
  47. Ai, X.; Yu, P.; Peng, L.; Luo, L.; Liu, J.; Li, S.; Lai, X.; Luan, F.; Meng, X. Berberine: A Review of its Pharmacokinetics Properties and Therapeutic Potentials in Diverse Vascular Diseases. Front. Pharmacol. 2021, 12, 762654. [Google Scholar] [CrossRef]
  48. Gour, A.; Manhas, D.; Bag, S.; Gorain, B.; Nandi, U. Flavonoids as potential phytotherapeutics to combat cytokine storm in SARS-CoV-2. Phytother. Res. 2021, 35, 4258–4283. [Google Scholar] [CrossRef] [PubMed]
  49. Fattori, V.; Hohmann, M.S.; Rossaneis, A.C.; Pinho-Ribeiro, F.A.; Verri, W.A. Capsaicin: Current Understanding of Its Mechanisms and Therapy of Pain and Other Pre-Clinical and Clinical Uses. Molecules 2016, 21, 844. [Google Scholar] [CrossRef] [PubMed]
  50. Gonçalves, L.M.; Valente, I.M.; Rodrigues, J.A. An overview on cardamonin. J. Med. Food 2014, 17, 633–640. [Google Scholar] [CrossRef] [PubMed]
  51. Sharma, P.; Kumari, A.; Gulati, A.; Krishnamurthy, S.; Hemalatha, S. Chrysin isolated from Pyrus pashia fruit ameliorates convulsions in experimental animals. Nutr. Neurosci. 2019, 22, 569–577. [Google Scholar] [CrossRef] [PubMed]
  52. Tay, K.C.; Tan, L.T.; Chan, C.K.; Hong, S.L.; Chan, K.G.; Yap, W.H.; Pusparajah, P.; Lee, L.H.; Goh, B.H. Formononetin: A Review of Its Anticancer Potentials and Mechanisms. Front. Pharmacol. 2019, 10, 820. [Google Scholar] [CrossRef]
  53. Ea, S.; Giacometti, S.; Ciccolini, J.; Akhmedjanova, V.; Aubert, C. Cytotoxic effects of haplamine and its major metabolites on human cancer cell lines. Planta Medica 2008, 74, 1265–1268. [Google Scholar] [CrossRef]
  54. Lawrence, P.F. Chapter 78—Pharmacologic Adjuncts to Endovascular Procedures. In Endovascular Surgery, 4th ed.; Moore, W.S., Ahn, S.S., Eds.; W.B. Saunders: Philadelphia, PA, USA, 2011; pp. 807–813. [Google Scholar]
  55. Mariadoss, A.V.A.; Vinyagam, R.; Rajamanickam, V.; Sankaran, V.; Venkatesan, S.; David, E. Pharmacological Aspects and Potential Use of Phloretin: A Systemic Review. Mini Rev. Med. Chem. 2019, 19, 1060–1067. [Google Scholar] [CrossRef]
  56. Roh, K.; Kim, S.; Kang, H.; Ku, J.M.; Park, K.W.; Lee, S. Sulfuretin has therapeutic activity against acquired lymphedema by reducing adipogenesis. Pharmacol. Res. 2017, 121, 230–239. [Google Scholar] [CrossRef]
  57. Pascarella, L.; Schmid-Schönbein, G.W.; Bergan, J. An animal model of venous hypertension: The role of inflammation in venous valve failure. J. Vasc. Surg. 2005, 41, 303–311. [Google Scholar] [CrossRef] [PubMed]
  58. Jones, G.T.; Grant, M.W.; Thomson, I.A.; Hill, B.G.; van Rij, A.M. Characterization of a porcine model of chronic superficial varicose veins. J. Vasc. Surg. 2009, 49, 1554–1561. [Google Scholar] [CrossRef] [PubMed]
  59. das Graças, C.d.S.M.; Cyrino, F.Z.; de Carvalho, J.J.; Blanc-Guillemaud, V.; Bouskela, E. Protective Effects of Micronized Purified Flavonoid Fraction (MPFF) on a Novel Experimental Model of Chronic Venous Hypertension. Eur. J. Vasc. Endovasc. Surg. Off. J. Eur. Soc. Vasc. Surg. 2018, 55, 694–702. [Google Scholar] [CrossRef] [PubMed]
  60. Prieto-Martínez, F.D.; Norinder, U.; Medina-Franco, J.L. Cheminformatics Explorations of Natural Products. Prog. Chem. Org. Nat. Prod. 2019, 110, 1–35. [Google Scholar] [CrossRef]
  61. Atta, H.M. Varicose veins: Role of mechanotransduction of venous hypertension. Int. J. Vasc. Med. 2012, 2012, 538627. [Google Scholar] [CrossRef]
  62. Grant, Y.; Onida, S.; Davies, A. Genetics in chronic venous disease. Phlebology 2017, 32, 3–5. [Google Scholar] [CrossRef]
  63. Serra, R.; Ssempijja, L.; Provenzano, M.; Andreucci, M. Genetic biomarkers in chronic venous disease. Biomark. Med. 2020, 14, 75–80. [Google Scholar] [CrossRef]
  64. Ligi, D.; Croce, L.; Mannello, F. Chronic Venous Disorders: The Dangerous, the Good, and the Diverse. Int. J. Mol. Sci. 2018, 19, 2544. [Google Scholar] [CrossRef]
  65. Castro-Ferreira, R.; Cardoso, R.; Leite-Moreira, A.; Mansilha, A. The Role of Endothelial Dysfunction and Inflammation in Chronic Venous Disease. Ann. Vasc. Surg. 2018, 46, 380–393. [Google Scholar] [CrossRef]
  66. Mansilha, A.; Sousa, J. Pathophysiological Mechanisms of Chronic Venous Disease and Implications for Venoactive Drug Therapy. Int. J. Mol. Sci. 2018, 19, 1669. [Google Scholar] [CrossRef]
  67. Filis, K.; Kavantzas, N.; Isopoulos, T.; Antonakis, P.; Sigalas, P.; Vavouranakis, E.; Sigala, F. Increased vein wall apoptosis in varicose vein disease is related to venous hypertension. Eur. J. Vasc. Endovasc. Surg. 2011, 41, 533–539. [Google Scholar] [CrossRef] [PubMed]
  68. Mikuła-Pietrasik, J.; Uruski, P.; Aniukiewicz, K.; Sosińska, P.; Krasiński, Z.; Tykarski, A.; Książek, K. Serum from Varicose Patients Induces Senescence-Related Dysfunction of Vascular Endothelium Generating Local and Systemic Proinflammatory Conditions. Oxid. Med. Cell Longev. 2016, 2016, 2069290. [Google Scholar] [CrossRef] [PubMed]
  69. Singh, S.; Genilloud, O.; Peláez, F. Terrestrial Microorganisms–Filamentous Bacteria. In Comprehensive Natural Products II; Liu, H.-W., Mander, L., Eds.; Elsevier: Oxford, UK, 2010; pp. 109–140. [Google Scholar]
  70. Moreira, D.M.; da Silva, R.L.; Vieira, J.L.; Fattah, T.; Lueneberg, M.E.; Gottschall, C.A. Role of vascular inflammation in coronary artery disease: Potential of anti-inflammatory drugs in the prevention of atherothrombosis. Inflammation and anti-inflammatory drugs in coronary artery disease. Am. J. Cardiovasc. Drugs Drugs Devices Other Interv. 2015, 15, 1–11. [Google Scholar] [CrossRef] [PubMed]
  71. Fujita, T.; Yoshimoto, T.; Kajiya, M.; Ouhara, K.; Matsuda, S.; Takemura, T.; Akutagawa, K.; Takeda, K.; Mizuno, N.; Kurihara, H. Regulation of defensive function on gingival epithelial cells can prevent periodontal disease. Jpn. Dent. Sci. Rev. 2018, 54, 66–75. [Google Scholar] [CrossRef] [PubMed]
  72. Vitiello, D.; Neagoe, P.E.; Sirois, M.G.; White, M. Effect of everolimus on the immunomodulation of the human neutrophil inflammatory response and activation. Cell Mol. Immunol. 2015, 12, 40–52. [Google Scholar] [CrossRef] [PubMed]
  73. Choi, W.Y.; Lim, H.W.; Lim, C.J. Anti-inflammatory, antioxidative and matrix metalloproteinase inhibitory properties of 20(R)-ginsenoside Rh2 in cultured macrophages and keratinocytes. J. Pharm. Pharmacol. 2013, 65, 310–316. [Google Scholar] [CrossRef]
  74. Belcaro, G.; Rosaria Cesarone, M.; Ledda, A.; Cacchio, M.; Ruffini, I.; Ricci, A.; Ippolito, E.; Di Renzo, A.; Dugall, M.; Corsi, M.; et al. O-(beta-hydroxyethyl)-rutosides systemic and local treatment in chronic venous disease and microangiopathy: An independent prospective comparative study. Angiology 2008, 59 (Suppl. S1), 7s–13s. [Google Scholar] [CrossRef]
  75. Wu, Y.; Xu, S.; Tian, X.Y. The Effect of Salvianolic Acid on Vascular Protection and Possible Mechanisms. Oxid. Med. Cell Longev. 2020, 2020, 5472096. [Google Scholar] [CrossRef]
  76. Lichota, A.; Gwozdzinski, L.; Gwozdzinski, K. Therapeutic potential of natural compounds in inflammation and chronic venous insufficiency. Eur. J. Med. Chem. 2019, 176, 68–91. [Google Scholar] [CrossRef]
  77. Bickerton, G.R.; Paolini, G.V.; Besnard, J.; Muresan, S.; Hopkins, A.L. Quantifying the chemical beauty of drugs. Nat. Chem. 2012, 4, 90–98. [Google Scholar] [CrossRef]
  78. Voicu, A.; Duteanu, N.; Voicu, M.; Vlad, D.; Dumitrascu, V. The rcdk and cluster R packages applied to drug candidate selection. J. Cheminform 2020, 12, 3. [Google Scholar] [CrossRef] [PubMed]
  79. Talevi, A. Multi-target pharmacology: Possibilities and limitations of the “skeleton key approach” from a medicinal chemist perspective. Front. Pharmacol. 2015, 6, 205. [Google Scholar] [CrossRef] [PubMed]
  80. Oklei, D. Systemic phlebotropic drugs in pharmacotherapy of chronic venous insufficiency of the lower extremities. Vìsnik Farm. 2015, 4, 74–77. [Google Scholar] [CrossRef]
  81. Le Dévéhat, C.; Khodabandehlou, T.; Vimeux, M.; Kempf, C. Evaluation of haemorheological and microcirculatory disturbances in chronic venous insufficiency: Activity of Daflon 500 mg. Int. J. Microcirc. Clin. Exp. 1997, 17 (Suppl. S1), 27–33. [Google Scholar] [CrossRef] [PubMed]
  82. Ramelet, A.-A. Chapter 14—Venoactive Drugs. In Sclerotherapy, 5th ed.; Goldman, M.P., Guex, J.-J., Weiss, R.A., Eds.; W.B. Saunders: Edinburgh, UK, 2011; pp. 369–377. [Google Scholar]
  83. May, J.M.; Harrison, F.E. Role of vitamin C in the function of the vascular endothelium. Antioxid. Redox Signal. 2013, 19, 2068–2083. [Google Scholar] [CrossRef] [PubMed]
  84. Berthet, P.; Farine, J.C.; Barras, J.P. Calcium dobesilate: Pharmacological profile related to its use in diabetic retinopathy. Int. J. Clin. Pract. 1999, 53, 631–636. [Google Scholar] [CrossRef] [PubMed]
  85. Dardi, P.; Dos Reis Costa, D.E.F.; Assunção, H.C.R.; Rossoni, L.V. Venous endothelial function in cardiovascular disease. Biosci. Rep. 2022, 42, BSR20220285. [Google Scholar] [CrossRef] [PubMed]
  86. Xia, M.L.; Zhou, X.M.; Yao, H.; Jiang, H.D.; Bruce, I.C.; Wei, E.Q.; Xia, Q. Rutin-induced endothelium-dependent vasorelaxation in rat aortic Rings and the underlying mechanism. In Proceedings of the 2005 IEEE Engineering in Medicine and Biology 27th Annual Conference, Shanghai, China, 17–18 January 2006; Volume 2005, pp. 5595–5597. [Google Scholar] [CrossRef]
  87. Zhang, M.; Feng, L.; Li, J.; Chen, L. Therapeutic Potential and Mechanisms of Berberine in Cardiovascular Disease. Curr. Pharmacol. Rep. 2016, 2, 281–292. [Google Scholar] [CrossRef]
  88. Oshiro, C.; Thorn, C.F.; Roden, D.M.; Klein, T.E.; Altman, R.B. KCNH2 pharmacogenomics summary. Pharmacogenetics Genom. 2010, 20, 775–777. [Google Scholar] [CrossRef]
  89. Kang, R.; Tang, D. Chapter 14—Ferroptosis, free radicals, and cancer. In Cancer, 2nd ed.; Preedy, V.R., Patel, V.B., Eds.; Academic Press: San Diego, CA, USA, 2021; pp. 149–158. [Google Scholar]
  90. Kim, S.; Chen, J.; Cheng, T.; Gindulyte, A.; He, J.; He, S.; Li, Q.; Shoemaker, B.A.; Thiessen, P.A.; Yu, B.; et al. PubChem 2019 update: Improved access to chemical data. Nucleic Acids Res. 2019, 47, D1102–D1109. [Google Scholar] [CrossRef]
  91. Daina, A.; Michielin, O.; Zoete, V. SwissADME: A free web tool to evaluate pharmacokinetics, drug-likeness and medicinal chemistry friendliness of small molecules. Sci. Rep. 2017, 7, 42717. [Google Scholar] [CrossRef] [PubMed]
  92. Sander, T.; Freyss, J.; von Korff, M.; Rufener, C. DataWarrior: An open-source program for chemistry aware data visualization and analysis. J. Chem. Inf. Model. 2015, 55, 460–473. [Google Scholar] [CrossRef] [PubMed]
  93. Gu, Z.; Eils, R.; Schlesner, M. Complex heatmaps reveal patterns and correlations in multidimensional genomic data. Bioinformatics 2016, 32, 2847–2849. [Google Scholar] [CrossRef] [PubMed]
  94. Gu, J.; Gui, Y.; Chen, L.; Yuan, G.; Lu, H.Z.; Xu, X. Use of natural products as chemical library for drug discovery and network pharmacology. PLoS ONE 2013, 8, e62839. [Google Scholar] [CrossRef] [PubMed]
  95. O’Hagan, S.; Kell, D.B. Structural Similarities between Some Common Fluorophores Used in Biology, Marketed Drugs, Endogenous Metabolites, and Natural Products. Mar. Drugs 2020, 18, 582. [Google Scholar] [CrossRef] [PubMed]
  96. Oselusi, S.O.; Egieyeh, S.A.; Christoffels, A. Cheminformatic Profiling and Hit Prioritization of Natural Products with Activities against Methicillin-Resistant Staphylococcus aureus (MRSA). Molecules 2021, 26, 3674. [Google Scholar] [CrossRef] [PubMed]
  97. Wenderski, T.A.; Stratton, C.F.; Bauer, R.A.; Kopp, F.; Tan, D.S. Principal component analysis as a tool for library design: A case study investigating natural products, brand-name drugs, natural product-like libraries, and drug-like libraries. Methods Mol. Biol. 2015, 1263, 225–242. [Google Scholar] [CrossRef]
  98. Fernández-de Gortari, E.; García-Jacas, C.R.; Martinez-Mayorga, K.; Medina-Franco, J.L. Database fingerprint (DFP): An approach to represent molecular databases. J. Cheminform. 2017, 9, 9. [Google Scholar] [CrossRef]
  99. Capecchi, A.; Probst, D.; Reymond, J.-L. One molecular fingerprint to rule them all: Drugs, biomolecules, and the metabolome. J. Cheminform. 2020, 12, 43. [Google Scholar] [CrossRef]
  100. Shannon, P.; Markiel, A.; Ozier, O.; Baliga, N.S.; Wang, J.T.; Ramage, D.; Amin, N.; Schwikowski, B.; Ideker, T. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 2003, 13, 2498–2504. [Google Scholar] [CrossRef]
  101. NOM-062-ZOO-1999; Especificaciones técnicas para la producción, cuidado y uso de los animales de laboratorio. Official Mexican Standards: Mexico City, Mexico, 1999. Available online: https://www.gob.mx/cms/uploads/attachment/file/203498/NOM-062-ZOO-1999_220801.pdf (accessed on 20 June 2023).
  102. NOM-087-SEMARNAT-SSA1-2002; Protección Ambiental-Salud Ambiental Residuos Peligrosos Biológico-Infecciosos-Clasificación y Especificaciones de Manejo. Official Mexican Standards: Mexico City, Mexico, 2002. Available online: https://www.cndh.org.mx/sites/default/files/doc/Programas/VIH/Leyes%20y%20normas%20y%20reglamentos/Norma%20Oficial%20Mexicana/NOM-087-SEMARNAT-SSA1-2002%20Proteccion%20ambiental-salud.pdf (accessed on 11 June 2023).
  103. Council, N.R. Committee for the Update of the Guide for the Care and Use of Laboratory Animals. In Guide for the Care and Use of Laboratory Animals, 8th ed.; National Academies Press (US): Washington, DC, USA, 2011. [Google Scholar]
  104. Kochi, T.; Imai, Y.; Takeda, A.; Watanabe, Y.; Mori, S.; Tachi, M.; Kodama, T. Characterization of the arterial anatomy of the murine hindlimb: Functional role in the design and understanding of ischemia models. PLoS ONE 2013, 8, e84047. [Google Scholar] [CrossRef] [PubMed]
  105. Lalka, S.G.; Unthank, J.L.; Nixon, J.C. Elevated cutaneous leukocyte concentration in a rodent model of acute venous hypertension. J. Surg. Res. 1998, 74, 59–63. [Google Scholar] [CrossRef] [PubMed]
  106. Hahn, T.L.; Whitfield, R.; Salter, J.; Granger, D.N.; Unthank, J.L.; Lalka, S.G. Evaluation of the role of intercellular adhesion molecule 1 in a rodent model of chronic venous hypertension. J. Surg. Res. 2000, 88, 150–154. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Hierarchical structural clustering of the previously reported reference compounds. Heatmap of hierarchical clustering was generated using the complexheatmap package in R-Studio separated by K-means. The histogram at the top represents the frequency distribution of each molecular descriptor (blue). Density plots at the top represent the distribution of each molecular descriptor (purple). Each compound is shown as a violin plot on the right, displaying its molecular descriptors’ distribution. On the right side, the boxplots show how the numerical data from molecular descriptors is spread out among reference compounds, the comparison circles display whether or not the mean values are significantly different from each other. As seen in the dendrogram at the left of the figure, there are two main clusters; we focus on the first framed in red lines because of its similarity coefficient.
Figure 1. Hierarchical structural clustering of the previously reported reference compounds. Heatmap of hierarchical clustering was generated using the complexheatmap package in R-Studio separated by K-means. The histogram at the top represents the frequency distribution of each molecular descriptor (blue). Density plots at the top represent the distribution of each molecular descriptor (purple). Each compound is shown as a violin plot on the right, displaying its molecular descriptors’ distribution. On the right side, the boxplots show how the numerical data from molecular descriptors is spread out among reference compounds, the comparison circles display whether or not the mean values are significantly different from each other. As seen in the dendrogram at the left of the figure, there are two main clusters; we focus on the first framed in red lines because of its similarity coefficient.
Molecules 28 07937 g001
Figure 2. Distribution of molecular descriptors for reference compounds, NP molecules and SMS molecules: molecular weight (Da), Ghose–Crippen–Viswanadhan octanol–water partition coefficient (AlogP), the number of H-acceptors (HBA), the number of H-donors (HBD), the number of rotatable bonds, total polar surface area (TPSA), aromatic bound count, and quantitative estimation of drug-likeness (QED). Vertical lines represent the frequency values, while horizontal lines represent the intervals of the data. The x-axis shows the calculated values for each compound and descriptor, while the y-axis shows the frequency of those values. The number of molecules is presented as NPs and SMSs (n = 41,655) and reference compounds (n = 145).
Figure 2. Distribution of molecular descriptors for reference compounds, NP molecules and SMS molecules: molecular weight (Da), Ghose–Crippen–Viswanadhan octanol–water partition coefficient (AlogP), the number of H-acceptors (HBA), the number of H-donors (HBD), the number of rotatable bonds, total polar surface area (TPSA), aromatic bound count, and quantitative estimation of drug-likeness (QED). Vertical lines represent the frequency values, while horizontal lines represent the intervals of the data. The x-axis shows the calculated values for each compound and descriptor, while the y-axis shows the frequency of those values. The number of molecules is presented as NPs and SMSs (n = 41,655) and reference compounds (n = 145).
Molecules 28 07937 g002
Figure 3. Boxplots showing the distribution and summary statistics of the molecular descriptors of reference compounds, NP molecules and SMS molecules: molecular weight range (Da), Ghose–Crippen–Viswanadhan octanol–water partition co-efficient (AlogP), the number of H-acceptors (HBA), the number of H-donors (HBD), the number of rotatable bonds, total polar surface area (TPSA), aromatic bound count of phlebotonic compounds and NPs. The horizontal axis displays calculated values for compounds and molecular descriptors, and the vertical axis shows their frequency. Dots represent outliers. The number of molecules is presented as NPs and SMSs (n = 41,655) and reference compounds (n = 145).
Figure 3. Boxplots showing the distribution and summary statistics of the molecular descriptors of reference compounds, NP molecules and SMS molecules: molecular weight range (Da), Ghose–Crippen–Viswanadhan octanol–water partition co-efficient (AlogP), the number of H-acceptors (HBA), the number of H-donors (HBD), the number of rotatable bonds, total polar surface area (TPSA), aromatic bound count of phlebotonic compounds and NPs. The horizontal axis displays calculated values for compounds and molecular descriptors, and the vertical axis shows their frequency. Dots represent outliers. The number of molecules is presented as NPs and SMSs (n = 41,655) and reference compounds (n = 145).
Molecules 28 07937 g003
Figure 4. The distribution in chemical space according to the PCA of NPs, SMSs, and reference compounds. The green color signifies molecules with a similar chemical distribution to the reference compounds. NP and SMS molecules (n = 41,655) and reference compounds (n = 145).
Figure 4. The distribution in chemical space according to the PCA of NPs, SMSs, and reference compounds. The green color signifies molecules with a similar chemical distribution to the reference compounds. NP and SMS molecules (n = 41,655) and reference compounds (n = 145).
Molecules 28 07937 g004
Figure 5. Steps of the selection of compounds with potential usefulness in VD. Comparison of 1846 drug-like NPs and SMSs with the reference compound diosmin. (A) The dashed line signs the suitable number of clusters determined using the Elbow method. (B) Clusterization of Ward’s method; the reference cluster is shaped in red, and the reference molecule is identified as 1. (C) Cluster plot using K-means shows the same molecules in the reference cluster marked in blue. (D) Silhouette cluster representation to corroborate the last clusterization methods. * means the matrix of distances used in clustering using Ward’s method in subfigure B.
Figure 5. Steps of the selection of compounds with potential usefulness in VD. Comparison of 1846 drug-like NPs and SMSs with the reference compound diosmin. (A) The dashed line signs the suitable number of clusters determined using the Elbow method. (B) Clusterization of Ward’s method; the reference cluster is shaped in red, and the reference molecule is identified as 1. (C) Cluster plot using K-means shows the same molecules in the reference cluster marked in blue. (D) Silhouette cluster representation to corroborate the last clusterization methods. * means the matrix of distances used in clustering using Ward’s method in subfigure B.
Molecules 28 07937 g005
Figure 6. Network of compounds with potential usefulness in VD and their biological targets. The compound–target network was generated with the organic layout of Cytoscape software, 3.8.2 version [20]. The network was constructed with VD pathways using colored nodes for compounds or targets and edges for their relationships.
Figure 6. Network of compounds with potential usefulness in VD and their biological targets. The compound–target network was generated with the organic layout of Cytoscape software, 3.8.2 version [20]. The network was constructed with VD pathways using colored nodes for compounds or targets and edges for their relationships.
Molecules 28 07937 g006
Figure 7. Effects of venous ligation on venous pressure (VP) of hind limbs before the administration of berberine (A) or β-escin (B). The abscissa corresponds to the limb where the registration is carried out and is the right hind limb (RHL) or sham control group and left hind limbs before (LHLB) and after (LHLA) venous ligation. In both treatments, the ordinates display VP variations highlighting the significant difference between control (RHL) and LHLB versus the group with venous hypertension (LHLA); the bars represent the mean ± SEM (standard error of the mean) of six experiments * (p < 0.05 ANOVA post hoc Dunnett’s test). ns = no significant differences.
Figure 7. Effects of venous ligation on venous pressure (VP) of hind limbs before the administration of berberine (A) or β-escin (B). The abscissa corresponds to the limb where the registration is carried out and is the right hind limb (RHL) or sham control group and left hind limbs before (LHLB) and after (LHLA) venous ligation. In both treatments, the ordinates display VP variations highlighting the significant difference between control (RHL) and LHLB versus the group with venous hypertension (LHLA); the bars represent the mean ± SEM (standard error of the mean) of six experiments * (p < 0.05 ANOVA post hoc Dunnett’s test). ns = no significant differences.
Molecules 28 07937 g007
Figure 8. Effect of increased berberine (A) and β-escin (B) doses on the MAP and HR of rats under anesthesia. Abscissae correspond to increased doses of the treatment; ordinates show changes in MAP or HR. The symbols represent the mean ± SEM of eight experiments using berberine and six using β-escin.
Figure 8. Effect of increased berberine (A) and β-escin (B) doses on the MAP and HR of rats under anesthesia. Abscissae correspond to increased doses of the treatment; ordinates show changes in MAP or HR. The symbols represent the mean ± SEM of eight experiments using berberine and six using β-escin.
Molecules 28 07937 g008
Figure 9. Effect of increased berberine (A) and β-escin (B) doses on the venous pressure (VP) of rats under anesthesia. Abscissae correspond to increased doses of the treatments. * indicates a significant difference between the control (RHL) and the values of VP of the limb with venous hypertension (LHLA) (p < 0.05 Student’s test for paired samples). The symbols represent the mean ± SEM (standard error of the mean) of eight experiments using berberine and six using β-escin.
Figure 9. Effect of increased berberine (A) and β-escin (B) doses on the venous pressure (VP) of rats under anesthesia. Abscissae correspond to increased doses of the treatments. * indicates a significant difference between the control (RHL) and the values of VP of the limb with venous hypertension (LHLA) (p < 0.05 Student’s test for paired samples). The symbols represent the mean ± SEM (standard error of the mean) of eight experiments using berberine and six using β-escin.
Molecules 28 07937 g009
Scheme 1. The workflow of the method for identifying compounds with potential usefulness in VD. (A) Phase 1: Creating a phlebotonic reference dataset from MEDLINE, EMBASE, and Scopus databases from the information obtained up to March 2023. The complexheatmap (K-means) in R-Studio was used to cluster the reference compounds based on their molecular descriptors. (B) Phase 2: Selecting drug-like NPs and SMSs from multiple datasets based on the Quantitative Estimate of Drug-likeness (QED). The molecules were searched on the PubChem server to obtain their structural information. Data filtering of NPs and SMSs to obtain drug-like molecules by QED was searched for the reported targets in PubChem. (C) Phase 3: Using the Tanimoto coefficient and Ward’s clustering method to identify compounds with potential usefulness in VD by comparing fingerprints of drug-like NPs, SMSs, and reference compounds. (D) Phase 4: Using Cytoscape software version 3.8 to identify the best candidates with potential usefulness in VD through network analysis. (E) Phase 5: Experimental evaluation of NPs and SMSs in an in vivo model of acute venous hypertension.
Scheme 1. The workflow of the method for identifying compounds with potential usefulness in VD. (A) Phase 1: Creating a phlebotonic reference dataset from MEDLINE, EMBASE, and Scopus databases from the information obtained up to March 2023. The complexheatmap (K-means) in R-Studio was used to cluster the reference compounds based on their molecular descriptors. (B) Phase 2: Selecting drug-like NPs and SMSs from multiple datasets based on the Quantitative Estimate of Drug-likeness (QED). The molecules were searched on the PubChem server to obtain their structural information. Data filtering of NPs and SMSs to obtain drug-like molecules by QED was searched for the reported targets in PubChem. (C) Phase 3: Using the Tanimoto coefficient and Ward’s clustering method to identify compounds with potential usefulness in VD by comparing fingerprints of drug-like NPs, SMSs, and reference compounds. (D) Phase 4: Using Cytoscape software version 3.8 to identify the best candidates with potential usefulness in VD through network analysis. (E) Phase 5: Experimental evaluation of NPs and SMSs in an in vivo model of acute venous hypertension.
Molecules 28 07937 sch001
Table 2. Number of drug-like molecules obtained using fingerprint analysis of each reference compound.
Table 2. Number of drug-like molecules obtained using fingerprint analysis of each reference compound.
Reference CompoundsNo. of Compounds
with Structural Similarity
Silhouette Coefficient (Sc)
Escin4680.5
Crocin4680.5
Echinomycin7640.35
Cyclosporin A4680.5
Amphotericin B4680.5
Everolimus4680.5
Rapamycin 4680.5
Chetomin 764 0.35
Astragaloside IV 418 0.5
20(R)-ginsenoside Rh2 418 0.5
Epigallocatechin 3-gallate 614 0.23
Troxerutin 614 0.23
O-(beta-hydroxyethyl)-rutoside 614 0.23
Salvianolic acid B 614 0.23
Keracyanin chloride 764 0.35
Rutin 614 0.23
Hidrosmin 614 0.23
Diosmin 614 0.23
Hesperidin 614 0.23
Linarin 614 0.23
Isorhoifolin 614 0.23
Naringin 614 0.23
Betanin4680.5
Table 3. Drug-like molecules obtained using fingerprint analysis, their number of pharmacological targets, and their interaction with the VD network.
Table 3. Drug-like molecules obtained using fingerprint analysis, their number of pharmacological targets, and their interaction with the VD network.
Compounds TypeNo. of Compounds
Obtained for Fingerprint Analysis
(Total n = 1846)
No. of Compounds with Targets
(Total n = 229)
No. of Compounds that Interact with the VD Network
NPs129320384
SMSs552260
Table 5. Baseline values of MAP and HR were recorded before (LHLB) and after (LHLA) venous ligation in anesthetized rats for the experiments of berberine and β-escin.
Table 5. Baseline values of MAP and HR were recorded before (LHLB) and after (LHLA) venous ligation in anesthetized rats for the experiments of berberine and β-escin.
LimbsTreatment
Berberineβ-Escin
MAP (mm Hg)HR (Beats/min)MAP (mm Hg)HR (Beats/min)
LHLB120.8 ± 6.0414.8 ± 13.4101.6 ± 14.9361.3 ± 29.3
LHLA121.9 ± 6.47410.0 ± 16.395.4 ± 4.2363.8 ± 32.6
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Barrera-Vázquez, O.S.; Escobar-Ramírez, J.L.; Santiago-Mejía, J.; Carrasco-Ortega, O.F.; Magos-Guerrero, G.A. Discovering Potential Compounds for Venous Disease Treatment through Virtual Screening and Network Pharmacology Approach. Molecules 2023, 28, 7937. https://doi.org/10.3390/molecules28247937

AMA Style

Barrera-Vázquez OS, Escobar-Ramírez JL, Santiago-Mejía J, Carrasco-Ortega OF, Magos-Guerrero GA. Discovering Potential Compounds for Venous Disease Treatment through Virtual Screening and Network Pharmacology Approach. Molecules. 2023; 28(24):7937. https://doi.org/10.3390/molecules28247937

Chicago/Turabian Style

Barrera-Vázquez, Oscar Salvador, Juan Luis Escobar-Ramírez, Jacinto Santiago-Mejía, Omar Francisco Carrasco-Ortega, and Gil Alfonso Magos-Guerrero. 2023. "Discovering Potential Compounds for Venous Disease Treatment through Virtual Screening and Network Pharmacology Approach" Molecules 28, no. 24: 7937. https://doi.org/10.3390/molecules28247937

APA Style

Barrera-Vázquez, O. S., Escobar-Ramírez, J. L., Santiago-Mejía, J., Carrasco-Ortega, O. F., & Magos-Guerrero, G. A. (2023). Discovering Potential Compounds for Venous Disease Treatment through Virtual Screening and Network Pharmacology Approach. Molecules, 28(24), 7937. https://doi.org/10.3390/molecules28247937

Article Metrics

Back to TopTop