Next Article in Journal
Quercetin and Ferulic Acid Elicit Estrogenic Activities In Vivo and In Silico
Previous Article in Journal
Hybrid Polymer-Silica Nanostructured Materials for Environmental Remediation
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 28
Issue 13
10.3390/molecules28135108
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Phytochemical Analysis, In Vitro Biological Activities, and Computer-Aided Analysis of Potentilla nepalensis Hook Compounds as Potential Melanoma Inhibitors Based on Molecular Docking, MD Simulations, and ADMET

by
Subhash Sharma
1,
Vikas Kumar
2,
Muhammad Yaseen
3,*,
Amr S. Abouzied
4,
Abgeena Arshad
5,
Mashooq Ahmad Bhat
6,
Ahmed M. Naglah
6,
Chirag N. Patel
7,8,
Prasanth Kumar Sivakumar
7,
Anuradha Sourirajan
1,
Adnan Shahzad
3 and
Kamal Dev
1,9,*
1
Faculty of Applied Sciences and Biotechnology, Shoolini University, P.O. Box 9, Head Post Office, Solan 173212, India
2
Department of Biotechnology, University Institute of Biotechnology, Chandigarh University, Gharuan, Mohali 140413, India
3
Institute of Chemical Sciences, University of Swat, Charbagh, Swat 19130, Pakistan
4
Department of Pharmaceutical Chemistry, National Organization for Drug Control and Research (NODCAR), Giza 12311, Egypt
5
Asian Medical Institute, Kant 725013, Kyrgyzstan
6
Department of Pharmaceutical Chemistry, College of Pharmacy, King Saud University, Riyadh 11451, Saudi Arabia
7
Department of Botany, Bioinformatics and Climate Change Impacts Management, University School of Science, Gujarat University, Ahmedabad 380009, India
8
Biotechnology Research Center, Technology Innovation Institute, Abu Dhabi 9639, United Arab Emirates
9
Department of Pharmacology and Toxicology, Wright State University, Dayton, OH 4543, USA
*
Authors to whom correspondence should be addressed.
Molecules 2023, 28(13), 5108; https://doi.org/10.3390/molecules28135108
Submission received: 15 May 2023 / Revised: 10 June 2023 / Accepted: 26 June 2023 / Published: 29 June 2023
Browse Figures
Versions Notes

Abstract

:
Potentilla nepalensis Hook is a perennial Himalayan medicinal herb of the Rosaceae family. The present study aimed to evaluate biological activities such as the antioxidant, antibacterial, and anticancer activities of roots and shoots of P. nepalensis and its synergistic antibacterial activity with antibacterial drugs. Folin–Ciocalteau and aluminium chloride methods were used for the calculation of total phenolic (TPC) and flavonoid content (TFC). A DPPH radical scavenging assay and broth dilution method were used for the determination of the antioxidant and antibacterial activity of the root and shoot extracts of P. nepalensis. Cytotoxic activity was determined using a colorimetric MTT assay. Further, phytochemical characterization of the root and shoot extracts was performed using the Gas chromatography–mass spectrophotometry (GC-MS) method. The TPC and TFC were found to be higher in the methanolic root extract of P. nepalensis. The methanolic shoot extract of P. nepalensis showed good antioxidant activity, while then-hexane root extract of P. nepalensis showed strong cytotoxic activity against tested SK-MEL-28 cells. Subsequently, in silico molecular docking studies of the identified bioactive compounds predicted potential anticancer properties. This study can lead to the production of new herbal medicines for various diseases employing P. nepalensis, leading to the creation of new medications.

1. Introduction

Throughout human history, medicinal plants have been employed in traditional medicine and are regarded as a source of healthy human habitation. Different plant sections, such as the roots, leaves, stems, bark, fruits, and seeds, have been used to boost immunity and prevent a number of ailments [1]. Plants produce secondary metabolites (small organic chemicals) that are mostly needed for reproduction and defense mechanisms against bacteria, fungi, viruses, vertebrates, etc., but are not essential for their regular growth or development. These goods have strong potential for use as drugs. The WHO has identified over 21,000 plants that are widely used for medical reasons around the world. About 2500 species have been found in India, and biopharmaceutical companies exploit over 150 of them commercially on a sizable scale as conventional medicine. India, which holds the title of “the botanical garden of the globe”, is the country that produces the most medicinal plants. Interesting options exist in traditional medicine to fight MDR (multidrug resistance). Herbal medications exhibit a diverse range of biological activity, making them effective tools for treating disease. Combining dietary and medicinal approaches could result in a potent method for managing a variety of disorders [2].
The Indian Himalayan region is one of the richest reservoirs of biodiversity in the world. This region is rich in medicinal herbs and plants, which are used by the local people for various medicinal purposes at home as well as being the basic constituents of the medicines and other products available on the market [3]. Potentilla nepalensis is a medicinal herb of the Rosaceae family found in the Indian Himalayan region, especially in the North-Western Himalayas region. The genus Potentilla is well known for its pharmacological activities and bioactive compounds. The existence of several phytochemicals originating from various plant parts of the Potentilla species can be used to explain their pharmacological effects. In many cultures around the world, Potentilla extracts have been used to cure a variety of ailments. Compounds extracted from various portions of Potentilla species plants have been shown to have anti-inflammatory, anti-hyperglycemic, anticancer, and anti-ulcerogenic activities [4]. Potentilla species plants have been found to be high in phenolics, flavonoids, and terpenoids, as well as having strong antioxidant and antibacterial effects. Some Potentilla species have been employed in traditional remedies for their anti-diabetic, anti-diarrheal, anti-viral, anti-inflammatory, wound-healing, and anticancer properties. Root extracts from some Potentilla species have also been used to treat viral infections in Tibetan traditional medicine [1]. The roots of P. nepalensis are traditionally used for thr treatment of headaches, cold, asthma, dysentery, skin diseases, and various other ailments [5].
Melanoma is a type of skin cancer that is highly aggressive and caused by the malignant proliferation of melanocytes [6,7]. The incidence of melanoma is increasing globally at a rate of approximately 3–7% per year, making it a significant public health concern [8]. Timely detection of melanoma usually results in successful surgical treatment [9]. However, advanced or metastatic melanoma does not respond well to current treatments, which include targeted therapy and immunotherapy due to their high cost and tumor resistance development [10,11]. Therefore, there is a need for novel, effective, and easily accessible therapeutic approaches for melanoma that overcome these limitations. Glycogen synthase kinase 3 (GSK3) is a serine/threonine protein kinase found in all cells, regulating multiple cellular processes such as glycogen metabolism, signal transduction, cell cycle regulation, and cell proliferation [12,13,14]. GSK3 plays a crucial role in regulating the oxidative stress response associated with cancer occurrence and progression. Recent studies [15,16] have demonstrated that GSK3 overexpression makes breast cancer cells more sensitive to chemotherapeutic drugs and facilitates elastin-induced ferroptosis. A report from John et al. [17] showed that low constitutive levels of GSK3 signaling control both N-cadherin expression and the formation of focal adhesion complexes, which in turn contribute to the oncogenic nature of melanoma. GSK-3 regulates melanoma proliferation and shape via phosphorylation and elevated PAX3 levels [18].
Molecular docking analysis is a significant tool that predicts molecule activity and affinity by depicting the binding position of molecules to protein targets [19,20]. This study aims to investigate the antioxidant potential, anticancer activities, and total phenolic and flavonoid content of the methanolic and n-hexane extracts of the roots and shoots of P. nepalensis, a medicinal plant from the North-Western Himalayas region, using MTT assay against SKMEL-28 cells. In silico molecular docking, MD simulations, and toxicity prediction were performed to identify potential bioactive compounds with anticancer potential to discover and formulate new drug formulations.

2. Results

2.1. Quantification of Total Phenolic and Flavonoid Content

The data of the TPC and TFC of the roots and shoots extracts of P. nepalensis are shown in Table 1. A higher TPC (21.21 ± 0.54 mg/g GAE) and TFC (4.24 ± 0.17 mg/g RE) was observed in the methanolic extract of P. nepalensis roots as compared to that of P. nepalensis shoots.

2.2. In Vitro Antioxidant Activity

The DPPH scavenging activity of the P. nepalensis roots and shoots was found to be concentration-dependent and expressed in terms of IC50 value. Among all extracts, the methanolic extract of the P. nepalensis shoots showed high antioxidant activity with the smallest IC50 value being 12.83 ± 0.35µg/mL, followed by the methanolic root extract (IC50-23.5 ± 0.92 µg/mL). The standard L-ascorbic acid was found to show strong antioxidant activity (IC50-5.86 ± 0.13µg/mL), while the n-hexane extract of both the roots and shoots of P. nepalensis exhibited low antioxidant activity (Table 1).

2.3. Cytotoxic Activity of Roots and Shoots of P. nepalensis

The results obtained for the MTT assay of the methanolic and n-hexane extracts of the roots and shoots of P. nepalensis are shown in Figure 1. Among all the extracts of roots and shoots, the n-hexane extract of the P. nepalensis roots showed comparatively high cytotoxic activity (IC50-65.37 ± 2.27 µg mL−1), followed by the n-hexane extract of the P. nepalensis shoots (IC50-81.97 ± 2.75 µg mL−1) and methanolic shoot extract (IC50-82.74 ± 1.54 µg mL−1) against tested SK-MEL-28 cells. The methanolic root extract was found to have the least cytotoxic effect (IC50-88.65 ± 3.71 µg mL−1). Doxorubicin was used as positive control showing IC50-4.55 ± 0.18 µg mL−1 against SK-MEL-28 cells. Both root and shoot samples elicited a concentration-dependent reduction in cell viability as indicated by the reduction in cell numbers from the Inverted phase-contrast microscopy images (Figure 2, Figure 3, Figure 4 and Figure 5). In addition to this, alterations in cell morphology, such as cell rounding up and cell fragmentation from the control group of cells, can be observed. All these together support the reduction in cell viability as indicated by the MTT assay.

2.4. Identification of Major Phytocompounds of Methanolic and n-Hexane Extracts of Roots and Shoots of P. nepalensis via GC-MS Profiling

The GC-MS chromatograms of the methanolic and n-hexane extracts of the roots and shoots from P. nepalensis revealed the presence of 10 compounds each (Figure 6A–D). The major phytocompounds reported in the GC-MS chromatograms of the methanolic root extract of P. nepalensis were tetradecanoic acid, 10,13-dimethyl-, methyl ester (12.64%), and Heptadecanoic acid, 16-methyl-, methyl ester (10.58%) (Table 2); while the methanolic shoot extracts of P. nepalensis showed the presence of Hexadecanoic acid, and methyl ester (11.82%),1,1,1,5,7,7,7-Heptamethyl-3,3-bis(trimethylsiloxy)tetrasiloxane (8.58%) (Table 3). GC-MS analysis of the n-hexane root extract of P. nepalensis showed the presence of Trichloromethyl 9-anthracenecarbodithioate (22.90%), Heptane,3,3-dimethyl-(17.71%), and Hexadecanene (12.28%) (Table 4); while the n-hexane shoot extract of P. nepalensis showed the presence of Benzene, 1,3,5-tri-tert-butyl-(23.19%), 1,1,1,3,5,5,5-Heptamethyltrisiloxane (23.19%), and 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester (17.32%) as major phytocompounds (Table 5).

2.5. Molecular Docking Analysis

The binding energies of selected phytocompounds of P. nepalensis with target proteins using the Glide (grid-based ligand docking) program are summarized in Table 6. Trichloromethyl 9-anthracenecarbodithioate (−8.9 kcal/mol) and 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester (−7.4 kcal/mol) have shown the highest binding affinity among all selected phytocompounds. However, encorafenib (drug) was found to a show binding affinity of -8.6 kcal/mol with the selected target protein. Binding interactions of selected phytocompounds of P. nepalensis with interacting amino acids were analyzed using the Discovery Studio (DS) visualizer. Figure 7A–D showed the binding interactions of the best docked phytocompounds, viz., trichloromethyl 9-anthracenecarbodithioate and 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester with amino acids of the selected target proteins (PDB ID: 5K5N).

2.6. MD Simulations Study

MD simulations were conducted for 100 ns using the Academic version of the Desmond programme version 2.0 (Schrödinger LLC, New York, NY, USA) to evaluate the overall stability and flexibility of the ligand–protein complexes. Two phytochemicals that exhibited strong binding interactions with the 5K5N receptor protein were selected for the MD simulations. RMSD (root-mean-square deviation) was used to assess the fluctuations of the ligand within the active site of the receptor protein based on the MD trajectories. The RMSD values over time for the Cα atoms of the protein–ligand complex are presented in Figure 8. Both complexes remained stable during the MD simulations, as indicated by the protein backbone RMSD values hovering around 4.0 Å and well under 8 Å for the trichloromethyl 9-anthracenecarbodithioate (Figure 8A), and 4H-1-benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester-5K5N (Figure 8B) complexes, respectively.
The root-mean-square fluctuation (RMSF) is a useful tool for identifying changes in protein structure, as shown in Figure 9A,B. The amino acids of the 5K5N proteins with trichloromethyl 9-anthracenecarbodithioate displayed fluctuations ranging between 0.5–5.0 Å (Figure 9A), while the fluctuations of the amino acids of the 5K5N proteins with 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester were represented in Figure 9B. In the RMSF plot, peaks correspond to the portions of the receptor proteins that exhibited maximum fluctuations during the MD simulation. Generally, protein tails, identified through the N and C terminals, display higher fluctuations than rigid structures such as α-helices and β-strands. These structured regions are stiffer than unstructured regions, hence showing minimal fluctuations. The highest RMSF values obtained during the MD simulation for both protein–ligand complexes were 4.6 Å in Glu249 (Figure 5A,B), likely due to its high flexibility resulting from the formation of a β-turn in the protein structure at that position [37].
Protein–ligand interactions offer valuable insights into simulation techniques, conformational stability, and correlated effects. These interactions can be classified into four main subtypes—Hydrogen Bonds, Hydrophobic, Ionic, and Water Bridges—each with explicit subtypes that can be studied using the ‘Simulation Interactions Diagram’ board. The stacked bar outlines used to represent the interactions are standardized, and an estimation of 0.8 indicates that during 80% of the simulation time the interaction is sustained. However, values over 1.0 are possible as some protein build-ups may create multiple contacts of the same subtype with the ligand. In the complex consisting of trichloromethyl 9-anthracenecarbodithioate-5K5N, hydrophobic interactions with the ligand were found to occur with Ile 62, Phe 67, Val 70, Ala 83, Leu 132, Tyr 134, and Leu 188 (Figure 10A). Amino acids Val 135 and Arg 141 were the most important for hydrogen bonding in the 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester-5K5N complex, while amino acids Ala 83, leu 132, and Leu 188 were the most significant for hydrophobic interactions (Figure 10B).

2.7. Assessment of Drug Likeness and Toxicity Prediction

Table 7 represents the drug likeness and toxicity prediction of the best docked phytocompounds of P. nepalensis. Lipinski’s rule of five and ADMET prediction were used to evaluate the pharmacokinetic and toxicity properties of the top-ranked compounds. The druggability and toxicity parameters were achieved for both top-ranked phytocompounds, except for carcinogenicity. Trichloromethyl 9-anthracenecarbodithioate was predicted to be carcinogenic, with an LD50 of 493 mg/kg (Class IV). In contrast, 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester had an LD50 of 100 mg/kg (Class III) (Table 7).

3. Discussion

Plant phenolics and flavonoids are currently a major study focus since they are assumed to be responsible for the bulk of biological activities of plants, such as their anti-inflammatory, antibacterial, antiviral, antioxidant, and anticancer capabilities. Medicinal plants have played an important role in the discovery of approximately 50% of anticancer drugs. In our study, we have observed high phenolic and flavonoid content in the methanolic extract of roots and shoots of P. nepalensis, whereas n-hexane extracts of roots and shoots of P. nepalensis showed lower amounts of phenolic and flavonoid content, which may be due to the higher solubility of phenolics and flavonoids in methanol extract as compared to that of n-hexane extract [38]. Tomczyk et al. [39] reported a higher TPC (73.9 ± 3.7 mg GAE/g dw) in the aqueous extract of areal parts of P. nepalensis, while the TFC in the aqueous extract (2.1 ± 0.5 mg QE/g dw) was found to be comparable with our results [39]. This variation in TPC content can be related to differences in geographical regions, solvent preference, plant age, and plant part selection [40,41,42]. Recently, a study by Sharma et al. [43] reported a higher phenolic (21.21 ± 0.54 mg g−1 GAE) and flavonoid content in wild grown plants (4.24 ± 0.17 mg g−1 RE) as compared to that of in vitro propagated plants (TPC-3.55 ± 0.72 mg g−1 GAE; TFC-0.33 ± 0.07 mg g−1 RE).
Free radicals, which can harm cell membranes and other structures, are unstable molecules produced by the oxidation process in the human body. These free radicals have been linked to several diseases, including heart disease and various cancers. Antioxidants, which are chemicals that scavenge and battle free radicals and are present in medicinal plants, may reduce the risk of a variety of diseases, including heart disease and several types of cancer. We examined the in vitro antioxidant potential of methanolic and n-hexane extracts of P. nepalensis roots and shoots while taking the positive effects of antioxidants into consideration. The antioxidant activity of methanolic root extract of P. nepalensis was also reported by Sharma et al. [43]. Low IC50 values of methanolic extracts indicating strong free radical scavenging or antioxidant capabilities were observed in our study. The strong antioxidant potential in methanolic root and shoot extracts as compared to n-hexane extracts can be attributed to the higher amount of phenolic and flavonoid content [44].
It is well recognized that medicinal plants contain a variety of chemicals with the ability to prevent or treat a wide range of illnesses. As cancer is a serious issue for public health, we preferred using an MTT assay for the methanolic and n-hexane extracts of P. nepalensis roots and shoots against SKMEL-28 cells. The n-hexane extract of the P. nepalensis roots was found to be most effective against SKMEL-28 cells with the smallest IC50 value. Additionally, there is very little analytical data on the chemical makeup of P. nepalensis, and we used GC-MS profiling to analyze the methanolic and n-hexane extracts of the plant’s roots and shoots to determine their primary phytocompounds. The majority of phytocompounds identified via GC-MS profiling of the methanolic and n-hexane extracts of P. nepalensis roots and shoots have been reported for biological activities such as antimicrobial, antioxidant, and anticancer activities and have also been discovered in the plant extracts of some other therapeutic plants. The presence of these phytocompounds may be responsible for the medicinal properties of P. nepalensis.
Computational approaches are extremely valuable in pharmaceutical research since they aid in the discovery and development of new, promising medicines, particularly when used in conjunction with molecular docking techniques. These techniques have been used by several research teams to screen possible new chemicals against a variety of ailments. In silico pharmacokinetic, pharmacological, and toxicological performance have also been predicted using them. In the present study, we have utilized a molecular docking system to analyze the binding affinities of selected phytocompounds of P. nepalensis with amino acids of glycogen synthase kinase 3β (GSK3β) protein to check the overall stability and/or flexibility of the ligand–protein complexes as well as the drug likeness and toxicity prediction of the best docked phytocompounds of P. nepalensis. To our knowledge, this is the first report on the anticancer activity and GC-MS characterization of the shoots and roots of P. nepalensis. Molecular docking and MD simulation studies also support the anticancer activity of P. nepalensis. Further, analysis of the drug-like effects of selected compounds in this study displayed the potential of these compounds as anticancer drugs. The selected potential drug candidates were also found to fulfill the toxicity criteria. Overall, the selected phytochemicals are good potential drug candidates to treat myeloma in the future.

4. Materials and Methods

4.1. Chemicals and Reagents

The aluminum chloride, ascorbic acid, dimethyl sulfoxide (DMSO), 2,2- diphenyl-2-picrylhydrazyl (DPPH), ferric chloride, Folin–Ciocalteau reagent, gallic acid, n-hexane, methanol, rutin, and sodium nitrite (NaNO2) were purchased from Loba Chemie Pvt. Ltd., Mumbai, India. 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2H-tetrazolium bromide (MTT), Dulbecco’s modified eagle medium (DMEM), nutrient broth, nutrient agar, ampicillin, chloramphenicol, and resazurin were purchased from Himedia Laboratories Pvt. Limited, Mumbai, India. All the chemicals and reagents utilized in this study were of analytical grade.

4.2. Collection and Identification of Plant Material

Potentilla nepalensis plants collected during July to September 2020 from Kurri, Shimla, Himachal Pradesh, at a height of 2600 m above sea level were used as experimental material. The sample was validated by the Department of Forest Products, Dr. Y.S. Parmar University of Horticulture and Forestry, Nauni, Solan, Himachal Pradesh under voucher no 13925.

4.3. Preparation of Extracts of Roots and Shoots of P. nepalensis

The roots and all the above ground parts were harvested from P. nepalensis plants and cut into pieces. The cut pieces were washed with running tap water followed by distilled water and dried at 40 °C and converted to fine powder using an electric grinder. About 10 g of fine powder was added to 100 mL of 70% methanol and n-hexane in separate culture tubes and macerated for 24 h at 30 °C for extract preparation [31,32,33]. The filtrate was collected and evaporated at 40 °C. The dried crude extracts were stored in airtight bottles at 4 °C for further use.

4.4. Spectrophotometric Quantification of Total Phenolic (TPC) and Flavonoid Content (TFC)

The TPC and TFC of the methanolic and n-hexane extracts of roots and shoots of P. nepalensis was performed using the Folin–Ciocalteau (FC) reagent method [45,46,47] and aluminum chloride (AlCl3) method [43,48]. Standard curve of Gallic acid (5–80 µg/mL) and Rutin (5–80 µg/mL) was used for calculation of the TPC and TFC, respectively.

4.5. In Vitro Antioxidant Potential Using DPPH Radical Scavenging Assay

The DPPH radical scavenging method was used to assess the in vitro antioxidant capacity of the roots and shoots of P. nepalensis [49]. To 100 µL of varied final concentrations of extracts (5–20 µg/mL), 900 µL of 0.004% DPPH solution (w/v, in methanol) was added and allowed to react in the dark for 15 min. In the control reaction, 100 µL of methanol was added to 900 µL of 0.004% DPPH solution. After incubation, absorbance was measured at 517 nm using a UV-visible spectrophotometer. Ascorbic acid (5–20 µg/mL) was used as standard. The following equation was used to calculate the DPPH scavenging activity:
DPPH   radical   scavenging = O . D .   C o n t r o l − O . D .   S t a n d a r d / T e s t O . D .   C o n t r o l × 100
where O.D. (control) is the absorbance of the control, and O.D. (standard/test) is the absorbance of the extract/standard. The IC50 value was calculated to express the free radical scavenging activity of extracts, showing the effective concentration of extract/standard used to scavenge 50% of DPPH radicals. The lower the IC50 value, the more the scavenging ability of the extract.

4.6. In Vitro Cytotoxic Activity of P. nepalensis Roots and Shoots

The cytotoxic properties of the extracts of the roots and shoots of P. nepalensis against myeloma (SK-MEL-28) cells were determined using the MTT assay. To begin with, the cells were obtained from culture flasks through trypsinization and transferred into 96-well microculture plates containing 100 μL aliquots of DMEM supplemented with 10% heat-inactivated fetal bovine serum, 1 mM sodium pyruvate, 2 mM L-glutamine, and 1% non-essential amino acids (100×). The cell density was adjusted to 2.5 × 103 cells per well to ensure their exponential growth. The cells were allowed to settle and grow in a drug-free complete culture medium for 24 h, after which they were exposed to various dilutions of the test samples (6.25, 12.5, 25, 50, and 100 µg/mL) dissolved in DMSO for 48 h. The control wells were left untreated. Doxorubicin was used as a positive control. Next, the medium was replaced with an RPMI 1640 medium along with MTT solution in phosphate-buffered saline (PBS), and the mixture was incubated for another 2 h. The viable cells produced formazan product that was dissolved in DMSO (150 μL/well), and the optical densities were measured at 570 nm with a microplate reader, using the reference wavelength of 690 nm to correct for non-specific absorption. The percentage of viable cells was determined, and the IC50 values were calculated [50]. The IC50 values were calculated using the equation for slope (y = mx + C) obtained by plotting the absorbance of the different concentrations of the test/drug sample (6.25–100 µg/mL) in Microsoft Excel.

4.7. GC-MS Profiling of Root and Shoot Extracts of P. nepalensis to Identify Major Phytocompounds

The GC-MS analysis of the roots and shoots of P. nepalensis was conducted using a Thermo Fisher Scientific Gas Chromatograph equipped with a Tri Plus RSH Autosampler, GC trace-1300, and MS-TSQ Duo. The Thermo Fisher Scientific TG-5MS column was utilized, which measured 40 m in length, 0.15 mm in film, and 0.15 m in internal diameter. The method involved setting the first oven temperature to 80 °C, with a temperature increase of 8 °C/min and a 1-min hold period, followed by increasing the temperature to 150 °C, with a rate of 10 °C/min and a 6-min hold period. The total run time was 32 min, with a 1 µL sample volume injected using helium at a flow rate of 0.7 mL/min as the carrier gas. The MS was operated within the electron ionization (EI) mode, scanning within a 40–450 amu range with a mass spectrometer source temperature and transfer line temperature set at 230 °C and 250 °C, respectively, and an electron multiplier voltage of 1 kV. Mass spectra were interpreted using the NIST/EPA/NIH Mass Spectral Library Version 2.2, 2014, and fragmentation patterns were compared with the instrument database data for all constituents detected.

4.8. Ligand Preparation

The three-dimensional structures of identified phytochemicals and Encorafenib (drug) were obtained from Pubchem (https://pubchem.ncbi.nlm.nih.gov/, accessed on 1 March 2023) (Figure S1). Extracts containing phytocompounds above 10% were selected as ligands for docking studies. The ligands and drug were minimized using Chem3D structure software and converted to pdbqt format using Open Babel (http://openbabel.org/, accessed on 1 March 2023).

4.9. Retrieval and Preparation of Target Proteins

The three-dimensional structure of GSK-3 beta complexed with PF-04802367 (PDB ID: 5K5N) [51] (Figure S2) was downloaded from the protein data bank (https://www.rcsb.org/, accessed on 5 March 2023) and used Auto-Dock tool 1.5.6 for the preparation of the protein structure. The active site of the target protein was predicted based on the previously bounded native ligand [38], and the grid box parameters for 5K5N were set to a size of 37.05 Å × 23.225 Å × 23.35 Å (x, y, and z) with a center at coordinates 4.056 Å × 2.085 Å × 29.143 Å (x, y, and z).

4.10. Molecular Docking

A total of seven major phytocompounds of P. nepalensis were docked against GSK-3 beta protein using the AutoDock vina tool [52]. The top scoring phytocompounds were selected based on their binding energy with target protein receptors. The best pose based on binding energies for each ligand–protein interaction was further analyzed using the Discovery Studio (DS) visualizer (Accelrys, San Diego, CA, USA).

4.11. Molecular Dynamics Simulations

Desmond program version 2.0 (academic version) was used to investigate the stability of the receptor–ligand complexes formed via molecular docking [53,54]. The system was prepared using the TIP3P water model with a cubic periodic box incorporating Simple Point Charge (SPC) (10 Å × 10 Å × 10 Å) with Optimized Potentials for the Liquid Simulations (OPLS) all-atom force field 2005 [55]. Sodium ions were added to neutralize the system. The receptor–ligand complexes were subjected to energy minimization and pre-equilibration through several restricted steps. MD simulations were conducted with periodic boundary conditions using the OPLS 2005 force field parameters with a relaxation time of 1 ps at a constant temperature of 300 K and constant volume in the NPT ensemble system [56,57]. The Smooth Particle Mesh Ewald (PME) approach was used to analyze protein structures every 1 ns. The stability was calculated from an average structure obtained from the production phase of the MD simulation. The authors examined the structural changes of the receptor–ligand complexes using the histogram for torsional bonds, the radius of gyration (Rg), the root-mean-square deviation (RMSD), and the root-mean-square fluctuation (RMSF) to understand the dynamic role played by the complexes [58,59,60,61].

4.12. Evaluation of Drug-Likeness and ADME/Toxicity Properties

Lipinski’s rule (rule of five, RO5) was considered the primary factor for screening of the molecules, and it was evaluated using the SWISS ADME web server (http://www.swissadme.ch/, accessed on 10 March 2023). Further, the toxicity of selected compounds was analyzed using the Protox-II tool to ascertain their risk of druggability [62,63].

5. Conclusions

The present study provides evidence that P. nepalensis shoots and roots have antioxidant and anticancer properties. Various phytocompounds with known biological activity were detected via GC-MS profiling of P. nepalensis, further supporting the therapeutic significance of this Himalayan plant. In addition, molecular docking, MD simulation studies, and druggability analyses were performed on GSK-3 beta protein with the selected phytochemicals, trichloromethyl 9-anthracenecarbodithioate and 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo, ethyl ester. Throughout the duration of the MD simulation, both phytocompounds stayed firmly connected to the target proteins, i.e., the GSK-3 beta protein. The drug-likeness and ADMET analyses of both of these drug candidates for the treatment of melanomas revealed their drug-like properties. Both phytochemicals were found to be non-toxic; therefore, their anticancer activities to treat or suppress melanoma could be the subject of future experimental studies.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/molecules28135108/s1, Figure S1: Chemical structures of identified compounds used for molecular docking.; Figure S2: Three-dimensional structure of target protein, Glycogen Synthase Kinase-3(PDB ID: 5K5N).

Author Contributions

Conceptualization, S.S. and V.K., M.A.B.; methodology, S.S., V.K., C.N.P., P.K.S., A.A., A.S. (Adnan Shahzad); software, C.N.P. and P.K.S.; validation, C.N.P., P.K.S., A.A., A.S. (Anuradha Sourirajan); investigation, S.S., V.K., M.A.B.; resources, A.S.A., A.S. (Adnan Shahzad), K.D., A.M.N.; writing—original draft preparation, S.S., V.K., A.S.A.; writing—review and editing, A.S. (Anuradha Sourirajan), K.D., A.M.N., M.Y.; visualization, V.K., C.N.P.; Funding execution, M.A.B., A.M.N.; Supervision, K.D., M.Y. All authors have read and agreed to the published version of the manuscript.

Funding

Research supporting project number (RSPD2023R740), King Saud University, Riyadh, Saudi Arabia.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

We are thankful to the Research supporting project number (RSPD2023R740), King Saud University, Riyadh, Saudi Arabia. The authors are also thankful to Shoolini University for providing infra-structural support.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Kumari, S.; Seth, A.; Sharma, S.; Attri, C. A holistic overview of different species of Potentilla a medicinally important plant along with their pharmaceutical significance: A review. J. Herb Med. 2021, 29, 100460. [Google Scholar] [CrossRef]
  2. Anand, U.; Jacobo-Herrera, N.; Altemimi, A.; Lakhssassi, N. A comprehensive review on medicinal plants as antimicrobial therapeutics: Potential avenues of biocompatible drug discovery. Metabolites 2019, 9, 258. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Mehta, P.; Bisht, K.; Sekar, K.C. Diversity of threatened medicinal plants of Indian Himalayan Region. Plant Biosyst. 2021, 155, 1121–1132. [Google Scholar] [CrossRef]
  4. Tomczyk, M.; Latté, K.P. Potentilla—A review of its phytochemical and pharmacological profile. J. Ethnopharmacol. 2009, 122, 184–204. [Google Scholar] [CrossRef] [PubMed]
  5. Jain, M.; Jain, V.K.; Sharma, S. Rosaceae of Solan District of Himachal Pradesh. IJSRST 2018, 4, 1580–1589. [Google Scholar]
  6. Eddy, K.; Chen, S. Overcoming immune evasion in melanoma. Int. J. Mol. Sci. 2020, 21, 8984. [Google Scholar] [CrossRef]
  7. Feng, S.; Zhou, Y.; Huang, H.; Lin, Y.; Zeng, Y.; Han, S.; Huang, K.; Liu, Q.; Zhu, W.; Yuan, Z.; et al. Nobiletin Induces Ferroptosis in Human Skin Melanoma Cells Through the GSK3β-Mediated Keap1/Nrf2/HO-1 Signalling Pathway. Front. Gen. 2022, 13, 865073. [Google Scholar] [CrossRef]
  8. Leiter, U.; Keim, U.; Garbe, C. Epidemiology of skin cancer: Update 2019. In Sunlight, Vitamin D and Skin Cancer; Springer: Cham, Switzerland, 2020; pp. 123–139. [Google Scholar]
  9. Scala, J.; Vojvodic, A.; Vojvodic, P.; Vlaskovic-Jovicevic, T.; Peric-Hajzler, Z.; Matovic, D.; Dimitrijevic, S.; Vojvodic, J.; Sijan, G.; Stepic, N.; et al. New Trends in Cutaneous Melanoma Surgery. Open Access Maced. J. Med. Sci. 2019, 7, 3090–3092. [Google Scholar] [CrossRef] [Green Version]
  10. Swayden, M.; Chhouri, H.; Anouar, Y.; Grumolato, L. Tolerant/persister cancer cells and the path to resistance to targeted therapy. Cells 2020, 9, 2601. [Google Scholar] [CrossRef]
  11. Matias, M.; Pinho, J.O.; Penetra, M.J.; Campos, G.; Reis, C.P.; Gaspar, M.M. The challenging melanoma landscape: From early drug discovery to clinical approval. Cells 2021, 10, 3088. [Google Scholar] [CrossRef]
  12. Theeuwes, W.F.; Gosker, H.R.; Langen, R.C.J.; Verhees, K.J.P.; Pansters, N.A.M.; Schols, A.M.W.J.; Remels, A.H.V. Inactivation of glycogen synthase kinase-3β (GSK-3β) enhances skeletal muscle oxidative metabolism. Biochim. Biophys. Acta BBA 2017, 1863, 3075–3086. [Google Scholar] [CrossRef] [PubMed]
  13. Domoto, T.; Uehara, M.; Bolidong, D.; Minamoto, T. Glycogen synthase kinase 3β in cancer biology and treatment. Cells 2020, 9, 1388. [Google Scholar] [CrossRef] [PubMed]
  14. He, R.; Du, S.; Lei, T.; Xie, X.; Wang, Y. Glycogen Synthase Kinase 3β in Tumorigenesis and Oncotherapy (Review). Oncol. Rep 2020, 44, 2373–2385. [Google Scholar] [CrossRef] [PubMed]
  15. Wu, X.; Liu, C.; Li, Z.; Gai, C.; Ding, D.; Chen, W.; Hao, F.; Li, W. Regulation of GSK3β/Nrf2 signaling pathway modulated erastin-induced ferroptosis in breast cancer. Mol. Cell. Biochem. 2020, 473, 217–228. [Google Scholar] [CrossRef] [PubMed]
  16. Wang, L.; Ouyang, S.; Li, B.; Wu, H.; Wang, F. GSK-3β manipulates ferroptosis sensitivity by dominating iron homeostasis. Cell Death Disc. 2021, 7, 334. [Google Scholar] [CrossRef] [PubMed]
  17. John, J.K.; Paraiso, K.H.; Rebecca, V.W.; Cantini, L.P.; Abel, E.V.; Pagano, N.; Meggers, E.; Mathew, R.; Krepler, C.; Izumi, V.; et al. GSK3β inhibition blocks melanoma cell/host interactions by downregulating N-cadherin expression and decreasing FAK phosphorylation. J. Invest. Dermat. 2012, 132, 2818–2827. [Google Scholar] [CrossRef] [Green Version]
  18. Kubic, J.D.; Mascarenhas, J.B.; Iizuka, T.; Wolfgeher, D.; Lang, D. GSK-3 promotes cell survival, growth, and PAX3 levels in human melanoma cells. Mol. Can. Res. 2012, 10, 1065–1076. [Google Scholar] [CrossRef] [Green Version]
  19. Sanghani, H.V.; Ganatra, S.H.; Pande, R. Molecular—Docking studies of potent anticancer agent. J. Comput. Sci. Syst. Biol. 2012, 5, 12–15. [Google Scholar] [CrossRef]
  20. Taghizadeh, M.S.; Niazi, A.; Moghadam, A.; Afsharifar, A. Experimental, molecular docking and molecular dynamic studies of natural products targeting overexpressed receptors in breast cancer. PLoS ONE 2022, 17, 0267961. [Google Scholar] [CrossRef]
  21. Kandasamy, S.; Sahu, S.K.; Kandasamy, K. In Silico studies on fungal metabolite against skin cancer protein (4, 5-Diarylisoxazole HSP90 Chaperone). Int. Sch. Res. Notices 2012, 2012, 626214. [Google Scholar]
  22. Ajayi, G.O.; Olagunju, J.A.; Ademuyiwa, O.; Martins, O.C. Gas chromatography-mass spectrometry analysis and phytochemical screening of ethanolic root extract of Plumbago zeylanica, Linn. J. Med. Plants Res. 2011, 5, 1756–1761. [Google Scholar]
  23. Musthafa, K.S.; Sahu, S.K.; Ravi, A.V.; Kathiresan, K. Anti-quorum sensing potential of the mangrove Rhizophora annamalayana. World J. Microbiol. Biot. 2013, 29, 1851–1858. [Google Scholar] [CrossRef] [PubMed]
  24. Patil, A.; Jadhav, V. GC-MS analysis of bioactive components from methanol leaf extract of Toddalia asiatica (L.). Int. J. Pharm. Sci. Rev. Res. 2014, 29, 18–20. [Google Scholar]
  25. Ferdosi, M.F.; Javaid, A.; Khan, I.H.; Khan, S.; Shad, N. Analysis of n-butanol flower extract of Cassia fistula through GC-MS and identification of antimicrobial compounds. Pak. J. Phytopathol. 2021, 33, 103–107. [Google Scholar] [CrossRef]
  26. Yakubu, O.E.; Udeh, C.S.; Asuk, A.A.; Moses, A.A.; Richard, B. Comparative determination of antioxidant activities and phytochemicals from fractions of ethanol extract of Senna occidentalis using GC-MS. J. Emerg. Technol. Innov. Res. 2021, 8, 195–211. [Google Scholar]
  27. Qureshi, M.Z.; Javed, S.; Javaid, A.; Al-Taie, A.H. Identification of antimicrobial compounds from n-hexane stem extract of Kochia indica by GC-MS analysis. MYCOPATH 2020, 16, 51–55. [Google Scholar]
  28. Abdel-Lateff, A.; Fisch, K.M.; Wright, A.D.; König, G.M. A new antioxidant isobenzofuranone derivative from the algicolous marine fungus Epicoccum sp. Planta Med. 2003, 69, 831–834. [Google Scholar]
  29. Ingole, S.N. Phytochemical analysis of leaf extract of Ocimum americanum L.(Lamiaceae) by GCMS method. WSN 2016, 37, 76–87. [Google Scholar]
  30. Selvamangai, G.; Bhaskar, A. GC–MS analysis of phytocomponents in the methanolic extract of Eupatorium triplinerve. Asian Pac. J. Trop. Biomed. 2012, 2, S1329–S1332. [Google Scholar] [CrossRef]
  31. Cottaz, A.; Bouarab, L.; De Clercq, J.; Oulahal, N.; Degraeve, P.; Joly, C. Potential of incorporation of antimicrobial plant phenolics into polyolefin-based food contact materials to produce active packaging by melt-blending: Proof of concept with isobutyl-4-hydroxybenzoate. Front. Chem. 2019, 7, 148. [Google Scholar] [CrossRef]
  32. Akpuaka, A.; Ekwenchi, M.M.; Dashak, D.A.; Dildar, A. Biological activities of characterized isolates of n-hexane extract of Azadirachta indica A. Juss (Neem) leaves. Nat. Sci. 2013, 11, 141–147. [Google Scholar]
  33. Kumaresan, S.; Senthilkumar, V.; Stephen, A.; Balakumar, B.S. GC-MS analysis and pass-assisted prediction of biological activity spectra of extract of Phomopsis sp. isolated from Andrographis paniculata. World J. Pharm. Res. 2015, 4, 1035–1053. [Google Scholar]
  34. Khan, M.S.; Yusufzai, S.K.; Ying, L.Y.; Zulnashriq, W. GC-MS based chemical profiling and evaluation of antioxidant potential of leaves and stems of Alternanthera sessilis red from Sabah, Malaysia. Int. J. Pharm. Sci. 2018, 10, 4–9. [Google Scholar] [CrossRef]
  35. Bagheri, G.; Ayatollahi, S.A.; Ramírez-Alarcón, K.; Fernández, M.; Salehi, B.; Forman, K.; Martorell, M.; Moghadam, M.H.; Sharifi-Rad, J. Phytochemical screening of Alstonia scholaris leaf and bark extracts and their antimicrobial activities. Cell. Mol. Biol. 2020, 66, 270–279. [Google Scholar] [CrossRef]
  36. Konappa, N.; Udayashankar, A.C.; Krishnamurthy, S.; Pradeep, C.K.; Chowdappa, S.; Jogaiah, S. GC–MS analysis of phytoconstituents from Amomum nilgiricum and molecular docking interactions of bioactive serverogenin acetate with target proteins. Sci. Rep. 2020, 10, 16438. [Google Scholar] [CrossRef] [PubMed]
  37. Huang, F.; Nau, W.M. A conformational flexibility scale for amino acids in peptides. Angew. Chem. Int. Ed. 2003, 42, 2269–2272. [Google Scholar] [CrossRef]
  38. Muhamad, N.; Muhmed, S.A.; Yusoff, M.M.; Gimbun, J. Influence of solvent polarity and conditions on extraction of antioxidant, flavonoids and phenolic content from Averrhoa bilimbi. J. Food Eng. 2014, 4, 255–260. [Google Scholar]
  39. Tomczyk, M.; Pleszczyńska, M.; Wiater, A. Variation in total polyphenolics contents of aerial parts of Potentilla species and their anticariogenic activity. Molecules 2010, 15, 4639–4651. [Google Scholar] [CrossRef]
  40. Dolkar, P.; Dolkar, D.; Angmo, S.; Kumar, B.; Stobdan, T. Variability in phenolics, flavonoids and antioxidants in seabuckthorn (Hippophae rhamnoides L.) seed from nine trans-Himalayan natural population. J. Berry Res. 2017, 7, 109–116. [Google Scholar] [CrossRef]
  41. Zhang, J.; Chen, G.; Zhao, P.; Zhou, Q.; Zhao, X. The abundance of certain metabolites responds to drought stress in the highly drought tolerant plant Caragana korshinskii. Acta Physiol. Plant 2017, 39, 116. [Google Scholar] [CrossRef]
  42. Kumar, V.; Sharma, N.; Sourirajan, A.; Khosla, P.K.; Dev, K. Correlation of biological activities of bark and leaves of Terminalia arjuna collected from different geographical regions of Himachal Pradesh, India. Biochem. Sys. Ecol. 2023, 106, 104563. [Google Scholar] [CrossRef]
  43. Sharma, S.; Kumar, V.; Attri, C.; Sourirajan, A.; Dev, K. Comparison of phytochemicals, antioxidant, and antimicrobial activities of in vitro propagated and wild grown Potentilla Nepalensis, an endemic medicinal plant from North Western Himalayas. J. Herbs Spices Med. Plants 2022, 28, 324–336. [Google Scholar] [CrossRef]
  44. Taruscio, T.G.; Barney, D.L.; Exon, J. Content and profile of flavanoid and phenolic acid compounds in conjunction with the antioxidant capacity for a variety of northwest Vaccinium berries. J. Agric. Food Chem. 2004, 52, 3169–3176. [Google Scholar] [CrossRef]
  45. Kumar, V.; Sharma, N.; Sourirajan, A.; Khosla, P.K.; Dev, K. Comparative evaluation of antimicrobial and antioxidant potential of ethanolic extract and its fractions of bark and leaves of Terminalia arjuna from north-western Himalayas, India. J. Tradit. Complement. Med. 2018, 8, 100–106. [Google Scholar] [CrossRef]
  46. Kumar, V.; Chandel, S.R.; Guleria, S.; Sharma, N.; Sourirajan, A.; Khosla, P.K.; Dev, K. Comparative analysis of phytochemicals, antimicrobial and antioxidant activity of different species of Terminalia from Himachal Pradesh, India. Vegetos 2021, 34, 528–539. [Google Scholar] [CrossRef]
  47. Singleton, V.L.; Orthofer, R.; Lamuela-Raventos, R.M. Analysis of total phenols and other oxidation substrates and antioxidants by means of Folin-Ciocalteu reagent. Meth. Enzymol. 1999, 299, 152–178. [Google Scholar]
  48. Zhishen, J.; Mengcheng, T.; Jianming, W. The determination of flavonoid contents in mulberry and their scavenging effects on superoxide radicals. Food Chem. 1999, 64, 555–559. [Google Scholar] [CrossRef]
  49. Barros, L.; Ferreira, M.J.; Queiros, B.; Ferreira, I.C.; Baptista, P. Total phenols, ascorbic acid, β-carotene and lycopene in Portuguese wild edible mushrooms and their antioxidant activities. Food Chem. 2007, 103, 413–419. [Google Scholar] [CrossRef]
  50. Sreelekha, T. Evaluation of antioxidant, antitumor and immunomodulatory properties of polysaccharide isolated from fruit rind of Punica granatum. Mol. Med. Rep. 2011, 5, 489–496. [Google Scholar] [CrossRef] [Green Version]
  51. Liang, S.H.; Chen, J.M.; Normandin, M.D.; Chang, J.S.; Chang, G.C.; Taylor, C.K.; Trapa, P.; Plummer, M.S.; Para, K.S.; Conn, E.L.; et al. Discovery of a highly selective glycogen synthase kinase-3 inhibitor (PF-04802367) that modulates tau phosphorylation in the brain: Translation for PET neuroimaging. Angew. Chem. Int. Ed. 2016, 55, 9601–9605. [Google Scholar] [CrossRef] [Green Version]
  52. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comp. Chem. 2010, 31, 455–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Banks, J.L.; Beard, H.S.; Cao, Y.; Cho, A.E.; Damm, W.; Farid, R.; Felts, A.K.; Halgren, T.A.; Mainz, D.T.; Maple, J.R.; et al. Integrated modeling program, applied chemical theory (IMPACT). J. Comput. Chem. 2005, 26, 1752–1780. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Shivakumar, D.; Williams, J.; Wu, Y.; Damm, W.; Shelley, J.; Sherman, W. Prediction of absolute solvation free energies using molecular dynamics free energy perturbation and the OPLS force field. J. Chem. Theory Comput. 2010, 6, 1509–1519. [Google Scholar] [CrossRef]
  55. Patel, C.N.; Georrge, J.J.; Modi, K.M.; Narechania, M.B.; Patel, D.P.; Gonzalez, F.J.; Pandya, H.A. Pharmacophore-based virtual screening of catechol-o-methyltransferase (COMT) inhibitors to combat Alzheimer’s disease. J. Biomol. Struct. Dyn. 2018, 36, 3938–3957. [Google Scholar] [CrossRef] [PubMed]
  56. Patel, C.N.; Kumar, S.P.; Modi, K.M.; Soni, M.N.; Modi, N.R.; Pandya, H.A. Cardiotonic steroids as potential Na+/K+-ATPase inhibitors–a computational study. J. Recept. Signal Transduct. 2019, 39, 226–234. [Google Scholar] [CrossRef]
  57. Reddy, S.V.G.; Reddy, K.T.; Kumari, V.V.; Basha, S.H. Molecular docking and dynamic simulation studies evidenced plausible immunotherapeutic anticancer property by Withaferin A targeting indoleamine 2, 3-dioxygenase. J. Biomol. Struct. Dyn. 2015, 33, 2695–2709. [Google Scholar] [CrossRef]
  58. Jorgensen, W.L.; Maxwell, D.S.; Tirado-Rives, J. Development and testing of the OPLS all-atom force field on conformational energetics and properties of organic liquids. J. Am. Chem. Soc. 1996, 118, 11225–11236. [Google Scholar] [CrossRef]
  59. Kaminski, G.A.; Friesner, R.A.; Tirado-Rives, J.; Jorgensen, W.L. Evaluation and reparametrization of the OPLS-AA force field for proteins via comparison with accurate quantum chemical calculations on peptides. J. Phys. Chem. B 2001, 105, 6474–6487. [Google Scholar] [CrossRef]
  60. Panchal, U.; Modi, K.; Dey, S.; Prajapati, U.; Patel, C.; Jain, V.K. A resorcinarene-based “turn-off” fluorescence sensor for 4-nitrotoluene: Insights from fluorescence and 1H NMR titration with computational approach. J. Lumin. 2017, 184, 74–82. [Google Scholar] [CrossRef]
  61. Mall, S.; Srivastava, R.; Sharma, N.; Patel, C.N.; Rolta, R.; Sourirajan, A.; Dev, K.; Ghosh, A.; Kumar, V. Antihypertensive activity of phytocompounds from selected medicinal plants via inhibition of angiotensin-converting enzyme (ACE) protein: An in-silico approach. Nat. Prod. Res. 2022, 36, 4526–4529. [Google Scholar] [CrossRef]
  62. Patel, C.N.; Goswami, D.; Jaiswal, D.G.; Parmar, R.M.; Solanki, H.A.; Pandya, H.A. Pinpointing the potential hits for hindering interaction of SARS-CoV-2 S-protein with ACE2 from the pool of antiviral phytochemicals utilizing molecular docking and molecular dynamics (MD) simulations. J. Mol. Graph. Model. 2021, 105, 107874. [Google Scholar] [CrossRef] [PubMed]
  63. Banerjee, P.; Eckert, A.O.; Schrey, A.K.; Preissner, R. ProTox-II: A webserver for the prediction of toxicity of chemicals. Nucleic Acids Res. 2018, 46, W257–W263. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Molecules 28 05108 g001 550
Figure 1. Concentration–effect of the roots and shoots of P. nepalensis in SK-MEL-28 cells obtained using the MTT assay (48 h exposure). Values were expressed as mean ± S.D. of three independent experiments.
Figure 1. Concentration–effect of the roots and shoots of P. nepalensis in SK-MEL-28 cells obtained using the MTT assay (48 h exposure). Values were expressed as mean ± S.D. of three independent experiments.
Molecules 28 05108 g001
Molecules 28 05108 g002 550
Figure 2. Cell morphology of SKMEL-28 cells after treatment with different concentrations of methanolic extract of roots of P. nepalensis. (A) Control; (B) 6.25 µg/mL; (C)12.5 µg/mL; (D) 25 µg/mL; (E) 50 µg/mL; (F) 100 µg/mL. Concentration-dependent reduction in cell number was observed in SKMEL-28 cells when administered with different concentrations of the extract. Cells showed alterations in morphology, such as rounding up and detachment from the surface, indicating the probable progression to apoptosis. The presence of cellular debris can also be seen, indicating the damage to cell structure.
Figure 2. Cell morphology of SKMEL-28 cells after treatment with different concentrations of methanolic extract of roots of P. nepalensis. (A) Control; (B) 6.25 µg/mL; (C)12.5 µg/mL; (D) 25 µg/mL; (E) 50 µg/mL; (F) 100 µg/mL. Concentration-dependent reduction in cell number was observed in SKMEL-28 cells when administered with different concentrations of the extract. Cells showed alterations in morphology, such as rounding up and detachment from the surface, indicating the probable progression to apoptosis. The presence of cellular debris can also be seen, indicating the damage to cell structure.
Molecules 28 05108 g002
Molecules 28 05108 g003 550
Figure 3. Cell morphology of SKMEL-28 cells after treatment with different concentrations of methanolic extract of shoots of P. nepalensis. (A) Control; (B) 6.25 µg/mL; (C)12.5 µg/mL; (D) 25 µg/mL; (E) 50 µg/mL; (F) 100 µg/mL. Concentration-dependent reduction in cell number was observed in SKMEL-28 cells when administered with different concentrations of the extract. Cells showed alterations in morphology, such as rounding up and detachment from the surface, indicating the probable progression to apoptosis. The presence of cellular debris can also be seen, indicating the damage to cell structure.
Figure 3. Cell morphology of SKMEL-28 cells after treatment with different concentrations of methanolic extract of shoots of P. nepalensis. (A) Control; (B) 6.25 µg/mL; (C)12.5 µg/mL; (D) 25 µg/mL; (E) 50 µg/mL; (F) 100 µg/mL. Concentration-dependent reduction in cell number was observed in SKMEL-28 cells when administered with different concentrations of the extract. Cells showed alterations in morphology, such as rounding up and detachment from the surface, indicating the probable progression to apoptosis. The presence of cellular debris can also be seen, indicating the damage to cell structure.
Molecules 28 05108 g003
Molecules 28 05108 g004 550
Figure 4. Cell morphology of SKMEL-28 cells after treatment with different concentrations of n-hexane extract of roots of P. nepalensis. (A) Control; (B) 6.25 µg/mL; (C)12.5 µg/mL; (D) 25 µg/mL; (E) 50 µg/mL; (F) 100 µg/mL. Concentration-dependent reduction in cell number was observed in SKMEL-28 cells when administered with different concentrations of the extract. Cells showed alterations in morphology, such as rounding up and detachment from the surface, indicating the probable progression to apoptosis. The presence of cellular debris can also be seen, indicating the damage to cell structure.
Figure 4. Cell morphology of SKMEL-28 cells after treatment with different concentrations of n-hexane extract of roots of P. nepalensis. (A) Control; (B) 6.25 µg/mL; (C)12.5 µg/mL; (D) 25 µg/mL; (E) 50 µg/mL; (F) 100 µg/mL. Concentration-dependent reduction in cell number was observed in SKMEL-28 cells when administered with different concentrations of the extract. Cells showed alterations in morphology, such as rounding up and detachment from the surface, indicating the probable progression to apoptosis. The presence of cellular debris can also be seen, indicating the damage to cell structure.
Molecules 28 05108 g004
Molecules 28 05108 g005 550
Figure 5. Cell morphology of SKMEL-28 cells after treatment with different concentrations of n-hexane extract of shoots of P. nepalensis. (A) Control; (B) 6.25 µg/mL; (C)12.5 µg/mL; (D) 25 µg/mL; (E) 50 µg/mL; (F) 100 µg/mL. Concentration-dependent reduction in cell number was observed in SKMEL-28 cells when administered with different concentrations of the extract. Cells showed alterations in morphology, such as rounding up and detachment from the surface, indicating the probable progression to apoptosis. The presence of cellular debris can also be seen, indicating the damage to cell structure.
Figure 5. Cell morphology of SKMEL-28 cells after treatment with different concentrations of n-hexane extract of shoots of P. nepalensis. (A) Control; (B) 6.25 µg/mL; (C)12.5 µg/mL; (D) 25 µg/mL; (E) 50 µg/mL; (F) 100 µg/mL. Concentration-dependent reduction in cell number was observed in SKMEL-28 cells when administered with different concentrations of the extract. Cells showed alterations in morphology, such as rounding up and detachment from the surface, indicating the probable progression to apoptosis. The presence of cellular debris can also be seen, indicating the damage to cell structure.
Molecules 28 05108 g005
Molecules 28 05108 g006 550
Figure 6. GC-MS chromatogram of roots and shoots of P. nepalensis: (A) methanolic root extract; (B) methanolic shoot extract; (C) n-hexane root extract; (D) n-hexane shoot extract. X-axis represents the time in minutes and, Y-axis represents the relative abundance of the phytocompound.
Figure 6. GC-MS chromatogram of roots and shoots of P. nepalensis: (A) methanolic root extract; (B) methanolic shoot extract; (C) n-hexane root extract; (D) n-hexane shoot extract. X-axis represents the time in minutes and, Y-axis represents the relative abundance of the phytocompound.
Molecules 28 05108 g006
Molecules 28 05108 g007 550
Figure 7. Docked pose of best docked phytocompounds from P. nepalensis against target protein (PDB ID: 5K5N): (A) 3D interactions of trichloromethyl 9-anthracenecarbodithioate with interacting amino acids of 5K5N; (B) 2D interactions of trichloromethyl 9-anthracenecarbodithioate with interacting amino acids of 5K5N; (C) 3D interactions of 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester with interacting amino acids of 5K5N; (D) 2D interactions of 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester with interacting amino acids of 5K5N.
Figure 7. Docked pose of best docked phytocompounds from P. nepalensis against target protein (PDB ID: 5K5N): (A) 3D interactions of trichloromethyl 9-anthracenecarbodithioate with interacting amino acids of 5K5N; (B) 2D interactions of trichloromethyl 9-anthracenecarbodithioate with interacting amino acids of 5K5N; (C) 3D interactions of 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester with interacting amino acids of 5K5N; (D) 2D interactions of 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester with interacting amino acids of 5K5N.
Molecules 28 05108 g007
Molecules 28 05108 g008 550
Figure 8. MD simulation of protein–ligand interaction root-mean-square deviation (RMSD) profiles of (A) trichloromethyl 9-anthracenecarbodithioate-5K5N complex and (B) 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester-5K5N complex. Color legends: Ca (blue color), side chains (green color), heavy atoms (yellow color), ligand with protein (dark pink color), ligand with ligand (pink color).
Figure 8. MD simulation of protein–ligand interaction root-mean-square deviation (RMSD) profiles of (A) trichloromethyl 9-anthracenecarbodithioate-5K5N complex and (B) 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester-5K5N complex. Color legends: Ca (blue color), side chains (green color), heavy atoms (yellow color), ligand with protein (dark pink color), ligand with ligand (pink color).
Molecules 28 05108 g008
Molecules 28 05108 g009 550
Figure 9. Analysis of RMSF trajectories over the 100 ns MD simulation for (A) trichloromethyl 9-anthracenecarbodithioate-5K5N complex and (B) 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester-5K5N complex. Protein residues that interact with ligand are marked with green-colored vertical bars.
Figure 9. Analysis of RMSF trajectories over the 100 ns MD simulation for (A) trichloromethyl 9-anthracenecarbodithioate-5K5N complex and (B) 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester-5K5N complex. Protein residues that interact with ligand are marked with green-colored vertical bars.
Molecules 28 05108 g009
Molecules 28 05108 g010 550
Figure 10. Protein–ligand contact histograms (hydrogen bonding, hydrophobic interactions, ionic interactions, and water bridge of (A) trichloromethyl 9-anthracenecarbodithioate-5K5N complex and (B) 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester-5K5N complex. Different color in bars indicates different interaction as hydrogen bond (green), hydrophobic contacts (purple) and water-bridge (blue).
Figure 10. Protein–ligand contact histograms (hydrogen bonding, hydrophobic interactions, ionic interactions, and water bridge of (A) trichloromethyl 9-anthracenecarbodithioate-5K5N complex and (B) 4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester-5K5N complex. Different color in bars indicates different interaction as hydrogen bond (green), hydrophobic contacts (purple) and water-bridge (blue).
Molecules 28 05108 g010
Table
Table 1. TPC, TFC, and DPPH free radical scavenging activity of the methanolic and n-hexane extracts of the roots and shoots of P. nepalensis.
Table 1. TPC, TFC, and DPPH free radical scavenging activity of the methanolic and n-hexane extracts of the roots and shoots of P. nepalensis.
Plant PartSolvent UsedTPC (mg g−1 GAE)TFC (mg g−1 RE)DPPH Activity
IC50 (µg mL−1)
RootsMethanol21.21 ± 0.544.24 ± 0.1723.5 ± 0.92
ShootsMethanol15.68 ± 2.792.58 ± 0.112.83 ± 0.35
Rootsn-Hexane0.90 ± 0.190.06 ± 0.0365.69 ± 0.77
Shootsn-Hexane1.59 ± 0.130.49 ± 0.0674.93 ± 1.01
Ascorbic AcidMethanol--5.86 ± 0.13
Table
Table 2. Major phytocompounds identified in GC-MS profiling of methanolic extracts of roots of P. nepalensis.
Table 2. Major phytocompounds identified in GC-MS profiling of methanolic extracts of roots of P. nepalensis.
RT (min)Area%Name of PhytocompoundMolecular FormulaBioactivity ReportedReferences
28.6112.64Tetradecanoic acid, 10,13-dimethyl-, methyl esterC17H34O2--
31.3210.58Heptadecanoic acid, 16-methyl-, methyl esterC19H38O2Effective against skin cancer [21]
27.156.441,1,1,5,7,7,7-Heptamethyl-3,3-bis (trimethyl siloxy) tetrasiloxaneC13H40O5Si6Anti-quorum sensing [22,23]
29.266.34Phthalic acid, butyl hept-4-yl esterC19H28O4--
29.355.871,1,1,3,5,5,7,7,7-Nonamethyl-3-(trimethyl siloxy) tetrasiloxaneC12H36O4Si5--
24.684.95Norepinephrine, (R)-, 4TMS derivativeC20H43NO3Si4--
31.024.632-Isopropyl-6-phenylnicotinonitrileC15H14N2--
33.354.56Cyclononasiloxane, octadecamethylC18H54O9Si9antifungal, antibacterial and/or antiviral properties [24,25]
39.034.23Oxalic acid, 2TMS derivativeC8H18O4Si2--
23.473.82-TetradecanolC14H30OIngredient in cosmetics such as cold creams for its emollient properties [26]
Table
Table 3. Major phytocompounds identified in GC-MS profiling of methanolic extracts of shoots of P. nepalensis.
Table 3. Major phytocompounds identified in GC-MS profiling of methanolic extracts of shoots of P. nepalensis.
RT (min)Area%Name of PhytocompoundMolecular FormulaBioactivity ReportedReference
28.6111.82Hexadecanoic acid, methyl esterC17H34O2Antimicrobial [27]
31.368.581,1,1,5,7,7,7-Heptamethyl-3,3-bis(trimethylsiloxy)tetrasiloxaneC13H40O5Si6Anti-quorum sensing [22,23]
29.355.54Oxalic acid, 2TMS derivativeC8H18O4Si2--
24.685.4Cyclooctasiloxane, hexadecamethyl-C16H48O8Si8--
36.115.113-Ethyl-7-hydroxyphthalideC10H10O3Antioxidant [28]
29.264.7Phthalic acid, butyl hex-3-yl esterC18H26O4Antimicrobial, insecticidal activity [29]
33.354.361,1,1,3,5,5,7,7,7-Nonamethyl-3-(trimethylsiloxy)tetrasiloxaneC12H36O4Si5--
31.024.091-Hexyl-1-nitrocyclohexaneC12H23NO2Antimicrobial and anti-inflammatory [30]
22.443.71Isobutyl 4-hydroxybenzoateC11H14O3Antimicrobial activity [31]
39.043.55Octanoic acid, 5-(acetyloxy)-, methyl esterC11H20O4--
Table
Table 4. Major phytocompounds identified in GC-MS profiling of n-hexane extracts of roots of P. nepalensis.
Table 4. Major phytocompounds identified in GC-MS profiling of n-hexane extracts of roots of P. nepalensis.
RT (min)Area%Name of PhytocompoundMolecular FormulaBioactivity ReportedReference
44.8522.9Trichloromethyl 9-anthracenecarbodithioateC16H9Cl3S2--
20.0117.71Heptane, 3,3-dimethyl-C9H20--
23.5912.28HexadecaneC16H34Antifungal, Antibacterial, antioxidant activity [32]
397.011,1,1,5,7,7,7-Heptamethyl-3, 3-bis(trimethylsiloxy)tetrasiloxaneC13H40O5Si6Anti-quorum sensing [22,23]
15.85.57DodecaneC12H26Antifungal activity [32]
26.764.99EicosaneC20H42Antimicrobial activity against clinical pathogens [33]
29.344.791,1,1,3,5,5,7,7,7-Nonamethyl-3-(trimethylsiloxy) tetrasiloxaneC12H36O4Si5- [34]
33.334.68Propanoic acid, 2-oxo-3-(trimethylsilyl)-, trimethylsilyl esterC9H20O3Si2--
43.524.31,1,1,3,5,7,7,7-Octamethyl-3,5-bis(trimethylsiloxy) tetrasiloxaneC14H42O5Si6- [35]
1122,4,6,9-DehydroadamantaneC10H12--
Table
Table 5. Major phytocompounds identified in GC-MS profiling of n-hexane extracts of shoots of P. nepalensis.
Table 5. Major phytocompounds identified in GC-MS profiling of n-hexane extracts of shoots of P. nepalensis.
RT (min)Area%Name of PhytocompoundMolecular FormulaBioactivity ReportedReference
25.0323.19Benzene, 1,3,5-tri-tert-butyl-C18H30--
32.7923.191,1,1,3,5,5,5-HeptamethyltrisiloxaneC7H22O2Si3Anti-inflammatory and antimicrobial properties [36]
46.4617.324H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl esterC12H11NO5--
42.1510.131-Propene, 3-methoxy-C4H8O--
29.248.05Phthalic acid, butyl hept-4-yl esterC19H28O4- [29]
203.14Heptane, 3,3-dimethyl-C9H20--
36.592.71N- (Methyl sulfonyl)-N, O-bis (trimethyl silyl) hydroxylamineC7H21NO3SSi2--
28.62.3Dodecanoic acid, 2-methyl-C13H26O2--
38.82.272-Acetyl-3-ethylpyrazineC8H10N2O--
26.761.87EicosaneC20H42Antimicrobial activity against food-borne pathogens [33]
Table
Table 6. Binding energy of docked phytocompounds from P. nepalensis against targeted amino acids. Binding energy was expressed in terms of Kcal/mol. The amino acids showing H-bonding with ligands were represented in red color fonts.
Table 6. Binding energy of docked phytocompounds from P. nepalensis against targeted amino acids. Binding energy was expressed in terms of Kcal/mol. The amino acids showing H-bonding with ligands were represented in red color fonts.
Sr. No.PhytocompoundsPubChem IDPDB ID: 5K5N
Binding Energy (Kcal/mol)Interacting Amino Acids
1Trichloromethyl 9-anthracenecarbodithioate613595−8.9Val A:70, Ala A:83, Lys A:85, Tyr A:134, Leu A:188, Cys A:199
24H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester619354−7.4Val A:70, Ala A:83, Met A:101, Val A:110, Leu A:130, Leu A:132, Val A:135, Leu A:188, Cys A:199, Asp A:200
3Benzene, 1,3,5-tri-tert-butyl-15089−7.1Ile A:62, Phe A:67, Val A:70, Leu A:132, Cys A:199
4Heptane, 3,3-dimethyl-520991−4.6Val A:70, Ala A:83, Lys A:85, A:132, Tyr A:134, Leu A:188, Cys A:199
5Hexadecane11006−5.1Ile A:62, Val A:70, Ala A:83, Lys A:85, Leu A:132, Tyr A:134, Leu A:188, Cys A:199
61,1,1,3,5,5,5-Heptamethyltrisiloxane6327366−3.6-
71-Propene, 3-methoxy-69392−2.8Tyr A71, Ile A: 84
8Encorafenib (drug)50922675−8.6Gly A:65, Phe A:67, Lys A:183, Asp A:200
Table
Table 7. Drug likeness and toxicity prediction of best docked phytocompounds of P. nepalensis.
Table 7. Drug likeness and toxicity prediction of best docked phytocompounds of P. nepalensis.
CompoundcLogPnrotMW
HBDHBALipinski RuleHepato-ToxicityImmuno-GenicityCarcino-GenicityCyto-
Toxicity		LD50 (mg/kg)
Trichloromethyl 9-anthracenecarbodithioate3.293371.7300YesNoNoYesNo493 mg/kg (Class IV)
4H-1-Benzopyran-2-carboxylic acid, 5-amino-6-hydroxy-4-oxo-, ethyl ester1.893249.2225YesNoNoNoNo100 mg/kg (Class III)
clogP (<5)—Measure of molecular hydrophobicity; nrot (<5)—Number of rotatable bonds; MW (<500Da)—Molecular weight; HBA (<10)—H-bond acceptor; HBD (<5)—H-bond donor; LD50—Lethal dose.
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

Sharma, S.; Kumar, V.; Yaseen, M.; S. Abouzied, A.; Arshad, A.; Bhat, M.A.; Naglah, A.M.; Patel, C.N.; Sivakumar, P.K.; Sourirajan, A.; et al. Phytochemical Analysis, In Vitro Biological Activities, and Computer-Aided Analysis of Potentilla nepalensis Hook Compounds as Potential Melanoma Inhibitors Based on Molecular Docking, MD Simulations, and ADMET. Molecules 2023, 28, 5108. https://doi.org/10.3390/molecules28135108

AMA Style

Sharma S, Kumar V, Yaseen M, S. Abouzied A, Arshad A, Bhat MA, Naglah AM, Patel CN, Sivakumar PK, Sourirajan A, et al. Phytochemical Analysis, In Vitro Biological Activities, and Computer-Aided Analysis of Potentilla nepalensis Hook Compounds as Potential Melanoma Inhibitors Based on Molecular Docking, MD Simulations, and ADMET. Molecules. 2023; 28(13):5108. https://doi.org/10.3390/molecules28135108

Chicago/Turabian Style

Sharma, Subhash, Vikas Kumar, Muhammad Yaseen, Amr S. Abouzied, Abgeena Arshad, Mashooq Ahmad Bhat, Ahmed M. Naglah, Chirag N. Patel, Prasanth Kumar Sivakumar, Anuradha Sourirajan, and et al. 2023. "Phytochemical Analysis, In Vitro Biological Activities, and Computer-Aided Analysis of Potentilla nepalensis Hook Compounds as Potential Melanoma Inhibitors Based on Molecular Docking, MD Simulations, and ADMET" Molecules 28, no. 13: 5108. https://doi.org/10.3390/molecules28135108

APA Style

Sharma, S., Kumar, V., Yaseen, M., S. Abouzied, A., Arshad, A., Bhat, M. A., Naglah, A. M., Patel, C. N., Sivakumar, P. K., Sourirajan, A., Shahzad, A., & Dev, K. (2023). Phytochemical Analysis, In Vitro Biological Activities, and Computer-Aided Analysis of Potentilla nepalensis Hook Compounds as Potential Melanoma Inhibitors Based on Molecular Docking, MD Simulations, and ADMET. Molecules, 28(13), 5108. https://doi.org/10.3390/molecules28135108

Article Metrics

Back to TopTop